38 CEED_QFUNCTION(SetupDiffGeo)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *…  in SetupDiffGeo()  argument
47   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  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()
76 CEED_QFUNCTION(SetupDiffRhs)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *…  in SetupDiffRhs()  argument
84   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupDiffRhs()
88 … sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) * sin(M_PI * (c[1] + k[1] * x[i + Q * 1])) * sin(M_PI * …  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()
105 CEED_QFUNCTION(Diff)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) {  in Diff()  argument
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()
115         {q_data[i + 1 * Q], q_data[i + 2 * Q], q_data[i + 3 * Q]},  in Diff()
116         {q_data[i + 2 * Q], q_data[i + 4 * Q], q_data[i + 5 * Q]},  in Diff()
117         {q_data[i + 3 * Q], q_data[i + 5 * Q], q_data[i + 6 * Q]}  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()