396   PetscInt     i, j, k, Mx, dof, xs, xm, len_slow1, len_slow2;  in FVRHSFunctionslow2()  local
414   PetscCall(ISGetSize(ctx->iss, &len_slow1));  in FVRHSFunctionslow2()
427     if (i < len_slow1 + len_slow2 + 1 && i > len_slow1 - 2) {  in FVRHSFunctionslow2()
457     if (i < len_slow1 + len_slow2 && i > len_slow1) { /* slow parts can be changed based on a */  in FVRHSFunctionslow2()
467         for (j = 0; j < dof; j++) f[(i - len_slow1 - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionslow2()
470         for (j = 0; j < dof; j++) f[(i - len_slow1) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionslow2()
473     if (i == len_slow1 + len_slow2) { /* slow parts can be changed based on a */  in FVRHSFunctionslow2()
483         for (j = 0; j < dof; j++) f[(i - len_slow1 - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionslow2()
486     if (i == len_slow1) { /* slow parts can be changed based on a */  in FVRHSFunctionslow2()
496         for (j = 0; j < dof; j++) f[(i - len_slow1) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionslow2()
511   PetscInt     i, j, k, Mx, dof, xs, xm, len_slow1, len_slow2;  in FVRHSFunctionfast2()  local
529   PetscCall(ISGetSize(ctx->iss, &len_slow1));  in FVRHSFunctionfast2()
543     if (i > len_slow1 + len_slow2 - 2) {  in FVRHSFunctionfast2()
573     if (i > len_slow1 + len_slow2) {  in FVRHSFunctionfast2()
583         for (j = 0; j < dof; j++) f[(i - len_slow1 - len_slow2 - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionfast2()
586         for (j = 0; j < dof; j++) f[(i - len_slow1 - len_slow2) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionfast2()
589     if (i == len_slow1 + len_slow2) {  in FVRHSFunctionfast2()
599         for (j = 0; j < dof; j++) f[(i - len_slow1 - len_slow2) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionfast2()