1a5b23f4aSJose E. Roman static char help[] = "Test MatrixMarket outputting.\n\n";
2c4762a1bSJed Brown
3c4762a1bSJed Brown /*
4c4762a1bSJed Brown Include "petscmat.h" so that we can use matrices.
5c4762a1bSJed Brown automatically includes:
6c4762a1bSJed Brown petscsys.h - base PETSc routines petscvec.h - vectors
7c4762a1bSJed Brown petscmat.h - matrices
8c4762a1bSJed Brown petscis.h - index sets petscviewer.h - viewers
9c4762a1bSJed Brown */
10c4762a1bSJed Brown
11c4762a1bSJed Brown #include <petscmat.h>
12c4762a1bSJed Brown
main(int argc,char ** args)13d71ae5a4SJacob Faibussowitsch int main(int argc, char **args)
14d71ae5a4SJacob Faibussowitsch {
15c4762a1bSJed Brown Mat A;
16c4762a1bSJed Brown PetscViewer fd; /* viewer */
17c4762a1bSJed Brown char file[PETSC_MAX_PATH_LEN]; /* input file name */
18c4762a1bSJed Brown PetscBool flg;
19c4762a1bSJed Brown
20327415f7SBarry Smith PetscFunctionBeginUser;
21*c8025a54SPierre Jolivet PetscCall(PetscInitialize(&argc, &args, NULL, help));
229566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetString(NULL, NULL, "-f0", file, sizeof(file), &flg));
2328b400f6SJacob Faibussowitsch PetscCheck(flg, PETSC_COMM_WORLD, PETSC_ERR_USER, "Must indicate binary file with the -f0 option");
249566063dSJacob Faibussowitsch PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, file, FILE_MODE_READ, &fd));
259566063dSJacob Faibussowitsch PetscCall(MatCreate(PETSC_COMM_WORLD, &A));
269566063dSJacob Faibussowitsch PetscCall(MatLoad(A, fd));
279566063dSJacob Faibussowitsch PetscCall(PetscViewerDestroy(&fd));
28c4762a1bSJed Brown
299566063dSJacob Faibussowitsch PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD, PETSC_VIEWER_ASCII_MATRIXMARKET));
309566063dSJacob Faibussowitsch PetscCall(MatView(A, PETSC_VIEWER_STDOUT_WORLD));
319566063dSJacob Faibussowitsch PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD));
329566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A));
339566063dSJacob Faibussowitsch PetscCall(PetscFinalize());
34b122ec5aSJacob Faibussowitsch return 0;
35c4762a1bSJed Brown }
36c4762a1bSJed Brown
37c4762a1bSJed Brown /*TEST
38c4762a1bSJed Brown
39c4762a1bSJed Brown test:
40c4762a1bSJed Brown args: -f0 ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64
41dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
42c4762a1bSJed Brown
43c4762a1bSJed Brown TEST*/
44