Lines Matching full:r
18 static PetscErrorCode PetscConvEstInitGuessTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u) in PetscConvEstInitGuessTS_Private() argument
25 static PetscErrorCode PetscConvEstComputeErrorTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u,… in PetscConvEstComputeErrorTS_Private() argument
55 PetscInt Ns, oNs, Nf = ce->Nf, f, Nr = ce->Nr, r; in PetscConvEstGetConvRateTS_Temporal_Private() local
64 for (r = 0; r <= Nr; ++r) { in PetscConvEstGetConvRateTS_Temporal_Private()
65 if (r > 0) { in PetscConvEstGetConvRateTS_Temporal_Private()
66 dt[r] = dt[r - 1] / ce->r; in PetscConvEstGetConvRateTS_Temporal_Private()
67 Ns = PetscCeilReal(Ns * ce->r); in PetscConvEstGetConvRateTS_Temporal_Private()
71 PetscCall(TSSetTimeStep(ts, dt[r])); in PetscConvEstGetConvRateTS_Temporal_Private()
74 PetscCall(PetscConvEstComputeInitialGuess(ce, r, NULL, u)); in PetscConvEstGetConvRateTS_Temporal_Private()
78 PetscCall(PetscConvEstComputeError(ce, r, ce->idm, u, &ce->errors[r * Nf])); in PetscConvEstGetConvRateTS_Temporal_Private()
81 ce->dofs[r * Nf + f] = 1.0 / dt[r]; in PetscConvEstGetConvRateTS_Temporal_Private()
82 PetscCall(PetscLogEventSetDof(ce->event, f, ce->dofs[r * Nf + f])); in PetscConvEstGetConvRateTS_Temporal_Private()
83 PetscCall(PetscLogEventSetError(ce->event, f, ce->errors[r * Nf + f])); in PetscConvEstGetConvRateTS_Temporal_Private()
86 PetscCall(PetscConvEstMonitorDefault(ce, r)); in PetscConvEstGetConvRateTS_Temporal_Private()
92 for (r = 0; r <= Nr; ++r) { in PetscConvEstGetConvRateTS_Temporal_Private()
93 x[r] = PetscLog10Real(dt[r]); in PetscConvEstGetConvRateTS_Temporal_Private()
94 y[r] = PetscLog10Real(ce->errors[r * Nf + f]); in PetscConvEstGetConvRateTS_Temporal_Private()
123 PetscInt Nr = ce->Nr, r, Nf = ce->Nf, f, dim, oldlevel, oldnlev; in PetscConvEstGetConvRateTS_Spatial_Private() local
130 …Check(ce->r == 2.0, PetscObjectComm((PetscObject)ce), PETSC_ERR_SUP, "Only refinement factor 2 is … in PetscConvEstGetConvRateTS_Spatial_Private()
141 for (r = 0; r <= Nr; ++r) { in PetscConvEstGetConvRateTS_Spatial_Private()
147 …ll(PetscSNPrintf(stageName, PETSC_MAX_PATH_LEN - 1, "ConvEst Refinement Level %" PetscInt_FMT, r)); in PetscConvEstGetConvRateTS_Spatial_Private()
151 if (r > 0) { in PetscConvEstGetConvRateTS_Spatial_Private()
153 PetscCall(DMRefine(dm[r - 1], MPI_COMM_NULL, &dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
154 PetscCall(DMSetCoarseDM(dm[r], dm[r - 1])); in PetscConvEstGetConvRateTS_Spatial_Private()
158 PetscCall(DMClone(dm[r - 1], &dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
159 PetscCall(DMCopyDisc(dm[r - 1], dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
160 PetscCall(DMGetCoordinateDM(dm[r - 1], &cdm)); in PetscConvEstGetConvRateTS_Spatial_Private()
161 PetscCall(DMGetCoordinateDM(dm[r], &rcdm)); in PetscConvEstGetConvRateTS_Spatial_Private()
164 PetscCall(DMCopyTransform(ce->idm, dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
165 PetscCall(PetscObjectGetName((PetscObject)dm[r - 1], &dmname)); in PetscConvEstGetConvRateTS_Spatial_Private()
166 PetscCall(PetscObjectSetName((PetscObject)dm[r], dmname)); in PetscConvEstGetConvRateTS_Spatial_Private()
170 PetscCall(DMGetNullSpaceConstructor(dm[r - 1], f, &nspconstr)); in PetscConvEstGetConvRateTS_Spatial_Private()
171 PetscCall(DMSetNullSpaceConstructor(dm[r], f, nspconstr)); in PetscConvEstGetConvRateTS_Spatial_Private()
174 PetscCall(DMViewFromOptions(dm[r], NULL, "-conv_dm_view")); in PetscConvEstGetConvRateTS_Spatial_Private()
176 PetscCall(DMCreateGlobalVector(dm[r], &u)); in PetscConvEstGetConvRateTS_Spatial_Private()
177 PetscCall(DMGetField(dm[r], 0, NULL, &disc)); in PetscConvEstGetConvRateTS_Spatial_Private()
182 PetscCall(TSSetDM(ts, dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
183 PetscCall(DMTSSetBoundaryLocal(dm[r], DMPlexTSComputeBoundary, ctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
184 if (r > 0) { in PetscConvEstGetConvRateTS_Spatial_Private()
185 PetscCall(DMTSGetIFunctionLocal(dm[r - 1], &ifunc, &fctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
186 PetscCall(DMTSGetIJacobianLocal(dm[r - 1], &ijac, &jctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
187 PetscCall(DMTSGetRHSFunctionLocal(dm[r - 1], &rhsfunc, &rctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
188 if (ifunc) PetscCall(DMTSSetIFunctionLocal(dm[r], ifunc, fctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
189 if (ijac) PetscCall(DMTSSetIJacobianLocal(dm[r], ijac, jctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
190 if (rhsfunc) PetscCall(DMTSSetRHSFunctionLocal(dm[r], rhsfunc, rctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
199 PetscCall(PetscConvEstComputeInitialGuess(ce, r, dm[r], u)); in PetscConvEstGetConvRateTS_Spatial_Private()
203 PetscCall(PetscConvEstComputeError(ce, r, dm[r], u, &ce->errors[r * Nf])); in PetscConvEstGetConvRateTS_Spatial_Private()
210 PetscCall(DMGetLocalSection(dm[r], &s)); in PetscConvEstGetConvRateTS_Spatial_Private()
213 …PetscCallMPI(MPIU_Allreduce(&lsize, &ce->dofs[r * Nf + f], 1, MPIU_INT, MPI_SUM, PetscObjectComm((… in PetscConvEstGetConvRateTS_Spatial_Private()
214 PetscCall(PetscLogEventSetDof(ce->event, f, ce->dofs[r * Nf + f])); in PetscConvEstGetConvRateTS_Spatial_Private()
215 PetscCall(PetscLogEventSetError(ce->event, f, ce->errors[r * Nf + f])); in PetscConvEstGetConvRateTS_Spatial_Private()
218 PetscCall(PetscConvEstMonitorDefault(ce, r)); in PetscConvEstGetConvRateTS_Spatial_Private()
219 if (!r) { in PetscConvEstGetConvRateTS_Spatial_Private()
233 PetscCall(DMTSGetIFunctionLocal(dm[r - 1], &ifunc, &fctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
234 PetscCall(DMTSGetIJacobianLocal(dm[r - 1], &ijac, &jctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
235 PetscCall(DMTSGetRHSFunctionLocal(dm[r - 1], &rhsfunc, &rctx)); in PetscConvEstGetConvRateTS_Spatial_Private()
236 for (r = 1; r <= Nr; ++r) PetscCall(DMDestroy(&dm[r])); in PetscConvEstGetConvRateTS_Spatial_Private()
240 for (r = 0; r <= Nr; ++r) { in PetscConvEstGetConvRateTS_Spatial_Private()
241 x[r] = PetscLog10Real(ce->dofs[r * Nf + f]); in PetscConvEstGetConvRateTS_Spatial_Private()
242 y[r] = PetscLog10Real(ce->errors[r * Nf + f]); in PetscConvEstGetConvRateTS_Spatial_Private()