| /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 198 if (XU) PetscValidHeaderSpecific(XU, VEC_CLASSID, 3); in TaoEstimateActiveBounds() 204 if (XU) PetscCheckSameType(X, 1, XU, 3); in TaoEstimateActiveBounds() 209 if (XU) PetscCheckSameComm(X, 1, XU, 3); in TaoEstimateActiveBounds() 214 if (XU) VecCheckSameSize(X, 1, XU, 3); in TaoEstimateActiveBounds() 222 PetscCall(TaoBoundSolution(W, XL, XU, 0.0, &nDiff, W)); in TaoEstimateActiveBounds() 236 if (!XL && !XU) { in TaoEstimateActiveBounds() 243 PetscCall(VecGetArrayRead(XU, &xu)); in TaoEstimateActiveBounds() 280 PetscCall(VecRestoreArrayRead(XU, &xu)); 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 22 if (XU) PetscValidHeaderSpecific(XU, VEC_CLASSID, 3); in TaoSetVariableBounds() 24 PetscCall(PetscObjectReference((PetscObject)XU)); in TaoSetVariableBounds() 26 PetscCall(VecDestroy(&tao->XU)); in TaoSetVariableBounds() 28 tao->XU = XU; in TaoSetVariableBounds() 29 tao->bounded = (PetscBool)(XL || XU); in TaoSetVariableBounds() 82 PetscErrorCode TaoGetVariableBounds(Tao tao, Vec *XL, Vec *XU) in TaoGetVariableBounds() argument 87 if (XU) *XU = tao->XU; in TaoGetVariableBounds() 113 if (!tao->XU) { in TaoComputeVariableBounds() 114 PetscCall(VecDuplicate(tao->solution, &tao->XU)); 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 766 if (XU) PetscValidHeaderSpecific(XU, VEC_CLASSID, 4); in VecBoundGradientProjection() 768 if (!XL && !XU) { in VecBoundGradientProjection() 777 PetscCall(VecGetArrayRead(XU, &xuptr)); in VecBoundGradientProjection() 793 PetscCall(VecRestoreArrayRead(XU, &xuptr)); in VecBoundGradientProjection() 816 PetscErrorCode VecStepMaxBounded(Vec X, Vec DX, Vec XL, Vec XU, PetscReal *stepmax) in VecStepMaxBounded() argument 826 PetscValidHeaderSpecific(XU, VEC_CLASSID, 4); in VecStepMaxBounded() 830 PetscCall(VecGetArrayRead(XU, &xu)); in VecStepMaxBounded() 842 PetscCall(VecRestoreArrayRead(XU, &xu)); 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 | 52 PetscCall(VecCopy(qp->XU, qp->T)); in QPIPSetInitialPoint() 104 PetscCall(VecAXPY(qp->R5, -1.0, qp->XU)); in QPIPSetInitialPoint() 212 PetscCall(VecDuplicate(tao->solution, &qp->XU)); in TaoSetup_BQPIP() 267 PetscCall(VecSet(qp->XU, 1.0e20)); in TaoSolve_BQPIP() 270 if (tao->XU) PetscCall(VecPointwiseMin(qp->XU, qp->XU, tao->XU)); in TaoSolve_BQPIP() 271 PetscCall(VecMedian(qp->XL, tao->solution, qp->XU, tao->solution)); in TaoSolve_BQPIP() 489 PetscCall(VecDestroy(&qp->XU)); 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 259 PetscCall(VecGetArrayRead(XU, &u)); in MatDFischer() 363 PetscCall(VecRestoreArrayRead(XU, &u)); 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() 409 PetscCall(VecGetArrayRead(XU, &u)); in MatDSFischer() 465 PetscCall(VecRestoreArrayRead(XU, &u)); in MatDSFischer()
|
| /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/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 …undGradientProjection(blmP->unprojected_gradient, tao->solution, tao->XL, tao->XU, tao->gradient)); in TaoSolve_BLMVM()
|
| /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 | 254 if (tao->XU) { in IPMInitializeBounds() 256 PetscCall(VecWhichLessThan(tao->XU, xtmp, &ipmP->isxu)); in IPMInitializeBounds() 365 PetscCall(VecScatterCreate(tao->XU, isx, ipmP->ci, isc, &ipmP->xu_scat)); in IPMInitializeBounds() 645 PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution)); in IPMPushInitialPoint() 774 if (tao->XU) { in IPMUpdateAi() 778 PetscCall(VecAXPY(ipmP->work, 1.0, tao->XU)); in IPMUpdateAi()
|
| /petsc/include/petsc/private/ |
| H A D | taoimpl.h | 78 Vec XU; member
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | pounders.c | 91 if (tao->XU && tao->XL) { in gqtwrap() 92 PetscCall(VecCopy(tao->XU, mfqP->subxu)); 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() 613 PetscCall(VecAXPY(mfqP->Xhist[0], -1.0, tao->XU)); in TaoSolve_POUNDERS() 627 PetscCall(VecAXPY(mfqP->Xhist[0], -1.0, tao->XU)); in TaoSolve_POUNDERS()
|
| /petsc/src/tao/constrained/impls/almm/ |
| H A D | almm.c | 269 PetscCall(TaoALMMCombinePrimal_Private(tao, tao->XU, SU, auglag->PU)); in TaoSetUp_ALMM() 276 PetscCall(VecCopy(tao->XU, auglag->PU)); in TaoSetUp_ALMM() 578 …etscCall(VecBoundGradientProjection(auglag->LgradX, auglag->Px, tao->XL, tao->XU, auglag->LgradX)); in TaoALMMComputeOptimalityNorms_Private()
|
| /petsc/src/tao/unconstrained/impls/lmvm/ |
| H A D | lmvm.c | 16 …if (tao->XL || tao->XU || tao->ops->computebounds) PetscCall(PetscInfo(tao, "WARNING: Variable bou… in TaoSolve_LMVM()
|
| /petsc/src/tao/unconstrained/impls/cg/ |
| H A D | taocg.c | 21 …if (tao->XL || tao->XU || tao->ops->computebounds) PetscCall(PetscInfo(tao, "WARNING: Variable bou… in TaoSolve_CG()
|