| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmvmimpl.c | 80 PetscCall(VecDestroy(&lmvm->Fprev)); in MatReset_LMVM()
95 …PetscCheck(lmvm->Xprev != NULL && lmvm->Fprev != NULL, PetscObjectComm((PetscObject)B), PETSC_ERR_… in MatLMVMAllocateBases()
97 …if (!lmvm->basis[LMBASIS_Y]) PetscCall(LMBasisCreate(lmvm->Fprev, lmvm->m, &lmvm->basis[LMBASIS_Y]… in MatLMVMAllocateBases()
107 if (!lmvm->Fprev) PetscCall(MatCreateVecs(B, NULL, &lmvm->Fprev)); in MatLMVMAllocateVecs()
162 PetscCall(VecAXPBY(lmvm->Fprev, 1.0, -1.0, F)); in MatUpdate_LMVM()
164 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVM()
169 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_LMVM()
249 PetscCall(MatLMVMAllocate(M, bctx->Xprev, bctx->Fprev)); in MatCopy_LMVM()
261 PetscCall(VecCopy(bctx->Fprev, mctx->Fprev)); 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 …LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmvm->m, &lmvm->b… in MatLMVMGetUpdatedBasis()
1071 …LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmvm->m, &lmvm->b… in MatLMVMGetUpdatedBasis()
|
| /petsc/src/ksp/ksp/utils/lmvm/diagbrdn/ |
| H A D | diagbrdn.c | 37 PetscCall(VecAYPX(lmvm->Fprev, -1.0, F)); 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()
59 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_DiagBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/sr1/ |
| H A D | sr1.c | 83 …minus_BkS_t]) PetscCall(LMBasisCreate(mode == MATLMVM_MODE_PRIMAL ? lmvm->Fprev : lmvm->Xprev, lmv… in SR1RecursiveBasisUpdate()
296 PetscCall(VecCopy(lmvm->Fprev, Fprev_old)); in MatUpdate_LMVMSR1()
311 PetscCall(VecAYPX(lmvm->Fprev, -1.0, F)); in MatUpdate_LMVMSR1()
331 PetscCall(VecAYPX(work, -1.0, lmvm->Fprev)); in MatUpdate_LMVMSR1()
338 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMSR1()
363 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_LMVMSR1()
|
| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | brdn.c | 359 PetscCall(VecCopy(lmvm->Fprev, Fprev_old)); in MatUpdate_LMVMBrdn()
364 PetscCall(VecAYPX(lmvm->Fprev, -1.0, F)); in MatUpdate_LMVMBrdn()
366 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMBrdn()
387 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_LMVMBrdn()
|
| H A D | badbrdn.c | 108 …if (!lbrdn->basis[Y_minus_BkS_t]) PetscCall(LMBasisCreate(Y_t == LMBASIS_Y ? lmvm->Fprev : lmvm->X… in BadBroydenRecursiveBasisUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 111 …MBasisCreate(MatLMVMBasisSizeOf(B0S_t) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmvm->m, &lsb->ba… in SymBroydenRecursiveBasisUpdate()
612 PetscCall(VecCopy(lmvm->Fprev, Fprev_old)); in MatUpdate_LMVMSymBrdn()
617 PetscCall(VecAYPX(lmvm->Fprev, -1.0, F)); in MatUpdate_LMVMSymBrdn()
621 Vec sy[2] = {lmvm->Xprev, lmvm->Fprev}; in MatUpdate_LMVMSymBrdn()
662 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMSymBrdn()
723 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_LMVMSymBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/dense/ |
| H A D | denseqn.c | 296 FX[0] = lmvm->Fprev; /* dotFX[0] = s^T Fprev */ in MatUpdate_LMVMDQN()
299 PetscCall(VecAYPX(lmvm->Fprev, -1.0, F)); in MatUpdate_LMVMDQN()
300 PetscCall(VecDot(lmvm->Fprev, lmvm->Fprev, &yTy)); in MatUpdate_LMVMDQN()
316 PetscCall(MatUpdateKernel_LMVM(B, lmvm->Xprev, lmvm->Fprev)); in MatUpdate_LMVMDQN()
432 PetscCall(VecCopy(F, lmvm->Fprev)); in MatUpdate_LMVMDQN()
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 61 …MBasisCreate(MatLMVMBasisSizeOf(B0S_t) == LMBASIS_S ? lmvm->Xprev : lmvm->Fprev, lmvm->m, &lsb->ba… 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.
|