30       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in build_mass_diff()  local
32         q_data[0][i] = w[i] * J[0][0][i];  in build_mass_diff()
35         q_data[1][i] = w[i] / J[0][0][i];  in build_mass_diff()
41       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in build_mass_diff()  local
44         const CeedScalar J00 = J[0][0][i];  in build_mass_diff()
45         const CeedScalar J10 = J[0][1][i];  in build_mass_diff()
46         const CeedScalar J01 = J[1][0][i];  in build_mass_diff()
47         const CeedScalar J11 = J[1][1][i];  in build_mass_diff()
48         const CeedScalar qw  = w[i] / (J00 * J11 - J10 * J01);  in build_mass_diff()
51         q_data[0][i] = w[i] * (J00 * J11 - J10 * J01);  in build_mass_diff()
54         q_data[1][i] = qw * (J01 * J01 + J11 * J11);  in build_mass_diff()
55         q_data[2][i] = qw * (J00 * J00 + J10 * J10);  in build_mass_diff()
56         q_data[3][i] = -qw * (J00 * J01 + J10 * J11);  in build_mass_diff()
62       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in build_mass_diff()  local
71 …(k + 1) % 3][(j + 1) % 3][i] * J[(k + 2) % 3][(j + 2) % 3][i] - J[(k + 2) % 3][(j + 1) % 3][i] * J…  in build_mass_diff()
76 … 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()
79         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()
86         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()
87         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()
88         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()
89         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()
90         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()
91         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()
112       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in apply_mass_diff()  local
114         v[i] = q_data[0][i] * u[i];  in apply_mass_diff()
117         vg[i] = q_data[1][i] * ug[i];  in apply_mass_diff()
126       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in apply_mass_diff()  local
128         v[i] = q_data[0][i] * u[i];  in apply_mass_diff()
136             {q_data[1][i], q_data[3][i]},  in apply_mass_diff()
137             {q_data[3][i], q_data[2][i]}  in apply_mass_diff()
141 …for (int j = 0; j < 2; j++) vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j]);  in apply_mass_diff()
150       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in apply_mass_diff()  local
152         v[i] = q_data[0][i] * u[i];  in apply_mass_diff()
161             {q_data[1][i], q_data[6][i], q_data[5][i]},  in apply_mass_diff()
162             {q_data[6][i], q_data[2][i], q_data[4][i]},  in apply_mass_diff()
163             {q_data[5][i], q_data[4][i], q_data[3][i]}  in apply_mass_diff()
167 …for (int j = 0; j < 3; j++) vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j] …  in apply_mass_diff()