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
750         for (j = 0; j < dof; j++) f[imedium * dof + j] += ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
751         imedium++;  in FVRHSFunctionmedium_3WaySplit()
761         for (j = 0; j < dof; j++) f[(imedium - 1) * dof + j] -= ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
764         for (j = 0; j < dof; j++) f[imedium * dof + j] += ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
765         imedium++;  in FVRHSFunctionmedium_3WaySplit()
775         for (j = 0; j < dof; j++) f[(imedium - 1) * dof + j] -= ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
785         for (j = 0; j < dof; j++) f[imedium * dof + j] += ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
786         imedium++;  in FVRHSFunctionmedium_3WaySplit()
796         for (j = 0; j < dof; j++) f[(imedium - 1) * dof + j] -= ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
799         for (j = 0; j < dof; j++) f[imedium * dof + j] += ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
800         imedium++;  in FVRHSFunctionmedium_3WaySplit()
810         for (j = 0; j < dof; j++) f[(imedium - 1) * dof + j] -= ctx->flux[j] / hxm;  in FVRHSFunctionmedium_3WaySplit()
1090 …, xs, xm, Mx, draw = 0, count_slow, count_medium, count_fast, islow = 0, imedium = 0, ifast = 0, *…  in main()  local
1199       for (k = 0; k < dof; k++) index_medium[imedium++] = i * dof + k;  in main()
1206   PetscCall(ISCreateGeneral(PETSC_COMM_WORLD, imedium, index_medium, PETSC_COPY_VALUES, &ctx.ism));  in main()