| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 79 Vec StB0X, YtX, u, v; in DFPKernel_CompactDense() local
91 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in DFPKernel_CompactDense()
99 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in DFPKernel_CompactDense()
100 PetscCall(LMProductsSolve(YtS, oldest, next, YtX, YtX, /* ^H */ PETSC_FALSE)); in DFPKernel_CompactDense()
102 PetscCall(VecAXPBY(u, -1.0, 0.0, YtX)); in DFPKernel_CompactDense()
103 PetscCall(LMProductsMult(D, oldest, next, 1.0, YtX, 0.0, v, /* ^H */ PETSC_FALSE)); in DFPKernel_CompactDense()
104 PetscCall(LMProductsMultHermitian(StB0S, oldest, next, 1.0, YtX, 1.0, v)); in DFPKernel_CompactDense()
113 PetscCall(MatLMVMRestoreWorkRow(B, &YtX)); in DFPKernel_CompactDense()
122 Vec YtX = NULL, u = NULL; in DFPKernel_Dense() local
140 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in DFPKernel_Dense()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 24 Vec StBkX, YtX; in BFGSKernel_Recursive_Inner() local
30 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in BFGSKernel_Recursive_Inner()
34 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in BFGSKernel_Recursive_Inner()
35 PetscCall(LMProductsSolve(YtS, oldest, next, YtX, YtX, /* ^H */ PETSC_FALSE)); in BFGSKernel_Recursive_Inner()
36 PetscCall(LMBasisGEMV(Y, oldest, next, 1.0, YtX, 1.0, B0X)); in BFGSKernel_Recursive_Inner()
37 PetscCall(MatLMVMRestoreWorkRow(B, &YtX)); in BFGSKernel_Recursive_Inner()
189 Vec YtX, StB0X, u, v; in BFGSKernel_CompactDense() local
199 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in BFGSKernel_CompactDense()
204 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in BFGSKernel_CompactDense()
209 PetscCall(LMProductsSolve(D, oldest, next, YtX, YtX, /* ^H */ PETSC_FALSE)); in BFGSKernel_CompactDense()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 276 Vec YtX, StBHX; in BroydenKernelHermitianTranspose_Dense() local
281 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in BroydenKernelHermitianTranspose_Dense()
284 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in BroydenKernelHermitianTranspose_Dense()
285 PetscCall(VecAXPY(YtX, -1.0, StBHX)); in BroydenKernelHermitianTranspose_Dense()
286 PetscCall(LMProductsSolve(StS, oldest, next, YtX, YtX, PETSC_TRUE)); in BroydenKernelHermitianTranspose_Dense()
287 PetscCall(LMBasisGEMV(S, oldest, next, 1.0, YtX, 1.0, BHX)); in BroydenKernelHermitianTranspose_Dense()
289 PetscCall(MatLMVMRestoreWorkRow(B, &YtX)); in BroydenKernelHermitianTranspose_Dense()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 38 Vec StBkX, YtX, U, V; in SymBroydenKernel_Recursive_Inner() local
43 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in SymBroydenKernel_Recursive_Inner()
47 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in SymBroydenKernel_Recursive_Inner()
49 PetscCall(LMProductsMult(M01, oldest, next, 1.0, YtX, 1.0, U, PETSC_FALSE)); in SymBroydenKernel_Recursive_Inner()
51 PetscCall(LMProductsMult(M11, oldest, next, 1.0, YtX, 1.0, V, PETSC_FALSE)); in SymBroydenKernel_Recursive_Inner()
56 PetscCall(MatLMVMRestoreWorkRow(B, &YtX)); in SymBroydenKernel_Recursive_Inner()
504 Vec StB0X, YtX, u, v; in SymBroydenKernel_CompactDense() local
525 PetscCall(MatLMVMGetWorkRow(B, &YtX)); in SymBroydenKernel_CompactDense()
533 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, X, 0.0, YtX)); in SymBroydenKernel_CompactDense()
536 PetscCall(LMProductsMult(M01, oldest, next, 1.0, YtX, 1.0, u, PETSC_FALSE)); in SymBroydenKernel_CompactDense()
[all …]
|