Home
last modified time | relevance | path

Searched refs:qdata (Results 1 – 25 of 51) sorted by relevance

123

/libCEED/examples/deal.II/
H A Dbps-qfunctions.h40 CeedScalar *qdata = out[0]; in f_build_mass() local
47 qdata[i] = J[i] * w[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()
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()
141 CeedScalar *qdata = out[0]; in f_build_poisson() local
148 qdata[i] = w[i] / J[i]; 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()
186 qdata[i + Q * 0] = qw * (A11 * A11 + A12 * A12 + A13 * A13); in f_build_poisson()
[all …]
/libCEED/examples/rust/ex3-volume/src/
H A Dmain.rs139 [qdata, ..]: QFunctionOutputs| { in example_3()
143 let q = qdata.len() / 2; in example_3()
146 qdata[i + q * 0] = weights[i] * jacobian[i]; in example_3()
148 qdata[i + q * 1] = weights[i] / jacobian[i]; in example_3()
152 let q = qdata.len() / 4; in example_3()
159 qdata[i + q * 0] = weights[i] * (j11 * j22 - j21 * j12); in example_3()
162 qdata[i + q * 1] = qw * (j12 * j12 + j22 * j22); in example_3()
163 qdata[i + q * 2] = qw * (j11 * j11 + j21 * j21); in example_3()
164 qdata[i + q * 3] = -qw * (j11 * j12 + j21 * j22); in example_3()
168 let q = qdata.len() / 7; in example_3()
[all …]
/libCEED/examples/mfem/
H A Dbp3.h23 CeedScalar *qdata = out[0]; in f_build_diff() local
28 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = w[i] / J[i]; } 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()
70 qdata[i + Q * 0] = qw * (A11 * A11 + A12 * A12 + A13 * A13); in f_build_diff()
71 qdata[i + Q * 1] = qw * (A21 * A21 + A22 * A22 + A23 * A23); in f_build_diff()
72 qdata[i + Q * 2] = qw * (A31 * A31 + A32 * A32 + A33 * A33); in f_build_diff()
73 qdata[i + Q * 3] = qw * (A21 * A31 + A22 * A32 + A23 * A33); in f_build_diff()
74 qdata[i + Q * 4] = qw * (A11 * A31 + A12 * A32 + A13 * A33); in f_build_diff()
[all …]
H A Dbp1.h21 CeedScalar *qdata = out[0]; in f_build_mass() local
26 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = J[i] * w[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()
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()
54 const CeedScalar *u = in[0], *qdata = in[1]; in f_apply_mass() local
58 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = qdata[i] * u[i]; } in f_apply_mass()
/libCEED/examples/rust/ex3-volume-vector/src/
H A Dmain.rs141 [qdata, ..]: QFunctionOutputs| { in example_3_vector()
145 let q = qdata.len() / 2; in example_3_vector()
148 qdata[i + q * 0] = weights[i] * jacobian[i]; in example_3_vector()
150 qdata[i + q * 1] = weights[i] / jacobian[i]; in example_3_vector()
154 let q = qdata.len() / 4; in example_3_vector()
161 qdata[i + q * 0] = weights[i] * (j11 * j22 - j21 * j12); in example_3_vector()
164 qdata[i + q * 1] = qw * (j12 * j12 + j22 * j22); in example_3_vector()
165 qdata[i + q * 2] = qw * (j11 * j11 + j21 * j21); in example_3_vector()
166 qdata[i + q * 3] = -qw * (j11 * j12 + j21 * j22); in example_3_vector()
170 let q = qdata.len() / 7; in example_3_vector()
[all …]
/libCEED/examples/nek/bps/
H A Dbps.h22 CeedScalar *qdata = out[0], *rhs = out[1]; in masssetupf() local
32 qdata[i] = det * w[i]; 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()
40 const CeedScalar *qdata = in[1]; in massf() local
44 for (CeedInt i = 0; i < Q; i++) v[i] = qdata[i] * u[i]; in massf()
55 CeedScalar *qdata = out[0], *rhs = out[1]; in diffsetupf() local
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()
[all …]
/libCEED/examples/rust/ex2-surface/src/
H A Dmain.rs144 [qdata, ..]: QFunctionOutputs| { in example_2()
147 1 => qdata in example_2()
150 .for_each(|(qdata, (j, weight))| *qdata = weight / j), in example_2()
152 let q = qdata.len() / 3; in example_2()
159 qdata[i + q * 0] = qw * (j12 * j12 + j22 * j22); in example_2()
160 qdata[i + q * 1] = qw * (j11 * j11 + j21 * j21); in example_2()
161 qdata[i + q * 2] = -qw * (j11 * j12 + j21 * j22); in example_2()
165 let q = qdata.len() / 6; in example_2()
180 qdata[i + q * 0] = qw in example_2()
184 qdata[i + q * 1] = qw in example_2()
[all …]
/libCEED/examples/rust/ex2-surface-vector/src/
H A Dmain.rs147 [qdata, ..]: QFunctionOutputs| { in example_2_vector()
150 1 => qdata in example_2_vector()
153 .for_each(|(qdata, (j, weight))| *qdata = weight / j), in example_2_vector()
155 let q = qdata.len() / 3; in example_2_vector()
162 qdata[i + q * 0] = qw * (j12 * j12 + j22 * j22); in example_2_vector()
163 qdata[i + q * 1] = qw * (j11 * j11 + j21 * j21); in example_2_vector()
164 qdata[i + q * 2] = -qw * (j11 * j12 + j21 * j22); in example_2_vector()
168 let q = qdata.len() / 6; in example_2_vector()
183 qdata[i + q * 0] = qw in example_2_vector()
187 qdata[i + q * 1] = qw in example_2_vector()
[all …]
/libCEED/python/tests/
H A Dtest-qfunctions.h17 CeedScalar *qdata = out[0]; in setup_mass() local
20 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = J[i] * w[i]; } in setup_mass()
34 CeedScalar *qdata = out[0]; in setup_mass_2d() local
37 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { qdata[i] = w[i] * (J[i + Q * 0] * J[i + Q * 3] - … in setup_mass_2d()
52 const CeedScalar *qdata = in[0], *u = in[1]; in apply_mass() local
58 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = scale * qdata[i] * u[i]; } in apply_mass()
69 const CeedScalar *qdata = in[0], *u = in[1]; in apply_mass_two() local
76 v[i] = qdata[i] * u[i]; in apply_mass_two()
77 v[Q + i] = qdata[i] * u[Q + i]; in apply_mass_two()
H A Dtest-4-qfunction.py77 qdata = ceed.Vector(q)
78 qdata.set_value(0)
81 outputs = [qdata]
84 inputs = [qdata, u]
139 qdata = ceed.Vector(q)
140 qdata.set_value(0)
143 outputs = [qdata]
146 inputs = [qdata, u]
226 qdata = ceed.Vector(q)
227 qdata.set_value(0)
[all …]
/libCEED/examples/ceed/
H A Dex1-volume-f.h3 qdata, v2, v3, v4, v5, v6, v7, v8, v9, v10, v11, v12, v13, v14, v15, v16, ierr)
10 ! qdata is quadrature data with shape [1, Q]
11 real*8 qdata(1)
20 qdata(i) = j(i) * w(i)
25 qdata(i) = (j(0*q + i)*j(3*q + i) - j(1*q + i)*j(2*q + i)) * w(i)
30 qdata(i) = (j(0*q + i) * (j(4*q + i)*j(8*q + i) - j(5*q + i)*j(7*q + i)) -&
39 subroutine apply_mass(ctx, q, u, qdata, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12, u13, u14, u15, u…
43 ! qdata is quadrature data with shape [1, Q]
45 real*8 qdata(1)
51 v(i) = qdata(i) * u(i)
/libCEED/examples/rust/ex1-volume/src/
H A Dmain.rs135 [qdata, ..]: QFunctionOutputs| { in example_1()
138 1 => qdata in example_1()
141 .for_each(|(qdata, (j, weight))| *qdata = j * weight), in example_1()
143 let q = qdata.len(); in example_1()
144 qdata.iter_mut().zip(weights.iter()).enumerate().for_each( in example_1()
145 |(i, (qdata, weight))| { in example_1()
146 *qdata = (jacobian[i + q * 0] * jacobian[i + q * 3] in example_1()
153 let q = qdata.len(); in example_1()
154 qdata.iter_mut().zip(weights.iter()).enumerate().for_each( in example_1()
155 |(i, (qdata, weight))| { in example_1()
[all …]
/libCEED/examples/rust/ex1-volume-vector/src/
H A Dmain.rs139 [qdata, ..]: QFunctionOutputs| { in example_1_vector()
142 1 => qdata in example_1_vector()
145 .for_each(|(qdata, (j, weight))| *qdata = j * weight), in example_1_vector()
147 let q = qdata.len(); in example_1_vector()
148 qdata.iter_mut().zip(weights.iter()).enumerate().for_each( in example_1_vector()
149 |(i, (qdata, weight))| { in example_1_vector()
150 *qdata = (jacobian[i + q * 0] * jacobian[i + q * 3] in example_1_vector()
157 let q = qdata.len(); in example_1_vector()
158 qdata.iter_mut().zip(weights.iter()).enumerate().for_each( in example_1_vector()
159 |(i, (qdata, weight))| { in example_1_vector()
[all …]
/libCEED/julia/LibCEED.jl/examples/
H A Dex3-volume.jl83 (qdata, :out, EVAL_NONE, num_q_comp), # ← qdata output
89 qdata[1] = weights*det_J
95 qdata[idx] = dx[:, i]'*dx[:, j]
109 (:qdata, qdata_rstr_i, BasisNone(), CeedVectorActive()),
116 qdata = CeedVector(ceed, num_elem*elem_qpts*num_q_comp)
119 apply!(build_oper, mesh_coords, qdata)
128 (qdata, :in, EVAL_NONE, num_q_comp),
132 # Apply mass: v = qdata[1] * u
133 v .= qdata[1].*u
135 # Apply diffusion: dv = (qdata[2:end]) * du
[all …]
H A Dex2-surface.jl78 (qdata, :out, EVAL_NONE, dim*(dim + 1)÷2),
81 qdata .= setvoigt(w*det(J)*Jinv*Jinv')
96 (:qdata, qdata_rstr_i, BasisNone(), CeedVectorActive()),
103 qdata = CeedVector(ceed, num_elem*elem_qpts*div(dim*(dim + 1), 2))
106 apply!(build_oper, mesh_coords, qdata)
115 (qdata, :in, EVAL_NONE, dim*(dim + 1)÷2),
118 dXdxdXdxT = getvoigt(qdata)
132 (:qdata, qdata_rstr_i, BasisNone(), qdata),
H A Dex1-volume.jl75 (qdata, :out, EVAL_NONE),
77 qdata .= w*det(J)
91 (:qdata, sol_rstr_i, BasisNone(), CeedVectorActive()),
98 qdata = CeedVector(ceed, num_elem*elem_qpts)
102 apply!(build_oper, mesh_coords, qdata)
110 (qdata, :in, EVAL_NONE),
113 v .= qdata*u
126 (:qdata, sol_rstr_i, BasisNone(), qdata),
/libCEED/tests/
H A Dt410-qfunction-f.f907 integer qdata,j,w,u,v variable
45 call ceedvectorcreate(ceed,q,qdata,err)
46 call ceedvectorsetvalue(qdata,0.d0,err)
51 &qdata,ceed_null,ceed_null,ceed_null,ceed_null,ceed_null,&
55 call ceedqfunctionapply(qf_mass,q,u,qdata,ceed_null,ceed_null,&
75 call ceedvectordestroy(qdata,err)
H A Dt400-qfunction-f.f9012 integer qdata,w,u,v
57 call ceedvectorcreate(ceed,q,qdata,err)
58 call ceedvectorsetvalue(qdata,0.d0,err)
63 &qdata,ceed_null,ceed_null,ceed_null,ceed_null,ceed_null,&
67 call ceedqfunctionapply(qf_mass,q,u,qdata,ceed_null,ceed_null,&
87 call ceedvectordestroy(qdata,err)
H A Dt401-qfunction-f.f9012 integer qdata,w,u,v
69 call ceedvectorcreate(ceed,q,qdata,err)
70 call ceedvectorsetvalue(qdata,0.d0,err)
75 &qdata,ceed_null,ceed_null,ceed_null,ceed_null,ceed_null,&
79 call ceedqfunctionapply(qf_mass,q,u,qdata,ceed_null,ceed_null,&
99 call ceedvectordestroy(qdata,err)
H A Dt530-operator-f.f9017 integer qdata,x,a,u,v
54 call ceedvectorcreate(ceed,nqpts,qdata,err)
114 ceed_basis_none,qdata,err)
121 call ceedoperatorapply(op_setup,x,qdata,ceed_request_immediate,err)
129 call ceedvectorgetarrayread(qdata,ceed_mem_host,qq,qoffset,err)
139 call ceedvectorrestorearrayread(qdata,qq,qoffset,err)
163 call ceedvectorsetarray(qdata,ceed_mem_host,ceed_copy_values,aa,&
199 call ceedvectordestroy(qdata,err)
H A Dt500-operator-f.f9017 integer qdata,x,u,v
87 call ceedvectorcreate(ceed,nelem*q,qdata,err)
96 ceed_basis_none,qdata,err)
102 call ceedoperatorapply(op_setup,x,qdata,ceed_request_immediate,err)
119 call ceedvectordestroy(qdata,err)
H A Dt503-operator-f.f9017 integer qdata,x,u,v
52 call ceedvectorcreate(ceed,nelem*q,qdata,err)
99 ceed_basis_none,qdata,err)
101 ceed_basis_none,qdata,err)
122 call ceedvectordestroy(qdata,err)
/libCEED/tests/output/
H A Dt413-qfunction-f.out13 Name: "qdata"
23 Name: "qdata"
38 Name: "qdata"
H A Dt402-qfunction-f.out9 Name: "qdata"
15 Name: "qdata"
32 Name: "qdata"
H A Dt413-qfunction.out13 Name: "qdata"
23 Name: "qdata"
38 Name: "qdata"

123