Home
last modified time | relevance | path

Searched refs:num_qpts (Results 1 – 25 of 90) sorted by relevance

1234

/libCEED/tests/
H A Dt414-qfunction.c20 CeedInt num_qpts = CeedIntPow(q, dim); in main() local
22 CeedVectorCreate(ceed, num_qpts * dim * dim, &dx); in main()
23 CeedVectorCreate(ceed, num_qpts, &w); in main()
24 CeedVectorCreate(ceed, num_qpts * num_comp, &u); in main()
26 CeedScalar dx_array[num_qpts * dim * dim], w_array[num_qpts], u_array[num_qpts * num_comp]; in main()
28 for (CeedInt i = 0; i < num_qpts; i++) { in main()
29 w_array[i] = 1.0 / num_qpts; in main()
33 for (CeedInt i = 0; i < num_qpts; i++) { in main()
34 dx_array[i + (g * dim + d) * num_qpts] = d == g; in main()
39 for (CeedInt i = 0; i < num_qpts; i++) { in main()
[all …]
H A Dt415-qfunction.c20 CeedInt num_qpts = CeedIntPow(q, dim); in main() local
22 CeedVectorCreate(ceed, num_qpts * dim * dim, &dx); in main()
23 CeedVectorCreate(ceed, num_qpts, &w); in main()
24 CeedVectorCreate(ceed, num_qpts * dim * num_comp, &du); in main()
26 … CeedScalar dx_array[num_qpts * dim * dim], w_array[num_qpts], du_array[num_qpts * dim * num_comp]; in main()
28 for (CeedInt i = 0; i < num_qpts; i++) { in main()
29 w_array[i] = 1.0 / num_qpts; in main()
33 for (CeedInt i = 0; i < num_qpts; i++) { in main()
34 dx_array[i + (g * dim + d) * num_qpts] = d == g; in main()
40 for (CeedInt i = 0; i < num_qpts; i++) { in main()
[all …]
H A Dt332-basis.c13 const CeedInt p = 8, q = 3, dim = 2, num_qpts = q * q; in main() local
15 CeedScalar q_ref[dim * num_qpts], q_weights[num_qpts]; in main()
16 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main()
21 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main()
28 for (CeedInt i = 0; i < p * num_qpts; i++) { in main()
35 CeedVectorCreate(ceed, num_qpts, &v); in main()
44 for (CeedInt i = 0; i < num_qpts; i++) { in main()
60 …if (fabs(num_qpts * 0.25 - u_array[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", 2.0, u_array[i]… in main()
H A Dt331-basis.c13 const CeedInt p = 8, q = 3, dim = 2, num_qpts = q * q; in main() local
15 CeedScalar q_ref[dim * num_qpts], q_weights[num_qpts]; in main()
16 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main()
21 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main()
28 for (CeedInt i = 0; i < dim * p * num_qpts; i++) { in main()
35 CeedVectorCreate(ceed, dim * num_qpts, &v); in main()
44 for (CeedInt i = 0; i < dim * num_qpts; i++) { in main()
H A Dt330-basis.c11 const CeedInt p = 8, q = 3, dim = 2, num_qpts = q * q; in main() local
13 CeedScalar q_ref[dim * num_qpts], q_weights[num_qpts]; in main()
14 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main()
22 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main()
H A Dt580-operator.c24 …CeedInt num_dofs_x = (n_x + 1) * (n_y + 1), num_dofs_u = num_faces * 2, num_qpts = q *… in main() local
28 CeedScalar q_ref[dim * num_qpts], q_weight[num_qpts]; in main()
29 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main()
78 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weight, &basis… in main()
/libCEED/examples/rust/ex1-volume/src/
H A Dmain.rs45 num_qpts, in example_1()
54 assert!(num_qpts >= 1); in example_1()
73 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_1()
89 num_qpts, in example_1()
96 num_qpts, in example_1()
116 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_1()
118 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_1()
207 let elem_qpts = num_qpts.pow(dim as u32); in example_1()
295 num_qpts: 6, in example_1_1d()
311 num_qpts: 6, in example_1_2d()
[all …]
/libCEED/include/ceed/jit-source/cuda/
H A Dcuda-ref-qfunction.h16 inline __device__ void readQuads(const CeedInt quad, const CeedInt num_qpts, const CeedScalar *d_u,… in readQuads() argument
18 r_u[comp] = d_u[quad + num_qpts * comp]; in readQuads()
26 inline __device__ void writeQuads(const CeedInt quad, const CeedInt num_qpts, const CeedScalar *r_v… in writeQuads() argument
28 d_v[quad + num_qpts * comp] = r_v[comp]; in writeQuads()
/libCEED/include/ceed/jit-source/hip/
H A Dhip-ref-qfunction.h16 inline __device__ void readQuads(const CeedInt quad, const CeedInt num_qpts, const CeedScalar *d_u,… in readQuads() argument
18 r_u[comp] = d_u[quad + num_qpts * comp]; in readQuads()
26 inline __device__ void writeQuads(const CeedInt quad, const CeedInt num_qpts, const CeedScalar *r_v… in writeQuads() argument
28 d_v[quad + num_qpts * comp] = r_v[comp]; in writeQuads()
/libCEED/examples/rust/ex1-volume-vector/src/
H A Dmain.rs46 num_qpts, in example_1_vector()
55 assert!(num_qpts >= 1); in example_1_vector()
75 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_1_vector()
91 num_qpts, in example_1_vector()
98 num_qpts, in example_1_vector()
118 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_1_vector()
120 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_1_vector()
122 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, ncomp_u, num_qpts)?; in example_1_vector()
211 let elem_qpts = num_qpts.pow(dim as u32); in example_1_vector()
313 num_qpts: 6, in example_1_vector_1d()
[all …]
/libCEED/examples/rust/ex2-surface/src/
H A Dmain.rs46 num_qpts, in example_2()
55 assert!(num_qpts >= 1); in example_2()
74 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_2()
90 num_qpts, in example_2()
97 num_qpts, in example_2()
117 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_2()
124 num_qpts, in example_2()
127 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_2()
244 let elem_qpts = num_qpts.pow(dim as u32); in example_2()
373 num_qpts: 6, in example_2_1d()
[all …]
/libCEED/examples/petsc/src/
H A Dlibceedsetup.c54 CeedInt num_qpts, q_data_size = bp_data.q_data_size; in SetupLibceedByDegree() local
71 CeedBasisGetNumQuadraturePoints(basis_u, &num_qpts); in SetupLibceedByDegree()
73 …CeedElemRestrictionCreateStrided(ceed, num_elem, num_qpts, num_comp_u, num_comp_u * num_elem * num… in SetupLibceedByDegree()
74 …eedElemRestrictionCreateStrided(ceed, num_elem, num_qpts, q_data_size, q_data_size * num_elem * nu… in SetupLibceedByDegree()
85 CeedVectorCreate(ceed, q_data_size * num_elem * num_qpts, &q_data); in SetupLibceedByDegree()
152 CeedVectorCreate(ceed, num_elem * num_qpts * num_comp_u, target); in SetupLibceedByDegree()
252 CeedInt num_qpts, q_data_size = bp_data.q_data_size; in SetupErrorOperator() local
269 CeedBasisGetNumQuadraturePoints(basis_u, &num_qpts); in SetupErrorOperator()
271 …CeedElemRestrictionCreateStrided(ceed, num_elem, num_qpts, num_comp_u, num_comp_u * num_elem * num… in SetupErrorOperator()
272 …eedElemRestrictionCreateStrided(ceed, num_elem, num_qpts, q_data_size, q_data_size * num_elem * nu… in SetupErrorOperator()
[all …]
/libCEED/examples/rust/ex2-surface-vector/src/
H A Dmain.rs47 num_qpts, in example_2_vector()
56 assert!(num_qpts >= 1); in example_2_vector()
76 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_2_vector()
92 num_qpts, in example_2_vector()
99 num_qpts, in example_2_vector()
119 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_2_vector()
126 num_qpts, in example_2_vector()
130 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, ncomp_u, num_qpts)?; in example_2_vector()
247 let elem_qpts = num_qpts.pow(dim as u32); in example_2_vector()
404 num_qpts: 6, in example_2_vector_1d()
[all …]
/libCEED/backends/sycl-ref/
H A Dceed-sycl-ref-basis.sycl.cpp47 const CeedInt work_group_size = CeedIntMin(impl->num_qpts, max_work_group_size); in CeedBasisApplyInterp_Sycl()
71 const CeedInt num_qpts = CeedIntPow(Q_1d, dim); in CeedBasisApplyInterp_Sycl() local
76 const CeedInt u_stride = is_transpose ? num_qpts : num_nodes; in CeedBasisApplyInterp_Sycl()
77 const CeedInt v_stride = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyInterp_Sycl()
174 const CeedInt num_qpts = CeedIntPow(Q_1d, dim); in CeedBasisApplyGrad_Sycl() local
179 const CeedInt u_stride = is_transpose ? num_qpts : num_nodes; in CeedBasisApplyGrad_Sycl()
180 const CeedInt v_stride = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyGrad_Sycl()
183 const CeedInt u_dim_stride = is_transpose ? num_elem * num_qpts * num_comp : 0; in CeedBasisApplyGrad_Sycl()
184 const CeedInt v_dim_stride = is_transpose ? 0 : num_elem * num_qpts * num_comp; in CeedBasisApplyGrad_Sycl()
203 CeedInt pre = is_transpose ? num_qpts : num_nodes; in CeedBasisApplyGrad_Sycl()
[all …]
H A Dceed-sycl-ref.hpp52 CeedInt num_qpts; member
65 CeedInt num_qpts; member
90 CeedInt num_qpts, num_comp; // Kernel parameters member
98 …CeedInt num_eval_mode_in, num_eval_mode_out, num_qpts, num_nodes, block_size, num_comp; // Ke… member
120 …ateH1_Sycl(CeedElemTopology topo, CeedInt dim, CeedInt num_dof, CeedInt num_qpts, const CeedScalar…
H A Dceed-sycl-ref-operator.sycl.cpp750 CeedInt num_nodes, num_qpts; in CeedOperatorAssembleDiagonalSetup_Sycl() local
752 CeedCallBackend(CeedBasisGetNumQuadraturePoints(basis_in, &num_qpts)); in CeedOperatorAssembleDiagonalSetup_Sycl()
755 diag->num_qpts = num_qpts; in CeedOperatorAssembleDiagonalSetup_Sycl()
759 const CeedInt i_len = num_qpts * num_nodes; in CeedOperatorAssembleDiagonalSetup_Sycl()
760 const CeedInt g_len = num_qpts * num_nodes * dim; in CeedOperatorAssembleDiagonalSetup_Sycl()
776 CeedCallBackend(CeedCalloc(num_qpts * num_nodes, &identity)); in CeedOperatorAssembleDiagonalSetup_Sycl()
777 …for (CeedSize i = 0; i < (num_nodes < num_qpts ? num_nodes : num_qpts); i++) identity[i * num_node… in CeedOperatorAssembleDiagonalSetup_Sycl()
835 const CeedSize num_qpts = diag->num_qpts; in CeedOperatorLinearDiagonal_Sycl() local
865 …ter_Sycl(&bt, eval_mode_out[e_out], identity, interp_out, &grad_out[d_out * num_qpts * num_nodes]); in CeedOperatorLinearDiagonal_Sycl()
872 …isPointer_Sycl(&b, eval_mode_in[e_in], identity, interp_in, &grad_in[d_in * num_qpts * num_nodes]); in CeedOperatorLinearDiagonal_Sycl()
[all …]
/libCEED/examples/rust/ex3-volume/src/
H A Dmain.rs45 num_qpts, in example_3()
53 assert!(num_qpts >= 1); in example_3()
72 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_3()
85 num_qpts, in example_3()
92 num_qpts, in example_3()
112 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_3()
119 num_qpts, in example_3()
122 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_3()
244 let elem_qpts = num_qpts.pow(dim as u32); in example_3()
376 num_qpts: 6, in example_3_1d()
[all …]
/libCEED/backends/magma/
H A Dceed-magma-basis.c281 CeedInt num_comp, num_nodes, num_qpts, P, Q, N; in CeedBasisApplyNonTensorCore_Magma() local
291 CeedCallBackend(CeedBasisGetNumQuadraturePoints(basis, &num_qpts)); in CeedBasisApplyNonTensorCore_Magma()
293 Q = num_qpts; in CeedBasisApplyNonTensorCore_Magma()
642 int CeedBasisCreateH1_Magma(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts… in CeedBasisCreateH1_Magma() argument
654 CeedCallBackend(magma_malloc((void **)&impl->d_q_weight, num_qpts * sizeof(q_weight[0]))); in CeedBasisCreateH1_Magma()
655 magma_setvector(num_qpts, sizeof(q_weight[0]), q_weight, 1, impl->d_q_weight, 1, data->queue); in CeedBasisCreateH1_Magma()
661 …CeedCallBackend(magma_malloc((void **)&impl->d_interp, num_qpts * num_nodes * q_comp_interp * size… in CeedBasisCreateH1_Magma()
662 …magma_setvector(num_qpts * num_nodes * q_comp_interp, sizeof(interp[0]), interp, 1, impl->d_interp… in CeedBasisCreateH1_Magma()
668 …CeedCallBackend(magma_malloc((void **)&impl->d_grad, num_qpts * num_nodes * q_comp_grad * sizeof(g… in CeedBasisCreateH1_Magma()
669 …magma_setvector(num_qpts * num_nodes * q_comp_grad, sizeof(grad[0]), grad, 1, impl->d_grad, 1, dat… in CeedBasisCreateH1_Magma()
[all …]
/libCEED/julia/LibCEED.jl/examples/
H A Dex2-surface.jl15 function run_ex2(; ceed_spec, dim, mesh_order, sol_order, num_qpts, prob_size, gallery) argument
24 create_tensor_h1_lagrange_basis(ceed, dim, ncompx, mesh_order + 1, num_qpts, GAUSS)
26 create_tensor_h1_lagrange_basis(ceed, dim, 1, sol_order + 1, num_qpts, GAUSS)
39 num_qpts,
48 num_qpts,
57 num_qpts,
101 elem_qpts = num_qpts^dim
168 num_qpts=6,
/libCEED/examples/rust/ex3-volume-vector/src/
H A Dmain.rs46 num_qpts, in example_3_vector()
54 assert!(num_qpts >= 1); in example_3_vector()
74 println!(" Num. 1D quadr. pts [-q] : {}", num_qpts); in example_3_vector()
87 num_qpts, in example_3_vector()
94 num_qpts, in example_3_vector()
114 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_3_vector()
121 num_qpts, in example_3_vector()
124 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, ncomp_u, num_qpts)?; in example_3_vector()
246 let elem_qpts = num_qpts.pow(dim as u32); in example_3_vector()
399 num_qpts: 6, in example_3_vector_1d()
[all …]
/libCEED/examples/ceed/
H A Dex1-volume-f.f9086 integer num_qpts local
106 elem_qpts = num_qpts**fe_dim
258 integer fe_dim, num_comp_x, mesh_degree, sol_degree, num_qpts variable
287 num_qpts = mesh_degree + 2
323 read(arg_value, '(I10)') num_qpts
358 write (*, *) ' Num. 1D quadrature pts [-q] : ', num_qpts
379 …call ceedbasiscreatetensorh1lagrange(ceed, fe_dim, num_comp_x, mesh_degree + 1, num_qpts, ceed_gau…
381 …call ceedbasiscreatetensorh1lagrange(ceed, fe_dim, 1, sol_degree + 1, num_qpts, ceed_gauss, sol_ba…
399 …all buildcartesianrestriction(ceed, fe_dim, num_xyz, mesh_degree, num_comp_x, mesh_size, num_qpts,&
401 …call buildcartesianrestriction(ceed, fe_dim, num_xyz, sol_degree, 1, sol_size, num_qpts, sol_restr…
[all …]
H A Dex2-surface.c49 …CeedInt dim, CeedInt num_xyz[3], CeedInt degree, CeedInt num_comp, CeedInt *size, CeedInt num_qpts,
61 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() local
81 parse_error = next_arg ? num_qpts = atoi(argv[++ia]), 0 : 1; in main()
111 printf(" Num. 1D quadrature pts [-q] : %" CeedInt_FMT "\n", num_qpts); in main()
130 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_x, mesh_degree + 1, num_qpts, CEED_GAUSS, &mes… in main()
131 CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, sol_degree + 1, num_qpts, CEED_GAUSS, &sol_basis); in main()
151 …BuildCartesianRestriction(ceed, dim, num_xyz, mesh_degree, num_comp_x, &mesh_size, num_qpts, &mesh… in main()
152 …riction(ceed, dim, num_xyz, sol_degree, dim * (dim + 1) / 2, &sol_size, num_qpts, NULL, &q_data_re… in main()
153 …BuildCartesianRestriction(ceed, dim, num_xyz, sol_degree, 1, &sol_size, num_qpts, &sol_restriction… in main()
205 CeedInt elem_qpts = CeedIntPow(num_qpts, dim); in main()
[all …]
H A Dex3-volume.c48 …edInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp, CeedInt *size, CeedInt num_qpts,
60 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() local
80 parse_error = next_arg ? num_qpts = atoi(argv[++ia]), 0 : 1; in main()
104 printf(" Num. 1D quadrature pts [-q] : %" CeedInt_FMT "\n", num_qpts); in main()
123 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_x, mesh_degree + 1, num_qpts, CEED_GAUSS, &mes… in main()
124 CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, sol_degree + 1, num_qpts, CEED_GAUSS, &sol_basis); in main()
143 …BuildCartesianRestriction(ceed, dim, num_xyz, mesh_degree, num_comp_x, &mesh_size, num_qpts, &mesh… in main()
144 …ion(ceed, dim, num_xyz, sol_degree, 1 + dim * (dim + 1) / 2, &sol_size, num_qpts, NULL, &q_data_re… in main()
145 …BuildCartesianRestriction(ceed, dim, num_xyz, sol_degree, 1, &sol_size, num_qpts, &sol_restriction… in main()
189 CeedInt elem_qpts = CeedIntPow(num_qpts, dim); in main()
[all …]
/libCEED/backends/ref/
H A Dceed-ref-basis.c22 CeedInt dim, num_comp, q_comp, num_nodes, num_qpts; in CeedBasisApplyCore_Ref() local
33 CeedCallBackend(CeedBasisGetNumQuadraturePoints(basis, &num_qpts)); in CeedBasisApplyCore_Ref()
105 … (t_mode == CEED_NOTRANSPOSE ? (d == 0 ? u : tmp[d % 2]) : &u[d * num_qpts * num_comp * num_elem]), in CeedBasisApplyCore_Ref()
122 … (t_mode == CEED_NOTRANSPOSE ? &v[d * num_qpts * num_comp * num_elem] in CeedBasisApplyCore_Ref()
137 … t_mode == CEED_NOTRANSPOSE ? u : &u[d * num_comp * num_qpts * num_elem], in CeedBasisApplyCore_Ref()
138 … t_mode == CEED_TRANSPOSE ? v : &v[d * num_comp * num_qpts * num_elem])); in CeedBasisApplyCore_Ref()
160 … (d == 0 ? (t_mode == CEED_NOTRANSPOSE ? u : &u[p * num_comp * num_qpts * num_elem]) : tmp[d % 2]), in CeedBasisApplyCore_Ref()
161 …(d == dim - 1 ? (t_mode == CEED_TRANSPOSE ? v : &v[p * num_comp * num_qpts * num_elem]) : tmp[(d +… in CeedBasisApplyCore_Ref()
199 CeedInt P = num_nodes, Q = num_qpts; in CeedBasisApplyCore_Ref()
236 for (CeedInt i = 0; i < num_qpts; i++) { in CeedBasisApplyCore_Ref()
[all …]
/libCEED/examples/rust-qfunctions/
H A Dex1-volume.c50 …edInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp, CeedInt *size, CeedInt num_qpts,
62 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() local
82 parse_error = next_arg ? num_qpts = atoi(argv[++ia]), 0 : 1; in main()
108 printf(" Num. 1D quadrature pts [-q] : %" CeedInt_FMT "\n", num_qpts); in main()
136 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_x, mesh_degree + 1, num_qpts, CEED_GAUSS, &mes… in main()
137 CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, sol_degree + 1, num_qpts, CEED_GAUSS, &sol_basis); in main()
156 …BuildCartesianRestriction(ceed, dim, num_xyz, mesh_degree, num_comp_x, &mesh_size, num_qpts, &mesh… in main()
157 …BuildCartesianRestriction(ceed, dim, num_xyz, sol_degree, 1, &sol_size, num_qpts, &sol_restriction… in main()
209 CeedInt elem_qpts = CeedIntPow(num_qpts, dim); in main()
328 …edInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp, CeedInt *size, CeedInt num_qpts, in BuildCartesianRestriction() argument
[all …]

1234