Actual source code: ex25.c
1: /*$Id: ex25.c,v 1.17 2001/08/07 03:03:07 balay Exp $*/
3: static char help[] = "Tests MatTranspose()\n\n";
5: #include petscmat.h
9: int main(int argc,char **args)
10: {
11: Mat C,A;
12: PetscScalar v;
13: int i,j,m = 4,n = 4,rank,size,I,J,ierr,Istart,Iend;
15: PetscInitialize(&argc,&args,(char *)0,help);
16: MPI_Comm_rank(PETSC_COMM_WORLD,&rank);
17: MPI_Comm_size(PETSC_COMM_WORLD,&size);
18:
19: MatCreateMPIAIJ(PETSC_COMM_WORLD,PETSC_DECIDE,PETSC_DECIDE,
20: m*n,m*n,5,PETSC_NULL,5,PETSC_NULL,&C);
22: /* create the matrix for the five point stencil, YET AGAIN*/
23: MatGetOwnershipRange(C,&Istart,&Iend);
24: for (I=Istart; I<Iend; I++) {
25: v = -1.0; i = I/n; j = I - i*n;
26: if (i>0) {J = I - n; MatSetValues(C,1,&I,1,&J,&v,INSERT_VALUES);}
27: if (i<m-1) {J = I + n; MatSetValues(C,1,&I,1,&J,&v,INSERT_VALUES);}
28: if (j>0) {J = I - 1; MatSetValues(C,1,&I,1,&J,&v,INSERT_VALUES);}
29: if (j<n-1) {J = I + 1; MatSetValues(C,1,&I,1,&J,&v,INSERT_VALUES);}
30: v = 4.0; MatSetValues(C,1,&I,1,&I,&v,INSERT_VALUES);
31: }
32: MatAssemblyBegin(C,MAT_FINAL_ASSEMBLY);
33: MatAssemblyEnd(C,MAT_FINAL_ASSEMBLY);
34: MatView(C,PETSC_VIEWER_STDOUT_WORLD);
36: MatTranspose(C,&A);
37: MatView(A,PETSC_VIEWER_STDOUT_WORLD);
39: MatDestroy(C);
40: MatDestroy(A);
41: PetscFinalize();
42: return 0;
43: }