Home
last modified time | relevance | path

Searched refs:xu (Results 1 – 25 of 43) sorted by relevance

12

/petsc/src/snes/impls/vi/
H A Dvi.c100 const PetscScalar *x, *xl, *xu, *f; in SNESMonitorVI() local
112 PetscCall(VecGetArrayRead(snes->xu, &xu)); in SNESMonitorVI()
118 … || (PetscRealPart(f[i]) <= 0.0)) && ((PetscRealPart(x[i]) < PetscRealPart(xu[i]) - zerotolerance)… in SNESMonitorVI()
120 …else if (PetscRealPart(x[i]) >= PetscRealPart(xu[i]) - zerotolerance && PetscRealPart(f[i]) < 0.0)… in SNESMonitorVI()
126 else if (PetscRealPart(x[i]) >= PetscRealPart(xu[i]) - zerotolerance) act_bound[1]++; in SNESMonitorVI()
130 PetscCall(VecRestoreArrayRead(snes->xu, &xu)); in SNESMonitorVI()
237 const PetscScalar *xl, *xu; in SNESVIProjectOntoBounds() local
245 PetscCall(VecGetArrayRead(snes->xu, &xu)); in SNESVIProjectOntoBounds()
249 else if (PetscRealPart(x[i]) > PetscRealPart(xu[i])) x[i] = xu[i]; in SNESVIProjectOntoBounds()
253 PetscCall(VecRestoreArrayRead(snes->xu, &xu)); in SNESVIProjectOntoBounds()
[all …]
/petsc/src/binding/petsc4py/test/
H A Dtest_tao.py194 xu = PETSc.Vec().create(tao.getComm())
195 xu.setType('standard')
196 xu.setSizes(2)
197 xu.set(2.0)
198 tao.setVariableBounds((xl, xu))
223 xu = PETSc.Vec().create(tao.getComm())
224 xu.setType('standard')
225 xu.setSizes(2)
226 xu.set(2.0)
227 tao.setVariableBounds((xl, xu))
/petsc/src/tao/bound/utils/
H A Disutil.c192 const PetscScalar *xl, *xu, *x, *g; in TaoEstimateActiveBounds() local
243 PetscCall(VecGetArrayRead(XU, &xu)); in TaoEstimateActiveBounds()
253 if (xl[i] == xu[i]) { in TaoEstimateActiveBounds()
265 } else if (xu[i] < PETSC_INFINITY && x[i] >= xu[i] - *bound_tol && g[i] < zero) { in TaoEstimateActiveBounds()
280 PetscCall(VecRestoreArrayRead(XU, &xu)); in TaoEstimateActiveBounds()
409 const PetscScalar *x, *xl, *xu; in TaoBoundSolution() local
436 PetscCall(VecGetArrayRead(XU, &xu)); in TaoBoundSolution()
443 } else if (xu[i] < PETSC_INFINITY && x[i] >= xu[i] - bound_tol) { in TaoBoundSolution()
444 xout[i] = xu[i]; in TaoBoundSolution()
451 PetscCall(VecRestoreArrayRead(XU, &xu)); in TaoBoundSolution()
/petsc/src/ts/tutorials/
H A Dex21.c79 Vec xl, xu; /* Lower and upper bounds on variables */ in main() local
133 PetscCall(VecDuplicate(u, &xu)); in main()
134 PetscCall(SetBounds(xl, xu, ul, uh, &appctx)); in main()
186 PetscCall(TSVISetVariableBounds(ts, xl, xu)); in main()
210 PetscCall(VecDestroy(&xu)); in main()
298 PetscErrorCode SetBounds(Vec xl, Vec xu, PetscScalar ul, PetscScalar uh, AppCtx *appctx) in SetBounds() argument
306 PetscCall(VecSet(xu, uh)); in SetBounds()
309 PetscCall(VecGetArray(xu, &u)); in SetBounds()
322 PetscCall(VecRestoreArray(xu, &u)); in SetBounds()
/petsc/src/tao/complementarity/tutorials/
H A Dblackscholes.c258 PetscErrorCode ComputeVariableBounds(Tao tao, Vec xl, Vec xu, PetscCtx ctx) in ComputeVariableBounds() argument
268 PetscCall(VecSet(xu, ub)); in ComputeVariableBounds()
279 PetscCall(VecGetArray(xu, &xl_array)); in ComputeVariableBounds()
281 PetscCall(VecRestoreArray(xu, &xl_array)); in ComputeVariableBounds()
284 PetscCall(VecGetArray(xu, &xl_array)); in ComputeVariableBounds()
286 PetscCall(VecRestoreArray(xu, &xl_array)); in ComputeVariableBounds()
H A Dminsurf1.c40 Vec xl, xu; /* Bounds on the variables */ in main() local
93 PetscCall(VecDuplicate(x, &xu)); in main()
95 PetscCall(VecSet(xu, ub)); in main()
96 PetscCall(TaoSetVariableBounds(tao, xl, xu)); in main()
110 PetscCall(VecDestroy(&xu)); in main()
/petsc/src/snes/impls/composite/
H A Dsnescomposite.c90 if (snes->xl && snes->xu) { in SNESCompositeApply_Multiplicative()
100 if (snes->xl && snes->xu) { in SNESCompositeApply_Multiplicative()
162 if (snes->xl && snes->xu) { in SNESCompositeApply_Additive()
282 if (snes->xl && snes->xu) { in SNESCompositeApply_AdditiveOptimal()
322 if (!snes->xu) PetscCall(VecDuplicate(snes->vec_sol, &snes->xu)); in SNESSetUp_Composite()
323 PetscUseTypeMethod(snes, computevariablebounds, snes->xl, snes->xu); in SNESSetUp_Composite()
331 if (snes->xl && snes->xu) { in SNESSetUp_Composite()
335 PetscCall(SNESVISetVariableBounds(next->snes, snes->xl, snes->xu)); in SNESSetUp_Composite()
713 if (snes->xl && snes->xu) { in SNESSolve_Composite()
758 if (snes->xl && snes->xu) { in SNESSolve_Composite()
/petsc/src/dm/impls/stag/tutorials/
H A Dex1.c294 PetscScalar val[3], xu[2]; in main() local
308 PetscCall(DMStagVecGetValuesStencil(dmCoordSol, coordSolLocal, 2, col, xu)); in main()
309 h = xu[0] - xu[1]; in main()
/petsc/src/vec/vec/utils/
H A Dprojection.c819 const PetscScalar *xx, *dx, *xl, *xu; in VecStepMaxBounded() local
830 PetscCall(VecGetArrayRead(XU, &xu)); in VecStepMaxBounded()
835 localmax = PetscMax(localmax, PetscRealPart((xu[i] - xx[i]) / dx[i])); in VecStepMaxBounded()
842 PetscCall(VecRestoreArrayRead(XU, &xu)); in VecStepMaxBounded()
874 const PetscScalar *x, *xl, *xu, *dx; in VecStepBoundInfo() local
887 PetscCall(VecGetArrayRead(XU, &xu)); in VecStepBoundInfo()
891 if (PetscRealPart(dx[i]) > 0 && PetscRealPart(xu[i]) < PETSC_INFINITY) { in VecStepBoundInfo()
892 t = PetscRealPart((xu[i] - x[i]) / dx[i]); in VecStepBoundInfo()
906 PetscCall(VecRestoreArrayRead(XU, &xu)); in VecStepBoundInfo()
/petsc/src/tao/constrained/tutorials/
H A Dex1.c40 Vec x, xl, xu; member
78 if (!user.noboundflag) PetscCall(TaoSetVariableBounds(tao, user.xl, user.xu)); in main()
189 PetscCall(VecDuplicate(user->x, &user->xu)); in InitializeProblem()
191 PetscCall(VecSet(user->xu, 2.0)); in InitializeProblem()
244 PetscCall(VecDestroy(&user->xu)); in DestroyProblem()
/petsc/src/tao/bound/tutorials/
H A Djbearing2.c53 Vec xl, xu; /* bounds vectors */ in main() local
133 PetscCall(VecDuplicate(x, &xu)); in main()
135 PetscCall(VecSet(xu, d1000)); in main()
136 PetscCall(TaoSetVariableBounds(tao, xl, xu)); in main()
155 PetscCall(VecDestroy(&xu)); in main()
H A Dplate2.c53 Vec x, xl, xu; /* solution vector and bounds*/ in main() local
108 PetscCall(VecDuplicate(x, &xu)); in main()
145 PetscCall(MSA_Plate(xl, xu, (void *)&user)); in main()
146 PetscCall(TaoSetVariableBounds(tao, xl, xu)); in main()
162 PetscCall(VecDestroy(&xu)); in main()
H A Dplate2f.F90664 subroutine MSA_Plate(ta, xl, xu, dummy, ierr) argument
667 Vec xl, xu
686 PetscCall(VecSet(xu, ub, ierr))
/petsc/src/tao/linesearch/interface/
H A Dtaolinesearch.c1090 PetscErrorCode TaoLineSearchSetVariableBounds(TaoLineSearch ls, Vec xl, Vec xu) in TaoLineSearchSetVariableBounds() argument
1095 if (xu) PetscValidHeaderSpecific(xu, VEC_CLASSID, 3); in TaoLineSearchSetVariableBounds()
1097 PetscCall(PetscObjectReference((PetscObject)xu)); in TaoLineSearchSetVariableBounds()
1101 ls->upper = xu; in TaoLineSearchSetVariableBounds()
1102 ls->bounded = (PetscBool)(xl || xu); in TaoLineSearchSetVariableBounds()
/petsc/src/snes/impls/vi/ss/
H A Dviss.c56 Vec Xl = snes->xl, Xu = snes->xu, F = snes->vec_func; in SNESVIComputeFunction()
105 PetscCall(VecGetArray(snes->xu, &u)); in SNESVIComputeBsubdifferentialVectors()
136 PetscCall(VecRestoreArray(snes->xu, &u)); in SNESVIComputeBsubdifferentialVectors()
/petsc/src/snes/impls/ksponly/
H A Dksponly.c15 …PetscCheck(!snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObje… in SNESSolve_KSPONLY()
/petsc/src/dm/impls/stag/tests/
H A Dex18.c84 Vec su, xu, fu; in main() local
155 PetscCall(MatCreateVecs(Auu, &xu, &fu)); in main()
165 PetscCall(MatCreateVecs(Auu, &xu, NULL)); in main()
216 PetscCall(KSPSolve(ksp, fu, xu)); in main()
235 PetscCall(VecDestroy(&xu)); in main()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DTAO.pyx517 cdef Vec xl = None, xu = None
520 xl = <Vec?> ol; xu = <Vec?> ou
521 CHKERR(TaoSetVariableBounds(self.tao, xl.vec, xu.vec))
525 xl = <Vec?> ol; xu = <Vec?> ou
526 CHKERR(TaoSetVariableBounds(self.tao, xl.vec, xu.vec))
1004 def computeDualVariables(self, Vec xl, Vec xu) -> None: argument
1014 CHKERR(TaoComputeDualVariables(self.tao, xl.vec, xu.vec))
1016 def computeVariableBounds(self, Vec xl, Vec xu) -> None: argument
1034 if xu.vec != NULL:
1036 CHKERR(VecCopy(Uvec, xu.vec))
[all …]
/petsc/src/snes/impls/richardson/
H A Dsnesrichardson.c34 …PetscCheck(!snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObje… in SNESSolve_NRichardson()
/petsc/src/tao/constrained/impls/ipm/
H A Dpdipm.c151 const PetscScalar *xl, *xu; in TaoPDIPMSetUpBounds() local
165 PetscCall(VecGetArrayRead(tao->XU, &xu)); in TaoPDIPMSetUpBounds()
168 if ((PetscRealPart(xl[i]) > PETSC_NINFINITY) && (PetscRealPart(xu[i]) < PETSC_INFINITY)) { in TaoPDIPMSetUpBounds()
169 if (PetscRealPart(xl[i]) == PetscRealPart(xu[i])) { in TaoPDIPMSetUpBounds()
173 if ((PetscRealPart(xl[i]) > PETSC_NINFINITY) && (PetscRealPart(xu[i]) >= PETSC_INFINITY)) { in TaoPDIPMSetUpBounds()
175 … } else if ((PetscRealPart(xl[i]) <= PETSC_NINFINITY) && (PetscRealPart(xu[i]) < PETSC_INFINITY)) { in TaoPDIPMSetUpBounds()
181 PetscCall(VecRestoreArrayRead(tao->XU, &xu)); in TaoPDIPMSetUpBounds()
/petsc/src/dm/impls/da/
H A Ddalocal.c335 PetscErrorCode DMDASetVertexCoordinates(DM dm, PetscReal xl, PetscReal xu, PetscReal yl, PetscReal … in DMDASetVertexCoordinates() argument
348 h[0] = (xu - xl) / M; in DMDASetVertexCoordinates()
/petsc/src/tao/interface/ftn-custom/
H A Dztaosolverf.c107 static PetscErrorCode ourtaoboundsroutine(Tao tao, Vec xl, Vec xu, PetscCtx ctx) in ourtaoboundsroutine() argument
109 …k(tao, _cb.bounds, (Tao *, Vec *, Vec *, void *, PetscErrorCode *), (&tao, &xl, &xu, _ctx, &ierr)); in ourtaoboundsroutine()
/petsc/src/snes/impls/gs/
H A Dsnesgs.c223 …PetscCheck(!snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObje… in SNESSolve_NGS()
/petsc/src/snes/impls/ngmres/
H A Danderson.c39 …PetscCheck(!snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObje… in SNESSolve_Anderson()
/petsc/src/snes/tutorials/
H A Dex58.c120 PetscErrorCode FormBounds(SNES snes, Vec xl, Vec xu) in FormBounds() argument
127 PetscCall(VecSet(xu, ctx->ub)); in FormBounds()

12