191   PetscReal    hx, cfl_idt = 0;  in FVRHSFunctionslow()  local
200   hx = (ctx->xmax - ctx->xmin) / Mx;  in FVRHSFunctionslow()
223 …eristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin + hx * i));  in FVRHSFunctionslow()
239       info.hx = hx;  in FVRHSFunctionslow()
241       for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */  in FVRHSFunctionslow()
257         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionslow()
258         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionslow()
260 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionslow()
261 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionslow()
263         for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionslow()
266         for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionslow()
273         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionslow()
274         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionslow()
276 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionslow()
277 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionslow()
279         for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionslow()
296   PetscReal    hx, cfl_idt = 0;  in FVRHSFunctionfast()  local
305   hx = (ctx->xmax - ctx->xmin) / Mx;  in FVRHSFunctionfast()
327 …eristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin + hx * i));  in FVRHSFunctionfast()
342       info.hx = hx;  in FVRHSFunctionfast()
344       for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */  in FVRHSFunctionfast()
360         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionfast()
361         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionfast()
363 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionfast()
364 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionfast()
366         for (j = 0; j < dof; j++) f[(i - len_slow - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionfast()
369         for (j = 0; j < dof; j++) f[(i - len_slow) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionfast()
376         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionfast()
377         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionfast()
379 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionfast()
380 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionfast()
382         for (j = 0; j < dof; j++) f[(i - len_slow) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionfast()
397   PetscReal    hx, cfl_idt = 0;  in FVRHSFunctionslow2()  local
406   hx = (ctx->xmax - ctx->xmin) / Mx;  in FVRHSFunctionslow2()
428 …eristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin + hx * i));  in FVRHSFunctionslow2()
443       info.hx = hx;  in FVRHSFunctionslow2()
445       for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */  in FVRHSFunctionslow2()
461         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
462         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
464 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionslow2()
465 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  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()
477         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
478         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
480 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionslow2()
481 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionslow2()
483         for (j = 0; j < dof; j++) f[(i - len_slow1 - 1) * dof + j] -= ctx->flux[j] / hx;  in FVRHSFunctionslow2()
490         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
491         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionslow2()
493 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionslow2()
494 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionslow2()
496         for (j = 0; j < dof; j++) f[(i - len_slow1) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionslow2()
512   PetscReal    hx, cfl_idt = 0;  in FVRHSFunctionfast2()  local
521   hx = (ctx->xmax - ctx->xmin) / Mx;  in FVRHSFunctionfast2()
544 …eristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin + hx * i));  in FVRHSFunctionfast2()
559       info.hx = hx;  in FVRHSFunctionfast2()
561       for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */  in FVRHSFunctionfast2()
577         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionfast2()
578         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionfast2()
580 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionfast2()
581 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  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()
593         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2;  in FVRHSFunctionfast2()
594         uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2;  in FVRHSFunctionfast2()
596 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-…  in FVRHSFunctionfast2()
597 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */  in FVRHSFunctionfast2()
599         for (j = 0; j < dof; j++) f[(i - len_slow1 - len_slow2) * dof + j] += ctx->flux[j] / hx;  in FVRHSFunctionfast2()