#include #if defined(PETSC_NEEDS_UTYPE_TYPEDEFS) /* Some systems have inconsistent include files that use but do not ensure that the following definitions are made */ typedef unsigned char u_char; typedef unsigned short u_short; typedef unsigned short ushort; typedef unsigned int u_int; typedef unsigned long u_long; #endif #include #include #if defined(PETSC_HAVE_MACHINE_ENDIAN_H) #include #endif #if defined(PETSC_HAVE_UNISTD_H) #include #endif #if defined(PETSC_HAVE_SYS_SOCKET_H) #include #endif #if defined(PETSC_HAVE_SYS_WAIT_H) #include #endif #if defined(PETSC_HAVE_NETINET_IN_H) #include #endif #if defined(PETSC_HAVE_NETDB_H) #include #endif #if defined(PETSC_HAVE_FCNTL_H) #include #endif #if defined(PETSC_HAVE_IO_H) #include #endif #if defined(PETSC_HAVE_WINSOCK2_H) #include #endif #include #include <../src/sys/classes/viewer/impls/socket/socket.h> #if defined(PETSC_NEED_CLOSE_PROTO) PETSC_EXTERN int close(int); #endif #if defined(PETSC_NEED_SOCKET_PROTO) PETSC_EXTERN int socket(int,int,int); #endif #if defined(PETSC_NEED_SLEEP_PROTO) PETSC_EXTERN int sleep(unsigned); #endif #if defined(PETSC_NEED_CONNECT_PROTO) PETSC_EXTERN int connect(int,struct sockaddr*,int); #endif /*--------------------------------------------------------------*/ #undef __FUNCT__ #define __FUNCT__ "PetscViewerDestroy_Socket" static PetscErrorCode PetscViewerDestroy_Socket(PetscViewer viewer) { PetscViewer_Socket *vmatlab = (PetscViewer_Socket*)viewer->data; PetscErrorCode ierr; PetscFunctionBegin; if (vmatlab->port) { #if defined(PETSC_HAVE_CLOSESOCKET) ierr = closesocket(vmatlab->port); #else ierr = close(vmatlab->port); #endif if (ierr) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"System error closing socket"); } ierr = PetscFree(vmatlab);CHKERRQ(ierr); PetscFunctionReturn(0); } /*--------------------------------------------------------------*/ #undef __FUNCT__ #define __FUNCT__ "PetscOpenSocket" /* PetscSocketOpen - handles connected to an open port where someone is waiting. .seealso: PetscSocketListen(), PetscSocketEstablish() */ PetscErrorCode PetscOpenSocket(char *hostname,int portnum,int *t) { struct sockaddr_in sa; struct hostent *hp; int s = 0; PetscErrorCode ierr; PetscBool flg = PETSC_TRUE; PetscFunctionBegin; if (!(hp=gethostbyname(hostname))) { perror("SEND: error gethostbyname: "); SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_SYS,"system error open connection to %s",hostname); } ierr = PetscMemzero(&sa,sizeof(sa));CHKERRQ(ierr); ierr = PetscMemcpy(&sa.sin_addr,hp->h_addr_list[0],hp->h_length);CHKERRQ(ierr); sa.sin_family = hp->h_addrtype; sa.sin_port = htons((u_short) portnum); while (flg) { if ((s=socket(hp->h_addrtype,SOCK_STREAM,0)) < 0) { perror("SEND: error socket"); SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"system error"); } if (connect(s,(struct sockaddr*)&sa,sizeof(sa)) < 0) { #if defined(PETSC_HAVE_WSAGETLASTERROR) ierr = WSAGetLastError(); if (ierr == WSAEADDRINUSE) (*PetscErrorPrintf)("SEND: address is in use\n"); else if (ierr == WSAEALREADY) (*PetscErrorPrintf)("SEND: socket is non-blocking \n"); else if (ierr == WSAEISCONN) { (*PetscErrorPrintf)("SEND: socket already connected\n"); Sleep((unsigned) 1); } else if (ierr == WSAECONNREFUSED) { /* (*PetscErrorPrintf)("SEND: forcefully rejected\n"); */ Sleep((unsigned) 1); } else { perror(NULL); SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"system error"); } #else if (errno == EADDRINUSE) (*PetscErrorPrintf)("SEND: address is in use\n"); else if (errno == EALREADY) (*PetscErrorPrintf)("SEND: socket is non-blocking \n"); else if (errno == EISCONN) { (*PetscErrorPrintf)("SEND: socket already connected\n"); sleep((unsigned) 1); } else if (errno == ECONNREFUSED) { /* (*PetscErrorPrintf)("SEND: forcefully rejected\n"); */ ierr = PetscInfo(0,"Connection refused in attaching socket, trying again");CHKERRQ(ierr); sleep((unsigned) 1); } else { perror(NULL); SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"system error"); } #endif flg = PETSC_TRUE; #if defined(PETSC_HAVE_CLOSESOCKET) closesocket(s); #else close(s); #endif } else flg = PETSC_FALSE; } *t = s; PetscFunctionReturn(0); } #define MAXHOSTNAME 100 #undef __FUNCT__ #define __FUNCT__ "PetscSocketEstablish" /* PetscSocketEstablish - starts a listener on a socket .seealso: PetscSocketListen() */ PetscErrorCode PetscSocketEstablish(int portnum,int *ss) { char myname[MAXHOSTNAME+1]; int s; PetscErrorCode ierr; struct sockaddr_in sa; struct hostent *hp; PetscFunctionBegin; ierr = PetscGetHostName(myname,MAXHOSTNAME);CHKERRQ(ierr); ierr = PetscMemzero(&sa,sizeof(struct sockaddr_in));CHKERRQ(ierr); hp = gethostbyname(myname); if (!hp) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"Unable to get hostent information from system"); sa.sin_family = hp->h_addrtype; sa.sin_port = htons((u_short)portnum); if ((s = socket(AF_INET,SOCK_STREAM,0)) < 0) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"Error running socket() command"); #if defined(PETSC_HAVE_SO_REUSEADDR) { int optval = 1; /* Turn on the option */ ierr = setsockopt(s,SOL_SOCKET,SO_REUSEADDR,(char*)&optval,sizeof(optval));CHKERRQ(ierr); } #endif while (bind(s,(struct sockaddr*)&sa,sizeof(sa)) < 0) { #if defined(PETSC_HAVE_WSAGETLASTERROR) ierr = WSAGetLastError(); if (ierr != WSAEADDRINUSE) { #else if (errno != EADDRINUSE) { #endif close(s); SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"Error from bind()"); } } listen(s,0); *ss = s; return(0); } #undef __FUNCT__ #define __FUNCT__ "PetscSocketListen" /* PetscSocketListens - Listens at a socket created with PetscSocketEstablish() .seealso: PetscSocketEstablish() */ PetscErrorCode PetscSocketListen(int listenport,int *t) { struct sockaddr_in isa; #if defined(PETSC_HAVE_ACCEPT_SIZE_T) size_t i; #else int i; #endif PetscFunctionBegin; /* wait for someone to try to connect */ i = sizeof(struct sockaddr_in); if ((*t = accept(listenport,(struct sockaddr*)&isa,(socklen_t*)&i)) < 0) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SYS,"error from accept()\n"); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscViewerSocketOpen" /*@C PetscViewerSocketOpen - Opens a connection to a MATLAB or other socket based server. Collective on MPI_Comm Input Parameters: + comm - the MPI communicator . machine - the machine the server is running on,, use NULL for the local machine, use "server" to passively wait for a connection from elsewhere - port - the port to connect to, use PETSC_DEFAULT for the default Output Parameter: . lab - a context to use when communicating with the server Level: intermediate Notes: Most users should employ the following commands to access the MATLAB PetscViewers $ $ PetscViewerSocketOpen(MPI_Comm comm, char *machine,int port,PetscViewer &viewer) $ MatView(Mat matrix,PetscViewer viewer) $ $ or $ $ PetscViewerSocketOpen(MPI_Comm comm,char *machine,int port,PetscViewer &viewer) $ VecView(Vec vector,PetscViewer viewer) Options Database Keys: For use with PETSC_VIEWER_SOCKET_WORLD, PETSC_VIEWER_SOCKET_SELF, PETSC_VIEWER_SOCKET_() or if NULL is passed for machine or PETSC_DEFAULT is passed for port $ -viewer_socket_machine $ -viewer_socket_port Environmental variables: + PETSC_VIEWER_SOCKET_PORT portnumber - PETSC_VIEWER_SOCKET_MACHINE machine name Currently the only socket client available is MATLAB. See src/dm/da/examples/tests/ex12.c and ex12.m for an example of usage. Notes: The socket viewer is in some sense a subclass of the binary viewer, to read and write to the socket use PetscViewerBinaryRead/Write/GetDescriptor(). Concepts: MATLAB^sending data Concepts: sockets^sending data .seealso: MatView(), VecView(), PetscViewerDestroy(), PetscViewerCreate(), PetscViewerSetType(), PetscViewerSocketSetConnection(), PETSC_VIEWER_SOCKET_, PETSC_VIEWER_SOCKET_WORLD, PETSC_VIEWER_SOCKET_SELF, PetscViewerBinaryWrite(), PetscViewerBinaryRead(), PetscViewerBinaryWriteStringArray(), PetscBinaryViewerGetDescriptor() @*/ PetscErrorCode PetscViewerSocketOpen(MPI_Comm comm,const char machine[],int port,PetscViewer *lab) { PetscErrorCode ierr; PetscFunctionBegin; ierr = PetscViewerCreate(comm,lab);CHKERRQ(ierr); ierr = PetscViewerSetType(*lab,PETSCVIEWERSOCKET);CHKERRQ(ierr); ierr = PetscViewerSocketSetConnection(*lab,machine,port);CHKERRQ(ierr); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscViewerSetFromOptions_Socket" static PetscErrorCode PetscViewerSetFromOptions_Socket(PetscViewer v) { PetscErrorCode ierr; PetscInt def = -1; char sdef[256]; PetscBool tflg; PetscFunctionBegin; /* These options are not processed here, they are processed in PetscViewerSocketSetConnection(), they are listed here for the GUI to display */ ierr = PetscOptionsHead("Socket PetscViewer Options");CHKERRQ(ierr); ierr = PetscOptionsGetenv(PetscObjectComm((PetscObject)v),"PETSC_VIEWER_SOCKET_PORT",sdef,16,&tflg);CHKERRQ(ierr); if (tflg) { ierr = PetscOptionsStringToInt(sdef,&def);CHKERRQ(ierr); } else def = PETSCSOCKETDEFAULTPORT; ierr = PetscOptionsInt("-viewer_socket_port","Port number to use for socket","PetscViewerSocketSetConnection",def,0,0);CHKERRQ(ierr); ierr = PetscOptionsString("-viewer_socket_machine","Machine to use for socket","PetscViewerSocketSetConnection",sdef,0,0,0);CHKERRQ(ierr); ierr = PetscOptionsGetenv(PetscObjectComm((PetscObject)v),"PETSC_VIEWER_SOCKET_MACHINE",sdef,256,&tflg);CHKERRQ(ierr); if (!tflg) { ierr = PetscGetHostName(sdef,256);CHKERRQ(ierr); } ierr = PetscOptionsTail();CHKERRQ(ierr); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscViewerCreate_Socket" PETSC_EXTERN PetscErrorCode PetscViewerCreate_Socket(PetscViewer v) { PetscViewer_Socket *vmatlab; PetscErrorCode ierr; PetscFunctionBegin; ierr = PetscNewLog(v,PetscViewer_Socket,&vmatlab);CHKERRQ(ierr); vmatlab->port = 0; v->data = (void*)vmatlab; v->ops->destroy = PetscViewerDestroy_Socket; v->ops->flush = 0; v->ops->setfromoptions = PetscViewerSetFromOptions_Socket; /* lie and say this is a binary viewer; then all the XXXView_Binary() methods will work correctly on it */ ierr = PetscObjectChangeTypeName((PetscObject)v,PETSCVIEWERBINARY);CHKERRQ(ierr); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscViewerSocketSetConnection" /*@C PetscViewerSocketSetConnection - Sets the machine and port that a PETSc socket viewer is to use Logically Collective on PetscViewer Input Parameters: + v - viewer to connect . machine - host to connect to, use NULL for the local machine,use "server" to passively wait for a connection from elsewhere - port - the port on the machine one is connecting to, use PETSC_DEFAULT for default Level: advanced .seealso: PetscViewerSocketOpen() @*/ PetscErrorCode PetscViewerSocketSetConnection(PetscViewer v,const char machine[],int port) { PetscErrorCode ierr; PetscMPIInt rank; char mach[256]; PetscBool tflg; PetscViewer_Socket *vmatlab = (PetscViewer_Socket*)v->data; PetscFunctionBegin; /* PetscValidLogicalCollectiveInt(v,port,3); not a PetscInt */ if (port <= 0) { char portn[16]; ierr = PetscOptionsGetenv(PetscObjectComm((PetscObject)v),"PETSC_VIEWER_SOCKET_PORT",portn,16,&tflg);CHKERRQ(ierr); if (tflg) { PetscInt pport; ierr = PetscOptionsStringToInt(portn,&pport);CHKERRQ(ierr); port = (int)pport; } else port = PETSCSOCKETDEFAULTPORT; } if (!machine) { ierr = PetscOptionsGetenv(PetscObjectComm((PetscObject)v),"PETSC_VIEWER_SOCKET_MACHINE",mach,256,&tflg);CHKERRQ(ierr); if (!tflg) { ierr = PetscGetHostName(mach,256);CHKERRQ(ierr); } } else { ierr = PetscStrncpy(mach,machine,256);CHKERRQ(ierr); } ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)v),&rank);CHKERRQ(ierr); if (!rank) { ierr = PetscStrcmp(mach,"server",&tflg);CHKERRQ(ierr); if (tflg) { int listenport; ierr = PetscInfo1(v,"Waiting for connection from socket process on port %D\n",port);CHKERRQ(ierr); ierr = PetscSocketEstablish(port,&listenport);CHKERRQ(ierr); ierr = PetscSocketListen(listenport,&vmatlab->port);CHKERRQ(ierr); close(listenport); } else { ierr = PetscInfo2(v,"Connecting to socket process on port %D machine %s\n",port,mach);CHKERRQ(ierr); ierr = PetscOpenSocket(mach,port,&vmatlab->port);CHKERRQ(ierr); } } PetscFunctionReturn(0); } /* ---------------------------------------------------------------------*/ /* The variable Petsc_Viewer_Socket_keyval is used to indicate an MPI attribute that is attached to a communicator, in this case the attribute is a PetscViewer. */ static PetscMPIInt Petsc_Viewer_Socket_keyval = MPI_KEYVAL_INVALID; #undef __FUNCT__ #define __FUNCT__ "PETSC_VIEWER_SOCKET_" /*@C PETSC_VIEWER_SOCKET_ - Creates a socket viewer shared by all processors in a communicator. Collective on MPI_Comm Input Parameter: . comm - the MPI communicator to share the socket PetscViewer Level: intermediate Options Database Keys: For use with the default PETSC_VIEWER_SOCKET_WORLD or if NULL is passed for machine or PETSC_DEFAULT is passed for port $ -viewer_socket_machine $ -viewer_socket_port Environmental variables: + PETSC_VIEWER_SOCKET_PORT portnumber - PETSC_VIEWER_SOCKET_MACHINE machine name Notes: Unlike almost all other PETSc routines, PetscViewer_SOCKET_ does not return an error code. The socket PetscViewer is usually used in the form $ XXXView(XXX object,PETSC_VIEWER_SOCKET_(comm)); Currently the only socket client available is MATLAB. See src/dm/da/examples/tests/ex12.c and ex12.m for an example of usage. Connects to a waiting socket and stays connected until PetscViewerDestroy() is called. Use this for communicating with an interactive MATLAB session, see PETSC_VIEWER_MATLAB_() for communicating with the MATLAB engine. .seealso: PETSC_VIEWER_SOCKET_WORLD, PETSC_VIEWER_SOCKET_SELF, PetscViewerSocketOpen(), PetscViewerCreate(), PetscViewerSocketSetConnection(), PetscViewerDestroy(), PETSC_VIEWER_SOCKET_(), PetscViewerBinaryWrite(), PetscViewerBinaryRead(), PetscViewerBinaryWriteStringArray(), PetscBinaryViewerGetDescriptor(), PETSC_VIEWER_MATLAB_() @*/ PetscViewer PETSC_VIEWER_SOCKET_(MPI_Comm comm) { PetscErrorCode ierr; PetscBool flg; PetscViewer viewer; MPI_Comm ncomm; PetscFunctionBegin; ierr = PetscCommDuplicate(comm,&ncomm,NULL);if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} if (Petsc_Viewer_Socket_keyval == MPI_KEYVAL_INVALID) { ierr = MPI_Keyval_create(MPI_NULL_COPY_FN,MPI_NULL_DELETE_FN,&Petsc_Viewer_Socket_keyval,0); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} } ierr = MPI_Attr_get(ncomm,Petsc_Viewer_Socket_keyval,(void**)&viewer,(int*)&flg); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} if (!flg) { /* PetscViewer not yet created */ ierr = PetscViewerSocketOpen(ncomm,0,0,&viewer); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} ierr = PetscObjectRegisterDestroy((PetscObject)viewer); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} ierr = MPI_Attr_put(ncomm,Petsc_Viewer_Socket_keyval,(void*)viewer); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} } ierr = PetscCommDestroy(&ncomm); if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_SOCKET_",__FILE__,__SDIR__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," ");PetscFunctionReturn(0);} PetscFunctionReturn(viewer); } /* ---------------------------------------------------------------------*/ #if defined(PETSC_USE_SERVER) /* Implements a crude webserver allowing the snooping on running application codes. Developer Notes: Most of this code, including the webserver, perhaps, belongs properly in the AMS with perhaps a few hooks for application/libraries like PETSc to interact with it. */ #include #include #define PROTOCOL "HTTP/1.1" #define RFC1123FMT "%a, %d %b %Y %H:%M:%S GMT" #undef __FUNCT__ #define __FUNCT__ "PetscWebSendHeader" PetscErrorCode PetscWebSendHeader(FILE *f, int status, const char *title, const char *extra, const char *mime, int length) { time_t now; char timebuf[128]; PetscFunctionBegin; fprintf(f, "%s %d %s\r\n", PROTOCOL, status, title); fprintf(f, "Server: %s\r\n", "petscserver/1.0"); now = time(NULL); strftime(timebuf, sizeof(timebuf), RFC1123FMT, gmtime(&now)); fprintf(f, "Date: %s\r\n", timebuf); if (extra) fprintf(f, "%s\r\n", extra); if (mime) fprintf(f, "Content-Type: %s\r\n", mime); if (length >= 0) fprintf(f, "Content-Length: %d\r\n", length); fprintf(f, "Connection: close\r\n"); fprintf(f, "\r\n"); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscWebSendFooter" PetscErrorCode PetscWebSendFooter(FILE *fd) { PetscFunctionBegin; fprintf(fd, "\r\n"); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscWebSendError" PetscErrorCode PetscWebSendError(FILE *f, int status, const char *title, const char *extra, const char *text) { PetscErrorCode ierr; PetscFunctionBegin; ierr = PetscWebSendHeader(f, status, title, extra, "text/html", -1);CHKERRQ(ierr); fprintf(f, "%d %s\r\n", status, title); fprintf(f, "

%d %s

\r\n", status, title); fprintf(f, "%s\r\n", text); ierr = PetscWebSendFooter(f);CHKERRQ(ierr); PetscFunctionReturn(0); } #if defined(PETSC_HAVE_AMS) #include #undef __FUNCT__ #define __FUNCT__ "PetscAMSObjectsDisplayList" /* Displays all the PETSc objects published with AMS in a simple HTML list Does NOT use Javascript or JSON-RPC */ static PetscErrorCode PetscAMSObjectsDisplayList(FILE *fd) { PetscErrorCode ierr; char host[256],**comm_list,**mem_list,**fld_list; AMS_Comm ams; PetscInt i = 0,j; AMS_Memory_type mtype; AMS_Data_type dtype; AMS_Shared_type stype; AMS_Reduction_type rtype; AMS_Memory memory; int len; void *addr; ierr = PetscGetHostName(host,256);CHKERRQ(ierr); PetscStackCallAMS(AMS_Connect,(host, -1, &comm_list)); ierr = PetscWebSendHeader(fd, 200, "OK", NULL, "text/html", -1);CHKERRQ(ierr); if (!comm_list || !comm_list[0]) fprintf(fd, "AMS Communicator not running

\r\n"); else { PetscStackCallAMS(AMS_Comm_attach,(comm_list[0],&ams)); PetscStackCallAMS(AMS_Comm_get_memory_list,(ams,&mem_list)); if (!mem_list[0]) fprintf(fd, "AMS Communicator %s has no published memories

