45   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Setup2d()  local
47     InvertMappingJacobian_2D(Q, i, J, dXdx, &detJ);  in Setup2d()
48     const CeedScalar wdetJ = w[i] * detJ;  in Setup2d()
50     StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data);  in Setup2d()
51     StoredValuesPack(Q, i, 1, 4, (const CeedScalar *)dXdx, q_data);  in Setup2d()
86   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in SetupBoundary2d()  local
88     NormalVectorFromdxdX_2D(Q, i, J, normal, &detJb);  in SetupBoundary2d()
89     const CeedScalar wdetJ = w[i] * detJb;  in SetupBoundary2d()
91     StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur);  in SetupBoundary2d()
92     StoredValuesPack(Q, i, 1, 2, normal, q_data_sur);  in SetupBoundary2d()
140   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Setup2D_3Dcoords()  local
143     NormalVectorFromdxdX_3D(Q, i, J, normal, &detJb);  in Setup2D_3Dcoords()
144     InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx);  in Setup2D_3Dcoords()
145     const CeedScalar wdetJ = w[i] * detJb;  in Setup2D_3Dcoords()
147     StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur);  in Setup2D_3Dcoords()
148     StoredValuesPack(Q, i, 1, 6, (const CeedScalar *)dXdx, q_data_sur);  in Setup2D_3Dcoords()
176   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Setup2DBoundaryGradient()  local
179     NormalVectorFromdxdX_2D(Q, i, J_face, normal, &detJ_face);  in Setup2DBoundaryGradient()
180     const CeedScalar wdetJ = w[i] * detJ_face;  in Setup2DBoundaryGradient()
181     InvertMappingJacobian_2D(Q, i, J_cell, dXdx, NULL);  in Setup2DBoundaryGradient()
183     StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur);  in Setup2DBoundaryGradient()
184     StoredValuesPack(Q, i, 1, 4, (CeedScalar *)dXdx, q_data_sur);  in Setup2DBoundaryGradient()
185     StoredValuesPack(Q, i, 5, 2, normal, q_data_sur);  in Setup2DBoundaryGradient()