xref: /petsc/src/sys/fileio/smatlab.c (revision 3ba1676111f5c958fe6c2729b46ca4d523958bb3)
1e5c89e4eSSatish Balay 
2c6db04a5SJed Brown #include <petscsys.h>
3e5c89e4eSSatish Balay 
4e5c89e4eSSatish Balay /*@C
5e3c5b3baSBarry Smith     PetscStartMatlab - starts up MATLAB with a MATLAB script
6e5c89e4eSSatish Balay 
7d083f849SBarry Smith     Logically Collective, but only processor zero in the communicator does anything
8e5c89e4eSSatish Balay 
9e5c89e4eSSatish Balay     Input Parameters:
10e5c89e4eSSatish Balay +     comm - MPI communicator
11e3c5b3baSBarry Smith .     machine - optional machine to run MATLAB on
12e5c89e4eSSatish Balay -     script - name of script (without the .m)
13e5c89e4eSSatish Balay 
14e5c89e4eSSatish Balay     Output Parameter:
15811af0c4SBarry Smith .     fp - a file pointer returned from `PetscPOpen()`
16e5c89e4eSSatish Balay 
17e5c89e4eSSatish Balay     Level: intermediate
18e5c89e4eSSatish Balay 
19e5c89e4eSSatish Balay     Notes:
20e5c89e4eSSatish Balay      This overwrites your matlab/startup.m file
21e5c89e4eSSatish Balay 
22e3c5b3baSBarry Smith      The script must be in your MATLAB path or current directory
23e5c89e4eSSatish Balay 
24e5c89e4eSSatish Balay      Assumes that all machines share a common file system
25e5c89e4eSSatish Balay 
26db781477SPatrick Sanan .seealso: `PetscPOpen()`, `PetscPClose()`
27e5c89e4eSSatish Balay @*/
28d71ae5a4SJacob Faibussowitsch PetscErrorCode PetscStartMatlab(MPI_Comm comm, const char machine[], const char script[], FILE **fp)
29d71ae5a4SJacob Faibussowitsch {
30e5c89e4eSSatish Balay   FILE *fd;
31e5c89e4eSSatish Balay   char  command[512];
32e5c89e4eSSatish Balay #if defined(PETSC_HAVE_UCBPS) && defined(PETSC_HAVE_POPEN)
33e5c89e4eSSatish Balay   char        buf[1024], *found;
34e5c89e4eSSatish Balay   PetscMPIInt rank;
35e5c89e4eSSatish Balay #endif
36e5c89e4eSSatish Balay 
37e5c89e4eSSatish Balay   PetscFunctionBegin;
38e5c89e4eSSatish Balay #if defined(PETSC_HAVE_UCBPS) && defined(PETSC_HAVE_POPEN)
39e3c5b3baSBarry Smith   /* check if MATLAB is not already running */
409566063dSJacob Faibussowitsch   PetscCall(PetscPOpen(comm, machine, "/usr/ucb/ps -ugxww | grep matlab | grep -v grep", "r", &fd));
419566063dSJacob Faibussowitsch   PetscCallMPI(MPI_Comm_rank(comm, &rank));
42dd400576SPatrick Sanan   if (rank == 0) found = fgets(buf, 1024, fd);
439566063dSJacob Faibussowitsch   PetscCallMPI(MPI_Bcast(&found, 1, MPI_CHAR, 0, comm));
449566063dSJacob Faibussowitsch   PetscCall(PetscPClose(comm, fd));
45*3ba16761SJacob Faibussowitsch   if (found) PetscFunctionReturn(PETSC_SUCCESS);
46e5c89e4eSSatish Balay #endif
47e5c89e4eSSatish Balay 
48e5c89e4eSSatish Balay   if (script) {
49e3c5b3baSBarry Smith     /* the remote machine won't know about current directory, so add it to MATLAB path */
50e5c89e4eSSatish Balay     /* the extra \" are to protect possible () in the script command from the shell */
51a364092eSJacob Faibussowitsch     PetscCall(PetscSNPrintf(command, PETSC_STATIC_ARRAY_LENGTH(command), "echo \"delete ${HOMEDIRECTORY}/matlab/startup.m ; path(path,'${WORKINGDIRECTORY}'); %s  \" > ${HOMEDIRECTORY}/matlab/startup.m", script));
52e5c89e4eSSatish Balay #if defined(PETSC_HAVE_POPEN)
539566063dSJacob Faibussowitsch     PetscCall(PetscPOpen(comm, machine, command, "r", &fd));
549566063dSJacob Faibussowitsch     PetscCall(PetscPClose(comm, fd));
55e5c89e4eSSatish Balay #endif
56e5c89e4eSSatish Balay   }
57e5c89e4eSSatish Balay #if defined(PETSC_HAVE_POPEN)
589566063dSJacob Faibussowitsch   PetscCall(PetscPOpen(comm, machine, "xterm -display ${DISPLAY} -e matlab -nosplash", "r", fp));
59e5c89e4eSSatish Balay #endif
60*3ba16761SJacob Faibussowitsch   PetscFunctionReturn(PETSC_SUCCESS);
61e5c89e4eSSatish Balay }
62