228 CEED_QFUNCTION(PhysicsRiemann_SW_Rusanov_CEED)(PetscCtx ctx, CeedInt Q, const CeedScalar *const in[…  in PhysicsRiemann_SW_Rusanov_CEED()  argument
239   CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i)  in PhysicsRiemann_SW_Rusanov_CEED()
241     const CeedScalar qL[3] = {xL[i + Q * 0], xL[i + Q * 1], xL[i + Q * 2]};  in PhysicsRiemann_SW_Rusanov_CEED()
242     const CeedScalar qR[3] = {xR[i + Q * 0], xR[i + Q * 1], xR[i + Q * 2]};  in PhysicsRiemann_SW_Rusanov_CEED()
243     const CeedScalar n[2]  = {geom[i + Q * 0], geom[i + Q * 1]};  in PhysicsRiemann_SW_Rusanov_CEED()
247     PetscPrintf(PETSC_COMM_SELF, "Face %d Normal\n", (int)info[i + Q * 0]);  in PhysicsRiemann_SW_Rusanov_CEED()
249     PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: left state\n", (int)info[i + Q * 1]);  in PhysicsRiemann_SW_Rusanov_CEED()
251     PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: right state\n", (int)info[i + Q * 2]);  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()
260     PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: left flux\n", (int)info[i + Q * 1]);  in PhysicsRiemann_SW_Rusanov_CEED()
261 … j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cL[i + Q * j], geom[i + Q * 2]);  in PhysicsRiemann_SW_Rusanov_CEED()
262     PetscPrintf(PETSC_COMM_SELF, "Cell %d Element Residual: right flux\n", (int)info[i + Q * 2]);  in PhysicsRiemann_SW_Rusanov_CEED()
263 … j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cR[i + Q * j], geom[i + Q * 3]);  in PhysicsRiemann_SW_Rusanov_CEED()
840 CEED_QFUNCTION(PhysicsRiemann_Euler_Godunov_CEED)(PetscCtx ctx, CeedInt Q, const CeedScalar *const …  in PhysicsRiemann_Euler_Godunov_CEED()  argument
848   CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i)  in PhysicsRiemann_Euler_Godunov_CEED()
850     const CeedScalar qL[DIM + 2] = {xL[i + Q * 0], xL[i + Q * 1], xL[i + Q * 2], xL[i + Q * 3]};  in PhysicsRiemann_Euler_Godunov_CEED()
851     const CeedScalar qR[DIM + 2] = {xR[i + Q * 0], xR[i + Q * 1], xR[i + Q * 2], xR[i + Q * 3]};  in PhysicsRiemann_Euler_Godunov_CEED()
852     const CeedScalar n[DIM]      = {geom[i + Q * 0], geom[i + Q * 1]};  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()
870 … j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cL[i + Q * j], geom[i + Q * 2]);  in PhysicsRiemann_Euler_Godunov_CEED()
872 … j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cR[i + Q * j], geom[i + Q * 3]);  in PhysicsRiemann_Euler_Godunov_CEED()