\r\n",comm_list[0]); else { fprintf(fd, "Petsc Application Server\r\n"); fprintf(fd,"
    \r\n"); while (mem_list[i]) { fprintf(fd,"
  • %s
    • \r\n",mem_list[i]); PetscStackCallAMS(AMS_Memory_attach,(ams,mem_list[i],&memory,NULL)); PetscStackCallAMS(AMS_Memory_get_field_list,(memory, &fld_list)); j = 0; fprintf(fd,"
      \r\n"); while (fld_list[j]) { fprintf(fd,"
    • %s",fld_list[j]); PetscStackCallAMS(AMS_Memory_get_field_info,(memory, fld_list[j], &addr, &len, &dtype, &mtype, &stype, &rtype)); if (len == 1) { if (dtype == AMS_INT) fprintf(fd," %d",*(int*)addr); else if (dtype == AMS_STRING) fprintf(fd," %s",*(char**)addr); } fprintf(fd,"
      • \r\n"); j++; } fprintf(fd,"
    \r\n"); i++; } fprintf(fd,"
\r\n"); } } ierr = PetscWebSendFooter(fd);CHKERRQ(ierr); PetscStackCallAMS(AMS_Disconnect,()); PetscFunctionReturn(0); } #undef __FUNCT__ #define __FUNCT__ "PetscAMSObjectsDisplayTree" /* Displays all the PETSc objects published with AMS in very crude HTML 5 graphics Does NOT use Javascript or JSON-RPC */ static PetscErrorCode PetscAMSObjectsDisplayTree(FILE *fd) { PetscErrorCode ierr; char host[256],**comm_list,**mem_list,**fld_list; AMS_Comm ams; PetscInt i = 0,j; AMS_Memory_type mtype; AMS_Data_type dtype; AMS_Shared_type stype; AMS_Reduction_type rtype; AMS_Memory memory; int len; void *addr2,*addr3,*addr,*addr4; ierr = PetscGetHostName(host,256);CHKERRQ(ierr); PetscStackCallAMS(AMS_Connect,(host, -1, &comm_list)); ierr = PetscWebSendHeader(fd, 200, "OK", NULL, "text/html", -1);CHKERRQ(ierr); if (!comm_list || !comm_list[0]) fprintf(fd, "AMS Communicator not running

