Lines Matching refs:Q
36 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_build_mass() argument
45 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_build_mass()
51 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_build_mass()
53 qdata[i] = (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]) * w[i]; in f_build_mass()
57 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_build_mass()
59 qdata[i] = (J[i + Q * 0] * (J[i + Q * 4] * J[i + Q * 8] - J[i + Q * 5] * J[i + Q * 7]) - in f_build_mass()
60 J[i + Q * 1] * (J[i + Q * 3] * J[i + Q * 8] - J[i + Q * 5] * J[i + Q * 6]) + in f_build_mass()
61 J[i + Q * 2] * (J[i + Q * 3] * J[i + Q * 7] - J[i + Q * 4] * J[i + Q * 6])) * in f_build_mass()
75 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_apply_mass() argument
83 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_apply_mass()
96 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_apply_mass_vec() argument
107 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_apply_mass_vec()
109 v[i + Q * 0] = JxW[i] * u[i + Q * 0]; in f_apply_mass_vec()
113 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_apply_mass_vec()
115 v[i + Q * 0] = JxW[i] * u[i + Q * 0]; in f_apply_mass_vec()
116 v[i + Q * 1] = JxW[i] * u[i + Q * 1]; in f_apply_mass_vec()
120 CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) in f_apply_mass_vec()
122 v[i + Q * 0] = JxW[i] * u[i + Q * 0]; in f_apply_mass_vec()
123 v[i + Q * 1] = JxW[i] * u[i + Q * 1]; in f_apply_mass_vec()
124 v[i + Q * 2] = JxW[i] * u[i + Q * 2]; in f_apply_mass_vec()
137 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_build_poisson() argument
146 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_build_poisson()
152 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_build_poisson()
154 const CeedScalar J11 = J[i + Q * 0]; in f_build_poisson()
155 const CeedScalar J21 = J[i + Q * 1]; in f_build_poisson()
156 const CeedScalar J12 = J[i + Q * 2]; in f_build_poisson()
157 const CeedScalar J22 = J[i + Q * 3]; in f_build_poisson()
159 qdata[i + Q * 0] = qw * (J12 * J12 + J22 * J22); in f_build_poisson()
160 qdata[i + Q * 1] = qw * (J11 * J11 + J21 * J21); in f_build_poisson()
161 qdata[i + Q * 2] = -qw * (J11 * J12 + J21 * J22); in f_build_poisson()
165 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_build_poisson()
167 const CeedScalar J11 = J[i + Q * 0]; in f_build_poisson()
168 const CeedScalar J21 = J[i + Q * 1]; in f_build_poisson()
169 const CeedScalar J31 = J[i + Q * 2]; in f_build_poisson()
170 const CeedScalar J12 = J[i + Q * 3]; in f_build_poisson()
171 const CeedScalar J22 = J[i + Q * 4]; in f_build_poisson()
172 const CeedScalar J32 = J[i + Q * 5]; in f_build_poisson()
173 const CeedScalar J13 = J[i + Q * 6]; in f_build_poisson()
174 const CeedScalar J23 = J[i + Q * 7]; in f_build_poisson()
175 const CeedScalar J33 = J[i + Q * 8]; in f_build_poisson()
186 qdata[i + Q * 0] = qw * (A11 * A11 + A12 * A12 + A13 * A13); in f_build_poisson()
187 qdata[i + Q * 1] = qw * (A21 * A21 + A22 * A22 + A23 * A23); in f_build_poisson()
188 qdata[i + Q * 2] = qw * (A31 * A31 + A32 * A32 + A33 * A33); in f_build_poisson()
189 qdata[i + Q * 3] = qw * (A21 * A31 + A22 * A32 + A23 * A33); in f_build_poisson()
190 qdata[i + Q * 4] = qw * (A11 * A31 + A12 * A32 + A13 * A33); in f_build_poisson()
191 qdata[i + Q * 5] = qw * (A11 * A21 + A12 * A22 + A13 * A23); in f_build_poisson()
204 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_apply_poisson() argument
213 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson()
219 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson()
221 const CeedScalar ug0 = ug[i + Q * 0]; in f_apply_poisson()
222 const CeedScalar ug1 = ug[i + Q * 1]; in f_apply_poisson()
223 vg[i + Q * 0] = qdata[i + Q * 0] * ug0 + qdata[i + Q * 2] * ug1; in f_apply_poisson()
224 vg[i + Q * 1] = qdata[i + Q * 2] * ug0 + qdata[i + Q * 1] * ug1; in f_apply_poisson()
228 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson()
230 const CeedScalar ug0 = ug[i + Q * 0]; in f_apply_poisson()
231 const CeedScalar ug1 = ug[i + Q * 1]; in f_apply_poisson()
232 const CeedScalar ug2 = ug[i + Q * 2]; in f_apply_poisson()
233 vg[i + Q * 0] = qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2; in f_apply_poisson()
234 vg[i + Q * 1] = qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2; in f_apply_poisson()
235 vg[i + Q * 2] = qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2; in f_apply_poisson()
248 (void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) in f_apply_poisson_vec() argument
257 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson_vec()
263 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson_vec()
266 const CeedScalar ug0 = ug[i + Q * 0 + Q * 2 * 0]; in f_apply_poisson_vec()
267 const CeedScalar ug1 = ug[i + Q * 1 + Q * 2 * 0]; in f_apply_poisson_vec()
268 vg[i + Q * 0 + Q * 2 * 0] = qdata[i + Q * 0] * ug0 + qdata[i + Q * 2] * ug1; in f_apply_poisson_vec()
269 vg[i + Q * 1 + Q * 2 * 0] = qdata[i + Q * 2] * ug0 + qdata[i + Q * 1] * ug1; in f_apply_poisson_vec()
272 const CeedScalar ug0 = ug[i + Q * 0 + Q * 2 * 1]; in f_apply_poisson_vec()
273 const CeedScalar ug1 = ug[i + Q * 1 + Q * 2 * 1]; in f_apply_poisson_vec()
274 vg[i + Q * 0 + Q * 2 * 1] = qdata[i + Q * 0] * ug0 + qdata[i + Q * 2] * ug1; in f_apply_poisson_vec()
275 vg[i + Q * 1 + Q * 2 * 1] = qdata[i + Q * 2] * ug0 + qdata[i + Q * 1] * ug1; in f_apply_poisson_vec()
280 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) in f_apply_poisson_vec()
283 const CeedScalar ug0 = ug[i + Q * 0 + Q * 3 * 0]; in f_apply_poisson_vec()
284 const CeedScalar ug1 = ug[i + Q * 1 + Q * 3 * 0]; in f_apply_poisson_vec()
285 const CeedScalar ug2 = ug[i + Q * 2 + Q * 3 * 0]; in f_apply_poisson_vec()
286 vg[i + Q * 0 + Q * 3 * 0] = in f_apply_poisson_vec()
287 qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2; in f_apply_poisson_vec()
288 vg[i + Q * 1 + Q * 3 * 0] = in f_apply_poisson_vec()
289 qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2; in f_apply_poisson_vec()
290 vg[i + Q * 2 + Q * 3 * 0] = in f_apply_poisson_vec()
291 qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2; in f_apply_poisson_vec()
294 const CeedScalar ug0 = ug[i + Q * 0 + Q * 3 * 1]; in f_apply_poisson_vec()
295 const CeedScalar ug1 = ug[i + Q * 1 + Q * 3 * 1]; in f_apply_poisson_vec()
296 const CeedScalar ug2 = ug[i + Q * 2 + Q * 3 * 1]; in f_apply_poisson_vec()
297 vg[i + Q * 0 + Q * 3 * 1] = in f_apply_poisson_vec()
298 qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2; in f_apply_poisson_vec()
299 vg[i + Q * 1 + Q * 3 * 1] = in f_apply_poisson_vec()
300 qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2; in f_apply_poisson_vec()
301 vg[i + Q * 2 + Q * 3 * 1] = in f_apply_poisson_vec()
302 qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2; in f_apply_poisson_vec()
305 const CeedScalar ug0 = ug[i + Q * 0 + Q * 3 * 2]; in f_apply_poisson_vec()
306 const CeedScalar ug1 = ug[i + Q * 1 + Q * 3 * 2]; in f_apply_poisson_vec()
307 const CeedScalar ug2 = ug[i + Q * 2 + Q * 3 * 2]; in f_apply_poisson_vec()
308 vg[i + Q * 0 + Q * 3 * 2] = in f_apply_poisson_vec()
309 qdata[i + Q * 0] * ug0 + qdata[i + Q * 5] * ug1 + qdata[i + Q * 4] * ug2; in f_apply_poisson_vec()
310 vg[i + Q * 1 + Q * 3 * 2] = in f_apply_poisson_vec()
311 qdata[i + Q * 5] * ug0 + qdata[i + Q * 1] * ug1 + qdata[i + Q * 3] * ug2; in f_apply_poisson_vec()
312 vg[i + Q * 2 + Q * 3 * 2] = in f_apply_poisson_vec()
313 qdata[i + Q * 4] * ug0 + qdata[i + Q * 3] * ug1 + qdata[i + Q * 2] * ug2; in f_apply_poisson_vec()