xref: /petsc/src/mat/tests/ex136.c (revision c4762a1b19cd2af06abeed90e8f9d34fb975dd94) 
1*c4762a1bSJed Brown 
2*c4762a1bSJed Brown static char help[] = "Tests MatLoad() MatView() for MPIBAIJ.\n\n";
3*c4762a1bSJed Brown 
4*c4762a1bSJed Brown #include <petscmat.h>
5*c4762a1bSJed Brown 
6*c4762a1bSJed Brown int main(int argc,char **args)
7*c4762a1bSJed Brown {
8*c4762a1bSJed Brown   Mat            A,B;
9*c4762a1bSJed Brown   PetscErrorCode ierr;
10*c4762a1bSJed Brown   char           file[PETSC_MAX_PATH_LEN];
11*c4762a1bSJed Brown   PetscBool      flg;
12*c4762a1bSJed Brown   PetscViewer    fd;
13*c4762a1bSJed Brown 
14*c4762a1bSJed Brown   ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr;
15*c4762a1bSJed Brown   ierr = PetscOptionsGetString(NULL,NULL,"-f",file,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
16*c4762a1bSJed Brown   if (!flg) SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_USER,"Must indicate binary file with the -f option");
17*c4762a1bSJed Brown 
18*c4762a1bSJed Brown   /*
19*c4762a1bSJed Brown      Open binary file.  Note that we use FILE_MODE_READ to indicate
20*c4762a1bSJed Brown      reading from this file.
21*c4762a1bSJed Brown   */
22*c4762a1bSJed Brown   ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,file,FILE_MODE_READ,&fd);CHKERRQ(ierr);
23*c4762a1bSJed Brown 
24*c4762a1bSJed Brown   /*
25*c4762a1bSJed Brown      Load the matrix; then destroy the viewer.
26*c4762a1bSJed Brown   */
27*c4762a1bSJed Brown   ierr = MatCreate(PETSC_COMM_WORLD,&A);CHKERRQ(ierr);
28*c4762a1bSJed Brown   ierr = MatSetFromOptions(A);CHKERRQ(ierr);
29*c4762a1bSJed Brown   ierr = MatLoad(A,fd);CHKERRQ(ierr);
30*c4762a1bSJed Brown   ierr = PetscViewerDestroy(&fd);CHKERRQ(ierr);
31*c4762a1bSJed Brown 
32*c4762a1bSJed Brown   /*
33*c4762a1bSJed Brown      Open another binary file.  Note that we use FILE_MODE_WRITE to indicate writing to the file
34*c4762a1bSJed Brown   */
35*c4762a1bSJed Brown   ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,"fileoutput",FILE_MODE_WRITE,&fd);CHKERRQ(ierr);
36*c4762a1bSJed Brown   ierr = PetscViewerBinarySetFlowControl(fd,3);CHKERRQ(ierr);
37*c4762a1bSJed Brown   /*
38*c4762a1bSJed Brown      Save the matrix and vector; then destroy the viewer.
39*c4762a1bSJed Brown   */
40*c4762a1bSJed Brown   ierr = MatView(A,fd);CHKERRQ(ierr);
41*c4762a1bSJed Brown   ierr = PetscViewerDestroy(&fd);CHKERRQ(ierr);
42*c4762a1bSJed Brown 
43*c4762a1bSJed Brown   /* load the new matrix */
44*c4762a1bSJed Brown   ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,"fileoutput",FILE_MODE_READ,&fd);CHKERRQ(ierr);
45*c4762a1bSJed Brown   ierr = MatCreate(PETSC_COMM_WORLD,&B);CHKERRQ(ierr);
46*c4762a1bSJed Brown   ierr = MatSetFromOptions(B);CHKERRQ(ierr);
47*c4762a1bSJed Brown   ierr = MatLoad(B,fd);CHKERRQ(ierr);
48*c4762a1bSJed Brown   ierr = PetscViewerDestroy(&fd);CHKERRQ(ierr);
49*c4762a1bSJed Brown 
50*c4762a1bSJed Brown   ierr = MatEqual(A,B,&flg);CHKERRQ(ierr);
51*c4762a1bSJed Brown   if (flg) {
52*c4762a1bSJed Brown     ierr = PetscPrintf(PETSC_COMM_WORLD,"Matrices are equal\n");CHKERRQ(ierr);
53*c4762a1bSJed Brown   } else {
54*c4762a1bSJed Brown     ierr = PetscPrintf(PETSC_COMM_WORLD,"Matrices are not equal\n");CHKERRQ(ierr);
55*c4762a1bSJed Brown   }
56*c4762a1bSJed Brown 
57*c4762a1bSJed Brown   ierr = MatDestroy(&A);CHKERRQ(ierr);
58*c4762a1bSJed Brown   ierr = MatDestroy(&B);CHKERRQ(ierr);
59*c4762a1bSJed Brown   ierr = PetscFinalize();
60*c4762a1bSJed Brown   return ierr;
61*c4762a1bSJed Brown }
62*c4762a1bSJed Brown 
63*c4762a1bSJed Brown 
64*c4762a1bSJed Brown /*TEST
65*c4762a1bSJed Brown 
66*c4762a1bSJed Brown    test:
67*c4762a1bSJed Brown       nsize: 3
68*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
69*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info
70*c4762a1bSJed Brown 
71*c4762a1bSJed Brown    test:
72*c4762a1bSJed Brown       suffix: 2
73*c4762a1bSJed Brown       nsize: 5
74*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
75*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info
76*c4762a1bSJed Brown 
77*c4762a1bSJed Brown    test:
78*c4762a1bSJed Brown       suffix: 3
79*c4762a1bSJed Brown       nsize: 7
80*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
81*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info
82*c4762a1bSJed Brown 
83*c4762a1bSJed Brown    test:
84*c4762a1bSJed Brown       suffix: 4
85*c4762a1bSJed Brown       nsize: 3
86*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
87*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info -mat_type baij
88*c4762a1bSJed Brown 
89*c4762a1bSJed Brown    test:
90*c4762a1bSJed Brown       suffix: 5
91*c4762a1bSJed Brown       nsize: 5
92*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
93*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info -mat_type baij
94*c4762a1bSJed Brown 
95*c4762a1bSJed Brown    test:
96*c4762a1bSJed Brown       suffix: 6
97*c4762a1bSJed Brown       nsize: 7
98*c4762a1bSJed Brown       requires: datafilespath !complex double !define(PETSC_USE_64BIT_INDICES)
99*c4762a1bSJed Brown       args: -f ${DATAFILESPATH}/matrices/cfd.2.100 -mat_view ascii::ascii_info -mat_type baij
100*c4762a1bSJed Brown 
101*c4762a1bSJed Brown TEST*/
102