115229ffcSPierre Jolivet static char help[] = "Tests MatLoad() with blocksize set in program\n\n";
2c4762a1bSJed Brown
3c4762a1bSJed Brown #include <petscmat.h>
4c4762a1bSJed Brown
main(int argc,char ** args)5d71ae5a4SJacob Faibussowitsch int main(int argc, char **args)
6d71ae5a4SJacob Faibussowitsch {
7c4762a1bSJed Brown Mat A;
8c4762a1bSJed Brown PetscViewer fd;
9c4762a1bSJed Brown char file[PETSC_MAX_PATH_LEN];
10c4762a1bSJed Brown PetscBool flg;
11c4762a1bSJed Brown
12327415f7SBarry Smith PetscFunctionBeginUser;
13c8025a54SPierre Jolivet PetscCall(PetscInitialize(&argc, &args, NULL, help));
14c4762a1bSJed Brown /* Determine files from which we read the matrix */
159566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetString(NULL, NULL, "-f", file, sizeof(file), &flg));
1628b400f6SJacob Faibussowitsch PetscCheck(flg, PETSC_COMM_WORLD, PETSC_ERR_USER, "Must indicate binary file with the -f");
17c4762a1bSJed Brown
18c4762a1bSJed Brown /* Load matrices */
199566063dSJacob Faibussowitsch PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, file, FILE_MODE_READ, &fd));
209566063dSJacob Faibussowitsch PetscCall(MatCreate(PETSC_COMM_WORLD, &A));
219566063dSJacob Faibussowitsch PetscCall(MatSetType(A, MATSBAIJ));
229566063dSJacob Faibussowitsch PetscCall(MatSetFromOptions(A));
239566063dSJacob Faibussowitsch PetscCall(MatSetBlockSize(A, 2));
249566063dSJacob Faibussowitsch PetscCall(MatLoad(A, fd));
259566063dSJacob Faibussowitsch PetscCall(PetscViewerDestroy(&fd));
269566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A));
279566063dSJacob Faibussowitsch PetscCall(PetscFinalize());
28b122ec5aSJacob Faibussowitsch return 0;
29c4762a1bSJed Brown }
30c4762a1bSJed Brown
31c4762a1bSJed Brown /*TEST
32c4762a1bSJed Brown
33c4762a1bSJed Brown test:
34c4762a1bSJed Brown args: -mat_type aij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
35*3886731fSPierre Jolivet output_file: output/empty.out
36dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
37c4762a1bSJed Brown
38c4762a1bSJed Brown test:
39c4762a1bSJed Brown suffix: 2
40c4762a1bSJed Brown nsize: 2
41c4762a1bSJed Brown args: -mat_type aij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
42*3886731fSPierre Jolivet output_file: output/empty.out
43dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
44c4762a1bSJed Brown
45c4762a1bSJed Brown test:
46c4762a1bSJed Brown suffix: 3
47c4762a1bSJed Brown args: -mat_type baij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
48*3886731fSPierre Jolivet output_file: output/empty.out
49dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
50c4762a1bSJed Brown
51c4762a1bSJed Brown test:
52c4762a1bSJed Brown suffix: 4more
53c4762a1bSJed Brown nsize: 2
54c4762a1bSJed Brown args: -mat_type baij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
55*3886731fSPierre Jolivet output_file: output/empty.out
56dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
57c4762a1bSJed Brown
58c4762a1bSJed Brown test:
59c4762a1bSJed Brown suffix: 5
60c4762a1bSJed Brown args: -mat_type sbaij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
61*3886731fSPierre Jolivet output_file: output/empty.out
62dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
63c4762a1bSJed Brown
64c4762a1bSJed Brown test:
65c4762a1bSJed Brown suffix: 6
66c4762a1bSJed Brown nsize: 2
67c4762a1bSJed Brown args: -mat_type sbaij -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
68*3886731fSPierre Jolivet output_file: output/empty.out
69dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
70c4762a1bSJed Brown
71c4762a1bSJed Brown test:
72c4762a1bSJed Brown suffix: 7
73c4762a1bSJed Brown args: -mat_type sbaij -matload_block_size 4 -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
74*3886731fSPierre Jolivet output_file: output/empty.out
75dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
76c4762a1bSJed Brown
77c4762a1bSJed Brown test:
78c4762a1bSJed Brown suffix: 8
79c4762a1bSJed Brown nsize: 2
80c4762a1bSJed Brown args: -mat_type sbaij -matload_block_size 4 -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
81*3886731fSPierre Jolivet output_file: output/empty.out
82dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
83c4762a1bSJed Brown
84c4762a1bSJed Brown test:
85c4762a1bSJed Brown suffix: 9
86c4762a1bSJed Brown args: -mat_type baij -matload_block_size 4 -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
87*3886731fSPierre Jolivet output_file: output/empty.out
88dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
89c4762a1bSJed Brown
90c4762a1bSJed Brown test:
91c4762a1bSJed Brown suffix: 10
92c4762a1bSJed Brown nsize: 2
93c4762a1bSJed Brown args: -mat_type baij -matload_block_size 4 -f ${wPETSC_DIR}/share/petsc/datafiles/matrices/ns-real-int32-float64 -malloc_dump
94*3886731fSPierre Jolivet output_file: output/empty.out
95dfd57a17SPierre Jolivet requires: !complex double !defined(PETSC_USE_64BIT_INDICES)
96c4762a1bSJed Brown
97c4762a1bSJed Brown TEST*/
98