43 static inline PetscReal RangeMod(PetscReal a, PetscReal xmin, PetscReal xmax)  in RangeMod()  argument
45   PetscReal range = xmax - xmin;  in RangeMod()
46   return xmin + PetscFmodReal(range + PetscFmodReal(a, range), range);  in RangeMod()
80 …icsSample_Advect(void *vctx, PetscInt initial, FVBCType bctype, PetscReal xmin, PetscReal xmax, Pe…  in PhysicsSample_Advect()  argument
91     x0 = RangeMod(x - a * t, xmin, xmax);  in PhysicsSample_Advect()
301 …csSample_Shallow(void *vctx, PetscInt initial, FVBCType bctype, PetscReal xmin, PetscReal xmax, Pe…  in PhysicsSample_Shallow()  argument
487   hs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVSample_2WaySplit()
488   hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in FVSample_2WaySplit()
491       xi = ctx->xmin + 0.5 * hs + i * hs;  in FVSample_2WaySplit()
496 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_2WaySplit()
500       xi = ctx->xmin + ctx->sf * hs + 0.5 * hf + (i - ctx->sf) * hf;  in FVSample_2WaySplit()
505 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_2WaySplit()
509       xi = ctx->xmin + ctx->sf * hs + (ctx->fs - ctx->sf) * hf + 0.5 * hs + (i - ctx->fs) * hs;  in FVSample_2WaySplit()
514 …PetscCall((*ctx->physics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->…  in FVSample_2WaySplit()
535   hs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in SolutionErrorNorms_2WaySplit()
536   hf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in SolutionErrorNorms_2WaySplit()
562   hxs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVRHSFunction_2WaySplit()
563   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in FVRHSFunction_2WaySplit()
587       PetscCall(ctx->physics2.inflow(ctx, time, ctx->xmin, ctx->ub));  in FVRHSFunction_2WaySplit()
738   hxs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVRHSFunctionslow_2WaySplit()
739   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in FVRHSFunctionslow_2WaySplit()
760       PetscCall(ctx->physics2.inflow(ctx, time, ctx->xmin, ctx->ub));  in FVRHSFunctionslow_2WaySplit()
889   hxs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVRHSFunctionslowbuffer_2WaySplit()
890   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in FVRHSFunctionslowbuffer_2WaySplit()
911       PetscCall(ctx->physics2.inflow(ctx, time, ctx->xmin, ctx->ub));  in FVRHSFunctionslowbuffer_2WaySplit()
1060   hxs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVRHSFunctionfast_2WaySplit()
1061   hxf = (ctx->xmax - ctx->xmin) / 4.0 / (ctx->fs - ctx->sf);  in FVRHSFunctionfast_2WaySplit()
1082       PetscCall(ctx->physics2.inflow(ctx, time, ctx->xmin, ctx->ub));  in FVRHSFunctionfast_2WaySplit()
1215   ctx.xmin       = -1.0;  in main()
1222   PetscCall(PetscOptionsReal("-xmin", "X min", "", ctx.xmin, &ctx.xmin, NULL));  in main()
1260 …ll(DMDASetUniformCoordinates(da, ctx.xmin + 0.5 * (ctx.xmax - ctx.xmin) / Mx, ctx.xmax + 0.5 * (ct…  in main()
1323     const PetscReal    hs = (ctx.xmax - ctx.xmin) * 3.0 / 4.0 / count_slow;  in main()
1324     const PetscReal    hf = (ctx.xmax - ctx.xmin) / 4.0 / count_fast;  in main()