15   // in[1] is quadrature weights, size (Q)  in build_mass_diff()
16   const CeedScalar *w             = in[1];  in build_mass_diff()
23     case 11: {  // dim = 1, space_dim = 1  in build_mass_diff()
24       const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0];  in build_mass_diff()
31         q_data[1][i] = w[i] / J[0][0][i];  in build_mass_diff()
39         //    1 3           2 1           -J10  J00  in build_mass_diff()
41         const CeedScalar J10 = J[0][1][i];  in build_mass_diff()
42         const CeedScalar J01 = J[1][0][i];  in build_mass_diff()
43         const CeedScalar J11 = J[1][1][i];  in build_mass_diff()
50         q_data[1][i] = qw * (J01 * J01 + J11 * J11);  in build_mass_diff()
64 …J[(k + 1) % 3][(j + 1) % 3][i] * J[(k + 2) % 3][(j + 2) % 3][i] - J[(k + 2) % 3][(j + 1) % 3][i] *…  in build_mass_diff()
69 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]);  in build_mass_diff()
72         q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]);  in build_mass_diff()
76         // 1 6 5  in build_mass_diff()
79         q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]);  in build_mass_diff()
80         q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]);  in build_mass_diff()
81         q_data[3][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]);  in build_mass_diff()
82         q_data[4][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]);  in build_mass_diff()
83         q_data[5][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]);  in build_mass_diff()
84         q_data[6][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]);  in build_mass_diff()
94   // in[0], out[0] solution values with shape [1, 1, Q]  in apply_mass_diff()
95   // in[1], out[1] solution gradients with shape [dim, 1, Q]  in apply_mass_diff()
100     case 1: {  in apply_mass_diff()
101       const CeedScalar *u = in[0], *ug = in[1];  in apply_mass_diff()
102       CeedScalar       *v = out[0], *vg = out[1];  in apply_mass_diff()
109         vg[i] = q_data[1][i] * ug[i];  in apply_mass_diff()
114       const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1];  in apply_mass_diff()
116       CeedScalar(*vg)[CEED_Q_VLA]       = (CeedScalar(*)[CEED_Q_VLA])out[1];  in apply_mass_diff()
125         // 1 3  in apply_mass_diff()
128             {q_data[1][i], q_data[3][i]},  in apply_mass_diff()
134           vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j]);  in apply_mass_diff()
140       const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1];  in apply_mass_diff()
142       CeedScalar(*vg)[CEED_Q_VLA]       = (CeedScalar(*)[CEED_Q_VLA])out[1];  in apply_mass_diff()
151         // 1 6 5  in apply_mass_diff()
155             {q_data[1][i], q_data[6][i], q_data[5][i]},  in apply_mass_diff()
162 …vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j] + ug[2][i] * dXdxdXdx_T[2][j…  in apply_mass_diff()