549     CeedScalar wdetJ, dXdx[9], F_diff[3] = {0.};  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()  local
557       MatVecNM(dXdx, grad_E_ref, dim, dim, CEED_NOTRANSPOSE, F_diff);  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
558       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
562     MatVecNM(dXdx, F_diff, dim, dim, CEED_NOTRANSPOSE, F_diff_dXdx);  in DivDiffusiveFluxVolumeRHS_AdvDif_Generic()
586     CeedScalar wdetJ, normal[3], dXdx[9], F_diff[3] = {0.};  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()  local
594       MatVecNM(dXdx, grad_E_ref, dim, dim, CEED_NOTRANSPOSE, F_diff);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
595       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
598     v[i] = wdetJ * DotN(F_diff, normal, dim);  in DivDiffusiveFluxBoundaryRHS_AdvDif_Generic()
621     CeedScalar wdetJ, dXdx[9], F_diff[3] = {0.};  in DiffusiveFluxRHS_AdvDif_Generic()  local
629       MatVecNM(dXdx, grad_E_ref, dim, dim, CEED_NOTRANSPOSE, F_diff);  in DiffusiveFluxRHS_AdvDif_Generic()
630       ScaleN(F_diff, -context->diffusion_coeff, dim);  in DiffusiveFluxRHS_AdvDif_Generic()
632     for (CeedInt k = 0; k < dim; k++) v[k][i] = wdetJ * F_diff[k];  in DiffusiveFluxRHS_AdvDif_Generic()