| /petsc/src/tao/bound/utils/ |
| H A D | isutil.c | 184 PetscErrorCode TaoEstimateActiveBounds(Vec X, Vec XL, Vec XU, Vec G, Vec S, Vec W, PetscReal steple… in TaoEstimateActiveBounds() argument 197 if (XL) PetscValidHeaderSpecific(XL, VEC_CLASSID, 2); in TaoEstimateActiveBounds() 203 if (XL) PetscCheckSameType(X, 1, XL, 2); in TaoEstimateActiveBounds() 208 if (XL) PetscCheckSameComm(X, 1, XL, 2); in TaoEstimateActiveBounds() 213 if (XL) VecCheckSameSize(X, 1, XL, 2); in TaoEstimateActiveBounds() 222 PetscCall(TaoBoundSolution(W, XL, XU, 0.0, &nDiff, W)); in TaoEstimateActiveBounds() 236 if (!XL && !XU) { in TaoEstimateActiveBounds() 242 PetscCall(VecGetArrayRead(XL, &xl)); in TaoEstimateActiveBounds() 279 PetscCall(VecRestoreArrayRead(XL, &xl)); in TaoEstimateActiveBounds() 344 PetscErrorCode TaoBoundStep(Vec X, Vec XL, Vec XU, IS active_lower, IS active_upper, IS active_fixe… in TaoBoundStep() argument [all …]
|
| /petsc/src/tao/interface/ |
| H A D | taosolver_bounds.c | 17 PetscErrorCode TaoSetVariableBounds(Tao tao, Vec XL, Vec XU) in TaoSetVariableBounds() argument 21 if (XL) PetscValidHeaderSpecific(XL, VEC_CLASSID, 2); in TaoSetVariableBounds() 23 PetscCall(PetscObjectReference((PetscObject)XL)); in TaoSetVariableBounds() 25 PetscCall(VecDestroy(&tao->XL)); in TaoSetVariableBounds() 27 tao->XL = XL; in TaoSetVariableBounds() 29 tao->bounded = (PetscBool)(XL || XU); in TaoSetVariableBounds() 82 PetscErrorCode TaoGetVariableBounds(Tao tao, Vec *XL, Vec *XU) in TaoGetVariableBounds() argument 86 if (XL) *XL = tao->XL; in TaoGetVariableBounds() 109 if (!tao->XL) { in TaoComputeVariableBounds() 110 PetscCall(VecDuplicate(tao->solution, &tao->XL)); in TaoComputeVariableBounds() [all …]
|
| /petsc/src/tao/quadratic/impls/gpcg/ |
| H A D | gpcg.c | 117 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_GPCG() 118 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSolve_GPCG() 130 …PetscCall(VecWhichInactive(tao->XL, tao->solution, tao->gradient, tao->XU, PETSC_TRUE, &gpcg->Free… in TaoSolve_GPCG() 134 PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, gpcg->PG)); in TaoSolve_GPCG() 195 … PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, gpcg->PG)); in TaoSolve_GPCG() 199 …PetscCall(VecWhichInactive(tao->XL, tao->solution, tao->gradient, tao->XU, PETSC_TRUE, &gpcg->Free… in TaoSolve_GPCG() 224 …Vec DX = tao->stepdirection, XL = tao->XL, XU = tao->XU, Work = gpcg->Wor… in GPCGGradProjections() local 234 PetscCall(VecBoundGradientProjection(G, X, XL, XU, DX)); in GPCGGradProjections() 259 PetscCall(VecWhichInactive(XL, X, tao->gradient, XU, PETSC_TRUE, &gpcg->Free_Local)); in GPCGGradProjections() 271 PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, gpcg->Work)); in TaoComputeDual_GPCG()
|
| /petsc/src/tao/complementarity/impls/ssls/ |
| H A D | ssls.c | 30 PetscCall(VecFischer(X, tao->constraints, tao->XL, tao->XU, ssls->ff)); in Tao_SSLS_Function() 43 PetscCall(VecFischer(X, tao->constraints, tao->XL, tao->XU, ssls->ff)); in Tao_SSLS_FunctionGradient() 49 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, ssls->t1, … in Tao_SSLS_FunctionGradient()
|
| H A D | ssfls.c | 17 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSetUp_SSFLS() 36 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_SSFLS() 67 PetscCall(VecBoundGradientProjection(ssls->w, tao->solution, tao->XL, tao->XU, ssls->w)); in TaoSolve_SSFLS()
|
| H A D | ssils.c | 49 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_SSILS()
|
| /petsc/src/tao/bound/impls/tron/ |
| H A D | tron.c | 80 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSolve_TRON() 83 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_TRON() 91 …PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, tao->gradient… in TaoSolve_TRON() 116 …PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, tao->gradient… in TaoSolve_TRON() 127 …PetscCall(VecWhichInactive(tao->XL, tao->solution, tao->gradient, tao->XU, PETSC_TRUE, &tron->Free… in TaoSolve_TRON() 207 … PetscCall(VecBoundGradientProjection(tron->G_New, tron->X_New, tao->XL, tao->XU, tao->gradient)); in TaoSolve_TRON() 209 …PetscCall(VecWhichInactive(tao->XL, tron->X_New, tao->gradient, tao->XU, PETSC_TRUE, &tron->Free_L… in TaoSolve_TRON() 260 …PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, tao->gradient… in TronGradientProjections() 281 PetscCall(VecBoundGradientProjection(tao->gradient, tao->solution, tao->XL, tao->XU, tron->Work)); in TaoComputeDual_TRON()
|
| /petsc/src/vec/vec/utils/ |
| H A D | projection.c | 754 PetscErrorCode VecBoundGradientProjection(Vec G, Vec X, Vec XL, Vec XU, Vec GP) in VecBoundGradientProjection() argument 765 if (XL) PetscValidHeaderSpecific(XL, VEC_CLASSID, 3); in VecBoundGradientProjection() 768 if (!XL && !XU) { in VecBoundGradientProjection() 776 PetscCall(VecGetArrayRead(XL, &xlptr)); in VecBoundGradientProjection() 792 PetscCall(VecRestoreArrayRead(XL, &xlptr)); in VecBoundGradientProjection() 816 PetscErrorCode VecStepMaxBounded(Vec X, Vec DX, Vec XL, Vec XU, PetscReal *stepmax) in VecStepMaxBounded() argument 825 PetscValidHeaderSpecific(XL, VEC_CLASSID, 3); in VecStepMaxBounded() 829 PetscCall(VecGetArrayRead(XL, &xl)); in VecStepMaxBounded() 841 PetscCall(VecRestoreArrayRead(XL, &xl)); in VecStepMaxBounded() 871 PetscErrorCode VecStepBoundInfo(Vec X, Vec DX, Vec XL, Vec XU, PetscReal *boundmin, PetscReal *wolf… in VecStepBoundInfo() argument [all …]
|
| /petsc/src/tao/quadratic/impls/bqpip/ |
| H A D | bqpip.c | 55 PetscCall(VecAXPY(qp->G, -1.0, qp->XL)); in QPIPSetInitialPoint() 100 PetscCall(VecAXPY(qp->R3, -1.0, qp->XL)); in QPIPSetInitialPoint() 213 PetscCall(VecDuplicate(tao->solution, &qp->XL)); in TaoSetup_BQPIP() 268 PetscCall(VecSet(qp->XL, -1.0e20)); in TaoSolve_BQPIP() 269 if (tao->XL) PetscCall(VecPointwiseMax(qp->XL, qp->XL, tao->XL)); in TaoSolve_BQPIP() 271 PetscCall(VecMedian(qp->XL, tao->solution, qp->XU, tao->solution)); in TaoSolve_BQPIP() 488 PetscCall(VecDestroy(&qp->XL)); in TaoDestroy_BQPIP()
|
| H A D | bqpipimpl.h | 41 Vec XL, XU; member
|
| /petsc/src/tao/util/ |
| H A D | tao_util.c | 247 PetscErrorCode MatDFischer(Mat jac, Vec X, Vec Con, Vec XL, Vec XU, Vec T1, Vec T2, Vec Da, Vec Db) in MatDFischer() argument 258 PetscCall(VecGetArrayRead(XL, &l)); in MatDFischer() 362 PetscCall(VecRestoreArrayRead(XL, &l)); in MatDFischer() 393 PetscErrorCode MatDSFischer(Mat jac, Vec X, Vec Con, Vec XL, Vec XU, PetscReal mu, Vec T1, Vec T2, … in MatDSFischer() argument 403 PetscCall(MatDFischer(jac, X, Con, XL, XU, T1, T2, Da, Db)); in MatDSFischer() 408 PetscCall(VecGetArrayRead(XL, &l)); in MatDSFischer() 464 PetscCall(VecRestoreArrayRead(XL, &l)); in MatDSFischer()
|
| /petsc/src/tao/bound/impls/blmvm/ |
| H A D | blmvm.c | 15 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_BLMVM() 16 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSolve_BLMVM() 20 …PetscCall(VecBoundGradientProjection(blmP->unprojected_gradient, tao->solution, tao->XL, tao->XU, … in TaoSolve_BLMVM() 56 …PetscCall(VecBoundGradientProjection(tao->stepdirection, tao->solution, tao->XL, tao->XU, tao->gra… in TaoSolve_BLMVM() 106 …PetscCall(VecBoundGradientProjection(blmP->unprojected_gradient, tao->solution, tao->XL, tao->XU, … in TaoSolve_BLMVM()
|
| /petsc/src/tao/complementarity/impls/asls/ |
| H A D | asfls.c | 73 PetscCall(VecFischer(X, tao->constraints, tao->XL, tao->XU, asls->ff)); in Tao_ASLS_FunctionGradient() 78 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, asls->t1, … in Tao_ASLS_FunctionGradient() 125 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSolve_ASFLS() 127 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in TaoSolve_ASFLS() 242 PetscCall(VecBoundGradientProjection(asls->w, tao->solution, tao->XL, tao->XU, asls->w)); in TaoSolve_ASFLS()
|
| H A D | asils.c | 73 PetscCall(VecFischer(X, tao->constraints, tao->XL, tao->XU, asls->ff)); in Tao_ASLS_FunctionGradient() 78 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, asls->t1, … in Tao_ASLS_FunctionGradient()
|
| /petsc/src/tao/bound/impls/bnk/ |
| H A D | bntr.c | 140 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoSolve_BNTR() 187 PetscCall(VecFischer(tao->solution, bnk->unprojected_gradient, tao->XL, tao->XU, bnk->W)); in TaoSolve_BNTR()
|
| H A D | bntl.c | 154 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoSolve_BNTL() 219 PetscCall(VecFischer(tao->solution, bnk->unprojected_gradient, tao->XL, tao->XU, bnk->W)); in TaoSolve_BNTL()
|
| H A D | bnk.c | 55 PetscCall(TaoSetVariableBounds(bnk->bncg, tao->XL, tao->XU)); in TaoBNKInitialize() 56 if (tao->bounded) PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoBNKInitialize() 59 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoBNKInitialize() 69 PetscCall(VecFischer(tao->solution, bnk->unprojected_gradient, tao->XL, tao->XU, bnk->W)); in TaoBNKInitialize() 149 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoBNKInitialize() 239 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoBNKInitialize() 250 PetscCall(VecFischer(tao->solution, bnk->unprojected_gradient, tao->XL, tao->XU, bnk->W)); in TaoBNKInitialize() 305 …PetscCall(VecWhichInactive(tao->XL, tao->solution, bnk->unprojected_gradient, tao->XU, PETSC_TRUE,… in TaoBNKEstimateActiveSet() 332 …PetscCall(TaoEstimateActiveBounds(tao->solution, tao->XL, tao->XU, bnk->unprojected_gradient, bnk-… in TaoBNKEstimateActiveSet() 354 …PetscCall(TaoBoundStep(tao->solution, tao->XL, tao->XU, bnk->active_lower, bnk->active_upper, bnk-… in TaoBNKBoundStep()
|
| H A D | bnls.c | 162 PetscCall(VecFischer(tao->solution, bnk->unprojected_gradient, tao->XL, tao->XU, bnk->W)); in TaoSolve_BNLS()
|
| /petsc/src/tao/bound/impls/bncg/ |
| H A D | bncg.c | 26 …PetscCall(VecWhichInactive(tao->XL, tao->solution, cg->unprojected_gradient, tao->XU, PETSC_TRUE, … in TaoBNCGEstimateActiveSet() 34 …PetscCall(TaoEstimateActiveBounds(tao->solution, tao->XL, tao->XU, cg->unprojected_gradient, cg->W… in TaoBNCGEstimateActiveSet() 52 …PetscCall(TaoBoundStep(tao->solution, tao->XL, tao->XU, cg->active_lower, cg->active_upper, cg->ac… in TaoBNCGBoundStep() 69 PetscCall(TaoLineSearchSetVariableBounds(tao->linesearch, tao->XL, tao->XU)); in TaoSolve_BNCG() 72 PetscCall(TaoBoundSolution(tao->solution, tao->XL, tao->XU, 0.0, &nDiff, tao->solution)); in TaoSolve_BNCG() 97 PetscCall(VecFischer(tao->solution, cg->unprojected_gradient, tao->XL, tao->XU, cg->W)); in TaoSolve_BNCG() 955 PetscCall(VecFischer(tao->solution, cg->unprojected_gradient, tao->XL, tao->XU, cg->W)); in TaoBNCGConductIteration()
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | ipm.c | 247 if (tao->XL) { in IPMInitializeBounds() 249 PetscCall(VecWhichGreaterThan(tao->XL, xtmp, &ipmP->isxl)); in IPMInitializeBounds() 338 PetscCall(VecScatterCreate(tao->XL, isx, ipmP->ci, isc, &ipmP->xl_scat)); in IPMInitializeBounds() 645 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in IPMPushInitialPoint() 765 if (tao->XL) { in IPMUpdateAi() 766 PetscCall(VecAXPY(ipmP->work, -1.0, tao->XL)); in IPMUpdateAi()
|
| H A D | pdipm.c | 63 PetscCall(VecGetArrayRead(tao->XL, &xlarr)); in TaoPDIPMUpdateConstraints() 132 PetscCall(VecRestoreArrayRead(tao->XL, &xlarr)); in TaoPDIPMUpdateConstraints() 160 PetscCall(VecGetLocalSize(tao->XL, &n)); in TaoPDIPMSetUpBounds() 163 PetscCall(VecGetOwnershipRange(tao->XL, &low, &high)); in TaoPDIPMSetUpBounds() 164 PetscCall(VecGetArrayRead(tao->XL, &xl)); in TaoPDIPMSetUpBounds() 180 PetscCall(VecRestoreArrayRead(tao->XL, &xl)); in TaoPDIPMSetUpBounds()
|
| /petsc/include/petsc/private/ |
| H A D | taoimpl.h | 77 Vec XL; member
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2.c | 722 static PetscErrorCode MSA_Plate(Vec XL, Vec XU, PetscCtx ctx) in MSA_Plate() argument 742 PetscCall(VecSet(XL, lb)); in MSA_Plate() 745 PetscCall(VecGetArray(XL, &xl)); in MSA_Plate() 768 PetscCall(VecRestoreArray(XL, &xl)); in MSA_Plate()
|
| /petsc/src/tao/constrained/impls/almm/ |
| H A D | almm.c | 268 PetscCall(TaoALMMCombinePrimal_Private(tao, tao->XL, SL, auglag->PL)); in TaoSetUp_ALMM() 275 PetscCall(VecCopy(tao->XL, auglag->PL)); in TaoSetUp_ALMM() 578 …if (tao->bounded) PetscCall(VecBoundGradientProjection(auglag->LgradX, auglag->Px, tao->XL, tao->X… in TaoALMMComputeOptimalityNorms_Private()
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | pounders.c | 91 if (tao->XU && tao->XL) { in gqtwrap() 95 PetscCall(VecCopy(tao->XL, mfqP->subxl)); in gqtwrap() 471 …if (tao->XU && tao->XL) PetscCall(VecMedian(mfqP->Xhist[mfqP->nHist], tao->XL, tao->XU, mfqP->Xhis… in addpoint() 608 if (tao->XL && tao->XU) { in TaoSolve_POUNDERS() 618 PetscCall(VecCopy(tao->XL, mfqP->Xhist[0])); in TaoSolve_POUNDERS()
|