21 CEED_QFUNCTION(masssetupf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in masssetupf()
22   CeedScalar       *qdata = out[0], *rhs = out[1];  in masssetupf()
23   const CeedScalar *x = in[0];  in masssetupf()
24   const CeedScalar *J = in[1];  in masssetupf()
25   const CeedScalar *w = in[2];  in masssetupf()
29     CeedScalar det = (J[i + Q * 0] * (J[i + Q * 4] * J[i + Q * 8] - J[i + Q * 5] * J[i + Q * 7]) -  in masssetupf()
38 CEED_QFUNCTION(massf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) {  in massf()
39   const CeedScalar *u     = in[0];  in massf()
40   const CeedScalar *qdata = in[1];  in massf()
41   CeedScalar       *v     = out[0];  in massf()
51 CEED_QFUNCTION(diffsetupf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in diffsetupf()
52   const CeedScalar *x     = in[0];  in diffsetupf()
53   const CeedScalar *J     = in[1];  in diffsetupf()
54   const CeedScalar *w     = in[2];  in diffsetupf()
55   CeedScalar       *qdata = out[0], *rhs = out[1];  in diffsetupf()
63     const CeedScalar J11  = J[i + Q * 0];  in diffsetupf()
64     const CeedScalar J21  = J[i + Q * 1];  in diffsetupf()
65     const CeedScalar J31  = J[i + Q * 2];  in diffsetupf()
66     const CeedScalar J12  = J[i + Q * 3];  in diffsetupf()
67     const CeedScalar J22  = J[i + Q * 4];  in diffsetupf()
68     const CeedScalar J32  = J[i + Q * 5];  in diffsetupf()
69     const CeedScalar J13  = J[i + Q * 6];  in diffsetupf()
70     const CeedScalar J23  = J[i + Q * 7];  in diffsetupf()
71     const CeedScalar J33  = J[i + Q * 8];  in diffsetupf()
72     const CeedScalar A11  = J22 * J33 - J23 * J32;  in diffsetupf()
73     const CeedScalar A12  = J13 * J32 - J12 * J33;  in diffsetupf()
74     const CeedScalar A13  = J12 * J23 - J13 * J22;  in diffsetupf()
75     const CeedScalar A21  = J23 * J31 - J21 * J33;  in diffsetupf()
76     const CeedScalar A22  = J11 * J33 - J13 * J31;  in diffsetupf()
77     const CeedScalar A23  = J13 * J21 - J11 * J23;  in diffsetupf()
78     const CeedScalar A31  = J21 * J32 - J22 * J31;  in diffsetupf()
79     const CeedScalar A32  = J12 * J31 - J11 * J32;  in diffsetupf()
80     const CeedScalar A33  = J11 * J22 - J12 * J21;  in diffsetupf()
81     const CeedScalar qw   = w[i] / (J11 * A11 + J21 * A12 + J31 * A13);  in diffsetupf()
88     const CeedScalar c[3] = {0, 1., 2.};  in diffsetupf()
89     const CeedScalar k[3] = {1., 2., 3.};  in diffsetupf()
90     const CeedScalar rho  = w[i] * (J11 * A11 + J21 * A12 + J31 * A13);  in diffsetupf()
97 CEED_QFUNCTION(diffusionf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in diffusionf()
98   const CeedScalar *ug    = in[0];  in diffusionf()
99   const CeedScalar *qdata = in[1];  in diffusionf()
100   CeedScalar       *vg    = out[0];  in diffusionf()
104     const CeedScalar ug0 = ug[i + Q * 0];  in diffusionf()
105     const CeedScalar ug1 = ug[i + Q * 1];  in diffusionf()
106     const CeedScalar ug2 = ug[i + Q * 2];  in diffusionf()