Lines Matching refs:flux
216 PetscScalar *flux; /* Flux across interface */ member
290 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Advect() argument
297 flux[0] = PetscMax(0, speed) * uL[0] + PetscMin(0, speed) * uR[0]; in PhysicsRiemann_Advect()
426 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Burgers_Exact() argument
430 flux[0] = (uL[0] * uR[0] < 0) ? 0 /* sonic rarefaction */ in PhysicsRiemann_Burgers_Exact()
433 flux[0] = 0.5 * PetscMax(PetscSqr(uL[0]), PetscSqr(uR[0])); in PhysicsRiemann_Burgers_Exact()
439 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Burgers_Roe() argument
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 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Burgers_LxF() argument
460 flux[0] = 0.5 * (fL + fR) - 0.5 * c * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_LxF()
465 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Burgers_Rusanov() argument
474 flux[0] = 0.5 * (fL + fR) - 0.5 * c * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_Rusanov()
554 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Traffic_Exact() argument
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()
568 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Traffic_Roe() argument
575 …flux[0] = 0.5 * (TrafficFlux(a, uL[0]) + TrafficFlux(a, uR[0])) - 0.5 * PetscAbs(speed) * (uR[0]… in PhysicsRiemann_Traffic_Roe()
580 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Traffic_LxF() argument
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 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Traffic_Rusanov() argument
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 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… 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 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_IsoGas_Roe() argument
849 IsoGasFlux(c, ufan, flux); in PhysicsRiemann_IsoGas_Roe()
854 IsoGasFlux(c, ufan, flux); in PhysicsRiemann_IsoGas_Roe()
860 flux[i] = 0.5 * (fL[i] + fR[i]) - 0.5 * absdu; in PhysicsRiemann_IsoGas_Roe()
867 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_IsoGas_Exact() argument
909 IsoGasFlux(c, ufan, flux); in PhysicsRiemann_IsoGas_Exact()
914 IsoGasFlux(c, ufan, flux); in PhysicsRiemann_IsoGas_Exact()
917 IsoGasFlux(c, uL, flux); in PhysicsRiemann_IsoGas_Exact()
920 IsoGasFlux(c, uR, flux); in PhysicsRiemann_IsoGas_Exact()
924 IsoGasFlux(c, ustar, flux); in PhysicsRiemann_IsoGas_Exact()
930 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_IsoGas_Rusanov() argument
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 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Shallow_Exact() argument
1065 ShallowFlux(phys, ufan, flux); in PhysicsRiemann_Shallow_Exact()
1070 ShallowFlux(phys, ufan, flux); in PhysicsRiemann_Shallow_Exact()
1073 ShallowFlux(phys, uL, flux); in PhysicsRiemann_Shallow_Exact()
1076 ShallowFlux(phys, uR, flux); in PhysicsRiemann_Shallow_Exact()
1080 ShallowFlux(phys, ustar, flux); in PhysicsRiemann_Shallow_Exact()
1086 …x, PetscInt m, const PetscScalar *uL, const PetscScalar *uR, PetscScalar *flux, PetscReal *maxspee… in PhysicsRiemann_Shallow_Rusanov() argument
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()
1232 PetscCall((*ctx->physics.riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed)); in FVRHSFunction()
1236 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hx; in FVRHSFunction()
1239 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hx; in FVRHSFunction()
1480 PetscCall(PetscMalloc3(2 * dof, &ctx.uLR, dof, &ctx.flux, dof, &ctx.speeds)); in main()
1545 PetscCall(PetscFree3(ctx.uLR, ctx.flux, ctx.speeds)); in main()