46   const CeedScalar(*q_data)            = in[2];  in MassFunction_Newtonian()
59     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in MassFunction_Newtonian()
193   const CeedScalar(*q_data)               = in[2];  in RHSFunction_Newtonian()
207     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in RHSFunction_Newtonian()
315   const CeedScalar(*q_data)               = in[3];  in IFunction_Newtonian()
332     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in IFunction_Newtonian()
425   const CeedScalar(*q_data)          = in[2];  in IJacobian_Newtonian()
443     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in IJacobian_Newtonian()
601   const CeedScalar(*q_data)          = in[2];  in DivDiffusiveFluxVolumeRHS_NS()
614     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in DivDiffusiveFluxVolumeRHS_NS()
654   const CeedScalar(*q_data)        = in[2];  in DivDiffusiveFluxBoundaryRHS_NS()
667     QdataBoundaryGradientUnpack_3D(Q, i, q_data, &wdetJ, dXdx, normal);  in DivDiffusiveFluxBoundaryRHS_NS()
703   const CeedScalar(*q_data)        = in[2];  in DiffusiveFluxRHS_NS()
716     QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in DiffusiveFluxRHS_NS()