\r\n"); else { PetscStackCallAMS(AMS_Comm_attach,(comm_list[0],&ams)); PetscStackCallAMS(AMS_Comm_get_memory_list,(ams,&mem_list)); if (!mem_list[0]) fprintf(fd, "AMS Communicator %s has no published memories

\r\n",comm_list[0]); else { PetscInt Nlevels,*Level,*Levelcnt,*Idbylevel,*Column,*parentid,*Id,maxId = 0,maxCol = 0,*parentId,id,cnt,Nlevelcnt = 0; PetscBool *mask; char **classes,*clas,**subclasses,*sclas; /* get maximum number of objects */ while (mem_list[i]) { PetscStackCallAMS(AMS_Memory_attach,(ams,mem_list[i],&memory,NULL)); PetscStackCallAMS(AMS_Memory_get_field_list,(memory, &fld_list)); PetscStackCallAMS(AMS_Memory_get_field_info,(memory, "Id", &addr2, &len, &dtype, &mtype, &stype, &rtype)); Id = (int*) addr2; maxId = PetscMax(maxId,*Id); i++; } maxId++; /* Gets everyone's parent ID and which nodes are masked */ ierr = PetscMalloc4(maxId,PetscInt,&parentid,maxId,PetscBool,&mask,maxId,char**,&classes,maxId,char**,&subclasses);CHKERRQ(ierr); ierr = PetscMemzero(classes,maxId*sizeof(char*));CHKERRQ(ierr); ierr = PetscMemzero(subclasses,maxId*sizeof(char*));CHKERRQ(ierr); for (i=0; i 0 && mask[parentid[i]]) parentid[i] = 0; } ierr = PetscProcessTree(maxId,mask,parentid,&Nlevels,&Level,&Levelcnt,&Idbylevel,&Column);CHKERRQ(ierr); for (i=0; iPetsc Application Server\r\n"); fprintf(fd, "\r\n"); fprintf(fd, "\r\n"); fprintf(fd, "\r\n"); fprintf(fd, "\r\n"); } } ierr = PetscWebSendFooter(fd);CHKERRQ(ierr); PetscFunctionReturn(0); } #endif #undef __FUNCT__ #define __FUNCT__ "PetscWebServeRequestGet" /*@C PetscWebServeRequestGet - serves a single web Get request Not collective Input Parameters: + port - the network file to read and write from - path - the command from the server Level: developer .seealso: PetscWebServe() @*/ static PetscErrorCode PetscWebServeRequestGet(FILE *fd,const char path[]) { PetscErrorCode ierr; FILE *fdo; char fullpath[PETSC_MAX_PATH_LEN],truefullpath[PETSC_MAX_PATH_LEN],*qmark; const char *type; PetscBool flg; PetscFunctionBegin; fseek(fd, 0, SEEK_CUR); /* Force change of stream direction */ ierr = PetscStrcmp(path,"/favicon.ico",&flg);CHKERRQ(ierr); if (flg) { /* should have cool PETSc icon */; PetscFunctionReturn(0); } ierr = PetscStrcmp(path,"/",&flg);CHKERRQ(ierr); if (flg) { char program[128]; PetscMPIInt size; PetscViewer viewer; ierr = MPI_Comm_size(PETSC_COMM_WORLD,&size);CHKERRQ(ierr); ierr = PetscGetProgramName(program,128);CHKERRQ(ierr); ierr = PetscWebSendHeader(fd, 200, "OK", NULL, "text/html", -1);CHKERRQ(ierr); fprintf(fd, "Petsc Application Server\r\n"); fprintf(fd, "

