10   Mat         mat, RHS, SOLU;  in createMatsAndVecs()  local
40   PetscCall(MatCreate(PETSC_COMM_WORLD, &mat));  in createMatsAndVecs()
41   PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n));  in createMatsAndVecs()
42   PetscCall(MatSetType(mat, MATDENSE));  in createMatsAndVecs()
43   PetscCall(MatSetFromOptions(mat));  in createMatsAndVecs()
44   PetscCall(MatSetUp(mat));  in createMatsAndVecs()
45   PetscCall(MatGetOwnershipRange(mat, &rstart, &rend));  in createMatsAndVecs()
46   PetscCall(MatGetOwnershipRangeColumn(mat, &cstart, &cend));  in createMatsAndVecs()
51         PetscCall(MatSetValues(mat, 1, &i, 1, &i, &value, INSERT_VALUES));  in createMatsAndVecs()
55           PetscCall(MatSetValues(mat, 1, &i, 1, &j, &value, INSERT_VALUES));  in createMatsAndVecs()
59     PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY));  in createMatsAndVecs()
60     PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY));  in createMatsAndVecs()
62     PetscCall(MatSetRandom(mat, rand));  in createMatsAndVecs()
66       PetscCall(MatMatTransposeMult(mat, mat, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &T));  in createMatsAndVecs()
67       PetscCall(MatDestroy(&mat));  in createMatsAndVecs()
68       mat = T;  in createMatsAndVecs()
73   PetscCall(MatCreateVecs(mat, &x, &b));  in createMatsAndVecs()
77   *_mat  = mat;  in createMatsAndVecs()
88   Mat         mat, F, RHS, SOLU;  in main()  local
111   PetscCall(createMatsAndVecs(n, n, nrhs, full, &mat, &RHS, &SOLU, &x, &y, &b));  in main()
115   PetscCall(PetscObjectBaseTypeCompare((PetscObject)mat, MATSEQDENSE, &inplace));  in main()
116   PetscCall(MatGetLocalSize(mat, &i, NULL));  in main()
117   PetscCall(MatGetOwnershipRange(mat, &j, NULL));  in main()
120   PetscCall(MatGetInfo(mat, MAT_LOCAL, &info));  in main()
122   PetscCall(MatMult(mat, x, b));  in main()
129     if (!ldl) PetscCall(MatSetOption(mat, MAT_SPD, PETSC_TRUE));  in main()
130     PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &F));  in main()
138     PetscCall(MatMatMult(mat, SOLU, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &RHS2));  in main()
147   if (!ldl) PetscCall(MatSetOption(mat, MAT_SPD, PETSC_TRUE));  in main()
148   PetscCall(MatGetFactor(mat, solver, MAT_FACTOR_CHOLESKY, &F));  in main()
149   PetscCall(MatCholeskyFactorSymbolic(F, mat, perm, 0));  in main()
150   PetscCall(MatCholeskyFactorNumeric(F, mat, 0));  in main()
159   PetscCall(MatZeroRows(mat, 1, &i, -1.0, NULL, NULL));  in main()
160   PetscCall(MatMult(mat, x, b));  in main()
166     PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &F));  in main()
173     PetscCall(MatMatMult(mat, SOLU, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &RHS2));  in main()
182   PetscCall(MatGetFactor(mat, solver, MAT_FACTOR_LU, &F));  in main()
183   PetscCall(MatLUFactorSymbolic(F, mat, perm, perm, 0));  in main()
184   PetscCall(MatLUFactorNumeric(F, mat, 0));  in main()
193   PetscCall(MatDestroy(&mat));  in main()
203     PetscCall(createMatsAndVecs(m, n, nrhs, full, &mat, &RHS, &SOLU, &x, &y, &b));  in main()
207     PetscCall(MatMult(mat, x, b));  in main()
213       PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &F));  in main()
219       PetscCall(MatMatMult(mat, SOLU, MAT_REUSE_MATRIX, PETSC_DETERMINE, &RHS));  in main()
230     PetscCall(MatGetFactor(mat, solver, MAT_FACTOR_QR, &F));  in main()
231     PetscCall(MatQRFactorSymbolic(F, mat, NULL, NULL));  in main()
232     PetscCall(MatQRFactorNumeric(F, mat, NULL));  in main()
240       PetscCall(MatMultTranspose(mat, x, b));  in main()
249     PetscCall(MatDestroy(&mat));  in main()