139 …N_HELPER CeedScalar Covariant_length_along_vector(CeedScalar vec[3], const CeedScalar dXdx[3][3]) {  in Covariant_length_along_vector()
144       vec_dot_jacobian[i] += dXdx[j][i] * vec[i];  in Covariant_length_along_vector()
168 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed…  in Tau_spatial()
172 …CeedScalar h = 2 / sqrt(dXdx[0][i] * dXdx[0][i] + dXdx[1][i] * dXdx[1][i] + dXdx[2][i] * dXdx[2][i…  in Tau_spatial()
248     CeedScalar       wdetJ, dXdx[3][3];  in EulerShockTube()  local
249     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in EulerShockTube()
259         drhodx[j] += drho[k] * dXdx[k][j];  in EulerShockTube()
260         dEdx[j] += dE[k] * dXdx[k][j];  in EulerShockTube()
262           dUdx[j][k] += dU[j][l] * dXdx[l][k];  in EulerShockTube()
263           dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l];  // dXdx_j,k * dXdx_k,j  in EulerShockTube()
280 … < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rho * u[2] …  in EulerShockTube()
285 …tJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[1] + (j == 1 ? P : 0))…  in EulerShockTube()
286                                     (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]);  in EulerShockTube()
291 … 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()
309         h_shock = Covariant_length_along_vector(j_vec, dXdx);  in EulerShockTube()
351     Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau);  in EulerShockTube()
366 … < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] …  in EulerShockTube()