Lines Matching refs:k
7 PetscInt i, rstart, rend, N = 10, num_numfac = 5, col[3], k; in main() local
25 for (k = 0; k < num_numfac; k++) { in main()
26 PetscCall(MatCreate(PETSC_COMM_WORLD, &A[k])); in main()
27 PetscCall(MatSetSizes(A[k], PETSC_DECIDE, PETSC_DECIDE, N, N)); in main()
28 PetscCall(MatSetFromOptions(A[k])); in main()
29 PetscCall(MatSetUp(A[k])); in main()
30 PetscCall(MatGetOwnershipRange(A[k], &rstart, &rend)); in main()
32 value[0] = -1.0 * (k + 1); in main()
33 value[1] = 2.0 * (k + 1); in main()
34 value[2] = -1.0 * (k + 1); in main()
40 PetscCall(MatSetValues(A[k], 1, &i, 2, col + 1, value + 1, INSERT_VALUES)); in main()
42 PetscCall(MatSetValues(A[k], 1, &i, 2, col, value, INSERT_VALUES)); in main()
44 PetscCall(MatSetValues(A[k], 1, &i, 3, col, value, INSERT_VALUES)); in main()
47 PetscCall(MatAssemblyBegin(A[k], MAT_FINAL_ASSEMBLY)); in main()
48 PetscCall(MatAssemblyEnd(A[k], MAT_FINAL_ASSEMBLY)); in main()
49 PetscCall(MatSetOption(A[k], MAT_NEW_NONZERO_LOCATION_ERR, PETSC_TRUE)); in main()
93 for (k = 0; k < num_numfac; k++) { in main()
97 PetscCall(MatLUFactorNumeric(F, A[k], &info)); in main()
101 PetscCall(MatCholeskyFactorNumeric(F, A[k], &info)); in main()
111 PetscCall(MatMult(A[k], x, u)); in main()
114 …SC_COMM_WORLD, "%" PetscInt_FMT "-the %s numfact and solve: residual %g\n", k, factortype, (double… in main()
118 for (k = 0; k < num_numfac; k++) PetscCall(MatDestroy(&A[k])); in main()