Lines Matching refs:ce
27 PetscErrorCode PetscConvEstDestroy(PetscConvEst *ce) in PetscConvEstDestroy() argument
30 if (!*ce) PetscFunctionReturn(PETSC_SUCCESS); in PetscConvEstDestroy()
31 PetscValidHeaderSpecific(*ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstDestroy()
32 if (--((PetscObject)*ce)->refct > 0) { in PetscConvEstDestroy()
33 *ce = NULL; in PetscConvEstDestroy()
36 PetscCall(PetscFree3((*ce)->initGuess, (*ce)->exactSol, (*ce)->ctxs)); in PetscConvEstDestroy()
37 PetscCall(PetscFree2((*ce)->dofs, (*ce)->errors)); in PetscConvEstDestroy()
38 PetscCall(PetscHeaderDestroy(ce)); in PetscConvEstDestroy()
54 PetscErrorCode PetscConvEstSetFromOptions(PetscConvEst ce) in PetscConvEstSetFromOptions() argument
57 …PetscOptionsBegin(PetscObjectComm((PetscObject)ce), "", "Convergence Estimator Options", "PetscCon… in PetscConvEstSetFromOptions()
58 …e", "The number of refinements for the convergence check", "PetscConvEst", ce->Nr, &ce->Nr, NULL)); in PetscConvEstSetFromOptions()
59 …ine_factor", "The increase in resolution in each dimension", "PetscConvEst", ce->r, &ce->r, NULL)); in PetscConvEstSetFromOptions()
60 …r", "Monitor the error for each convergence check", "PetscConvEst", ce->monitor, &ce->monitor, NUL… in PetscConvEstSetFromOptions()
61 …"Debugging flag to run on the same mesh each time", "PetscConvEst", ce->noRefine, &ce->noRefine, N… in PetscConvEstSetFromOptions()
79 PetscErrorCode PetscConvEstView(PetscConvEst ce, PetscViewer viewer) in PetscConvEstView() argument
82 PetscCall(PetscObjectPrintClassNamePrefixType((PetscObject)ce, viewer)); in PetscConvEstView()
83 PetscCall(PetscViewerASCIIPrintf(viewer, "ConvEst with %" PetscInt_FMT " levels\n", ce->Nr + 1)); in PetscConvEstView()
102 PetscErrorCode PetscConvEstGetSolver(PetscConvEst ce, PetscObject *solver) in PetscConvEstGetSolver() argument
105 PetscValidHeaderSpecific(ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstGetSolver()
107 *solver = ce->solver; in PetscConvEstGetSolver()
124 PetscErrorCode PetscConvEstSetSolver(PetscConvEst ce, PetscObject solver) in PetscConvEstSetSolver() argument
127 PetscValidHeaderSpecific(ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstSetSolver()
129 ce->solver = solver; in PetscConvEstSetSolver()
130 PetscUseTypeMethod(ce, setsolver, solver); in PetscConvEstSetSolver()
146 PetscErrorCode PetscConvEstSetUp(PetscConvEst ce) in PetscConvEstSetUp() argument
151 PetscCall(DMGetNumFields(ce->idm, &Nf)); in PetscConvEstSetUp()
152 ce->Nf = PetscMax(Nf, 1); in PetscConvEstSetUp()
153 PetscCall(PetscMalloc2((ce->Nr + 1) * ce->Nf, &ce->dofs, (ce->Nr + 1) * ce->Nf, &ce->errors)); in PetscConvEstSetUp()
154 PetscCall(PetscCalloc3(ce->Nf, &ce->initGuess, ce->Nf, &ce->exactSol, ce->Nf, &ce->ctxs)); in PetscConvEstSetUp()
155 for (f = 0; f < Nf; ++f) ce->initGuess[f] = zero_private; in PetscConvEstSetUp()
156 PetscCall(DMGetNumDS(ce->idm, &Nds)); in PetscConvEstSetUp()
164 PetscCall(DMGetRegionNumDS(ce->idm, s, &label, &fieldIS, &ds, NULL)); in PetscConvEstSetUp()
169 PetscCall(PetscDSGetExactSolution(ds, field, &ce->exactSol[field], &ce->ctxs[field])); in PetscConvEstSetUp()
173 …for (f = 0; f < Nf; ++f) PetscCheck(ce->exactSol[f], PetscObjectComm((PetscObject)ce), PETSC_ERR_A… in PetscConvEstSetUp()
177 PetscErrorCode PetscConvEstComputeInitialGuess(PetscConvEst ce, PetscInt r, DM dm, Vec u) in PetscConvEstComputeInitialGuess() argument
180 PetscValidHeaderSpecific(ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstComputeInitialGuess()
183 PetscUseTypeMethod(ce, initguess, r, dm, u); in PetscConvEstComputeInitialGuess()
187 PetscErrorCode PetscConvEstComputeError(PetscConvEst ce, PetscInt r, DM dm, Vec u, PetscReal errors… in PetscConvEstComputeError() argument
190 PetscValidHeaderSpecific(ce, PETSC_OBJECT_CLASSID, 1); in PetscConvEstComputeError()
194 PetscUseTypeMethod(ce, computeerror, r, dm, u, errors); in PetscConvEstComputeError()
214 PetscErrorCode PetscConvEstMonitorDefault(PetscConvEst ce, PetscInt r) in PetscConvEstMonitorDefault() argument
220 if (ce->monitor) { in PetscConvEstMonitorDefault()
221 PetscInt *dofs = &ce->dofs[r * ce->Nf]; in PetscConvEstMonitorDefault()
222 PetscReal *errors = &ce->errors[r * ce->Nf]; in PetscConvEstMonitorDefault()
224 PetscCall(PetscObjectGetComm((PetscObject)ce, &comm)); in PetscConvEstMonitorDefault()
226 if (ce->Nf > 1) PetscCall(PetscPrintf(comm, "[")); in PetscConvEstMonitorDefault()
227 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstMonitorDefault()
231 if (ce->Nf > 1) PetscCall(PetscPrintf(comm, "]")); in PetscConvEstMonitorDefault()
234 if (ce->Nf > 1) PetscCall(PetscPrintf(comm, "[")); in PetscConvEstMonitorDefault()
235 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstMonitorDefault()
240 if (ce->Nf > 1) PetscCall(PetscPrintf(comm, "]")); in PetscConvEstMonitorDefault()
246 static PetscErrorCode PetscConvEstSetSNES_Private(PetscConvEst ce, PetscObject solver) in PetscConvEstSetSNES_Private() argument
251 PetscCall(PetscObjectGetClassId(ce->solver, &id)); in PetscConvEstSetSNES_Private()
252 …PetscCheck(id == SNES_CLASSID, PetscObjectComm((PetscObject)ce), PETSC_ERR_ARG_WRONG, "Solver was … in PetscConvEstSetSNES_Private()
253 PetscCall(SNESGetDM((SNES)ce->solver, &ce->idm)); in PetscConvEstSetSNES_Private()
257 static PetscErrorCode PetscConvEstInitGuessSNES_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u) in PetscConvEstInitGuessSNES_Private() argument
260 PetscCall(DMProjectFunction(dm, 0.0, ce->initGuess, ce->ctxs, INSERT_VALUES, u)); in PetscConvEstInitGuessSNES_Private()
264 static PetscErrorCode PetscConvEstComputeErrorSNES_Private(PetscConvEst ce, PetscInt r, DM dm, Vec … in PetscConvEstComputeErrorSNES_Private() argument
270 PetscCall(DMComputeL2FieldDiff(dm, 0.0, ce->exactSol, ce->ctxs, u, errors)); in PetscConvEstComputeErrorSNES_Private()
271 PetscCall(PetscObjectGetOptionsPrefix((PetscObject)ce, &prefix)); in PetscConvEstComputeErrorSNES_Private()
307 static PetscErrorCode PetscConvEstSetJacobianNullSpace_Private(PetscConvEst ce, SNES snes) in PetscConvEstSetJacobianNullSpace_Private() argument
314 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstSetJacobianNullSpace_Private()
333 static PetscErrorCode PetscConvEstGetConvRateSNES_Private(PetscConvEst ce, PetscReal alpha[]) in PetscConvEstGetConvRateSNES_Private() argument
335 SNES snes = (SNES)ce->solver; in PetscConvEstGetConvRateSNES_Private()
339 PetscInt Nr = ce->Nr, r, f, dim, oldlevel, oldnlev; in PetscConvEstGetConvRateSNES_Private()
343 …etscCheck(ce->r == 2.0, PetscObjectComm((PetscObject)ce), PETSC_ERR_SUP, "Only refinement factor 2… in PetscConvEstGetConvRateSNES_Private()
344 PetscCall(DMGetDimension(ce->idm, &dim)); in PetscConvEstGetConvRateSNES_Private()
345 PetscCall(DMGetApplicationContext(ce->idm, &ctx)); in PetscConvEstGetConvRateSNES_Private()
346 PetscCall(DMPlexSetRefinementUniform(ce->idm, PETSC_TRUE)); in PetscConvEstGetConvRateSNES_Private()
347 PetscCall(DMGetRefineLevel(ce->idm, &oldlevel)); in PetscConvEstGetConvRateSNES_Private()
350 dm[0] = ce->idm; in PetscConvEstGetConvRateSNES_Private()
362 if (!ce->noRefine) { in PetscConvEstGetConvRateSNES_Private()
374 PetscCall(DMCopyTransform(ce->idm, dm[r])); in PetscConvEstGetConvRateSNES_Private()
377 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstGetConvRateSNES_Private()
397 PetscCall(PetscConvEstSetJacobianNullSpace_Private(ce, snes)); in PetscConvEstGetConvRateSNES_Private()
399 PetscCall(PetscConvEstComputeInitialGuess(ce, r, dm[r], u)); in PetscConvEstGetConvRateSNES_Private()
401 PetscCall(PetscLogEventBegin(ce->event, ce, 0, 0, 0)); in PetscConvEstGetConvRateSNES_Private()
402 PetscCall(PetscConvEstComputeError(ce, r, dm[r], u, &ce->errors[r * ce->Nf])); in PetscConvEstGetConvRateSNES_Private()
403 PetscCall(PetscLogEventEnd(ce->event, ce, 0, 0, 0)); in PetscConvEstGetConvRateSNES_Private()
404 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstGetConvRateSNES_Private()
412 …PetscCallMPI(MPIU_Allreduce(&lsize, &ce->dofs[r * ce->Nf + f], 1, MPIU_INT, MPI_SUM, PetscObjectCo… in PetscConvEstGetConvRateSNES_Private()
413 PetscCall(PetscLogEventSetDof(ce->event, f, ce->dofs[r * ce->Nf + f])); in PetscConvEstGetConvRateSNES_Private()
414 PetscCall(PetscLogEventSetError(ce->event, f, ce->errors[r * ce->Nf + f])); in PetscConvEstGetConvRateSNES_Private()
417 PetscCall(PetscConvEstMonitorDefault(ce, r)); in PetscConvEstGetConvRateSNES_Private()
434 for (f = 0; f < ce->Nf; ++f) { in PetscConvEstGetConvRateSNES_Private()
436 x[r] = PetscLog10Real(ce->dofs[r * ce->Nf + f]); in PetscConvEstGetConvRateSNES_Private()
437 y[r] = PetscLog10Real(ce->errors[r * ce->Nf + f]); in PetscConvEstGetConvRateSNES_Private()
455 …PetscCall(DMSetRefineLevel(ce->idm, oldlevel)); /* The damn DMCoarsen() calls in PCMG can reset th… in PetscConvEstGetConvRateSNES_Private()
457 PetscCall(SNESSetDM(snes, ce->idm)); in PetscConvEstGetConvRateSNES_Private()
458 PetscCall(DMPlexSetSNESLocalFEM(ce->idm, PETSC_FALSE, ctx)); in PetscConvEstGetConvRateSNES_Private()
459 PetscCall(DMPlexSetSNESVariableBounds(ce->idm, snes)); in PetscConvEstGetConvRateSNES_Private()
461 PetscCall(PetscConvEstSetJacobianNullSpace_Private(ce, snes)); in PetscConvEstGetConvRateSNES_Private()
497 PetscErrorCode PetscConvEstGetConvRate(PetscConvEst ce, PetscReal alpha[]) in PetscConvEstGetConvRate() argument
502 …if (ce->event < 0) PetscCall(PetscLogEventRegister("ConvEst Error", PETSC_OBJECT_CLASSID, &ce->eve… in PetscConvEstGetConvRate()
503 for (f = 0; f < ce->Nf; ++f) alpha[f] = 0.0; in PetscConvEstGetConvRate()
504 PetscUseTypeMethod(ce, getconvrate, alpha); in PetscConvEstGetConvRate()
525 PetscErrorCode PetscConvEstRateView(PetscConvEst ce, const PetscReal alpha[], PetscViewer viewer) in PetscConvEstRateView() argument
532 PetscInt Nf = ce->Nf, f; in PetscConvEstRateView()
534 PetscCall(PetscViewerASCIIAddTab(viewer, ((PetscObject)ce)->tablevel)); in PetscConvEstRateView()
543 PetscCall(PetscViewerASCIISubtractTab(viewer, ((PetscObject)ce)->tablevel)); in PetscConvEstRateView()
563 PetscErrorCode PetscConvEstCreate(MPI_Comm comm, PetscConvEst *ce) in PetscConvEstCreate() argument
566 PetscAssertPointer(ce, 2); in PetscConvEstCreate()
568 …PetscCall(PetscHeaderCreate(*ce, PETSC_OBJECT_CLASSID, "PetscConvEst", "ConvergenceEstimator", "SN… in PetscConvEstCreate()
569 (*ce)->monitor = PETSC_FALSE; in PetscConvEstCreate()
570 (*ce)->r = 2.0; in PetscConvEstCreate()
571 (*ce)->Nr = 4; in PetscConvEstCreate()
572 (*ce)->event = -1; in PetscConvEstCreate()
573 (*ce)->ops->setsolver = PetscConvEstSetSNES_Private; in PetscConvEstCreate()
574 (*ce)->ops->initguess = PetscConvEstInitGuessSNES_Private; in PetscConvEstCreate()
575 (*ce)->ops->computeerror = PetscConvEstComputeErrorSNES_Private; in PetscConvEstCreate()
576 (*ce)->ops->getconvrate = PetscConvEstGetConvRateSNES_Private; in PetscConvEstCreate()