21 CEED_QFUNCTION(masssetupf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in masssetupf()  argument
28   for (CeedInt i = 0; i < Q; i++) {  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()
30                       J[i + Q * 1] * (J[i + Q * 3] * J[i + Q * 8] - J[i + Q * 5] * J[i + Q * 6]) +  in masssetupf()
31                       J[i + Q * 2] * (J[i + Q * 3] * J[i + Q * 7] - J[i + Q * 4] * J[i + Q * 6]));  in masssetupf()
33 … rhs[i]         = qdata[i] * sqrt(x[i] * x[i] + x[i + Q] * x[i + Q] + x[i + 2 * Q] * x[i + 2 * Q]);  in masssetupf()
38 CEED_QFUNCTION(massf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) {  in massf()  argument
44   for (CeedInt i = 0; i < Q; i++) v[i] = qdata[i] * u[i];  in massf()
51 CEED_QFUNCTION(diffsetupf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in diffsetupf()  argument
58   for (CeedInt i = 0; i < Q; i++) {  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()
82     qdata[i + Q * 0]      = qw * (A11 * A11 + A12 * A12 + A13 * A13);  in diffsetupf()
83     qdata[i + Q * 1]      = qw * (A21 * A21 + A22 * A22 + A23 * A23);  in diffsetupf()
84     qdata[i + Q * 2]      = qw * (A31 * A31 + A32 * A32 + A33 * A33);  in diffsetupf()
85     qdata[i + Q * 3]      = qw * (A21 * A31 + A22 * A32 + A23 * A33);  in diffsetupf()
86     qdata[i + Q * 4]      = qw * (A11 * A31 + A12 * A32 + A13 * A33);  in diffsetupf()
87     qdata[i + Q * 5]      = qw * (A11 * A21 + A12 * A22 + A13 * A23);  in diffsetupf()
91 …_PI * M_PI * (k[0] * k[0] + k[1] * k[1] + k[2] * k[2]) * sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) *  in diffsetupf()
92              sin(M_PI * (c[1] + k[1] * x[i + Q * 1])) * sin(M_PI * (c[2] + k[2] * x[i + Q * 2]));  in diffsetupf()
97 CEED_QFUNCTION(diffusionf)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou…  in diffusionf()  argument
103   for (CeedInt i = 0; i < Q; i++) {  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()
107     vg[i + Q * 0]        = qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2;  in diffusionf()
108     vg[i + Q * 1]        = qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2;  in diffusionf()
109     vg[i + Q * 2]        = qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2;  in diffusionf()