41   const CeedScalar(*w)                = in[2];  // Note: *X = in[0]  in SetupDiffGeo()
43   CeedScalar(*qd) = out[0];  in SetupDiffGeo()
47   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupDiffGeo()
50     for (CeedInt j = 0; j < dim; j++) {  in SetupDiffGeo()
51       for (CeedInt k = 0; k < dim; k++) {  in SetupDiffGeo()
58     const CeedScalar detJ = J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2];  in SetupDiffGeo()
61     qd[i + Q * 0]       = w[i] * detJ;  in SetupDiffGeo()
62     qd[i + Q * 1]       = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]);  in SetupDiffGeo()
63     qd[i + Q * 2]       = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]);  in SetupDiffGeo()
64     qd[i + Q * 3]       = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]);  in SetupDiffGeo()
65     qd[i + Q * 4]       = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]);  in SetupDiffGeo()
66     qd[i + Q * 5]       = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]);  in SetupDiffGeo()
67     qd[i + Q * 6]       = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]);  in SetupDiffGeo()
70   return 0;  in SetupDiffGeo()
80   const CeedScalar *x = in[0], *w = in[1];  in SetupDiffRhs()
81   CeedScalar       *true_soln = out[0], *rhs = out[1];  in SetupDiffRhs()
84   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupDiffRhs()
85     const CeedScalar c[3] = {0, 1., 2.};  in SetupDiffRhs()
88 …true_soln[i] = sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) * sin(M_PI * (c[1] + k[1] * x[i + Q * 1]))…  in SetupDiffRhs()
90     rhs[i] = w[i + Q * 0] * M_PI * M_PI * (k[0] * k[0] + k[1] * k[1] + k[2] * k[2]) * true_soln[i];  in SetupDiffRhs()
92   return 0;  in SetupDiffRhs()
106   const CeedScalar *ug = in[0], *q_data = in[1];  in Diff()
107   CeedScalar       *vg = out[0];  in Diff()
110   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Diff()
112     const CeedScalar du[3] = {ug[i + Q * 0], ug[i + Q * 1], ug[i + Q * 2]};  in Diff()
120     for (int j = 0; j < 3; j++) {  // j = direction of vg  in Diff()
121 …  vg[i + j * Q] = (du[0] * dXdxdXdx_T[0][j] + du[1] * dXdxdXdx_T[1][j] + du[2] * dXdxdXdx_T[2][j]);  in Diff()
124   return 0;  in Diff()