13   for (CeedInt i = 0; i < Q; i++) {  in setup_mass()  local
14     rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]);  in setup_mass()
31   for (CeedInt i = 0; i < Q; i++) {  in setup_diff()  local
34     const CeedScalar J11 = J[i + Q * 0];  in setup_diff()
35     const CeedScalar J21 = J[i + Q * 1];  in setup_diff()
36     const CeedScalar J12 = J[i + Q * 2];  in setup_diff()
37     const CeedScalar J22 = J[i + Q * 3];  in setup_diff()
38     const CeedScalar w   = qw[i] / (J11 * J22 - J21 * J12);  in setup_diff()
39     qd[i + Q * 0]        = w * (J12 * J12 + J22 * J22);  in setup_diff()
40     qd[i + Q * 1]        = w * (J11 * J11 + J21 * J21);  in setup_diff()
41     qd[i + Q * 2]        = -w * (J11 * J12 + J21 * J22);  in setup_diff()
59   for (CeedInt i = 0; i < Q; i++) {  in apply()  local
61     v[i] = qd_mass[i] * u[i];  in apply()
63     const CeedScalar du0 = du[i + Q * 0];  in apply()
64     const CeedScalar du1 = du[i + Q * 1];  in apply()
65     dv[i + Q * 0]        = qd_diff[i + Q * 0] * du0 + qd_diff[i + Q * 2] * du1;  in apply()
66     dv[i + Q * 1]        = qd_diff[i + Q * 2] * du0 + qd_diff[i + Q * 1] * du1;  in apply()
83   for (CeedInt i = 0; i < Q; i++) {  in apply_lin()  local
84     const CeedScalar du0 = du[i + Q * 0];  in apply_lin()
85     const CeedScalar du1 = du[i + Q * 1];  in apply_lin()
86     v[i + Q * 0]         = qd[i + Q * 0] * du0 + qd[i + Q * 3] * du1 + qd[i + Q * 6] * u[i];  in apply_lin()
87     dv[i + Q * 0]        = qd[i + Q * 1] * du0 + qd[i + Q * 4] * du1 + qd[i + Q * 7] * u[i];  in apply_lin()
88     dv[i + Q * 1]        = qd[i + Q * 2] * du0 + qd[i + Q * 5] * du1 + qd[i + Q * 8] * u[i];  in apply_lin()