| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmvmimpl.c | 79 PetscCall(VecDestroy(&lmvm->Xprev)); in MatReset_LMVM()
95 …PetscCheck(lmvm->Xprev != NULL && lmvm->Fprev != NULL, PetscObjectComm((PetscObject)B), PETSC_ERR_… in MatLMVMAllocateBases()
96 …if (!lmvm->basis[LMBASIS_S]) PetscCall(LMBasisCreate(lmvm->Xprev, lmvm->m, &lmvm->basis[LMBASIS_S]… in MatLMVMAllocateBases()
106 if (!lmvm->Xprev) PetscCall(MatCreateVecs(B, &lmvm->Xprev, NULL)); in MatLMVMAllocateVecs()
161 PetscCall(VecAXPBY(lmvm->Xprev, 1.0, -1.0, X)); in MatUpdate_LMVM()
164 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVM()
168 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_LMVM()
249 PetscCall(MatLMVMAllocate(M, bctx->Xprev, bctx->Fprev)); in MatCopy_LMVM()
260 PetscCall(VecCopy(bctx->Xprev, mctx->Xprev)); in MatCopy_LMVM()
297 PetscCall(MatLMVMAllocate(*mat, bctx->Xprev, bctx->Fprev)); in MatDuplicate_LMVM()
[all …]
|
| H A D | lmvm.h | 89 Vec Xprev, Fprev; member
|
| H A D | lmvmutils.c | 1058 …e]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmv… in MatLMVMGetUpdatedBasis()
1071 …e]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmv… in MatLMVMGetUpdatedBasis()
|
| /petsc/src/ksp/ksp/utils/lmvm/diagbrdn/ |
| H A D | diagbrdn.c | 36 PetscCall(VecAYPX(lmvm->Xprev, -1.0, X)); in MatUpdate_DiagBrdn()
40 PetscCall(VecDotNorm2(lmvm->Fprev, lmvm->Xprev, &curvature, &ststmp)); in MatUpdate_DiagBrdn()
47 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_DiagBrdn()
58 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_DiagBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/sr1/ |
| H A D | sr1.c | 83 …PetscCall(LMBasisCreate(mode == MATLMVM_MODE_PRIMAL ? lmvm->Fprev : lmvm->Xprev, lmvm->m, &lsr1->b… in SR1RecursiveBasisUpdate()
310 PetscCall(VecAYPX(lmvm->Xprev, -1.0, X)); in MatUpdate_LMVMSR1()
330 PetscCall(MatMult(B, lmvm->Xprev, work)); in MatUpdate_LMVMSR1()
332 PetscCall(VecDot(lmvm->Xprev, work, &sktw)); in MatUpdate_LMVMSR1()
333 PetscCall(VecNorm(lmvm->Xprev, NORM_2, &snorm)); in MatUpdate_LMVMSR1()
338 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMSR1()
362 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_LMVMSR1()
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 439 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->invDnew)); in SymBroydenRescaleAllocate()
440 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->BFGS)); in SymBroydenRescaleAllocate()
441 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->DFP)); in SymBroydenRescaleAllocate()
442 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->U)); in SymBroydenRescaleAllocate()
443 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->V)); in SymBroydenRescaleAllocate()
444 PetscCall(VecDuplicate(lmvm->Xprev, &ldb->W)); in SymBroydenRescaleAllocate()
|
| /petsc/src/ksp/ksp/utils/lmvm/dense/ |
| H A D | denseqn.c | 165 PetscCall(VecGetType(lmvm->Xprev, &vec_type)); in MatAllocate_LMVMDQN_Internal()
194 PetscCall(VecDuplicateVecs(lmvm->Xprev, lmvm->m, &lqn->PQ)); in MatAllocate_LMVMDQN_Internal()
195 PetscCall(VecDuplicate(lmvm->Xprev, &lqn->column_work2)); in MatAllocate_LMVMDQN_Internal()
209 PetscCall(VecDuplicate(lmvm->Xprev, &lqn->column_work)); in MatAllocate_LMVMDQN_Internal()
294 PetscCall(VecAYPX(lmvm->Xprev, -1.0, X)); in MatUpdate_LMVMDQN()
298 PetscCall(VecMDot(lmvm->Xprev, 2, FX, dotFX)); in MatUpdate_LMVMDQN()
316 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMDQN()
385 … PetscCall(MatMultHermitianTransposeColumnRange(Yfull, lmvm->Xprev, this_ys_col, 0, h_new)); in MatUpdate_LMVMDQN()
431 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_LMVMDQN()
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 363 PetscCall(VecAYPX(lmvm->Xprev, -1.0, X)); in MatUpdate_LMVMBrdn()
366 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMBrdn()
386 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_LMVMBrdn()
|
| H A D | badbrdn.c | 108 …us_BkS_t]) PetscCall(LMBasisCreate(Y_t == LMBASIS_Y ? lmvm->Fprev : lmvm->Xprev, lmvm->m, &lbrdn->… in BadBroydenRecursiveBasisUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 111 …]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(B0S_t) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmv… in SymBroydenRecursiveBasisUpdate()
616 PetscCall(VecAYPX(lmvm->Xprev, -1.0, X)); in MatUpdate_LMVMSymBrdn()
621 Vec sy[2] = {lmvm->Xprev, lmvm->Fprev}; in MatUpdate_LMVMSymBrdn()
624 PetscCall(VecMDot(lmvm->Xprev, 2, sy, stsy)); in MatUpdate_LMVMSymBrdn()
662 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMSymBrdn()
722 PetscCall(VecCopy(X, lmvm->Xprev)); in MatUpdate_LMVMSymBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 61 …]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(B0S_t) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmv… in BFGSRecursiveBasisUpdate()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | Mat.pyx | 5254 """Adds (X-Xprev) and (F-Fprev) updates to LMVM matrix.
|