59 CEED_QFUNCTION(SetupDiffGeo)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *…  in SetupDiffGeo()  argument
64   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupDiffGeo()
66     const CeedScalar xx[3] = {X[i + 0 * Q], X[i + 1 * Q], X[i + 2 * Q]};  in SetupDiffGeo()
73         {J[i + Q * 0], J[i + Q * 3]},  in SetupDiffGeo()
74         {J[i + Q * 1], J[i + Q * 4]},  in SetupDiffGeo()
75         {J[i + Q * 2], J[i + Q * 5]}  in SetupDiffGeo()
109     q_data[i + Q * 0] = mod_J * w[i];  in SetupDiffGeo()
130     q_data[i + Q * 1] = dxdXTdxdX_inv[0][0];  in SetupDiffGeo()
131     q_data[i + Q * 2] = dxdXTdxdX_inv[1][1];  in SetupDiffGeo()
132     q_data[i + Q * 3] = dxdXTdxdX_inv[0][1];  in SetupDiffGeo()
142 CEED_QFUNCTION(SetupDiffRhs)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *…  in SetupDiffRhs()  argument
153   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupDiffRhs()
155     CeedScalar x = X[i + Q * 0], y = X[i + Q * 1], z = X[i + Q * 2];  in SetupDiffRhs()
165     true_soln[i + Q * 0] = sin(lambda) * cos(theta);  in SetupDiffRhs()
167     rhs[i + Q * 0] = q_data[i + Q * 0] * 2 * sin(lambda) * cos(theta) / (R * R);  in SetupDiffRhs()
183 CEED_QFUNCTION(Diff)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) {  in Diff()  argument
190   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Diff()
192     const CeedScalar du[2] = {ug[i + Q * 0], ug[i + Q * 1]};  in Diff()
194     const CeedScalar w_det_J = q_data[i + Q * 0];  in Diff()
198         {q_data[i + Q * 1], q_data[i + Q * 3]},  in Diff()
199         {q_data[i + Q * 3], q_data[i + Q * 2]}  in Diff()
203       vg[i + j * Q] = w_det_J * (du[0] * dXdxdXdx_T[0][j] + du[1] * dXdxdXdx_T[1][j]);  in Diff()