| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 13 PETSC_INTERN PetscErrorCode DFPKernel_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec BX) in DFPKernel_Recursive() argument
46 PetscCall(MatLMVMApplyJ0Mode(mode)(B, G, BX)); in DFPKernel_Recursive()
52 PetscCall(VecDot(BX, s_i, &sitbx)); in DFPKernel_Recursive()
58 PetscCall(VecAXPY(BX, a - b, y_i)); in DFPKernel_Recursive()
68 PETSC_INTERN PetscErrorCode DFPKernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in DFPKernel_CompactDense() argument
73 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in DFPKernel_CompactDense()
97 else PetscCall(LMBasisGEMVH(S, oldest, next, 1.0, BX, 0.0, StB0X)); in DFPKernel_CompactDense()
108 PetscCall(LMBasisGEMV(Y, oldest, next, 1.0, v, 1.0, BX)); in DFPKernel_CompactDense()
109 PetscCall(MatLMVMBasisGEMV(B, B0S_t, oldest, next, 1.0, u, 1.0, BX)); in DFPKernel_CompactDense()
119 PETSC_INTERN PetscErrorCode DFPKernel_Dense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in DFPKernel_Dense() argument
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/sr1/ |
| H A D | sr1.c | 45 …e SR1Kernel_Recursive_Inner(Mat B, MatLMVMMode mode, PetscInt oldest, PetscInt next, Vec X, Vec BX) in SR1Kernel_Recursive_Inner() argument
59 PetscCall(LMBasisGEMV(Y_minus_BkS, oldest, next, 1.0, YmBkStX, 1.0, BX)); in SR1Kernel_Recursive_Inner()
126 static PetscErrorCode SR1Kernel_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec BX) in SR1Kernel_Recursive() argument
131 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in SR1Kernel_Recursive()
135 PetscCall(SR1Kernel_Recursive_Inner(B, mode, oldest, next, X, BX)); in SR1Kernel_Recursive()
203 static PetscErrorCode SR1Kernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in SR1Kernel_CompactDense() argument
208 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in SR1Kernel_CompactDense()
232 PetscCall(LMBasisGEMVH(S, oldest, next, -1.0, BX, 1.0, YmB0StX)); in SR1Kernel_CompactDense()
235 PetscCall(MatLMVMBasisGEMV(B, Y_minus_B0S_t, oldest, next, 1.0, v, 1.0, BX)); in SR1Kernel_CompactDense()
242 static PetscErrorCode MatMult_LMVMSR1_CompactDense(Mat B, Vec X, Vec BX) in MatMult_LMVMSR1_CompactDense() argument
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 47 PETSC_INTERN PetscErrorCode BroydenKernel_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BroydenKernel_Recursive() argument
50 Vec W = BX; in BroydenKernel_Recursive()
67 PetscCall(VecZeroEntries(BX)); in BroydenKernel_Recursive()
79 PetscCall(VecAXPY(BX, alphai, y_i)); in BroydenKernel_Recursive()
90 PetscCall(VecAXPY(BX, 1.0, W)); in BroydenKernel_Recursive()
153 PETSC_INTERN PetscErrorCode BroydenKernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BroydenKernel_CompactDense() argument
158 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in BroydenKernel_CompactDense()
172 PetscCall(MatLMVMBasisGEMV(B, Y_minus_B0S_t, oldest, next, 1.0, StX, 1.0, BX)); in BroydenKernel_CompactDense()
219 PETSC_INTERN PetscErrorCode BroydenKernel_Dense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BroydenKernel_Dense() argument
222 Vec W = BX; in BroydenKernel_Dense()
[all …]
|
| H A D | badbrdn.c | 203 …PetscErrorCode BadBroydenKernelHermitianTranspose_Recursive(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BadBroydenKernelHermitianTranspose_Recursive() argument
219 PetscCall(MatLMVMApplyJ0HermitianTransposeMode(mode)(B, G, BX)); in BadBroydenKernelHermitianTranspose_Recursive()
281 PETSC_INTERN PetscErrorCode BadBroydenKernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BadBroydenKernel_CompactDense() argument
286 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in BadBroydenKernel_CompactDense()
303 PetscCall(LMBasisGEMVH(Y, oldest, next, 1.0, BX, 0.0, YtB0X)); in BadBroydenKernel_CompactDense()
305 PetscCall(MatLMVMBasisGEMV(B, Y_minus_B0S_t, oldest, next, 1.0, v, 1.0, BX)); in BadBroydenKernel_CompactDense()
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 173 PETSC_INTERN PetscErrorCode BFGSKernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX) in BFGSKernel_CompactDense() argument
178 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in BFGSKernel_CompactDense()
207 else PetscCall(LMBasisGEMVH(S, oldest, next, 1.0, BX, 0.0, StB0X)); in BFGSKernel_CompactDense()
217 PetscCall(LMBasisGEMV(Y, oldest, next, 1.0, v, 1.0, BX)); in BFGSKernel_CompactDense()
218 PetscCall(MatLMVMBasisGEMV(B, B0S_t, oldest, next, 1.0, u, 1.0, BX)); in BFGSKernel_CompactDense()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 494 PETSC_INTERN PetscErrorCode SymBroydenKernel_CompactDense(Mat B, MatLMVMMode mode, Vec X, Vec BX, P… in SymBroydenKernel_CompactDense() argument
499 PetscCall(MatLMVMApplyJ0Mode(mode)(B, X, BX)); in SymBroydenKernel_CompactDense()
531 else PetscCall(LMBasisGEMVH(S, oldest, next, 1.0, BX, 0.0, StB0X)); in SymBroydenKernel_CompactDense()
540 PetscCall(LMBasisGEMV(Y, oldest, next, 1.0, v, 1.0, BX)); in SymBroydenKernel_CompactDense()
541 PetscCall(MatLMVMBasisGEMV(B, B0S_t, oldest, next, 1.0, u, 1.0, BX)); in SymBroydenKernel_CompactDense()
551 static PetscErrorCode MatMult_LMVMSymBrdn_CompactDense(Mat B, Vec X, Vec BX) in MatMult_LMVMSymBrdn_CompactDense() argument
559 PetscCall(BFGSKernel_CompactDense(B, MATLMVM_MODE_PRIMAL, X, BX)); in MatMult_LMVMSymBrdn_CompactDense()
561 PetscCall(DFPKernel_CompactDense(B, MATLMVM_MODE_PRIMAL, X, BX)); in MatMult_LMVMSymBrdn_CompactDense()
563 PetscCall(SymBroydenKernel_CompactDense(B, MATLMVM_MODE_PRIMAL, X, BX, PETSC_FALSE)); in MatMult_LMVMSymBrdn_CompactDense()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | surfacesphere_bin.msh | 21 …�|�AKm���Կ �B�X��?����g`*3տl���Z�����:��?�����䨟K��]��ƿ"�!�@�?����#�c�BX�� ��X"Ŀ���o��?����J…
|