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