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()
141     if (i < ctx->sm) {  in FVSample_3WaySplit()
151       xi = ctx->xmin + ctx->sm * hs + 0.5 * hm + (i - ctx->sm) * hm;  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()
169 …xi = ctx->xmin + ctx->sm * hs + (ctx->mf - ctx->sm) * hm + (ctx->fm - ctx->mf) * hf + 0.5 * hm + (…  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()
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()
210     if (i < ctx->sm || i > ctx->ms - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]);  in SolutionErrorNorms_3WaySplit()
223   PetscInt     i, j, k, Mx, dof, xs, xm, sm = ctx->sm, mf = ctx->mf, fm = ctx->fm, ms = ctx->ms;  in FVRHSFunction_3WaySplit()  local
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()
278     (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  in FVRHSFunction_3WaySplit()
291     if (i < sm || i > ms) { /* slow region */  in FVRHSFunction_3WaySplit()
303     } else if (i == sm) { /* interface between slow and medium component */  in FVRHSFunction_3WaySplit()
397 …PetscInt     i, j, k, Mx, dof, xs, xm, islow = 0, sm = ctx->sm, ms = ctx->ms, lsbwidth = ctx->lsbw…  in FVRHSFunctionslow_3WaySplit()  local
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()
429 …if (i < sm - lsbwidth + 1 || i > ms + rsbwidth - 2) { /* slow components and the first and last fa…  in FVRHSFunctionslow_3WaySplit()
450       (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  in FVRHSFunctionslow_3WaySplit()
464     if (i < sm - lsbwidth) { /* slow region */  in FVRHSFunctionslow_3WaySplit()
479     if (i == sm - lsbwidth) { /* interface between slow and medium regions */  in FVRHSFunctionslow_3WaySplit()
527 …PetscInt     i, j, k, Mx, dof, xs, xm, islowbuffer = 0, sm = ctx->sm, ms = ctx->ms, lsbwidth = ctx…  in FVRHSFunctionslowbuffer_3WaySplit()  local
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()
559     if ((i > sm - lsbwidth - 2 && i < sm + 1) || (i > ms - 2 && i < ms + rsbwidth + 1)) {  in FVRHSFunctionslowbuffer_3WaySplit()
580       (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  in FVRHSFunctionslowbuffer_3WaySplit()
594     if (i == sm - lsbwidth) {  in FVRHSFunctionslowbuffer_3WaySplit()
605     if (i > sm - lsbwidth && i < sm) {  in FVRHSFunctionslowbuffer_3WaySplit()
619     if (i == sm) { /* interface between the slow region and the medium region */  in FVRHSFunctionslowbuffer_3WaySplit()
676 …PetscInt     i, j, k, Mx, dof, xs, xm, imedium = 0, sm = ctx->sm, mf = ctx->mf, fm = ctx->fm, ms =…  in FVRHSFunctionmedium_3WaySplit()  local
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()
708 …if ((i > sm - 2 && i < mf - lmbwidth + 1) || (i > fm + rmbwidth - 2 && i < ms + 1)) { /* slow comp…  in FVRHSFunctionmedium_3WaySplit()
729       (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  in FVRHSFunctionmedium_3WaySplit()
743     if (i == sm) { /* interface between slow and medium regions */  in FVRHSFunctionmedium_3WaySplit()
754     if (i > sm && i < mf - lmbwidth) { /* medium region */  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()
877       (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  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()
1024       (*ctx->limit3)(&info, cjmpL, cjmpR, ctx->sm, ctx->mf, ctx->fm, ctx->ms, i, ctx->cslope);  in FVRHSFunctionfast_3WaySplit()
1172   ctx.sm     = count_slow / 2;  in main()
1173   ctx.mf     = ctx.sm + count_medium / 2;  in main()
1194     if (i < ctx.sm - ctx.lsbwidth || i > ctx.ms + ctx.rsbwidth - 1)  in main()
1196 …else if ((i >= ctx.sm - ctx.lsbwidth && i < ctx.sm) || (i > ctx.ms - 1 && i <= ctx.ms + ctx.rsbwid…  in main()
1255       if (i < ctx.sm || i > ctx.ms - 1)  in main()
1300         if (i < ctx.sm || i > ctx.ms - 1)  in main()