338   @param[in]   dim       Dimension of the element
343   @param[out]  dXdx      Inverse of the mapping Jacobian (shape [dim][dim]), or `NULL`
347 CEED_QFUNCTION_HELPER int QdataUnpack_ND(CeedInt dim, CeedInt Q, CeedInt i, const CeedScalar *q_dat…  in QdataUnpack_ND()  argument
348   switch (dim) {  in QdataUnpack_ND()
364   @param[in]   dim       Dimension of the element
369   @param[out]  dXdx      Inverse of the mapping Jacobian (shape [dim - 1][dim]), or `NULL`
370   @param[out]  normal    Components of the normal vector (shape [dim]), or `NULL`
374 CEED_QFUNCTION_HELPER int QdataBoundaryUnpack_ND(CeedInt dim, CeedInt Q, CeedInt i, const CeedScala…  in QdataBoundaryUnpack_ND()  argument
376   switch (dim) {  in QdataBoundaryUnpack_ND()
393   @param[in]   dim       Dimension of the element
398   @param[out]  dXdx      Inverse of the mapping Jacobian (shape [dim][dim]), or `NULL`
399   @param[out]  normal    Components of the normal vector (shape [dim]), or `NULL`
403 CEED_QFUNCTION_HELPER int QdataBoundaryGradientUnpack_ND(CeedInt dim, CeedInt Q, CeedInt i, const C…  in QdataBoundaryGradientUnpack_ND()  argument
405   switch (dim) {  in QdataBoundaryGradientUnpack_ND()
524 …@param[in]  dim        Topological dimension of the element (ie. number of derivative terms per co…
525   @param[in]  grad       QF gradient input, shape `[dim][num_comp][Q]`
526   @param[out] grad_local Gradient array at quadrature point Q, shape `[num_comp][dim]`
528 CEED_QFUNCTION_HELPER void GradUnpackND(CeedInt Q, CeedInt i, CeedInt num_comp, CeedInt dim, const …  in GradUnpackND()  argument
529   for (CeedInt d = 0; d < dim; d++) {  in GradUnpackND()
531       grad_local[dim * c + d] = grad[(Q * num_comp) * d + Q * c + i];  in GradUnpackND()
570   @param[in]  dXdx       Inverse of the mapping Jacobian (shape [dim][dim])
571   @param[in]  dim        Dimension of the problem
574 …ergenceND(const CeedScalar *grad_qn, const CeedScalar *dXdx, const CeedInt dim, CeedScalar *diverg…  in DivergenceND()  argument
575   for (CeedInt i = 0; i < dim; i++) {  in DivergenceND()
576     for (CeedInt j = 0; j < dim; j++) {  in DivergenceND()
577       *divergence += grad_qn[i * dim + j] * dXdx[j * dim + i];  in DivergenceND()