38 #define DIM 2 /* Geometric dimension */  macro
72   PetscReal uh[DIM];
77   PetscReal vals[DIM + 1];
104   PetscReal ru[DIM];
110   PetscReal vals[DIM + 2];
138   for (i = 0; i < DIM; i++) prod += x[i] * y[i];  in DotDIMReal()
159   PetscReal   uhn, u[DIM];  in SWFlux()
166   for (i = 0; i < DIM; i++) f->uh[i] = u[i] * uhn + sw->gravity * PetscSqr(x->h) * n[i];  in SWFlux()
174   PetscReal   nn[DIM];  in PhysicsRiemann_SW_Rusanov()
248     for (CeedInt j = 0; j < DIM; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", n[j]);  in PhysicsRiemann_SW_Rusanov_CEED()
250     for (CeedInt j = 0; j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", qL[j]);  in PhysicsRiemann_SW_Rusanov_CEED()
252     for (CeedInt j = 0; j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", qR[j]);  in PhysicsRiemann_SW_Rusanov_CEED()
254     PhysicsRiemann_SW_Rusanov(DIM, DIM + 1, NULL, n, qL, qR, 0, NULL, flux, &phys);  in PhysicsRiemann_SW_Rusanov_CEED()
261 …for (CeedInt j = 0; j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cL[i + Q * j…  in PhysicsRiemann_SW_Rusanov_CEED()
263 …for (CeedInt j = 0; j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cR[i + Q * j…  in PhysicsRiemann_SW_Rusanov_CEED()
274   PetscReal   nn[DIM];  in PhysicsRiemann_SW_HLL()
372   for (i = 0; i < DIM; i++) f->ru[i] = nu * x->ru[i] + n[i] * p; /* r u^2 + p */  in EulerFlux()
789   PetscReal       cL, cR, speed, velL, velR, nn[DIM], s2;  in PhysicsRiemann_Euler_Godunov()
794   for (i = 0, s2 = 0.; i < DIM; i++) {  in PhysicsRiemann_Euler_Godunov()
799   for (i = 0.; i < DIM; i++) nn[i] /= s2;  in PhysicsRiemann_Euler_Godunov()
833     godunovflux(xL, xR, flux, nn, DIM, gamma);  in PhysicsRiemann_Euler_Godunov()
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()
853     CeedScalar       flux[DIM + 2];  in PhysicsRiemann_Euler_Godunov_CEED()
857     for (CeedInt j = 0; j < DIM; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", n[j]);  in PhysicsRiemann_Euler_Godunov_CEED()
859     for (CeedInt j = 0; j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", qL[j]);  in PhysicsRiemann_Euler_Godunov_CEED()
861     for (CeedInt j = 0; j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g |\n", qR[j]);  in PhysicsRiemann_Euler_Godunov_CEED()
863     PhysicsRiemann_Euler_Godunov(DIM, DIM + 2, NULL, n, qL, qR, 0, NULL, flux, &phys);  in PhysicsRiemann_Euler_Godunov_CEED()
864     for (CeedInt j = 0; j < DIM + 2; ++j) {  in PhysicsRiemann_Euler_Godunov_CEED()
870 …for (CeedInt j = 0; j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cL[i + Q * j…  in PhysicsRiemann_Euler_Godunov_CEED()
872 …for (CeedInt j = 0; j < DIM + 2; ++j) PetscPrintf(PETSC_COMM_SELF, "  | %g | (%g)\n", cR[i + Q * j…  in PhysicsRiemann_Euler_Godunov_CEED()