Lines Matching refs:flux
170 …cScalar *xR, PetscInt numConstants, const PetscScalar constants[], PetscScalar *flux, Physics phys) in PhysicsRiemann_SW_Rusanov() argument
197 for (i = 0; i < 1 + dim; ++i) flux[i] = zero / zero; in PhysicsRiemann_SW_Rusanov()
220 …for (i = 0; i < 1 + dim; i++) flux[i] = (0.5 * (fL.vals[i] + fR.vals[i]) + 0.5 * speed * (xL[i] - … in PhysicsRiemann_SW_Rusanov()
223 for (PetscInt j = 0; j < 3; ++j) PetscPrintf(PETSC_COMM_SELF, " | %g |\n", flux[j]); in PhysicsRiemann_SW_Rusanov()
244 CeedScalar flux[3]; in PhysicsRiemann_SW_Rusanov_CEED() local
254 PhysicsRiemann_SW_Rusanov(DIM, DIM + 1, NULL, n, qL, qR, 0, NULL, flux, &phys); in PhysicsRiemann_SW_Rusanov_CEED()
256 cL[i + Q * j] = -flux[j] / geom[i + Q * 2]; in PhysicsRiemann_SW_Rusanov_CEED()
257 cR[i + Q * j] = flux[j] / geom[i + Q * 3]; in PhysicsRiemann_SW_Rusanov_CEED()
270 …cScalar *xR, PetscInt numConstants, const PetscScalar constants[], PetscScalar *flux, Physics phys) in PhysicsRiemann_SW_HLL() argument
295 for (i = 0; i < 1 + dim; i++) flux[i] = zero; in PhysicsRiemann_SW_HLL()
325 for (i = 0; i < dim + 1; i++) flux[i] = fL.vals[i] * Norm2Real(n); in PhysicsRiemann_SW_HLL()
327 for (i = 0; i < dim + 1; i++) flux[i] = fR.vals[i] * Norm2Real(n); in PhysicsRiemann_SW_HLL()
329 …for (i = 0; i < dim + 1; i++) flux[i] = ((sR * fL.vals[i] - sL * fR.vals[i] + sR * sL * (xR[i] - x… in PhysicsRiemann_SW_HLL()
675 static int godunovflux(const PetscScalar *ul, const PetscScalar *ur, PetscScalar *flux, const Petsc… in godunovflux() argument
769 flux[0] = rhom * unm; in godunovflux()
774 flux[1] = fn * nn[0] + ft * tg[0]; in godunovflux()
775 flux[2] = fn * nn[1] + ft * tg[1]; in godunovflux()
778 if (ndim == 3) flux[3] = rhom * unm * ubm; in godunovflux()
779 …flux[ndim + 1] = (rhom * .5 * (unm * unm + utm * utm + ubm * ubm) + gamm / (gamm - 1.) * pm) * unm; in godunovflux()
784 …cScalar *xR, PetscInt numConstants, const PetscScalar constants[], PetscScalar *flux, Physics phys) in PhysicsRiemann_Euler_Godunov() argument
830 …for (i = 0; i < 2 + dim; i++) flux[i] = 0.5 * ((fL.vals[i] + fR.vals[i]) + speed * (xL[i] - xR[i])… in PhysicsRiemann_Euler_Godunov()
833 godunovflux(xL, xR, flux, nn, DIM, gamma); in PhysicsRiemann_Euler_Godunov()
834 for (i = 0; i < 2 + dim; i++) flux[i] *= s2; in PhysicsRiemann_Euler_Godunov()
853 CeedScalar flux[DIM + 2]; in PhysicsRiemann_Euler_Godunov_CEED() local
863 PhysicsRiemann_Euler_Godunov(DIM, DIM + 2, NULL, n, qL, qR, 0, NULL, flux, &phys); in PhysicsRiemann_Euler_Godunov_CEED()
865 cL[i + Q * j] = -flux[j] / geom[i + Q * 2]; in PhysicsRiemann_Euler_Godunov_CEED()
866 cR[i + Q * j] = flux[j] / geom[i + Q * 3]; in PhysicsRiemann_Euler_Godunov_CEED()