Actual source code: ex173.c
petsc-3.14.6 2021-03-30
1: static char help[] = "Test MatrixMarket outputing.\n\n";
3: /*
4: Include "petscmat.h" so that we can use matrices.
5: automatically includes:
6: petscsys.h - base PETSc routines petscvec.h - vectors
7: petscmat.h - matrices
8: petscis.h - index sets petscviewer.h - viewers
9: */
11: #include <petscmat.h>
13: int main(int argc,char **args)
14: {
15: Mat A;
16: PetscViewer fd; /* viewer */
17: char file[PETSC_MAX_PATH_LEN]; /* input file name */
19: PetscBool flg;
21: PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr;
22: PetscOptionsGetString(NULL,NULL,"-f0",file,sizeof(file),&flg);
23: if (!flg) SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_USER,"Must indicate binary file with the -f0 option");
24: PetscViewerBinaryOpen(PETSC_COMM_WORLD,file,FILE_MODE_READ,&fd);
25: MatCreate(PETSC_COMM_WORLD,&A);
26: MatLoad(A,fd);
27: PetscViewerDestroy(&fd);
29: PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,PETSC_VIEWER_ASCII_MATRIXMARKET);
30: MatView(A,PETSC_VIEWER_STDOUT_WORLD);
31: PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD);
32: MatDestroy(&A);
33: PetscFinalize();
34: return ierr;
35: }
37: /*TEST
39: test:
40: args: -f0 ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64
41: requires: !complex double !define(PETSC_USE_64BIT_INDICES)
43: TEST*/