| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 34 Vec s_i, y_i; in DFPKernel_Recursive() local
37 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in DFPKernel_Recursive()
38 PetscCall(VecDot(G, y_i, &yitx)); in DFPKernel_Recursive()
39 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in DFPKernel_Recursive()
48 Vec s_i, y_i; in DFPKernel_Recursive() local
57 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in DFPKernel_Recursive()
58 PetscCall(VecAXPY(BX, a - b, y_i)); in DFPKernel_Recursive()
59 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in DFPKernel_Recursive()
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 70 Vec s_i, y_i; in BroydenKernel_Recursive() local
74 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in BroydenKernel_Recursive()
79 PetscCall(VecAXPY(BX, alphai, y_i)); in BroydenKernel_Recursive()
82 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in BroydenKernel_Recursive()
125 Vec s_i, y_i; in BroydenKernelHermitianTranspose_Recursive() local
129 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in BroydenKernelHermitianTranspose_Recursive()
133 PetscCall(VecDotBegin(X, y_i, &yitx)); in BroydenKernelHermitianTranspose_Recursive()
135 PetscCall(VecDotEnd(X, y_i, &yitx)); in BroydenKernelHermitianTranspose_Recursive()
137 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in BroydenKernelHermitianTranspose_Recursive()
|
| H A D | badbrdn.c | 73 Vec y_i, yimbisi; in BadBroydenKernel_Recursive_Inner() local
77 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in BadBroydenKernel_Recursive_Inner()
78 PetscCall(VecDot(B0X, y_i, &yitbix)); in BadBroydenKernel_Recursive_Inner()
79 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in BadBroydenKernel_Recursive_Inner()
188 Vec yimbisi, y_i; in BadBroydenKernelHermitianTranspose_Recursive_Inner() local
196 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in BadBroydenKernelHermitianTranspose_Recursive_Inner()
197 PetscCall(VecAXPY(X, yimbisitx / PetscConj(yitbisi), y_i)); in BadBroydenKernelHermitianTranspose_Recursive_Inner()
198 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in BadBroydenKernelHermitianTranspose_Recursive_Inner()
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 207 Vec s_i, y_i; in SymBroydenRescaleUpdateDiagonal() local
210 PetscCall(LMBasisGetVecRead(Y, oldest + i, &y_i)); in SymBroydenRescaleUpdateDiagonal()
214 PetscCall(VecPointwiseMult(ldb->V, y_i, D_minus_beta_squared)); in SymBroydenRescaleUpdateDiagonal()
217 PetscCall(VecDotRealPart(ldb->V, y_i, &ytDinvy)); in SymBroydenRescaleUpdateDiagonal()
226 PetscCall(VecDotNorm2(y_i, ldb->W, &ytDs_scalar, &stDsr)); in SymBroydenRescaleUpdateDiagonal()
234 PetscCall(VecPointwiseMult(ldb->V, y_i, D_minus_beta)); in SymBroydenRescaleUpdateDiagonal()
243 PetscCall(VecPointwiseMult(ldb->V, y_i, D_minus_beta)); in SymBroydenRescaleUpdateDiagonal()
251 PetscCall(LMBasisRestoreVecRead(Y, oldest + i, &y_i)); in SymBroydenRescaleUpdateDiagonal()
|
| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmvmutils.c | 1013 Vec y_i, op_x_i, y_m_op_x_i; in MatLMVMUpdateOpDiffVecs() local
1015 PetscCall(LMBasisGetVecRead(Y, i, &y_i)); in MatLMVMUpdateOpDiffVecs()
1018 PetscCall(VecAXPBYPCZ(y_m_op_x_i, 1.0, -alpha, 0.0, y_i, op_x_i)); in MatLMVMUpdateOpDiffVecs()
1021 PetscCall(LMBasisRestoreVecRead(Y, i, &y_i)); in MatLMVMUpdateOpDiffVecs()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexgeometry.c | 499 PetscReal y_i = PetscRealPart(coords[faces[2 * f + 0] * embedDim + 1]); in DMPlexLocatePoint_Quad_2D_Linear_Internal() local
508 if ((y_j > y) != (y_i > y)) { in DMPlexLocatePoint_Quad_2D_Linear_Internal()
509 PetscReal slope = (x - x_j) * (y_i - y_j) - (x_i - x_j) * (y - y_j); in DMPlexLocatePoint_Quad_2D_Linear_Internal()
515 if ((slope < 0) != (y_i < y_j)) ++crossings; in DMPlexLocatePoint_Quad_2D_Linear_Internal()
|
| /petsc/doc/manual/ |
| H A D | tao.md | 2417 $\beta(x) = ||Dx||_1 \approx \sum_{i} \sqrt{y_i^2 + \epsilon^2}-\epsilon$
|