Serving PETSc application code %s

\r\n\n",program); fprintf(fd, "Number of processes %d\r\n\n",size); fprintf(fd, "
\r\n"); ierr = PetscViewerASCIIOpenWithFILE(PETSC_COMM_WORLD,fd,&viewer);CHKERRQ(ierr); ierr = PetscOptionsView(viewer);CHKERRQ(ierr); ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr); fprintf(fd, "
\r\n"); #if defined(PETSC_HAVE_AMS) fprintf(fd, "Connect to Memory Snooper--Tree Display

\r\n\r\n"); fprintf(fd, "Connect to Memory Snooper--List Display

\r\n\r\n"); fprintf(fd, "Connect to Memory Snooper--Interactive Javascript

\r\n\r\n"); #endif ierr = PetscWebSendFooter(fd);CHKERRQ(ierr); PetscFunctionReturn(0); } #if defined(PETSC_HAVE_AMS) ierr = PetscStrcmp(path,"/ams-list",&flg);CHKERRQ(ierr); if (flg) { ierr = PetscAMSObjectsDisplayList(fd);CHKERRQ(ierr); PetscFunctionReturn(0); } ierr = PetscStrcmp(path,"/ams-tree",&flg);CHKERRQ(ierr); if (flg) { ierr = PetscAMSObjectsDisplayTree(fd);CHKERRQ(ierr); PetscFunctionReturn(0); } #endif ierr = PetscStrcpy(fullpath,"${PETSC_DIR}/include/web");CHKERRQ(ierr); ierr = PetscStrcat(fullpath,path);CHKERRQ(ierr); ierr = PetscInfo1(NULL,"Checking for file %s\n",fullpath);CHKERRQ(ierr); ierr = PetscStrstr(fullpath,"?",&qmark);CHKERRQ(ierr); if (qmark) *qmark = 0; ierr = PetscStrreplace(PETSC_COMM_SELF,fullpath,truefullpath,PETSC_MAX_PATH_LEN);CHKERRQ(ierr); fdo = fopen(truefullpath,"r"); if (fdo) { PetscInt length,index; char data[4096]; struct stat statbuf; int n; const char *suffixes[] = {".html",".js",".gif",0}, *mimes[] = {"text/html","text/javascript","image/gif","text/unknown"}; ierr = PetscStrendswithwhich(fullpath,suffixes,&index);CHKERRQ(ierr); type = mimes[index]; if (!stat(truefullpath, &statbuf)) length = -1; else length = S_ISREG(statbuf.st_mode) ? statbuf.st_size : -1; ierr = PetscWebSendHeader(fd, 200, "OK", NULL, type, length);CHKERRQ(ierr); while ((n = fread(data, 1, sizeof(data), fdo)) > 0) fwrite(data, 1, n, fd); fclose(fdo); ierr = PetscInfo2(NULL,"Sent file %s to browser using format %s\n",fullpath,type);CHKERRQ(ierr); PetscFunctionReturn(0); } ierr = PetscWebSendError(fd, 501, "Not supported", NULL, "Unknown request.");CHKERRQ(ierr); PetscFunctionReturn(0); } #if defined(PETSC_HAVE_YAML) /* Toy YAML/JSON-RPC function that returns all the arguments it is passed */ #undef __FUNCT__ #define __FUNCT__ "YAML_echo" PETSC_UNUSED static PetscErrorCode YAML_echo(PetscInt argc,char **args,PetscInt *argco,char ***argso) { PetscErrorCode ierr; PetscInt i; ierr = PetscPrintf(PETSC_COMM_SELF,"Number of arguments to function %d\n",argc);CHKERRQ(ierr); for (i=0; i 1) {ierr = PetscStrcat(*result,"[");CHKERRQ(ierr);} for (i=0; i 1) {ierr = PetscStrcat(*result,"]");CHKERRQ(ierr);} ierr = PetscStrcat(*result,"}");CHKERRQ(ierr); ierr = PetscInfo1(NULL,"YAML/JSON result of function %s\n",*result);CHKERRQ(ierr); /* free work space */ ierr = PetscFree(idname);CHKERRQ(ierr); for (i=0; i - start PETSc webserver (default port is 8080) - -xxx_view ams - publish object xxx to be accessible in the server Notes: Point your browser to http://hostname:8080 to access the PETSc web server, where hostname is the name of your machine. If you are running PETSc on your local machine you can use http://localhost:8080 If the PETSc program completes before you connect with the browser you will not be able to connect to the PETSc webserver. Read the top of $PETSC_DIR/include/web/AMSJavascript.py before running. Level: intermediate .seealso: PetscViewerSocketOpen() @*/ PetscErrorCode PetscWebServe(MPI_Comm comm,int port) { PetscErrorCode ierr; PetscMPIInt rank; pthread_t thread; int *trueport; PetscBool flg; PetscFunctionBegin; if (port < 1 && port != PETSC_DEFAULT && port != PETSC_DECIDE) SETERRQ1(PETSC_COMM_WORLD,PETSC_ERR_ARG_WRONG,"Cannot use negative port number %d",port); ierr = PetscMallocGetDebug(&flg);CHKERRQ(ierr); if (flg) SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"Cannot use PetscWebServe() (-server) with any malloc debugging, run with -malloc off -malloc_test off -malloc_debug off"); ierr = MPI_Comm_rank(comm,&rank);CHKERRQ(ierr); if (rank) PetscFunctionReturn(0); if (port == PETSC_DECIDE || port == PETSC_DEFAULT) port = 8080; ierr = PetscMalloc(1*sizeof(int),&trueport);CHKERRQ(ierr); /* malloc this so it still exists in thread */ *trueport = port; ierr = pthread_create(&thread, NULL, (void *(*)(void*))PetscWebServeWait, trueport);CHKERRQ(ierr); PetscFunctionReturn(0); } #endif