| /petsc/src/ksp/ksp/utils/lmvm/brdn/ |
| H A D | badbrdn.c | 64 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenKernel_Recursive_Inner() 96 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenRecursiveBasisUpdate() 108 …if (!lbrdn->basis[Y_minus_BkS_t]) PetscCall(LMBasisCreate(Y_t == LMBASIS_Y ? lmvm->Fprev : lmvm->X… in BadBroydenRecursiveBasisUpdate() 175 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenKernelHermitianTranspose_Recursive_Inner() 205 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenKernelHermitianTranspose_Recursive() 236 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenCompactProductsUpdate() 291 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenKernel_CompactDense() 321 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BadBroydenKernelHermitianTranspose_CompactDense()
|
| H A D | brdn.c | 52 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BroydenKernel_Recursive() 116 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BroydenKernelHermitianTranspose_Recursive() 224 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BroydenKernel_Dense() 273 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BroydenKernelHermitianTranspose_Dense() 352 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &Y, NULL, NULL)); in MatUpdate_LMVMBrdn() 357 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_Y, LMBLOCK_UPPER_TRIANGLE, &YtY)); in MatUpdate_LMVMBrdn() 382 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &Y, NULL, NULL)); in MatUpdate_LMVMBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/bfgs/ |
| H A D | bfgs.c | 19 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BFGSKernel_Recursive_Inner() 28 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_S, LMBLOCK_DIAGONAL, &YtS)); in BFGSKernel_Recursive_Inner() 76 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_S, LMBLOCK_DIAGONAL, &YtS)); in BFGSRecursiveBasisUpdate() 120 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BFGSCompactDenseProductsUpdate() 136 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_S, LMBLOCK_DIAGONAL, &D)); in BFGSCompactDenseProductsUpdate() 184 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in BFGSKernel_CompactDense() 196 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, LMBASIS_S, LMBLOCK_DIAGONAL, &D)); in BFGSKernel_CompactDense()
|
| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmvmutils.c | 1038 case LMBASIS_Y: in MatLMVMGetUpdatedBasis() 1050 *basis_p = lmvm->basis[type == LMBASIS_B0S ? LMBASIS_S : LMBASIS_Y]; in MatLMVMGetUpdatedBasis() 1051 *returned_type = (type == LMBASIS_B0S) ? LMBASIS_S : LMBASIS_Y; in MatLMVMGetUpdatedBasis() 1054 LMBasis orig_basis = (type == LMBASIS_B0S) ? lmvm->basis[LMBASIS_S] : lmvm->basis[LMBASIS_Y]; in MatLMVMGetUpdatedBasis() 1098 PetscCall(LMBasisGetVecRead(lmvm->basis[type == LMBASIS_B0S ? LMBASIS_S : LMBASIS_Y], i, y)); in MatLMVMBasisGetVecRead() 1102 LMBasis orig_basis = (type == LMBASIS_B0S) ? lmvm->basis[LMBASIS_S] : lmvm->basis[LMBASIS_Y]; in MatLMVMBasisGetVecRead() 1143 … PetscCall(LMBasisRestoreVecRead(lmvm->basis[type == LMBASIS_B0S ? LMBASIS_S : LMBASIS_Y], i, y)); in MatLMVMBasisRestoreVecRead() 1172 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &basis, NULL, NULL)); in MatLMVMGetWorkRow() 1182 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &basis, NULL, NULL)); in MatLMVMRestoreWorkRow() 1229 if (cache_J0_products || type == LMBASIS_S || type == LMBASIS_Y) { in MatLMVMBasisGEMVH() [all …]
|
| H A D | lmvmimpl.c | 97 …if (!lmvm->basis[LMBASIS_Y]) PetscCall(LMBasisCreate(lmvm->Fprev, lmvm->m, &lmvm->basis[LMBASIS_Y]… in MatLMVMAllocateBases() 143 PetscCall(LMBasisGetNextVec(lmvm->basis[LMBASIS_Y], &y_k)); in MatUpdateKernel_LMVM() 145 PetscCall(LMBasisRestoreNextVec(lmvm->basis[LMBASIS_Y], &y_k)); in MatUpdateKernel_LMVM() 149 …PetscAssert(lmvm->k == lmvm->basis[LMBASIS_Y]->k, PetscObjectComm((PetscObject)B), PETSC_ERR_PLIB,… in MatUpdateKernel_LMVM() 264 PetscCall(LMBasisCopy(bctx->basis[LMBASIS_Y], mctx->basis[LMBASIS_Y])); in MatCopy_LMVM()
|
| H A D | lmvm.h | 62 LMBASIS_Y = 1, // differences in function values, Y_i = (F_{i+1} - F_i) enumerator 80 #define MatLMVMBasisSizeOf(type) ((type) & LMBASIS_Y)
|
| /petsc/src/ksp/ksp/utils/lmvm/dfp/ |
| H A D | dfp.c | 16 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in DFPKernel_Recursive() 77 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in DFPKernel_CompactDense() 131 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in DFPKernel_Dense()
|
| /petsc/src/ksp/ksp/utils/lmvm/sr1/ |
| H A D | sr1.c | 72 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in SR1RecursiveBasisUpdate() 226 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in SR1Kernel_CompactDense() 294 …PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_S_MINUS_H0Y, LMBASIS_Y, LMBLOCK_UPPER_TRIANGLE, &Sm… in MatUpdate_LMVMSR1() 328 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &Y, NULL, NULL)); in MatUpdate_LMVMSR1()
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 134 PetscCall(MatLMVMProductsGetDiagonalValue(B, LMBASIS_Y, LMBASIS_S, next - 1, &yts)); in SymBroydenRescaleUpdateDiagonal() 136 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &Y, NULL, NULL)); in SymBroydenRescaleUpdateDiagonal() 338 PetscCall(MatLMVMProductsGetDiagonalValue(B, LMBASIS_Y, LMBASIS_Y, i, &yty)); in SymBroydenRescaleUpdate() 339 PetscCall(MatLMVMProductsGetDiagonalValue(B, LMBASIS_Y, LMBASIS_S, i, &yts)); in SymBroydenRescaleUpdate()
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 32 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in SymBroydenKernel_Recursive_Inner() 156 PetscCall(MatLMVMProductsGetDiagonalValue(B, LMBASIS_Y, LMBASIS_S, j, &yjtsj)); in SymBroydenRecursiveBasisUpdate() 346 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in SymBroydenCompactProductsUpdate() 506 MatLMVMBasisType Y_t = LMVMModeMap(LMBASIS_Y, mode); in SymBroydenKernel_CompactDense() 610 PetscCall(MatLMVMGetUpdatedBasis(B, LMBASIS_Y, &Y, NULL, NULL)); in MatUpdate_LMVMSymBrdn() 641 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_S, LMBASIS_Y, LMBLOCK_UPPER_TRIANGLE, &StY)); in MatUpdate_LMVMSymBrdn() 648 PetscCall(MatLMVMGetUpdatedProducts(B, LMBASIS_Y, H0Y_t, LMBLOCK_UPPER_TRIANGLE, &YtH0Y)); in MatUpdate_LMVMSymBrdn() 682 …PetscCall(MatLMVMProductsInsertDiagonalValue(B, LMBASIS_Y, LMBASIS_S, next, PetscRealPart(curvatur… in MatUpdate_LMVMSymBrdn() 683 …PetscCall(MatLMVMProductsInsertDiagonalValue(B, LMBASIS_S, LMBASIS_Y, next, PetscRealPart(curvatur… in MatUpdate_LMVMSymBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/diagbrdn/ |
| H A D | diagbrdn.c | 48 …PetscCall(MatLMVMProductsInsertDiagonalValue(B, LMBASIS_Y, LMBASIS_S, next, PetscRealPart(curvatur… in MatUpdate_DiagBrdn()
|
| /petsc/src/ksp/ksp/utils/lmvm/dense/ |
| H A D | denseqn.c | 275 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatUpdate_LMVMDQN() 672 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatLMVMDBFGSUpdateMultData() 767 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatSolve_LMVMDBFGS() 856 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatMult_LMVMDBFGS() 1097 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatLMVMDDFPUpdateSolveData() 1220 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatSolve_LMVMDDFP() 1376 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatMult_LMVMDDFP()
|