134 …N_HELPER CeedScalar Covariant_length_along_vector(CeedScalar vec[3], const CeedScalar dXdx[3][3]) {  in Covariant_length_along_vector()
137   MatVec3(dXdx, vec, CEED_TRANSPOSE, vec_dot_jacobian);  in Covariant_length_along_vector()
156 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed…  in Tau_spatial()
160     CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i]));  in Tau_spatial()
236     CeedScalar       wdetJ, dXdx[3][3];  in EulerShockTube()  local
237     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in EulerShockTube()
247         drhodx[j] += drho[k] * dXdx[k][j];  in EulerShockTube()
248         dEdx[j] += dE[k] * dXdx[k][j];  in EulerShockTube()
250           dUdx[j][k] += dU[j][l] * dXdx[l][k];  in EulerShockTube()
251           dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l];  // dXdx_j,k * dXdx_k,j  in EulerShockTube()
268 … < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rho * u[2] …  in EulerShockTube()
273 …tJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[1] + (j == 1 ? P : 0))…  in EulerShockTube()
274                                     (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]);  in EulerShockTube()
279 … j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[j][1] + u[2] * dXdx…  in EulerShockTube()
297         h_shock = Covariant_length_along_vector(j_vec, dXdx);  in EulerShockTube()
339     Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau);  in EulerShockTube()
354 … < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] …  in EulerShockTube()