Lines Matching refs:xu

100   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()
274 Vec Xl = snes->xl, Xu = snes->xu; in SNESVIGetActiveSetIS()
275 const PetscScalar *x, *f, *xl, *xu; in SNESVIGetActiveSetIS() local
284 PetscCall(VecGetArrayRead(Xu, &xu)); in SNESVIGetActiveSetIS()
288 … || (PetscRealPart(f[i]) <= 0.0)) && ((PetscRealPart(x[i]) < PetscRealPart(xu[i]) - zerotolerance)… in SNESVIGetActiveSetIS()
295 … || (PetscRealPart(f[i]) <= 0.0)) && ((PetscRealPart(x[i]) < PetscRealPart(xu[i]) - zerotolerance)… in SNESVIGetActiveSetIS()
303 PetscCall(VecRestoreArrayRead(Xu, &xu)); in SNESVIGetActiveSetIS()
325 const PetscScalar *x, *xl, *xu, *f; in SNESVIComputeInactiveSetFnorm() local
334 PetscCall(VecGetArrayRead(snes->xu, &xu)); in SNESVIComputeInactiveSetFnorm()
339 … || (PetscRealPart(f[i]) <= 0.0)) && ((PetscRealPart(x[i]) < PetscRealPart(xu[i]) - zerotolerance)… in SNESVIComputeInactiveSetFnorm()
343 PetscCall(VecRestoreArrayRead(snes->xu, &xu)); in SNESVIComputeInactiveSetFnorm()
369 const PetscScalar *x, *xl, *xu, *y, *f; in SNESVIComputeInactiveSetFtY() local
380 PetscCall(VecGetArrayRead(snes->xu, &xu)); in SNESVIComputeInactiveSetFtY()
384 … || (PetscRealPart(f[i]) <= 0.0)) && ((PetscRealPart(x[i]) < PetscRealPart(xu[i]) - zerotolerance)… in SNESVIComputeInactiveSetFtY()
389 PetscCall(VecRestoreArrayRead(snes->xu, &xu)); in SNESVIComputeInactiveSetFtY()
395 static PetscErrorCode SNESVIDMComputeVariableBounds(SNES snes, Vec xl, Vec xu) in SNESVIDMComputeVariableBounds() argument
398 PetscCall(DMComputeVariableBounds(snes->dm, xl, xu)); in SNESVIDMComputeVariableBounds()
432 if (!snes->xu) PetscCall(VecDuplicate(snes->work[0], &snes->xu)); in SNESSetUp_VI()
433 PetscUseTypeMethod(snes, computevariablebounds, snes->xl, snes->xu); in SNESSetUp_VI()
434 } else if (!snes->xl && !snes->xu) { in SNESSetUp_VI()
438 PetscCall(VecDuplicate(snes->work[0], &snes->xu)); in SNESSetUp_VI()
439 PetscCall(VecSet(snes->xu, PETSC_INFINITY)); in SNESSetUp_VI()
444 PetscCall(VecGetOwnershipRange(snes->xu, i_start + 2, i_end + 2)); in SNESSetUp_VI()
455 PetscCall(VecDestroy(&snes->xu)); in SNESReset_VI()
504 PetscErrorCode SNESVISetVariableBounds(SNES snes, Vec xl, Vec xu) in SNESVISetVariableBounds() argument
511 PetscValidHeaderSpecific(xu, VEC_CLASSID, 3); in SNESVISetVariableBounds()
513 if (f) PetscUseMethod(snes, "SNESVISetVariableBounds_C", (SNES, Vec, Vec), (snes, xl, xu)); in SNESVISetVariableBounds()
514 else PetscCall(SNESVISetVariableBounds_VI(snes, xl, xu)); in SNESVISetVariableBounds()
519 PetscErrorCode SNESVISetVariableBounds_VI(SNES snes, Vec xl, Vec xu) in SNESVISetVariableBounds_VI() argument
530 PetscCall(VecGetSize(xu, &xuN)); in SNESVISetVariableBounds_VI()
536 PetscCall(PetscObjectReference((PetscObject)xu)); in SNESVISetVariableBounds_VI()
538 PetscCall(VecDestroy(&snes->xu)); in SNESVISetVariableBounds_VI()
540 snes->xu = xu; in SNESVISetVariableBounds_VI()
543 PetscCall(VecGetArrayRead(xu, &xxu)); in SNESVISetVariableBounds_VI()
548 PetscCall(VecRestoreArrayRead(xu, &xxu)); in SNESVISetVariableBounds_VI()
568 PetscErrorCode SNESVIGetVariableBounds(SNES snes, Vec *xl, Vec *xu) in SNESVIGetVariableBounds() argument
573 if (xu) *xu = snes->xu; in SNESVIGetVariableBounds()