Lines Matching refs:PetscCall
11 PetscCall(MatGetType(mat, &mtype)); in TransposeAXPY()
13 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\nMatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZ… in TransposeAXPY()
15 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\nMatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZ… in TransposeAXPY()
17 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
18 PetscCall(f(C, &D)); in TransposeAXPY()
19 PetscCall(f(D, &E)); in TransposeAXPY()
20 PetscCall(MatAXPY(E, alpha, mat, SAME_NONZERO_PATTERN)); in TransposeAXPY()
21 PetscCall(MatConvert(E, mtype, MAT_INPLACE_MATRIX, &E)); in TransposeAXPY()
22 PetscCall(MatView(E, PETSC_VIEWER_STDOUT_WORLD)); in TransposeAXPY()
23 PetscCall(MatDestroy(&E)); in TransposeAXPY()
24 PetscCall(MatDestroy(&D)); in TransposeAXPY()
25 PetscCall(MatDestroy(&C)); in TransposeAXPY()
27 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATT… in TransposeAXPY()
29 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATT… in TransposeAXPY()
32 PetscCall(MatTranspose(mat, MAT_INITIAL_MATRIX, &D)); in TransposeAXPY()
34 PetscCall(MatHermitianTranspose(mat, MAT_INITIAL_MATRIX, &D)); in TransposeAXPY()
36 PetscCall(f(D, &E)); in TransposeAXPY()
37 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
38 PetscCall(MatAXPY(C, alpha, E, SAME_NONZERO_PATTERN)); in TransposeAXPY()
39 PetscCall(MatView(C, PETSC_VIEWER_STDOUT_WORLD)); in TransposeAXPY()
40 PetscCall(MatDestroy(&E)); in TransposeAXPY()
41 PetscCall(MatDestroy(&D)); in TransposeAXPY()
42 PetscCall(MatDestroy(&C)); in TransposeAXPY()
44 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_… in TransposeAXPY()
46 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_… in TransposeAXPY()
48 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
49 PetscCall(f(C, &D)); in TransposeAXPY()
50 PetscCall(f(D, &E)); in TransposeAXPY()
51 PetscCall(f(mat, &F)); in TransposeAXPY()
52 PetscCall(f(F, &G)); in TransposeAXPY()
53 PetscCall(MatAXPY(E, alpha, G, SAME_NONZERO_PATTERN)); in TransposeAXPY()
54 PetscCall(MatConvert(E, mtype, MAT_INPLACE_MATRIX, &E)); in TransposeAXPY()
55 PetscCall(MatView(E, PETSC_VIEWER_STDOUT_WORLD)); in TransposeAXPY()
56 PetscCall(MatDestroy(&G)); in TransposeAXPY()
57 PetscCall(MatDestroy(&F)); in TransposeAXPY()
58 PetscCall(MatDestroy(&E)); in TransposeAXPY()
59 PetscCall(MatDestroy(&D)); in TransposeAXPY()
60 PetscCall(MatDestroy(&C)); in TransposeAXPY()
63 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: Now without the transposition operation\n")); in TransposeAXPY()
64 PetscCall(MatConvert(mat, MATSHELL, MAT_INITIAL_MATRIX, &C)); in TransposeAXPY()
65 PetscCall(f(C, &D)); in TransposeAXPY()
66 PetscCall(f(D, &E)); in TransposeAXPY()
68 PetscCall(MatConvert(E, mtype, MAT_INITIAL_MATRIX, &F)); in TransposeAXPY()
69 PetscCall(MatAXPY(F, alpha, mat, SAME_NONZERO_PATTERN)); in TransposeAXPY()
70 PetscCall(MatView(F, PETSC_VIEWER_STDOUT_WORLD)); in TransposeAXPY()
71 PetscCall(MatDestroy(&F)); in TransposeAXPY()
72 PetscCall(MatDestroy(&E)); in TransposeAXPY()
73 PetscCall(MatDestroy(&D)); in TransposeAXPY()
74 PetscCall(MatDestroy(&C)); in TransposeAXPY()
88 PetscCall(PetscInitialize(&argc, &argv, NULL, help)); in main()
89 PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD, PETSC_VIEWER_ASCII_COMMON)); in main()
90 PetscCall(PetscOptionsGetInt(NULL, NULL, "-m", &m, NULL)); in main()
94 PetscCall(PetscOptionsHasName(NULL, NULL, "-rectA", &flg)); in main()
99 PetscCall(PetscOptionsHasName(NULL, NULL, "-rectB", &flg)); in main()
106 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
107 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n)); in main()
108 PetscCall(MatSetFromOptions(mat)); in main()
109 PetscCall(MatSetUp(mat)); in main()
110 PetscCall(MatGetOwnershipRange(mat, &rstart, &rend)); in main()
114 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
117 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
118 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in main()
121 PetscCall(MatNorm(mat, NORM_FROBENIUS, &normf)); in main()
122 PetscCall(MatNorm(mat, NORM_1, &norm1)); in main()
123 PetscCall(MatNorm(mat, NORM_INFINITY, &normi)); in main()
124 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "original A: Frobenius norm = %g, one norm = %g, infinity … in main()
125 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
128 PetscCall(PetscOptionsHasName(NULL, NULL, "-in_place", &flg)); in main()
130 PetscCall(MatTranspose(mat, MAT_REUSE_MATRIX, &mat)); /* in-place transpose */ in main()
134 PetscCall(MatTranspose(mat, MAT_INITIAL_MATRIX, &tmat)); in main()
139 PetscCall(MatNorm(tmat, NORM_FROBENIUS, &normf)); in main()
140 PetscCall(MatNorm(tmat, NORM_1, &norm1)); in main()
141 PetscCall(MatNorm(tmat, NORM_INFINITY, &normi)); in main()
142 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "B = A^T: Frobenius norm = %g, one norm = %g, infinity nor… in main()
143 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
148 PetscCall(PetscOptionsGetScalar(NULL, NULL, "-alpha", &alpha, NULL)); in main()
149 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: B = B + alpha * A\n")); in main()
150 PetscCall(MatAXPY(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
151 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
153 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAYPX: B = alpha*B + A\n")); in main()
154 PetscCall(MatAYPX(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
155 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
161 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN\n"… in main()
162 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in main()
163 PetscCall(MatAXPY(C, alpha, mat, SAME_NONZERO_PATTERN)); in main()
164 PetscCall(MatView(C, PETSC_VIEWER_STDOUT_WORLD)); in main()
165 PetscCall(MatDestroy(&C)); in main()
166 PetscCall(TransposeAXPY(C, alpha, mat, MatCreateTranspose)); in main()
167 PetscCall(TransposeAXPY(C, alpha, mat, MatCreateHermitianTranspose)); in main()
173 PetscCall(MatCreate(PETSC_COMM_WORLD, &matB)); in main()
174 PetscCall(MatSetSizes(matB, PETSC_DECIDE, PETSC_DECIDE, m, n)); in main()
175 PetscCall(MatSetFromOptions(matB)); in main()
176 PetscCall(MatSetUp(matB)); in main()
177 PetscCall(MatGetOwnershipRange(matB, &rstart, &rend)); in main()
182 PetscCall(MatSetValues(matB, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
185 PetscCall(MatAssemblyBegin(matB, MAT_FINAL_ASSEMBLY)); in main()
186 PetscCall(MatAssemblyEnd(matB, MAT_FINAL_ASSEMBLY)); in main()
187 PetscCall(PetscPrintf(PETSC_COMM_WORLD, " A: original matrix:\n")); in main()
188 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
189 PetscCall(PetscPrintf(PETSC_COMM_WORLD, " B(a subset of A):\n")); in main()
190 PetscCall(MatView(matB, PETSC_VIEWER_STDOUT_WORLD)); in main()
191 … PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN\n")); in main()
192 PetscCall(MatAXPY(mat, alpha, matB, SUBSET_NONZERO_PATTERN)); in main()
193 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
194 PetscCall(MatDestroy(&matB)); in main()
200 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatZeroRows:\n")); in main()
201 PetscCall(MatZeroRows(mat, 1, &j, 0.0, NULL, NULL)); in main()
202 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
205 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "MatShift: B = B - 2*I\n")); in main()
206 PetscCall(MatShift(mat, -2.0)); in main()
207 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
209 PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD)); in main()
211 PetscCall(MatDestroy(&mat)); in main()
212 PetscCall(MatDestroy(&tmat)); in main()
213 PetscCall(PetscFinalize()); in main()