| /libCEED/examples/deal.II/ |
| H A D | bps-qfunctions.h | 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()
[all …]
|
| /libCEED/examples/petsc/qfunctions/bps/ |
| H A D | bp24.h | 19 CEED_QFUNCTION(SetupMassDiffRhs3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co… in SetupMassDiffRhs3() argument
27 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupMassDiffRhs3()
32 true_soln[i + 0 * Q] = in SetupMassDiffRhs3()
33 …sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) * sin(M_PI * (c[1] + k[1] * x[i + Q * 1])) * sin(M_PI * (… in SetupMassDiffRhs3()
35 true_soln[i + 1 * Q] = 2 * true_soln[i + 0 * Q]; in SetupMassDiffRhs3()
37 true_soln[i + 2 * Q] = 3 * true_soln[i + 0 * Q]; in SetupMassDiffRhs3()
40 …rhs[i + 0 * Q] = w[i + Q * 0] * (M_PI * M_PI * (k[0] * k[0] + k[1] * k[1] + k[2] * k[2]) + 1.0) * … in SetupMassDiffRhs3()
42 rhs[i + 1 * Q] = 2 * rhs[i + 0 * Q]; in SetupMassDiffRhs3()
44 rhs[i + 2 * Q] = 3 * rhs[i + 0 * Q]; in SetupMassDiffRhs3()
61 CEED_QFUNCTION(MassDiff3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out… in MassDiff3() argument
[all …]
|
| H A D | bp4.h | 19 CEED_QFUNCTION(SetupDiffRhs3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const … in SetupDiffRhs3() argument
27 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupDiffRhs3()
32 true_soln[i + 0 * Q] = in SetupDiffRhs3()
33 …sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) * sin(M_PI * (c[1] + k[1] * x[i + Q * 1])) * sin(M_PI * (… in SetupDiffRhs3()
35 true_soln[i + 1 * Q] = 2 * true_soln[i + 0 * Q]; in SetupDiffRhs3()
37 true_soln[i + 2 * Q] = 3 * true_soln[i + 0 * Q]; in SetupDiffRhs3()
40 …rhs[i + 0 * Q] = w[i + Q * 0] * M_PI * M_PI * (k[0] * k[0] + k[1] * k[1] + k[2] * k[2]) * true_sol… in SetupDiffRhs3()
42 rhs[i + 1 * Q] = 2 * rhs[i + 0 * Q]; in SetupDiffRhs3()
44 rhs[i + 2 * Q] = 3 * rhs[i + 0 * Q]; in SetupDiffRhs3()
59 CEED_QFUNCTION(Diff3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { in Diff3() argument
[all …]
|
| H A D | bp2.h | 19 CEED_QFUNCTION(SetupMassRhs3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *… in SetupMassRhs3() argument
24 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupMassRhs3()
26 true_soln[i + 0 * Q] = sqrt(x[i] * x[i] + x[i + Q] * x[i + Q] + x[i + 2 * Q] * x[i + 2 * Q]); in SetupMassRhs3()
28 true_soln[i + 1 * Q] = 2 * true_soln[i + 0 * Q]; in SetupMassRhs3()
30 true_soln[i + 2 * Q] = 3 * true_soln[i + 0 * Q]; in SetupMassRhs3()
33 rhs[i + 0 * Q] = w[i] * true_soln[i + 0 * Q]; in SetupMassRhs3()
35 rhs[i + 1 * Q] = 2 * rhs[i + 0 * Q]; in SetupMassRhs3()
37 rhs[i + 2 * Q] = 3 * rhs[i + 0 * Q]; in SetupMassRhs3()
52 CEED_QFUNCTION(Mass3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in Mass3() argument
57 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Mass3()
[all …]
|
| H A D | bp4sphere.h | 19 CEED_QFUNCTION(SetupDiffRhs3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *… in SetupDiffRhs3() argument
30 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupDiffRhs3()
32 CeedScalar x = X[i + Q * 0], y = X[i + Q * 1], z = X[i + Q * 2]; in SetupDiffRhs3()
44 true_soln[i + 0 * Q] = sin(lambda) * cos(theta); in SetupDiffRhs3()
46 true_soln[i + 1 * Q] = 2 * true_soln[i + 0 * Q]; in SetupDiffRhs3()
48 true_soln[i + 2 * Q] = 3 * true_soln[i + 0 * Q]; in SetupDiffRhs3()
51 rhs[i + 0 * Q] = q_data[i + Q * 0] * 2 * sin(lambda) * cos(theta) / (R * R); in SetupDiffRhs3()
53 rhs[i + 1 * Q] = 2 * rhs[i + 0 * Q]; in SetupDiffRhs3()
55 rhs[i + 2 * Q] = 3 * rhs[i + 0 * Q]; in SetupDiffRhs3()
71 CEED_QFUNCTION(Diff3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in Diff3() argument
[all …]
|
| H A D | bp2sphere.h | 19 CEED_QFUNCTION(SetupMassRhs3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *… in SetupMassRhs3() argument
30 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupMassRhs3()
32 const CeedScalar theta = asin(X[i + 2 * Q] / R); in SetupMassRhs3()
36 true_soln[i + 0 * Q] = fabs(theta); in SetupMassRhs3()
38 true_soln[i + 1 * Q] = 2 * true_soln[i + 0 * Q]; in SetupMassRhs3()
40 true_soln[i + 2 * Q] = 3 * true_soln[i + 0 * Q]; in SetupMassRhs3()
43 rhs[i + 0 * Q] = q_data[i] * true_soln[i]; in SetupMassRhs3()
45 rhs[i + 1 * Q] = 2 * rhs[i + 0 * Q]; in SetupMassRhs3()
47 rhs[i + 2 * Q] = 3 * rhs[i + 0 * Q]; in SetupMassRhs3()
63 CEED_QFUNCTION(Mass3)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in Mass3() argument
[all …]
|
| H A D | common.h | 14 CEED_QFUNCTION(Error)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { in Error() argument
17 for (CeedInt i = 0; i < Q; i++) { in Error()
24 CEED_QFUNCTION(Error3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { in Error3() argument
27 for (CeedInt i = 0; i < Q; i++) { in Error3()
28 …error[i + 0 * Q] = (u[i + 0 * Q] - target[i + 0 * Q]) * (u[i + 0 * Q] - target[i + 0 * Q]) * q_dat… in Error3()
29 …error[i + 1 * Q] = (u[i + 1 * Q] - target[i + 1 * Q]) * (u[i + 1 * Q] - target[i + 1 * Q]) * q_dat… in Error3()
30 …error[i + 2 * Q] = (u[i + 2 * Q] - target[i + 2 * Q]) * (u[i + 2 * Q] - target[i + 2 * Q]) * q_dat… in Error3()
|
| H A D | bp13.h | 19 CEED_QFUNCTION(SetupMassDiffRhs)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScala… in SetupMassDiffRhs() argument
27 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupMassDiffRhs()
31 … sin(M_PI * (c[0] + k[0] * x[i + Q * 0])) * sin(M_PI * (c[1] + k[1] * x[i + Q * 1])) * sin(M_PI * … in SetupMassDiffRhs()
33 …rhs[i] = w[i + Q * 0] * (M_PI * M_PI * (k[0] * k[0] + k[1] * k[1] + k[2] * k[2]) + 1.0) * true_sol… in SetupMassDiffRhs()
50 CEED_QFUNCTION(MassDiff)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out)… in MassDiff() argument
55 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in MassDiff()
57 const CeedScalar du[3] = {ug[i + Q * 0], ug[i + Q * 1], ug[i + Q * 2]}; in MassDiff()
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()
[all …]
|
| H A D | bp3.h | 38 CEED_QFUNCTION(SetupDiffGeo)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *… in SetupDiffGeo() argument
47 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupDiffGeo()
61 qd[i + Q * 0] = w[i] * detJ; in SetupDiffGeo()
62 qd[i + Q * 1] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in SetupDiffGeo()
63 qd[i + Q * 2] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in SetupDiffGeo()
64 qd[i + Q * 3] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in SetupDiffGeo()
65 qd[i + Q * 4] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in SetupDiffGeo()
66 qd[i + Q * 5] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in SetupDiffGeo()
67 qd[i + Q * 6] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in SetupDiffGeo()
76 CEED_QFUNCTION(SetupDiffRhs)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *… in SetupDiffRhs() argument
[all …]
|
| H A D | bp3sphere.h | 59 CEED_QFUNCTION(SetupDiffGeo)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *… in SetupDiffGeo() argument
64 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupDiffGeo()
66 const CeedScalar xx[3] = {X[i + 0 * Q], X[i + 1 * Q], X[i + 2 * Q]}; in SetupDiffGeo()
73 {J[i + Q * 0], J[i + Q * 3]}, in SetupDiffGeo()
74 {J[i + Q * 1], J[i + Q * 4]}, in SetupDiffGeo()
75 {J[i + Q * 2], J[i + Q * 5]} in SetupDiffGeo()
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()
[all …]
|
| /libCEED/tests/ |
| H A D | t532-operator.h | 10 CEED_QFUNCTION(setup_mass)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in setup_mass() argument
13 for (CeedInt i = 0; i < Q; i++) { in setup_mass()
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup_mass()
19 CEED_QFUNCTION(setup_diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in setup_diff() argument
31 for (CeedInt i = 0; i < Q; i++) { in setup_diff()
34 const CeedScalar J11 = J[i + Q * 0]; in setup_diff()
35 const CeedScalar J21 = J[i + Q * 1]; in setup_diff()
36 const CeedScalar J12 = J[i + Q * 2]; in setup_diff()
37 const CeedScalar J22 = J[i + Q * 3]; in setup_diff()
39 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup_diff()
[all …]
|
| H A D | t531-operator.h | 10 CEED_QFUNCTION(setup)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in setup() argument
22 for (CeedInt i = 0; i < Q; i++) { in setup()
25 const CeedScalar J11 = J[i + Q * 0]; in setup()
26 const CeedScalar J21 = J[i + Q * 1]; in setup()
27 const CeedScalar J12 = J[i + Q * 2]; in setup()
28 const CeedScalar J22 = J[i + Q * 3]; in setup()
30 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup()
31 qd[i + Q * 1] = w * (J11 * J11 + J21 * J21); in setup()
32 qd[i + Q * 2] = -w * (J11 * J12 + J21 * J22); in setup()
38 CEED_QFUNCTION(diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in diff() argument
[all …]
|
| H A D | t568-operator.h | 10 CEED_QFUNCTION(setup)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in setup() argument
22 for (CeedInt i = 0; i < Q; i++) { in setup()
25 const CeedScalar J11 = J[i + Q * 0]; in setup()
26 const CeedScalar J21 = J[i + Q * 1]; in setup()
27 const CeedScalar J12 = J[i + Q * 2]; in setup()
28 const CeedScalar J22 = J[i + Q * 3]; in setup()
30 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup()
31 qd[i + Q * 2] = w * (J11 * J11 + J21 * J21); in setup()
32 qd[i + Q * 1] = -w * (J11 * J12 + J21 * J22); in setup()
38 CEED_QFUNCTION(diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in diff() argument
[all …]
|
| H A D | t535-operator.h | 10 CEED_QFUNCTION(setup_mass)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in setup_mass() argument
13 for (CeedInt i = 0; i < Q; i++) { in setup_mass()
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup_mass()
19 CEED_QFUNCTION(setup_diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in setup_diff() argument
31 for (CeedInt i = 0; i < Q; i++) { in setup_diff()
34 const CeedScalar J11 = J[i + Q * 0]; in setup_diff()
35 const CeedScalar J21 = J[i + Q * 1]; in setup_diff()
36 const CeedScalar J12 = J[i + Q * 2]; in setup_diff()
37 const CeedScalar J22 = J[i + Q * 3]; in setup_diff()
39 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup_diff()
[all …]
|
| H A D | t534-operator.h | 10 CEED_QFUNCTION(setup)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in setup() argument
22 for (CeedInt i = 0; i < Q; i++) { in setup()
25 const CeedScalar J11 = J[i + Q * 0]; in setup()
26 const CeedScalar J21 = J[i + Q * 1]; in setup()
27 const CeedScalar J12 = J[i + Q * 2]; in setup()
28 const CeedScalar J22 = J[i + Q * 3]; in setup()
30 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup()
31 qd[i + Q * 2] = w * (J11 * J11 + J21 * J21); in setup()
32 qd[i + Q * 1] = -w * (J11 * J12 + J21 * J22); in setup()
38 CEED_QFUNCTION(diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in diff() argument
[all …]
|
| H A D | t597-operator.h | 10 CEED_QFUNCTION(setup)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in setup() argument
22 for (CeedInt i = 0; i < Q; i++) { in setup()
25 const CeedScalar J11 = J[i + Q * 0]; in setup()
26 const CeedScalar J21 = J[i + Q * 1]; in setup()
27 const CeedScalar J12 = J[i + Q * 2]; in setup()
28 const CeedScalar J22 = J[i + Q * 3]; in setup()
30 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup()
31 qd[i + Q * 2] = w * (J11 * J11 + J21 * J21); in setup()
32 qd[i + Q * 1] = -w * (J11 * J12 + J21 * J22); in setup()
38 CEED_QFUNCTION(diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in diff() argument
[all …]
|
| H A D | t541-operator.h | 10 CEED_QFUNCTION(setup_diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in setup_diff() argument
19 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in setup_diff()
23 const CeedScalar J11 = J[i + Q * 0]; in setup_diff()
24 const CeedScalar J21 = J[i + Q * 1]; in setup_diff()
25 const CeedScalar J12 = J[i + Q * 2]; in setup_diff()
26 const CeedScalar J22 = J[i + Q * 3]; in setup_diff()
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()
35 CEED_QFUNCTION(apply)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in apply() argument
[all …]
|
| H A D | t522-operator.h | 10 CEED_QFUNCTION(setup)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in setup() argument
14 for (CeedInt i = 0; i < Q; i++) { in setup()
18 const CeedScalar J11 = J[i + Q * 0]; in setup()
19 const CeedScalar J21 = J[i + Q * 1]; in setup()
20 const CeedScalar J12 = J[i + Q * 2]; in setup()
21 const CeedScalar J22 = J[i + Q * 3]; in setup()
23 qd[i + Q * 0] = w * (J12 * J12 + J22 * J22); in setup()
24 qd[i + Q * 1] = w * (J11 * J11 + J21 * J21); in setup()
25 qd[i + Q * 2] = -w * (J11 * J12 + J21 * J22); in setup()
30 CEED_QFUNCTION(diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in diff() argument
[all …]
|
| H A D | t539-operator.h | 10 CEED_QFUNCTION(apply)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *o… in apply() argument
26 for (CeedInt i = 0; i < Q; i++) { in apply()
29 v_0[i + Q * c] = qd_mass[i] * (c + 1) * u_0[i + Q * c]; in apply()
31 dv_0[i + Q * (0 * num_comp_0 + c)] = in apply()
32 …qd_diff[i + Q * 0] * (c + 1) * du_0[i + Q * (0 * num_comp_0 + c)] + qd_diff[i + Q * 2] * du_0[i + … in apply()
33 dv_0[i + Q * (1 * num_comp_0 + c)] = in apply()
34 …qd_diff[i + Q * 2] * (c + 1) * du_0[i + Q * (0 * num_comp_0 + c)] + qd_diff[i + Q * 1] * du_0[i + … in apply()
|
| /libCEED/examples/nek/bps/ |
| H A D | bps.h | 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()
[all …]
|
| /libCEED/examples/mfem/ |
| H A D | bp3.h | 16 CEED_QFUNCTION(f_build_diff)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *c… in f_build_diff() argument
28 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = w[i] / J[i]; } in f_build_diff()
32 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in f_build_diff()
35 const CeedScalar J11 = J[i + Q * 0]; in f_build_diff()
36 const CeedScalar J21 = J[i + Q * 1]; in f_build_diff()
37 const CeedScalar J12 = J[i + Q * 2]; in f_build_diff()
38 const CeedScalar J22 = J[i + Q * 3]; in f_build_diff()
40 qdata[i + Q * 0] = qw * (J12 * J12 + J22 * J22); in f_build_diff()
41 qdata[i + Q * 1] = qw * (J11 * J11 + J21 * J21); in f_build_diff()
42 qdata[i + Q * 2] = -qw * (J11 * J12 + J21 * J22); in f_build_diff()
[all …]
|
| H A D | bp1.h | 16 CEED_QFUNCTION(f_build_mass)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *c… in f_build_mass() argument
26 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = J[i] * w[i]; } in f_build_mass()
30 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in f_build_mass()
33 qdata[i] = (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]) * w[i]; in f_build_mass()
38 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in f_build_mass()
42 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()
43 J[i + Q * 1] * (J[i + Q * 3] * J[i + Q * 8] - J[i + Q * 5] * J[i + Q * 6]) + in f_build_mass()
44 J[i + Q * 2] * (J[i + Q * 3] * J[i + Q * 7] - J[i + Q * 4] * J[i + Q * 6])) * in f_build_mass()
53 CEED_QFUNCTION(f_apply_mass)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *c… in f_apply_mass() argument
58 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = qdata[i] * u[i]; } in f_apply_mass()
|
| /libCEED/python/tests/ |
| H A D | buildmats.py | 5 P, Q, dim = 6, 4, 2
6 interp = np.empty(P * Q, dtype=mat_dtype)
7 grad = np.empty(dim * P * Q, dtype=mat_dtype)
23 for i in range(Q):
24 x1 = qref[0 * Q + i]
25 x2 = qref[1 * Q + i]
36 grad[(i + Q) * P + 0] = 2. * \
39 grad[(i + Q) * P + 1] = -4. * (x1 * 1.)
41 grad[(i + Q) * P + 2] = 2. * 0.
43 grad[(i + Q) * P + 3] = -4. * (1. * (x1 + x2 - 1.) + x2 * 1.)
[all …]
|
| /libCEED/include/ceed/jit-source/magma/ |
| H A D | magma-common-nontensor.h | 19 template <typename T, int P, int Q, int BY>
20 …line__ void read_A_notrans_g2r_1D_nosync(const int tx, const int ty, const T *dA, T *sA, T rA[Q]) { in read_A_notrans_g2r_1D_nosync() argument
25 for (i = 0; i < P * Q - P * BY; i += P * BY) { in read_A_notrans_g2r_1D_nosync()
28 if (i + tid < P * Q) { in read_A_notrans_g2r_1D_nosync()
34 for (int j = 0; j < Q; j++) { in read_A_notrans_g2r_1D_nosync()
44 template <typename T, int P, int Q, int BY>
45 …inline__ void read_A_trans_g2r_1D_nosync(const int tx, const int ty, const T *dA, T *sA, T rA[Q]) { in read_A_trans_g2r_1D_nosync() argument
50 for (i = 0; i < P * Q - P * BY; i += P * BY) { in read_A_trans_g2r_1D_nosync()
53 if (i + tid < P * Q) { in read_A_trans_g2r_1D_nosync()
59 for (int j = 0; j < Q; j++) { in read_A_trans_g2r_1D_nosync()
[all …]
|
| /libCEED/examples/fluids/qfunctions/ |
| H A D | bc_freestream.h | 22 CEED_QFUNCTION_HELPER int Freestream(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar … in Freestream() argument
33 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Freestream()
38 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, normal); in Freestream()
56 StoredValuesPack(Q, i, 0, 5, qi, jac_data_sur); in Freestream()
57 StoredValuesPack(Q, i, 5, 6, zeros, jac_data_sur); // Every output value must be set in Freestream()
63 CEED_QFUNCTION(Freestream_Conserv_HLL)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScala… in Freestream_Conserv_HLL() argument
64 return Freestream(ctx, Q, in, out, STATEVAR_CONSERVATIVE, RIEMANN_HLL); in Freestream_Conserv_HLL()
67 CEED_QFUNCTION(Freestream_Prim_HLL)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in Freestream_Prim_HLL() argument
68 return Freestream(ctx, Q, in, out, STATEVAR_PRIMITIVE, RIEMANN_HLL); in Freestream_Prim_HLL()
71 CEED_QFUNCTION(Freestream_Entropy_HLL)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScala… in Freestream_Entropy_HLL() argument
[all …]
|