| /libCEED/examples/fluids/qfunctions/ |
| H A D | grid_anisotropy_tensor.h | 39 const CeedScalar(*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in AnisotropyTensorProjection() local
43 const CeedScalar wdetJ = q_data[0][i]; in AnisotropyTensorProjection()
45 {q_data[1][i], q_data[2][i], q_data[3][i]}, in AnisotropyTensorProjection()
46 {q_data[4][i], q_data[5][i], q_data[6][i]}, in AnisotropyTensorProjection()
47 {q_data[7][i], q_data[8][i], q_data[9][i]} in AnisotropyTensorProjection()
63 const CeedScalar(*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in AnisotropyTensorCollocate() local
68 {q_data[1][i], q_data[2][i], q_data[3][i]}, in AnisotropyTensorCollocate()
69 {q_data[4][i], q_data[5][i], q_data[6][i]}, in AnisotropyTensorCollocate()
70 {q_data[7][i], q_data[8][i], q_data[9][i]} in AnisotropyTensorCollocate()
|
| /libCEED/examples/python/qfunctions/ |
| H A D | ex3-volume.h | 17 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_mass_diff() local
28 q_data[0][i] = w[i] * J[0][0][i]; in build_mass_diff()
31 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
47 q_data[0][i] = w[i] * (J00 * J11 - J10 * J01); in build_mass_diff()
50 q_data[1][i] = qw * (J01 * J01 + J11 * J11); in build_mass_diff()
51 q_data[2][i] = qw * (J00 * J00 + J10 * J10); in build_mass_diff()
52 q_data[3][i] = -qw * (J00 * J01 + J10 * J11); in build_mass_diff()
72 q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
79 q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_mass_diff()
80 q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_mass_diff()
[all …]
|
| H A D | ex2-surface.h | 17 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_diff() local
26 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { q_data[0][i] = w[i] / J[0][0][i]; } // End of Qu… in build_diff()
40 q_data[0][i] = qw * (J01 * J01 + J11 * J11); in build_diff()
41 q_data[1][i] = qw * (J00 * J00 + J10 * J10); in build_diff()
42 q_data[2][i] = -qw * (J00 * J01 + J10 * J11); in build_diff()
69 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_diff()
70 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_diff()
71 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_diff()
72 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_diff()
73 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_diff()
[all …]
|
| H A D | ex1-volume.h | 17 CeedScalar *q_data = out[0]; in build_mass() local
25 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { q_data[i] = J[0][0][i] * w[i]; } // End of Quadr… in build_mass()
32 q_data[i] = (J[0][0][i] * J[1][1][i] - J[0][1][i] * J[1][0][i]) * w[i]; in build_mass()
40 q_data[i] = in build_mass()
54 const CeedScalar *u = in[0], *q_data = in[1]; in apply_mass() local
58 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = q_data[i] * u[i]; } // End of Quadrature … in apply_mass()
|
| /libCEED/examples/ceed/ |
| H A D | ex3-volume.h | 22 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_mass_diff() local
32 q_data[0][i] = w[i] * J[0][0][i]; in build_mass_diff()
35 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
51 q_data[0][i] = w[i] * (J00 * J11 - J10 * J01); in build_mass_diff()
54 q_data[1][i] = qw * (J01 * J01 + J11 * J11); in build_mass_diff()
55 q_data[2][i] = qw * (J00 * J00 + J10 * J10); in build_mass_diff()
56 q_data[3][i] = -qw * (J00 * J01 + J10 * J11); in build_mass_diff()
79 q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
86 q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_mass_diff()
87 q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_mass_diff()
[all …]
|
| H A D | ex2-surface.h | 22 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_diff() local
30 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { q_data[0][i] = w[i] / J[0][0][i]; } // End of Qu… in build_diff()
44 q_data[0][i] = qw * (J01 * J01 + J11 * J11); in build_diff()
45 q_data[1][i] = qw * (J00 * J00 + J10 * J10); in build_diff()
46 q_data[2][i] = -qw * (J00 * J01 + J10 * J11); in build_diff()
73 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_diff()
74 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_diff()
75 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_diff()
76 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_diff()
77 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_diff()
[all …]
|
| /libCEED/examples/solids/qfunctions/ |
| H A D | linear.h | 30 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasResidual_Linear() local
51 const CeedScalar wdetJ = q_data[0][i]; in ElasResidual_Linear()
53 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasResidual_Linear()
54 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasResidual_Linear()
55 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasResidual_Linear()
127 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasJacobian_Linear() local
148 const CeedScalar wdetJ = q_data[0][i]; in ElasJacobian_Linear()
150 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasJacobian_Linear()
151 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasJacobian_Linear()
152 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasJacobian_Linear()
[all …]
|
| H A D | common.h | 39 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in SetupGeo() local
67 q_data[0][i] = w[i] * detJ; in SetupGeo()
71 q_data[1][i] = A11 / detJ; in SetupGeo()
72 q_data[2][i] = A12 / detJ; in SetupGeo()
73 q_data[3][i] = A13 / detJ; in SetupGeo()
74 q_data[4][i] = A21 / detJ; in SetupGeo()
75 q_data[5][i] = A22 / detJ; in SetupGeo()
76 q_data[6][i] = A23 / detJ; in SetupGeo()
77 q_data[7][i] = A31 / detJ; in SetupGeo()
78 q_data[8][i] = A32 / detJ; in SetupGeo()
[all …]
|
| H A D | finite-strain-neo-hookean.h | 149 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSResidual_NH() local
188 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_NH()
190 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasFSResidual_NH()
191 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasFSResidual_NH()
192 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasFSResidual_NH()
256 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasFSJacobian_NH() local
283 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_NH()
285 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasFSJacobian_NH()
286 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasFSJacobian_NH()
287 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasFSJacobian_NH()
[all …]
|
| H A D | finite-strain-mooney-rivlin.h | 160 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSResidual_MR() local
190 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_MR()
192 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasFSResidual_MR()
193 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasFSResidual_MR()
194 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasFSResidual_MR()
259 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasFSJacobian_MR() local
281 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_MR()
283 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasFSJacobian_MR()
284 {q_data[4][i], q_data[5][i], q_data[6][i]}, in ElasFSJacobian_MR()
285 {q_data[7][i], q_data[8][i], q_data[9][i]} in ElasFSJacobian_MR()
[all …]
|
| /libCEED/include/ceed/jit-source/gallery/ |
| H A D | ceed-poisson3dapply.h | 16 …const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*q_data)[CEED_Q_VLA] … in Poisson3DApply() local
30 {q_data[0][i], q_data[5][i], q_data[4][i]}, in Poisson3DApply()
31 {q_data[5][i], q_data[1][i], q_data[3][i]}, in Poisson3DApply()
32 {q_data[4][i], q_data[3][i], q_data[2][i]} in Poisson3DApply()
|
| H A D | ceed-vectorpoisson3dapply.h | 16 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in Vector3Poisson3DApply() local
30 {q_data[0][i], q_data[5][i], q_data[4][i]}, in Vector3Poisson3DApply()
31 {q_data[5][i], q_data[1][i], q_data[3][i]}, in Vector3Poisson3DApply()
32 {q_data[4][i], q_data[3][i], q_data[2][i]} in Vector3Poisson3DApply()
|
| H A D | ceed-poisson3dbuild.h | 19 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Poisson3DBuild() local
42 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in Poisson3DBuild()
43 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in Poisson3DBuild()
44 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in Poisson3DBuild()
45 q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in Poisson3DBuild()
46 q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in Poisson3DBuild()
47 q_data[5][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in Poisson3DBuild()
|
| H A D | ceed-poisson2dapply.h | 16 …const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*q_data)[CEED_Q_VLA] … in Poisson2DApply() local
29 {q_data[0][i], q_data[2][i]}, in Poisson2DApply()
30 {q_data[2][i], q_data[1][i]} in Poisson2DApply()
|
| H A D | ceed-vectorpoisson2dapply.h | 16 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in Vector3Poisson2DApply() local
29 {q_data[0][i], q_data[2][i]}, in Vector3Poisson2DApply()
30 {q_data[2][i], q_data[1][i]} in Vector3Poisson2DApply()
|
| /libCEED/examples/petsc/qfunctions/bps/ |
| H A D | bp13.h | 51 const CeedScalar *u = in[0], *ug = in[1], *q_data = in[2]; in MassDiff() local
60 {q_data[i + 1 * Q], q_data[i + 2 * Q], q_data[i + 3 * Q]}, in MassDiff()
61 {q_data[i + 2 * Q], q_data[i + 4 * Q], q_data[i + 5 * Q]}, in MassDiff()
62 {q_data[i + 3 * Q], q_data[i + 5 * Q], q_data[i + 6 * Q]} in MassDiff()
66 v[i] = q_data[i + 0 * Q] * u[i]; in MassDiff()
|
| H A D | bp24.h | 62 const CeedScalar *u = in[0], *ug = in[1], *q_data = in[2]; in MassDiff3() local
75 {q_data[i + 1 * Q], q_data[i + 2 * Q], q_data[i + 3 * Q]}, in MassDiff3()
76 {q_data[i + 2 * Q], q_data[i + 4 * Q], q_data[i + 5 * Q]}, in MassDiff3()
77 {q_data[i + 3 * Q], q_data[i + 5 * Q], q_data[i + 6 * Q]} in MassDiff3()
82 v[i + k * Q] = q_data[i + 0 * Q] * u[i + k * Q]; in MassDiff3()
|
| H A D | bp4.h | 60 const CeedScalar *ug = in[0], *q_data = in[1]; in Diff3() local
73 {q_data[i + 1 * Q], q_data[i + 2 * Q], q_data[i + 3 * Q]}, in Diff3()
74 {q_data[i + 2 * Q], q_data[i + 4 * Q], q_data[i + 5 * Q]}, in Diff3()
75 {q_data[i + 3 * Q], q_data[i + 5 * Q], q_data[i + 6 * Q]} in Diff3()
|
| H A D | bp3sphere.h | 61 CeedScalar *q_data = out[0]; in SetupDiffGeo() local
109 q_data[i + Q * 0] = mod_J * w[i]; in SetupDiffGeo()
130 q_data[i + Q * 1] = dxdXTdxdX_inv[0][0]; in SetupDiffGeo()
131 q_data[i + Q * 2] = dxdXTdxdX_inv[1][1]; in SetupDiffGeo()
132 q_data[i + Q * 3] = dxdXTdxdX_inv[0][1]; in SetupDiffGeo()
144 const CeedScalar *X = in[0], *q_data = in[1]; in SetupDiffRhs() local
167 rhs[i + Q * 0] = q_data[i + Q * 0] * 2 * sin(lambda) * cos(theta) / (R * R); in SetupDiffRhs()
185 const CeedScalar *ug = in[0], *q_data = in[1]; in Diff() local
194 const CeedScalar w_det_J = q_data[i + Q * 0]; in Diff()
198 {q_data[i + Q * 1], q_data[i + Q * 3]}, in Diff()
[all …]
|
| H A D | bp4sphere.h | 21 const CeedScalar *X = in[0], *q_data = in[1]; in SetupDiffRhs3() local
51 rhs[i + 0 * Q] = q_data[i + Q * 0] * 2 * sin(lambda) * cos(theta) / (R * R); in SetupDiffRhs3()
72 const CeedScalar *ug = in[0], *q_data = in[1]; in Diff3() local
84 const CeedScalar w_det_J = q_data[i + Q * 0]; in Diff3()
88 {q_data[i + Q * 1], q_data[i + Q * 3]}, in Diff3()
89 {q_data[i + Q * 3], q_data[i + Q * 2]} in Diff3()
|
| H A D | common.h | 15 const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; in Error() local
18 error[i] = (u[i] - target[i]) * (u[i] - target[i]) * q_data[i]; in Error()
25 const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; in Error3() local
28 …r[i + 0 * Q] = (u[i + 0 * Q] - target[i + 0 * Q]) * (u[i + 0 * Q] - target[i + 0 * Q]) * q_data[i]; in Error3()
29 …r[i + 1 * Q] = (u[i + 1 * Q] - target[i + 1 * Q]) * (u[i + 1 * Q] - target[i + 1 * Q]) * q_data[i]; in Error3()
30 …r[i + 2 * Q] = (u[i + 2 * Q] - target[i + 2 * Q]) * (u[i + 2 * Q] - target[i + 2 * Q]) * q_data[i]; in Error3()
|
| H A D | bp3.h | 106 const CeedScalar *ug = in[0], *q_data = in[1]; in Diff() local
115 {q_data[i + 1 * Q], q_data[i + 2 * Q], q_data[i + 3 * Q]}, in Diff()
116 {q_data[i + 2 * Q], q_data[i + 4 * Q], q_data[i + 5 * Q]}, in Diff()
117 {q_data[i + 3 * Q], q_data[i + 5 * Q], q_data[i + 6 * Q]} in Diff()
|
| H A D | bp2sphere.h | 21 const CeedScalar *X = in[0], *q_data = in[1]; in SetupMassRhs3() local
43 rhs[i + 0 * Q] = q_data[i] * true_soln[i]; in SetupMassRhs3()
64 const CeedScalar *u = in[0], *q_data = in[1]; in Mass3() local
70 v[i + 0 * Q] = q_data[i] * u[i + 0 * Q]; in Mass3()
72 v[i + 1 * Q] = q_data[i] * u[i + 1 * Q]; in Mass3()
74 v[i + 2 * Q] = q_data[i] * u[i + 2 * Q]; in Mass3()
|
| /libCEED/tests/ |
| H A D | t541-operator.h | 16 CeedScalar *q_data = out[0]; in setup_diff() local
28 q_data[i + Q * 0] = qw * (J12 * J12 + J22 * J22); in setup_diff()
29 q_data[i + Q * 1] = qw * (J11 * J11 + J21 * J21); in setup_diff()
30 q_data[i + Q * 2] = -qw * (J11 * J12 + J21 * J22); in setup_diff()
38 const CeedScalar *ug = in[0], *q_data = in[1]; in apply() local
53 {q_data[i + 0 * Q], q_data[i + 2 * Q]}, in apply()
54 {q_data[i + 2 * Q], q_data[i + 1 * Q]} in apply()
|
| H A D | t567-operator.h | 12 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in setup() local
24 q_data[0][i] = qw * (J01 * J01 + J11 * J11); in setup()
25 q_data[1][i] = qw * (J00 * J00 + J10 * J10); in setup()
26 q_data[2][i] = -qw * (J00 * J01 + J10 * J11); in setup()
33 …const CeedScalar(*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[2][CEED_Q_VL… in diff() local
50 {q_data[0][i], q_data[2][i]}, in diff()
51 {q_data[2][i], q_data[1][i]} in diff()
|