Actual source code: ex31.c
1: /*$Id: ex31.c,v 1.6 2001/03/23 23:21:30 balay Exp $*/
3: /*
4: Demonstrates PetscMatlabEngineXXX()
5: */
7: #include petscvec.h
9: int main(int argc,char **argv)
10: {
11: int ierr,rank,n = 5;
12: char *output;
13: Vec x;
15: PetscInitialize(&argc,&argv,(char *)0,0);
17: VecCreate(PETSC_COMM_WORLD,PETSC_DECIDE,n,&x);
18: VecSetFromOptions(x);
20: MPI_Comm_rank(PETSC_COMM_WORLD,&rank);
21: PetscMatlabEngineGetOutput(PETSC_COMM_WORLD,&output);
22: PetscMatlabEngineEvaluate(PETSC_COMM_WORLD,"MPI_Comm_rank");
23: PetscSynchronizedPrintf(PETSC_COMM_WORLD,"[%d]Processor rank is %s",rank,output);
24: PetscSynchronizedFlush(PETSC_COMM_WORLD);
26: PetscObjectSetName((PetscObject)x,"x");
27: PetscMatlabEnginePut((PetscObject)x);
28: PetscMatlabEngineEvaluate(PETSC_COMM_WORLD,"x = x + MPI_Comm_rank;n");
29: PetscMatlabEngineGet((PetscObject)x);
31: PetscMatlabEngineEvaluate(PETSC_COMM_WORLD,"whosn");
32: PetscSynchronizedPrintf(PETSC_COMM_WORLD,"[%d]The result is %s",rank,output);
33: PetscSynchronizedFlush(PETSC_COMM_WORLD);
35: VecView(x,PETSC_VIEWER_STDOUT_WORLD);
36: VecDestroy(x);
37: PetscFinalize();
38: return 0;
39: }
40: