Actual source code: ex65.c
1: /*$Id: ex65.c,v 1.15 2001/08/07 21:30:08 bsmith Exp $*/
3: static char help[] = "Saves a rectangular sparse matrix to disk.nn";
5: #include petscmat.h
7: int main(int argc,char **args)
8: {
9: Mat A;
10: int m = 100,n = 11,ierr,js[11],i,j,cnt;
11: PetscScalar values[11];
12: PetscViewer view;
14: PetscInitialize(&argc,&args,(char *)0,help);
16: MatCreateSeqAIJ(PETSC_COMM_WORLD,m,n,20,0,&A);
18: for (i=0; i<n; i++) values[i] = (PetscReal)i;
20: for (i=0; i<m; i++) {
21: cnt = 0;
22: if (i % 2) {
23: for (j=0; j<n; j += 2) {
24: js[cnt++] = j;
25: }
26: } else {
27: ;
28: }
29: MatSetValues(A,1,&i,cnt,js,values,INSERT_VALUES);
30: }
31: MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY);
32: MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY);
34: PetscViewerBinaryOpen(PETSC_COMM_WORLD,"rect",PETSC_BINARY_CREATE,&view);
35: MatView(A,view);
36: PetscViewerDestroy(view);
38: MatDestroy(A);
40: PetscFinalize();
41: return 0;
42: }