| /petsc/src/snes/impls/vi/ss/ |
| H A D | viss.c | 55 SNES_VINEWTONSSLS *vi = (SNES_VINEWTONSSLS *)snes->data; in SNESVIComputeFunction() local
62 PetscCall((*vi->computeuserfunction)(snes, X, F, functx)); in SNESVIComputeFunction()
197 SNES_VINEWTONSSLS *vi = (SNES_VINEWTONSSLS *)snes->data; in SNESSolve_VINEWTONSSLS() local
210 vi->computeuserfunction = sdm->ops->computefunction; in SNESSolve_VINEWTONSSLS()
228 PetscCall(SNESComputeFunction(snes, X, vi->phi)); in SNESSolve_VINEWTONSSLS()
232 sdm->ops->computefunction = vi->computeuserfunction; in SNESSolve_VINEWTONSSLS()
236 PetscCall(SNESVIComputeMeritFunction(vi->phi, &vi->merit, &vi->phinorm)); in SNESSolve_VINEWTONSSLS()
240 SNESCheckFunctionDomainError(snes, vi->merit); in SNESSolve_VINEWTONSSLS()
243 snes->norm = vi->phinorm; in SNESSolve_VINEWTONSSLS()
245 PetscCall(SNESLogConvergenceHistory(snes, vi->phinorm, 0)); in SNESSolve_VINEWTONSSLS()
[all …]
|
| /petsc/src/mat/tutorials/ |
| H A D | ex18.c | 60 for (PetscInt vi = 0; vi < 3; vi++) { in CreateMatrix() local
62 oor[cnt] = fe->vertices[3 * e + vi]; in CreateMatrix()
85 for (PetscInt vi = 0; vi < 3; vi++) { in FillMatrixCPU() local
86 for (PetscInt vj = 0; vj < 3; vj++) s[3 * vi + vj] = vi + 2 * vj; in FillMatrixCPU()
108 for (PetscInt vi = 0; vi < 3; vi++) { in FillMatrixCPUCOO() local
109 for (PetscInt vj = 0; vj < 3; vj++) s[3 * vi + vj] = vi + 2 * vj; in FillMatrixCPUCOO()
129 for (PetscInt vi = 0; vi < 3; vi++) { in FillMatrixCPUCOO3d() local
130 for (PetscInt vj = 0; vj < 3; vj++) s[e][vi][vj] = vi + 2 * vj; in FillMatrixCPUCOO3d()
|
| H A D | ex18k.kokkos.cxx | 16 for (PetscInt vi = 0; vi < 3; vi++) { in FillMatrixKokkosCOO()
17 for (PetscInt vj = 0; vj < 3; vj++) s[vi * 3 + vj] = vi + 2 * vj; in FillMatrixKokkosCOO()
|
| H A D | ex18cu.cu | 10 for (PetscInt vi = 0; vi < 3; vi++) { in FillValues() local
11 for (PetscInt vj = 0; vj < 3; vj++) s[vi * 3 + vj] = vi + 2 * vj; in FillValues()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvnat1.c | 11 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() local
27 vi = aj + ai[i]; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
32 jdx = *vi++; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
42 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
47 idx = *vi++; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
64 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *vi; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering() local
79 vi = aj; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
83 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
85 vi += nz; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
97 const PetscInt n = a->mbs, *aj = a->j, *adiag = a->diag, *vi; in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering() local
[all …]
|
| H A D | baijsolvnat2.c | 15 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
31 vi = aj + ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
37 jdx = 2 * (*vi++); in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
50 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
56 idx = 2 * (*vi++); in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
77 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_2_NaturalOrdering() local
92 vi = aj + ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
97 PetscPrefetchBlock(vi + nz, nz, 0, PETSC_PREFETCH_HINT_NTA); in MatSolve_SeqBAIJ_2_NaturalOrdering()
100 jdx = 2 * vi[k]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
114 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_2_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat1.c | 6 const PetscInt *adiag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering() local
24 vi = aj + adiag[i + 1] + 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
28 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
35 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
38 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
55 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace() local
68 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
70 while (nz--) x[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
76 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
79 while (nz--) x[*vi--] -= (*v--) * s1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
|
| H A D | baijsolvtran1.c | 8 const PetscInt *rout, *cout, *r, *c, *adiag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_1() local
31 vi = aj + adiag[i + 1] + 1; in MatSolveTranspose_SeqBAIJ_1()
35 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
42 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_1()
45 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
65 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_1_inplace() local
89 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_1_inplace()
91 while (nz--) t[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_inplace()
97 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_inplace()
100 while (nz--) t[*vi--] -= (*v--) * s1; in MatSolveTranspose_SeqBAIJ_1_inplace()
|
| H A D | baijsolvtrannat2.c | 7 const PetscInt *diag = a->diag, *vi, n = a->mbs, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace() local
26 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
29 oidx = 2 * (*vi++); in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
41 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
47 idx = 2 * (*vi--); in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
61 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering() local
82 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
85 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
97 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
103 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
|
| H A D | baijsolvtrannat3.c | 6 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
28 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
31 oidx = 3 * (*vi++); in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
45 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
52 idx = 3 * (*vi--); in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
67 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering() local
90 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
93 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
107 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
114 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
|
| H A D | baijsolvnat3.c | 12 const PetscInt *diag = a->diag, *vi; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
29 vi = aj + ai[i]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
36 jdx = 3 * (*vi++); in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
52 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
59 idx = 3 * (*vi++); in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
83 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_3_NaturalOrdering() local
101 vi = aj + ai[i]; in MatSolve_SeqBAIJ_3_NaturalOrdering()
108 jdx = bs * vi[k]; in MatSolve_SeqBAIJ_3_NaturalOrdering()
127 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_3_NaturalOrdering()
135 idx = bs * vi[k]; in MatSolve_SeqBAIJ_3_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrann.c | 8 const PetscInt *r, *c, *rout, *cout, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_N_inplace() local
37 vi = aj + a->diag[i] + 1; in MatSolveTranspose_SeqBAIJ_N_inplace()
40 PetscKernel_v_gets_v_minus_transA_times_w(bs, t + bs * (*vi++), v, t + i * bs); in MatSolveTranspose_SeqBAIJ_N_inplace()
48 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_N_inplace()
51 PetscKernel_v_gets_v_minus_transA_times_w(bs, t + bs * (*vi++), v, t + i * bs); in MatSolveTranspose_SeqBAIJ_N_inplace()
74 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *vi, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_N() local
103 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_N()
106 PetscKernel_v_gets_v_minus_transA_times_w(bs, t + bs * (vi[j]), v, t + i * bs); in MatSolveTranspose_SeqBAIJ_N()
114 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_N()
117 PetscKernel_v_gets_v_minus_transA_times_w(bs, t + bs * (vi[j]), v, t + i * bs); in MatSolveTranspose_SeqBAIJ_N()
|
| H A D | baijsolvtran2.c | 9 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_2_inplace() local
45 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_2_inplace()
48 oidx = 2 * (*vi++); in MatSolveTranspose_SeqBAIJ_2_inplace()
60 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_inplace()
66 idx = 2 * (*vi--); in MatSolveTranspose_SeqBAIJ_2_inplace()
94 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_2() local
131 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2()
134 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_2()
146 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_2()
152 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_2()
|
| H A D | baijsolv.c | 9 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *vi; in MatSolve_SeqBAIJ_N_inplace() local
31 vi = aj + ai[i]; in MatSolve_SeqBAIJ_N_inplace()
36 PetscKernel_v_gets_v_minus_A_times_w(bs, s, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
44 vi = aj + a->diag[i] + 1; in MatSolve_SeqBAIJ_N_inplace()
48 PetscKernel_v_gets_v_minus_A_times_w(bs, ls, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
68 const PetscInt *rout, *cout, *diag = a->diag, *vi, n = a->mbs; in MatSolve_SeqBAIJ_7_inplace() local
96 vi = aj + ai[i]; in MatSolve_SeqBAIJ_7_inplace()
107 idx = 7 * (*vi++); in MatSolve_SeqBAIJ_7_inplace()
136 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_7_inplace()
147 idx = 7 * (*vi++); in MatSolve_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtrannat4.c | 6 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace() local
30 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
33 oidx = 4 * (*vi++); in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
49 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
57 idx = 4 * (*vi--); in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
73 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering() local
98 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
101 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
117 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
125 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
|
| H A D | baijsolvtran3.c | 9 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_3_inplace() local
48 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_3_inplace()
51 oidx = 3 * (*vi++); in MatSolveTranspose_SeqBAIJ_3_inplace()
65 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3_inplace()
72 idx = 3 * (*vi--); in MatSolveTranspose_SeqBAIJ_3_inplace()
102 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_3() local
142 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3()
145 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_3()
159 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_3()
166 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_3()
|
| H A D | baijsolvnat4.c | 34 const PetscInt *vi; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
44 vi = aj + ai[i]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
52 jdx = 4 * (*vi++); in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
73 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
80 idx = 4 * (*vi++); in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
110 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_4_NaturalOrdering() local
129 vi = aj + ai[i]; in MatSolve_SeqBAIJ_4_NaturalOrdering()
137 jdx = bs * vi[k]; in MatSolve_SeqBAIJ_4_NaturalOrdering()
159 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_4_NaturalOrdering()
168 idx = bs * vi[k]; in MatSolve_SeqBAIJ_4_NaturalOrdering()
|
| H A D | baijsolvtrannat5.c | 6 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace() local
32 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
35 oidx = 5 * (*vi++); in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
53 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
62 idx = 5 * (*vi--); in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
79 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering() local
106 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
109 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
127 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
136 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
|
| H A D | baijsolvnat14.c | 10 …const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *adiag = a->diag, *vi, bs = A->rmap->bs, bs… in MatSolve_SeqBAIJ_14_NaturalOrdering() local
24 vi = aj + ai[i]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
42 idx = bs * vi[m]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
66 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_14_NaturalOrdering()
85 idx = bs * vi[m]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
136 …const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *adiag = a->diag, *vi, bs = A->rmap->bs, bs… in MatSolve_SeqBAIJ_13_NaturalOrdering() local
150 vi = aj + ai[i]; in MatSolve_SeqBAIJ_13_NaturalOrdering()
167 idx = bs * vi[m]; in MatSolve_SeqBAIJ_13_NaturalOrdering()
190 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_13_NaturalOrdering()
208 idx = bs * vi[m]; in MatSolve_SeqBAIJ_13_NaturalOrdering()
[all …]
|
| H A D | baijsolvnat5.c | 7 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() local
25 vi = aj + ai[i]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
34 jdx = 5 * (*vi++); in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
56 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
65 idx = 5 * (*vi++); in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
95 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_5_NaturalOrdering() local
113 vi = aj + ai[i]; in MatSolve_SeqBAIJ_5_NaturalOrdering()
122 jdx = 5 * vi[k]; in MatSolve_SeqBAIJ_5_NaturalOrdering()
145 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_5_NaturalOrdering()
154 idx = 5 * vi[k]; in MatSolve_SeqBAIJ_5_NaturalOrdering()
|
| H A D | baijsolvtran4.c | 9 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_4_inplace() local
51 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_4_inplace()
54 oidx = 4 * (*vi++); in MatSolveTranspose_SeqBAIJ_4_inplace()
70 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_4_inplace()
78 idx = 4 * (*vi--); in MatSolveTranspose_SeqBAIJ_4_inplace()
110 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_4() local
153 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_4()
156 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_4()
172 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_4()
180 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_4()
|
| H A D | baijsolvtrannat6.c | 6 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace() local
34 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
37 oidx = 6 * (*vi++); in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
57 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
67 idx = 6 * (*vi--); in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
85 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering() local
114 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering()
117 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering()
137 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering()
147 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering()
|
| H A D | baijsolvtran5.c | 9 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_5_inplace() local
54 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_5_inplace()
57 oidx = 5 * (*vi++); in MatSolveTranspose_SeqBAIJ_5_inplace()
75 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_5_inplace()
84 idx = 5 * (*vi--); in MatSolveTranspose_SeqBAIJ_5_inplace()
118 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_5() local
164 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_5()
167 oidx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_5()
185 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_5()
194 idx = bs * vi[j]; in MatSolveTranspose_SeqBAIJ_5()
|
| /petsc/src/snes/impls/vi/rs/ |
| H A D | virs.c | 21 SNES_VINEWTONRSLS *vi = (SNES_VINEWTONRSLS *)snes->data; in SNESVIGetInactiveSet() local
24 *inact = vi->IS_inact; in SNESVIGetInactiveSet()
298 SNES_VINEWTONRSLS *vi = (SNES_VINEWTONRSLS *)snes->data; in SNESSolve_VINEWTONRSLS() local
367 if (vi->checkredundancy) { in SNESSolve_VINEWTONRSLS()
368 PetscCall((*vi->checkredundancy)(snes, IS_act, &IS_redact, vi->ctxP)); in SNESSolve_VINEWTONRSLS()
371 PetscCall(ISComplement(IS_redact, X->map->rstart, X->map->rend, &vi->IS_inact)); in SNESSolve_VINEWTONRSLS()
374 PetscCall(ISComplement(IS_act, X->map->rstart, X->map->rend, &vi->IS_inact)); in SNESSolve_VINEWTONRSLS()
377 PetscCall(ISComplement(IS_act, X->map->rstart, X->map->rend, &vi->IS_inact)); in SNESSolve_VINEWTONRSLS()
381 …PetscCall(MatCreateSubMatrix(snes->jacobian, vi->IS_inact, vi->IS_inact, MAT_INITIAL_MATRIX, &jac_… in SNESSolve_VINEWTONRSLS()
397 PetscCall(ISGetIndices(vi->IS_inact, &inact)); in SNESSolve_VINEWTONRSLS()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 745 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatSolve_SeqAIJ_inplace() local
772 vi = aj + ai[i]; in MatSolve_SeqAIJ_inplace()
775 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_inplace()
782 vi = aj + adiag[i] + 1; in MatSolve_SeqAIJ_inplace()
785 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_inplace()
801 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatMatSolve_SeqAIJ_inplace() local
833 vi = aj + ai[i]; in MatMatSolve_SeqAIJ_inplace()
836 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatMatSolve_SeqAIJ_inplace()
842 vi = aj + adiag[i] + 1; in MatMatSolve_SeqAIJ_inplace()
845 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatMatSolve_SeqAIJ_inplace()
[all …]
|