Lines Matching refs:A
13 static PetscErrorCode CheckValuesAIJ(Mat A) in CheckValuesAIJ() argument
20 PetscCall(MatGetSize(A, &M, &N)); in CheckValuesAIJ()
21 PetscCall(MatGetOwnershipRange(A, &rstart, &rend)); in CheckValuesAIJ()
24 PetscCall(MatGetValue(A, i, j, &val)); in CheckValuesAIJ()
35 Mat A; in main() local
45 …TSC_COMM_WORLD, bs, PETSC_DECIDE, PETSC_DECIDE, M, N, PETSC_DECIDE, NULL, PETSC_DECIDE, NULL, &A)); in main()
50 PetscCall(MatGetSize(A, &M, &N)); in main()
51 PetscCall(MatGetOwnershipRange(A, &rstart, &rend)); in main()
55 if (PetscAbsReal(v) > 0) PetscCall(MatSetValue(A, i, j, v, INSERT_VALUES)); in main()
58 PetscCall(MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY)); in main()
59 PetscCall(MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY)); in main()
60 PetscCall(MatViewFromOptions(A, NULL, "-mat_base_view")); in main()
66 for (i = 0; i < 3; i++) PetscCall(MatView(A, view)); in main()
68 PetscCall(MatDestroy(&A)); in main()
74 PetscCall(MatCreate(PETSC_COMM_WORLD, &A)); in main()
75 PetscCall(MatSetType(A, MATBAIJ)); in main()
77 if (i > 0) PetscCall(MatZeroEntries(A)); in main()
78 PetscCall(MatLoad(A, view)); in main()
79 PetscCall(CheckValuesAIJ(A)); in main()
82 PetscCall(MatViewFromOptions(A, NULL, "-mat_load_view")); in main()
83 PetscCall(MatDestroy(&A)); in main()
89 PetscCall(MatCreate(PETSC_COMM_SELF, &A)); in main()
90 PetscCall(MatSetType(A, MATSEQBAIJ)); in main()
92 if (i > 0) PetscCall(MatZeroEntries(A)); in main()
93 PetscCall(MatLoad(A, view)); in main()
94 PetscCall(CheckValuesAIJ(A)); in main()
97 PetscCall(MatDestroy(&A)); in main()
103 PetscCall(MatCreate(PETSC_COMM_WORLD, &A)); in main()
104 PetscCall(MatSetType(A, MATMPIBAIJ)); in main()
106 if (i > 0) PetscCall(MatZeroEntries(A)); in main()
107 PetscCall(MatLoad(A, view)); in main()
108 PetscCall(CheckValuesAIJ(A)); in main()
111 PetscCall(MatDestroy(&A)); in main()