Lines Matching refs:p
12 PetscInt m, n, M, N, p = 1, q = 1, i, j; in main() local
30 PetscCall(PetscOptionsGetInt(NULL, NULL, "-p", &p, NULL)); in main()
32 PetscCall(PetscMalloc2(p * q, &S, p * q, &T)); in main()
33 for (i = 0; i < p * q; i++) S[i] = 0; in main()
35 for (i = 0; i < p; i++) { in main()
38 S[i + p * j] = ((PetscReal)((i + 1) * (j + 1))) / ((PetscReal)(p + q)); in main()
39 T[i + p * j] = ((PetscReal)((p - i) + j)) / ((PetscReal)(p * q)); in main()
46 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main()
68 PetscCall(MatCreateKAIJ(A, p, q, NULL, T, &TA)); in main()
90 PetscCall(MatCreateKAIJ(A, p, q, S, NULL, &TA)); in main()
111 if (p == q) { in main()
112 for (i = 0; i < p; i++) { in main()
114 if (i == j) T[i + j * p] = 1.0; in main()
115 else T[i + j * p] = 0.0; in main()
120 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main()
139 PetscCall(MatCreateKAIJ(A, p, q, NULL, T, &TA)); in main()
145 for (i = 0; i < p; i++) { in main()
147 if (i == j) S[i + j * p] = T[i + j * p] = 2.0; in main()
148 else S[i + j * p] = T[i + j * p] = 0.0; in main()
151 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main()