Lines Matching refs:uR
290 …n_Advect(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Advect() argument
297 flux[0] = PetscMax(0, speed) * uL[0] + PetscMin(0, speed) * uR[0]; in PhysicsRiemann_Advect()
426 …rs_Exact(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Burgers_Exact() argument
429 if (uL[0] < uR[0]) { /* rarefaction */ in PhysicsRiemann_Burgers_Exact()
430 flux[0] = (uL[0] * uR[0] < 0) ? 0 /* sonic rarefaction */ in PhysicsRiemann_Burgers_Exact()
431 : 0.5 * PetscMin(PetscSqr(uL[0]), PetscSqr(uR[0])); in PhysicsRiemann_Burgers_Exact()
433 flux[0] = 0.5 * PetscMax(PetscSqr(uL[0]), PetscSqr(uR[0])); in PhysicsRiemann_Burgers_Exact()
435 *maxspeed = (PetscAbs(uL[0]) > PetscAbs(uR[0])) ? uL[0] : uR[0]; in PhysicsRiemann_Burgers_Exact()
439 …gers_Roe(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Burgers_Roe() argument
444 speed = 0.5 * (uL[0] + uR[0]); in PhysicsRiemann_Burgers_Roe()
445 flux[0] = 0.25 * (PetscSqr(uL[0]) + PetscSqr(uR[0])) - 0.5 * PetscAbs(speed) * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_Roe()
446 if (uL[0] <= 0 && 0 <= uR[0]) flux[0] = 0; /* Entropy fix for sonic rarefaction */ in PhysicsRiemann_Burgers_Roe()
451 …gers_LxF(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Burgers_LxF() argument
459 fR = 0.5 * PetscSqr(uR[0]); in PhysicsRiemann_Burgers_LxF()
460 flux[0] = 0.5 * (fL + fR) - 0.5 * c * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_LxF()
465 …_Rusanov(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Burgers_Rusanov() argument
471 c = PetscMax(PetscAbs(uL[0]), PetscAbs(uR[0])); in PhysicsRiemann_Burgers_Rusanov()
473 fR = 0.5 * PetscSqr(uR[0]); in PhysicsRiemann_Burgers_Rusanov()
474 flux[0] = 0.5 * (fL + fR) - 0.5 * c * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_Rusanov()
554 …ic_Exact(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Traffic_Exact() argument
559 if (uL[0] < uR[0]) { in PhysicsRiemann_Traffic_Exact()
560 flux[0] = PetscMin(TrafficFlux(a, uL[0]), TrafficFlux(a, uR[0])); in PhysicsRiemann_Traffic_Exact()
562 …flux[0] = (uR[0] < 0.5 && 0.5 < uL[0]) ? TrafficFlux(a, 0.5) : PetscMax(TrafficFlux(a, uL[0]), Tra… in PhysicsRiemann_Traffic_Exact()
564 *maxspeed = a * MaxAbs(1 - 2 * uL[0], 1 - 2 * uR[0]); in PhysicsRiemann_Traffic_Exact()
568 …ffic_Roe(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Traffic_Roe() argument
574 speed = a * (1 - (uL[0] + uR[0])); in PhysicsRiemann_Traffic_Roe()
575 …flux[0] = 0.5 * (TrafficFlux(a, uL[0]) + TrafficFlux(a, uR[0])) - 0.5 * PetscAbs(speed) * (uR[0]… in PhysicsRiemann_Traffic_Roe()
580 …ffic_LxF(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Traffic_LxF() argument
587 speed = a * (1 - (uL[0] + uR[0])); in PhysicsRiemann_Traffic_LxF()
588 …flux[0] = 0.5 * (TrafficFlux(a, uL[0]) + TrafficFlux(a, uR[0])) - 0.5 * phys->lxf_speed * (uR[0]… in PhysicsRiemann_Traffic_LxF()
593 …_Rusanov(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Traffic_Rusanov() argument
599 speed = a * PetscMax(PetscAbs(1 - 2 * uL[0]), PetscAbs(1 - 2 * uR[0])); in PhysicsRiemann_Traffic_Rusanov()
600 flux[0] = 0.5 * (TrafficFlux(a, uL[0]) + TrafficFlux(a, uR[0])) - 0.5 * speed * (uR[0] - uL[0]); in PhysicsRiemann_Traffic_Rusanov()
735 …cs_Exact(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Acoustics_Exact() argument
747 flux[0] = Al[0][0] * uR[0] + Al[0][1] * uR[1] + Ar[0][0] * uL[0] + Ar[0][1] * uL[1]; in PhysicsRiemann_Acoustics_Exact()
748 flux[1] = Al[1][0] * uR[0] + Al[1][1] * uR[1] + Ar[1][0] * uL[0] + Ar[1][1] * uL[1]; in PhysicsRiemann_Acoustics_Exact()
821 …oGas_Roe(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_IsoGas_Roe() argument
829 … = (uL[1] / PetscSqrtScalar(uL[0]) + uR[1] / PetscSqrtScalar(uR[0])) / (PetscSqrtScalar(uL[0]) + P… in PhysicsRiemann_IsoGas_Roe()
831 du[0] = uR[0] - uL[0]; in PhysicsRiemann_IsoGas_Roe()
832 du[1] = uR[1] - uL[1]; in PhysicsRiemann_IsoGas_Roe()
850 } else if (ustar[1] / ustar[0] < -c && -c < uR[1] / uR[0]) { /* 2-wave is sonic rarefaction */ in PhysicsRiemann_IsoGas_Roe()
852 ufan[0] = uR[0] * PetscExpScalar(-uR[1] / (uR[0] * c) - 1); in PhysicsRiemann_IsoGas_Roe()
857 IsoGasFlux(c, uR, fR); in PhysicsRiemann_IsoGas_Roe()
867 …as_Exact(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_IsoGas_Exact() argument
874 } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}, star; in PhysicsRiemann_IsoGas_Exact()
920 IsoGasFlux(c, uR, flux); in PhysicsRiemann_IsoGas_Exact()
930 …_Rusanov(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_IsoGas_Rusanov() argument
936 } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}; in PhysicsRiemann_IsoGas_Rusanov()
941 IsoGasFlux(c, uR, fR); in PhysicsRiemann_IsoGas_Rusanov()
943 flux[0] = 0.5 * (fL[0] + fR[0]) + 0.5 * s * (uL[0] - uR[0]); in PhysicsRiemann_IsoGas_Rusanov()
944 flux[1] = 0.5 * (fL[1] + fR[1]) + 0.5 * s * (uL[1] - uR[1]); in PhysicsRiemann_IsoGas_Rusanov()
1012 …ow_Exact(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Shallow_Exact() argument
1018 } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}, star; in PhysicsRiemann_Shallow_Exact()
1076 ShallowFlux(phys, uR, flux); in PhysicsRiemann_Shallow_Exact()
1086 …_Rusanov(void *vctx, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, … in PhysicsRiemann_Shallow_Rusanov() argument
1092 } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}; in PhysicsRiemann_Shallow_Rusanov()
1097 ShallowFlux(phys, uR, fR); in PhysicsRiemann_Shallow_Rusanov()
1099 flux[0] = 0.5 * (fL[0] + fR[0]) + 0.5 * s * (uL[0] - uR[0]); in PhysicsRiemann_Shallow_Rusanov()
1100 flux[1] = 0.5 * (fL[1] + fR[1]) + 0.5 * s * (uL[1] - uR[1]); in PhysicsRiemann_Shallow_Rusanov()
1225 PetscScalar *uL, *uR; in FVRHSFunction() local
1227 uR = &ctx->uLR[dof]; in FVRHSFunction()
1230 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunction()
1232 PetscCall((*ctx->physics.riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed)); in FVRHSFunction()