28       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = w[i] / J[i]; }  in f_build_diff()  local
32       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in f_build_diff()  local
35         const CeedScalar J11 = J[i + Q * 0];  in f_build_diff()
36         const CeedScalar J21 = J[i + Q * 1];  in f_build_diff()
37         const CeedScalar J12 = J[i + Q * 2];  in f_build_diff()
38         const CeedScalar J22 = J[i + Q * 3];  in f_build_diff()
39         const CeedScalar qw  = w[i] / (J11 * J22 - J21 * J12);  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()
47       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in f_build_diff()  local
51         const CeedScalar J11 = J[i + Q * 0];  in f_build_diff()
52         const CeedScalar J21 = J[i + Q * 1];  in f_build_diff()
53         const CeedScalar J31 = J[i + Q * 2];  in f_build_diff()
54         const CeedScalar J12 = J[i + Q * 3];  in f_build_diff()
55         const CeedScalar J22 = J[i + Q * 4];  in f_build_diff()
56         const CeedScalar J32 = J[i + Q * 5];  in f_build_diff()
57         const CeedScalar J13 = J[i + Q * 6];  in f_build_diff()
58         const CeedScalar J23 = J[i + Q * 7];  in f_build_diff()
59         const CeedScalar J33 = J[i + Q * 8];  in f_build_diff()
69         const CeedScalar qw  = w[i] / (J11 * A11 + J21 * A12 + J31 * A13);  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()
92       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { vg[i] = ug[i] * qdata[i]; }  in f_apply_diff()  local
96       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in f_apply_diff()  local
97         const CeedScalar ug0 = ug[i + Q * 0];  in f_apply_diff()
98         const CeedScalar ug1 = ug[i + Q * 1];  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()
105       CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in f_apply_diff()  local
106         const CeedScalar ug0 = ug[i + Q * 0];  in f_apply_diff()
107         const CeedScalar ug1 = ug[i + Q * 1];  in f_apply_diff()
108         const CeedScalar ug2 = ug[i + Q * 2];  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()