479   PetscReal      hs, hf;  in FVSample_2WaySplit()  local
487   hs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in FVSample_2WaySplit()
491       xi = ctx->xmin + 0.5 * hs + i * hs;  in FVSample_2WaySplit()
495         xj = xi + hs * (j - N / 2) / (PetscReal)N;  in FVSample_2WaySplit()
500       xi = ctx->xmin + ctx->sf * hs + 0.5 * hf + (i - ctx->sf) * hf;  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()
513         xj = xi + hs * (j - N / 2) / (PetscReal)N;  in FVSample_2WaySplit()
529   PetscReal          hs, hf;  in SolutionErrorNorms_2WaySplit()  local
535   hs = (ctx->xmax - ctx->xmin) * 3.0 / 8.0 / ctx->sf;  in SolutionErrorNorms_2WaySplit()
540     if (i < ctx->sf || i > ctx->fs - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]);  in SolutionErrorNorms_2WaySplit()
1323     const PetscReal    hs = (ctx.xmax - ctx.xmin) * 3.0 / 4.0 / count_slow;  in main()  local
1337           mass_initial = mass_initial + hs * ptr_X0[i * dof + k];  in main()
1338           mass_final   = mass_final + hs * ptr_X[i * dof + k];  in main()
1376 …   for (k = 0; k < dof; k++) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]);  in main()