| /petsc/src/dm/dt/tests/ |
| H A D | ex13.c | 30 const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints]; in constructTabulationAndMass() local
32 for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt]; in constructTabulationAndMass()
109 const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints]; in test() local
111 for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt]; in test()
178 const PetscReal *p_j; in test() local
189 p_j = &p_trimmed1[(b * Nf1 + j) * npoints]; in test()
190 for (PetscInt pt = 0; pt < npoints; pt++) p_i[pt] += p_j[pt] * points[pt * dim + k] * sign; in test()
202 const PetscReal *p_j = &p_koszul[(j * Nf + v) * npoints]; in test() local
204 for (PetscInt pt = 0; pt < npoints; pt++) val += p_i[pt] * p_j[pt] * weights[pt]; in test()
218 const PetscReal *p_j = &p_trimmed[(j * Nf + v) * Nk * npoints]; in test() local
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | badbrdn.c | 124 Vec p_j, y_j, B0s_j; in BadBroydenRecursiveBasisUpdate() local
127 PetscCall(LMBasisGetWorkVec(Y_minus_BkS, &p_j)); in BadBroydenRecursiveBasisUpdate()
131 PetscCall(VecAXPBY(p_j, alpha, 0.0, B0s_j)); in BadBroydenRecursiveBasisUpdate()
136 if (j > oldest) PetscCall(BadBroydenKernel_Recursive_Inner(B, mode, oldest, j, p_j)); in BadBroydenRecursiveBasisUpdate()
138 PetscCall(VecDot(p_j, y_j, &yjtbjsj)); in BadBroydenRecursiveBasisUpdate()
140 PetscCall(VecAYPX(p_j, -1.0, y_j)); in BadBroydenRecursiveBasisUpdate()
142 PetscCall(LMBasisSetNextVec(Y_minus_BkS, p_j)); in BadBroydenRecursiveBasisUpdate()
143 PetscCall(LMBasisRestoreWorkVec(Y_minus_BkS, &p_j)); in BadBroydenRecursiveBasisUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/sr1/ |
| H A D | sr1.c | 100 Vec s_j, B0s_j, p_j, y_j; in SR1RecursiveBasisUpdate() local
103 PetscCall(LMBasisGetWorkVec(Y_minus_BkS, &p_j)); in SR1RecursiveBasisUpdate()
107 PetscCall(VecAXPBY(p_j, alpha, 0.0, B0s_j)); in SR1RecursiveBasisUpdate()
113 if (j > oldest) PetscCall(SR1Kernel_Recursive_Inner(B, mode, oldest, j, s_j, p_j)); in SR1RecursiveBasisUpdate()
115 PetscCall(VecAYPX(p_j, -1.0, y_j)); in SR1RecursiveBasisUpdate()
116 PetscCall(VecDot(s_j, p_j, &ymbksts)); in SR1RecursiveBasisUpdate()
120 PetscCall(LMBasisSetNextVec(Y_minus_BkS, p_j)); in SR1RecursiveBasisUpdate()
121 PetscCall(LMBasisRestoreWorkVec(Y_minus_BkS, &p_j)); in SR1RecursiveBasisUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 78 Vec p_j, s_j, B0s_j; in BFGSRecursiveBasisUpdate() local
81 PetscCall(LMBasisGetWorkVec(BkS, &p_j)); in BFGSRecursiveBasisUpdate()
84 PetscCall(VecAXPBY(p_j, alpha, 0.0, B0s_j)); in BFGSRecursiveBasisUpdate()
89 if (j > oldest) PetscCall(BFGSKernel_Recursive_Inner(B, mode, oldest, j, s_j, p_j)); in BFGSRecursiveBasisUpdate()
90 PetscCall(VecDot(p_j, s_j, &sjtbjsj)); in BFGSRecursiveBasisUpdate()
93 PetscCall(LMBasisSetNextVec(BkS, p_j)); in BFGSRecursiveBasisUpdate()
94 PetscCall(LMBasisRestoreWorkVec(BkS, &p_j)); in BFGSRecursiveBasisUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 136 Vec p_j, s_j, B0s_j; in SymBroydenRecursiveBasisUpdate() local
141 PetscCall(LMBasisGetWorkVec(BkS, &p_j)); in SymBroydenRecursiveBasisUpdate()
144 PetscCall(VecAXPBY(p_j, alpha, 0.0, B0s_j)); in SymBroydenRecursiveBasisUpdate()
149 if (j > oldest) PetscCall(SymBroydenKernel_Recursive_Inner(B, mode, oldest, j, s_j, p_j)); in SymBroydenRecursiveBasisUpdate()
150 PetscCall(VecDot(p_j, s_j, &sjtbjsj)); in SymBroydenRecursiveBasisUpdate()
153 PetscCall(LMBasisSetNextVec(BkS, p_j)); in SymBroydenRecursiveBasisUpdate()
154 PetscCall(LMBasisRestoreWorkVec(BkS, &p_j)); in SymBroydenRecursiveBasisUpdate()
|