xref: /petsc/src/sys/tests/ex16.c (revision 9566063d113dddea24716c546802770db7481bc0)
1c4762a1bSJed Brown 
2c4762a1bSJed Brown static char help[] = "Demonstrates PetscMatlabEngineXXX()\n";
3c4762a1bSJed Brown 
4c4762a1bSJed Brown #include <petscsys.h>
5c4762a1bSJed Brown #include <petscmatlab.h>
6c4762a1bSJed Brown 
7c4762a1bSJed Brown int main(int argc,char **argv)
8c4762a1bSJed Brown {
9c4762a1bSJed Brown   PetscErrorCode ierr;
10c4762a1bSJed Brown   PetscMPIInt    rank;
11c4762a1bSJed Brown   char           buffer[256],*output,user[256];
12c4762a1bSJed Brown   PetscBool      userhappy = PETSC_FALSE;
13c4762a1bSJed Brown 
14*9566063dSJacob Faibussowitsch   PetscCall(PetscInitialize(&argc,&argv,(char*)0,help));
15c4762a1bSJed Brown   ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);
16c4762a1bSJed Brown 
17*9566063dSJacob Faibussowitsch   PetscCall(PetscMatlabEngineGetOutput(PETSC_MATLAB_ENGINE_(PETSC_COMM_WORLD),&output));
18c4762a1bSJed Brown 
19*9566063dSJacob Faibussowitsch   PetscCall(PetscMatlabEngineEvaluate(PETSC_MATLAB_ENGINE_(PETSC_COMM_WORLD),"MPI_Comm_rank"));
20*9566063dSJacob Faibussowitsch   PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD,"[%d]Processor rank is %s",rank,output));
21*9566063dSJacob Faibussowitsch   PetscCall(PetscSynchronizedFlush(PETSC_COMM_WORLD,PETSC_STDOUT));
22*9566063dSJacob Faibussowitsch   PetscCall(PetscPrintf(PETSC_COMM_WORLD,">>"));
23*9566063dSJacob Faibussowitsch   PetscCall(PetscSynchronizedFGets(PETSC_COMM_WORLD,stdin,256,user));
24*9566063dSJacob Faibussowitsch   PetscCall(PetscStrncmp(user,"exit",4,&userhappy));
25c4762a1bSJed Brown   while (!userhappy) {
26*9566063dSJacob Faibussowitsch     PetscCall(PetscMatlabEngineEvaluate(PETSC_MATLAB_ENGINE_(PETSC_COMM_WORLD),user));
27*9566063dSJacob Faibussowitsch     PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD,"[%d]The result is %s",rank,output));
28*9566063dSJacob Faibussowitsch     PetscCall(PetscSynchronizedFlush(PETSC_COMM_WORLD,PETSC_STDOUT));
29*9566063dSJacob Faibussowitsch     PetscCall(PetscPrintf(PETSC_COMM_WORLD,">>"));
30*9566063dSJacob Faibussowitsch     PetscCall(PetscSynchronizedFGets(PETSC_COMM_WORLD,stdin,256,user));
31*9566063dSJacob Faibussowitsch     PetscCall(PetscStrncmp(user,"exit",4,&userhappy));
32c4762a1bSJed Brown   }
33*9566063dSJacob Faibussowitsch   PetscCall(PetscFinalize());
34b122ec5aSJacob Faibussowitsch   return 0;
35c4762a1bSJed Brown }
36