118     for (CeedInt j = 0; j < 3; j++) {                         // Rows of each Jacobian matrix  in ConvectiveFluxJacobian_Euler()  local
119       dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j];  in ConvectiveFluxJacobian_Euler()
122 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (…  in ConvectiveFluxJacobian_Euler()
125       dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.);  in ConvectiveFluxJacobian_Euler()
143     for (CeedInt j = 0; j < 3; j++) {  in Covariant_length_along_vector()  local
144       vec_dot_jacobian[i] += dXdx[j][i] * vec[i];  in Covariant_length_along_vector()
192     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsShockTube()  local
256     for (CeedInt j = 0; j < 3; j++) {  in EulerShockTube()  local
258         du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho;  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()
273     for (CeedInt j = 0; j < 5; j++) {  in EulerShockTube()  local
274       v[j][i] = 0;  in EulerShockTube()
275       for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0;  in EulerShockTube()
280 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd…  in EulerShockTube()  local
283     for (CeedInt j = 0; j < 3; j++) {  in EulerShockTube()  local
285 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[…  in EulerShockTube()
286                                     (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]);  in EulerShockTube()
291 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[…  in EulerShockTube()  local
304       for (CeedInt j = 0; j < 3; j++) j_vec[j] = drhodx[j] / (drho_norm + 1e-20);  in EulerShockTube()  local
316       for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * nu_shock * drhodx[j];  in EulerShockTube()  local
319         for (CeedInt j = 0; j < 3; j++) dv[j][k][i] -= wdetJ * nu_shock * du[k][j];  in EulerShockTube()  local
322       for (CeedInt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * nu_shock * dEdx[j];  in EulerShockTube()  local
333     for (CeedInt j = 0; j < 3; j++) {  in EulerShockTube()  local
334       dqdx[0][j] = drhodx[j];  in EulerShockTube()
335       dqdx[4][j] = dEdx[j];  in EulerShockTube()
336       for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j];  in EulerShockTube()
341     for (CeedInt j = 0; j < 3; j++) {  in EulerShockTube()  local
343         for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j];  in EulerShockTube()
358         for (CeedInt j = 0; j < 3; j++) {  in EulerShockTube()  local
361               stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l];  in EulerShockTube()
365         for (CeedInt j = 0; j < 5; j++) {  in EulerShockTube()  local
366 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd…  in EulerShockTube()