38 static PetscErrorCode PhysicsRiemann_Advect(void *vctx, PetscInt m, const PetscScalar *uL, const Pe…  in PhysicsRiemann_Advect()  argument
45   flux[0]   = PetscMax(0, speed) * uL[0] + PetscMin(0, speed) * uR[0];  in PhysicsRiemann_Advect()
271     PetscScalar *uL, *uR;  in FVRHSFunction_2WaySplit()  local
272     uL = &ctx->uLR[0];  in FVRHSFunction_2WaySplit()
276         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunction_2WaySplit()
279       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunction_2WaySplit()
288         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunction_2WaySplit()
291       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunction_2WaySplit()
300         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2;  in FVRHSFunction_2WaySplit()
303       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunction_2WaySplit()
312         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2;  in FVRHSFunction_2WaySplit()
315       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunction_2WaySplit()
324         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunction_2WaySplit()
327       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunction_2WaySplit()
421     PetscScalar *uL, *uR;  in FVRHSFunctionslow_2WaySplit()  local
422     uL = &ctx->uLR[0];  in FVRHSFunctionslow_2WaySplit()
426         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslow_2WaySplit()
429       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslow_2WaySplit()
441         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslow_2WaySplit()
444       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslow_2WaySplit()
451         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslow_2WaySplit()
454       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslow_2WaySplit()
462         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslow_2WaySplit()
465       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslow_2WaySplit()
548     PetscScalar *uL, *uR;  in FVRHSFunctionslowbuffer_2WaySplit()  local
549     uL = &ctx->uLR[0];  in FVRHSFunctionslowbuffer_2WaySplit()
553         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
556       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
564         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
567       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
578         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
581       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
588         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
591       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
599         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
602       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
613         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionslowbuffer_2WaySplit()
616       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionslowbuffer_2WaySplit()
693     PetscScalar *uL, *uR;  in FVRHSFunctionfast_2WaySplit()  local
694     uL = &ctx->uLR[0];  in FVRHSFunctionfast_2WaySplit()
698         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2;  in FVRHSFunctionfast_2WaySplit()
701       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionfast_2WaySplit()
709         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2;  in FVRHSFunctionfast_2WaySplit()
712       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionfast_2WaySplit()
723         uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2;  in FVRHSFunctionfast_2WaySplit()
726       PetscCall((*ctx->physics2.riemann2)(ctx->physics2.user, dof, uL, uR, ctx->flux, &maxspeed));  in FVRHSFunctionfast_2WaySplit()