23   CeedScalar       *qdata = out[0];  in f_build_diff()  local
28       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = w[i] / J[i]; }  in f_build_diff()
40         qdata[i + Q * 0]     = qw * (J12 * J12 + J22 * J22);  in f_build_diff()
41         qdata[i + Q * 1]     = qw * (J11 * J11 + J21 * J21);  in f_build_diff()
42         qdata[i + Q * 2]     = -qw * (J11 * J12 + J21 * J22);  in f_build_diff()
70         qdata[i + Q * 0]     = qw * (A11 * A11 + A12 * A12 + A13 * A13);  in f_build_diff()
71         qdata[i + Q * 1]     = qw * (A21 * A21 + A22 * A22 + A23 * A23);  in f_build_diff()
72         qdata[i + Q * 2]     = qw * (A31 * A31 + A32 * A32 + A33 * A33);  in f_build_diff()
73         qdata[i + Q * 3]     = qw * (A21 * A31 + A22 * A32 + A23 * A33);  in f_build_diff()
74         qdata[i + Q * 4]     = qw * (A11 * A31 + A12 * A32 + A13 * A33);  in f_build_diff()
75         qdata[i + Q * 5]     = qw * (A11 * A21 + A12 * A22 + A13 * A23);  in f_build_diff()
86   const CeedScalar *ug = in[0], *qdata = in[1];  in f_apply_diff()  local
92       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { vg[i] = ug[i] * qdata[i]; }  in f_apply_diff()
99         vg[i + Q * 0]        = qdata[i + Q * 0] * ug0 + qdata[i + Q * 2] * ug1;  in f_apply_diff()
100         vg[i + Q * 1]        = qdata[i + Q * 2] * ug0 + qdata[i + Q * 1] * ug1;  in f_apply_diff()
109 …   vg[i + Q * 0]        = qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2;  in f_apply_diff()
110 …   vg[i + Q * 1]        = qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2;  in f_apply_diff()
111 …   vg[i + Q * 2]        = qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2;  in f_apply_diff()