20 static inline PetscReal RangeMod(PetscReal a, PetscReal xmin, PetscReal xmax)  in RangeMod()  argument
22   PetscReal range = xmax - xmin;  in RangeMod()
23   return xmin + PetscFmodReal(range + PetscFmodReal(a, range), range);  in RangeMod()
54 …icsSample_Advect(void *vctx, PetscInt initial, FVBCType bctype, PetscReal xmin, PetscReal xmax, Pe…  in PhysicsSample_Advect()  argument
65     x0 = RangeMod(x - a * t, xmin, xmax);  in PhysicsSample_Advect()
137   hs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVSample_3WaySplit()
138   hm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVSample_3WaySplit()
139   hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVSample_3WaySplit()
142       xi = ctx->xmin + 0.5 * hs + i * hs;  in FVSample_3WaySplit()
147 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_3WaySplit()
151       xi = ctx->xmin + ctx->sm * hs + 0.5 * hm + (i - ctx->sm) * hm;  in FVSample_3WaySplit()
156 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_3WaySplit()
160       xi = ctx->xmin + ctx->sm * hs + (ctx->mf - ctx->sm) * hm + 0.5 * hf + (i - ctx->mf) * hf;  in FVSample_3WaySplit()
165 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_3WaySplit()
169 …xi = ctx->xmin + ctx->sm * hs + (ctx->mf - ctx->sm) * hm + (ctx->fm - ctx->mf) * hf + 0.5 * hm + (…  in FVSample_3WaySplit()
174 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_3WaySplit()
178 …xi = ctx->xmin + ctx->sm * hs + (ctx->mf - ctx->sm) * hm + (ctx->fm - ctx->mf) * hf + (ctx->ms - c…  in FVSample_3WaySplit()
183 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_3WaySplit()
202   hs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in SolutionErrorNorms_3WaySplit()
203   hm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in SolutionErrorNorms_3WaySplit()
204   hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in SolutionErrorNorms_3WaySplit()
233   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunction_3WaySplit()
234   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunction_3WaySplit()
235   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunction_3WaySplit()
407   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunctionslow_3WaySplit()
408   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunctionslow_3WaySplit()
409   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunctionslow_3WaySplit()
537   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunctionslowbuffer_3WaySplit()
538   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunctionslowbuffer_3WaySplit()
539   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunctionslowbuffer_3WaySplit()
686   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunctionmedium_3WaySplit()
687   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunctionmedium_3WaySplit()
688   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunctionmedium_3WaySplit()
834   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunctionmediumbuffer_3WaySplit()
835   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunctionmediumbuffer_3WaySplit()
836   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunctionmediumbuffer_3WaySplit()
983   hxs = (ctx->xmax - ctx->xmin) / 8.0 / ctx->sm;  in FVRHSFunctionfast_3WaySplit()
984   hxm = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->mf - ctx->sm);  in FVRHSFunctionfast_3WaySplit()
985   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fm - ctx->mf);  in FVRHSFunctionfast_3WaySplit()
1115   ctx.xmin   = -1.0;  in main()
1118   PetscCall(PetscOptionsReal("-xmin", "X min", "", ctx.xmin, &ctx.xmin, NULL));  in main()
1155 …ll(DMDASetUniformCoordinates(da, ctx.xmin + 0.5 * (ctx.xmax - ctx.xmin) / Mx, ctx.xmax + 0.5 * (ct…  in main()
1242     const PetscReal    hs = (ctx.xmax - ctx.xmin) / 4.0 / count_slow;  in main()
1243     const PetscReal    hm = (ctx.xmax - ctx.xmin) / 2.0 / count_medium;  in main()
1244     const PetscReal    hf = (ctx.xmax - ctx.xmin) / 4.0 / count_fast;  in main()