| /libCEED/tests/ |
| H A D | t570-operator.c | 17 CeedInt num_nodes = num_elem * (p - 1) + 1; in main() local
19 CeedScalar assembled_values[num_nodes * q * num_elem]; in main()
20 CeedScalar assembled_true[num_nodes * q * num_elem]; in main()
24 CeedVectorCreate(ceed, num_nodes, &u); in main()
33 …CeedElemRestrictionCreate(ceed, num_elem, p, 1, 1, num_nodes, CEED_MEM_HOST, CEED_USE_POINTER, ind… in main()
55 for (CeedInt k = 0; k < num_nodes * q * num_elem; ++k) { in main()
67 assembled_values[rows[k] * num_nodes + cols[k]] += assembled_array[k]; in main()
74 for (CeedInt j = 0; j < num_nodes; j++) { in main()
88 for (CeedInt i = 0; i < q * num_elem; i++) assembled_true[i * num_nodes + j] = v_array[i]; in main()
94 for (CeedInt j = 0; j < num_nodes; j++) { in main()
[all …]
|
| H A D | t509-operator.c | 16 CeedInt num_nodes = num_elem * p; in main() local
20 CeedVectorCreate(ceed, num_nodes, &u); in main()
21 CeedVectorCreate(ceed, num_nodes, &v); in main()
43 for (CeedInt i = 0; i < num_nodes; i++) { in main()
|
| H A D | t591-operator.c | 14 …CeedInt num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem… in main() local
54 CeedInt p = 2, num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1); in main() local
82 …CeedElemRestrictionCreate(ceed, num_elem, p * p, dim, 1, dim * num_nodes, CEED_MEM_HOST, CEED_COPY… in main()
83 CeedVectorCreate(ceed, dim * num_nodes, &x_elem); in main()
85 CeedScalar x_array[dim * num_nodes]; in main()
128 …CeedElemRestrictionCreate(ceed, num_elem, p * p, 1, 1, num_nodes, CEED_MEM_HOST, CEED_COPY_VALUES,… in main()
161 CeedVectorCreate(ceed, num_nodes, &u); in main()
163 CeedVectorCreate(ceed, num_nodes, &v); in main()
171 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
|
| H A D | t592-operator.c | 14 …CeedInt num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem… in main() local
54 CeedInt p = 2, num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1); in main() local
82 …CeedElemRestrictionCreate(ceed, num_elem, p * p, dim, 1, dim * num_nodes, CEED_MEM_HOST, CEED_COPY… in main()
83 CeedVectorCreate(ceed, dim * num_nodes, &x_elem); in main()
85 CeedScalar x_array[dim * num_nodes]; in main()
128 …CeedElemRestrictionCreate(ceed, num_elem, p * p, 1, 1, num_nodes, CEED_MEM_HOST, CEED_COPY_VALUES,… in main()
161 CeedVectorCreate(ceed, num_nodes, &u); in main()
163 CeedVectorCreate(ceed, num_nodes, &v); in main()
197 for (CeedInt i = 0; i < num_nodes; i++) area += v_array[i]; in main()
220 for (CeedInt i = 0; i < num_nodes; i++) area += v_array[i]; in main()
|
| H A D | t590-operator.c | 14 …CeedInt num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem… in main() local
75 …CeedElemRestrictionCreate(ceed, num_elem, p * p, 1, 1, num_nodes, CEED_MEM_HOST, CEED_COPY_VALUES,… in main()
91 CeedVectorCreate(ceed, num_nodes, &u); in main()
93 CeedVectorCreate(ceed, num_nodes, &v); in main()
101 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
|
| H A D | t599-operator.c | 13 …CeedInt num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem… in main() local
77 …CeedElemRestrictionCreate(ceed, num_elem, p * p, 1, 1, num_nodes, CEED_MEM_HOST, CEED_COPY_VALUES,… in main()
110 CeedVectorCreate(ceed, num_nodes, &u); in main()
112 CeedVectorCreate(ceed, num_nodes, &v); in main()
123 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
|
| /libCEED/examples/python/ |
| H A D | ex_common.py | 108 num_nodes = p ** dim
123 elem_nodes = np.zeros(num_elem * num_nodes, dtype=np.int32)
133 for n in range(num_nodes):
141 elem_nodes[e * num_nodes + n] = g_node
145 num_elem, num_nodes, num_comp, scalar_size, size, elem_nodes)
211 num_nodes = mesh_size // dim
214 for i in range(num_nodes):
219 for i in range(num_nodes):
221 v = np.pi / 2. * coords[i + num_nodes]
223 coords[i + num_nodes] = u * np.sin(v)
[all …]
|
| /libCEED/examples/rust/ex3-volume/src/ |
| H A D | transform.rs | 28 let num_nodes = mesh_size / dim; in transform_mesh_coordinates() localVariable
30 for i in 0..num_nodes { in transform_mesh_coordinates()
34 let v = coords[i + num_nodes] * std::f64::consts::PI / 2.; in transform_mesh_coordinates()
36 coords[i + num_nodes] = u * v.sin(); in transform_mesh_coordinates()
|
| /libCEED/examples/rust/ex3-volume-vector/src/ |
| H A D | transform.rs | 28 let num_nodes = mesh_size / dim; in transform_mesh_coordinates() localVariable
30 for i in 0..num_nodes { in transform_mesh_coordinates()
34 let v = coords[i + num_nodes] * std::f64::consts::PI / 2.; in transform_mesh_coordinates()
36 coords[i + num_nodes] = u * v.sin(); in transform_mesh_coordinates()
|
| /libCEED/examples/rust/ex1-volume/src/ |
| H A D | transform.rs | 26 let num_nodes = mesh_size / dim; in transform_mesh_coordinates() localVariable
27 for i in 0..num_nodes { in transform_mesh_coordinates()
31 let v = std::f64::consts::PI as libceed::Scalar / 2.0 * coords[i + num_nodes]; in transform_mesh_coordinates()
33 coords[i + num_nodes] = u * v.sin(); in transform_mesh_coordinates()
|
| /libCEED/examples/rust/ex1-volume-vector/src/ |
| H A D | transform.rs | 26 let num_nodes = mesh_size / dim; in transform_mesh_coordinates() localVariable
27 for i in 0..num_nodes { in transform_mesh_coordinates()
31 let v = std::f64::consts::PI as libceed::Scalar / 2.0 * coords[i + num_nodes]; in transform_mesh_coordinates()
33 coords[i + num_nodes] = u * v.sin(); in transform_mesh_coordinates()
|
| /libCEED/benchmarks/ |
| H A D | postprocess_plot.py | 81 num_nodes = num_procs / num_procs_node variable
89 backend, num_nodes, num_procs))
102 d = [[run['degree'], run['num_elem'], 1. * run['num_unknowns'] / num_nodes / vdim,
103 run['cg_iteration_dps'] / num_nodes]
126 d = [[run['degree'], run['num_elem'], 1. * run['num_unknowns'] / num_nodes / vdim,
127 run['cg_iteration_dps'] / num_nodes]
157 num_nodes, '' if num_nodes == 1 else 's',
181 code, test_short_save, short_backend, backend_memtype, num_nodes, num_procs_node)
|
| /libCEED/backends/sycl-ref/ |
| H A D | ceed-sycl-ref.hpp | 29 CeedInt num_nodes; member
51 CeedInt num_nodes; member
64 CeedInt num_nodes; member
89 CeedInt num_eval_mode_in, num_eval_mode_out, num_nodes; member
98 …CeedInt num_eval_mode_in, num_eval_mode_out, num_qpts, num_nodes, block_size, num_comp; // Ke… member
|
| H A D | ceed-sycl-ref-basis.sycl.cpp | 70 const CeedInt num_nodes = CeedIntPow(P_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()
173 const CeedInt num_nodes = CeedIntPow(P_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()
203 CeedInt pre = is_transpose ? num_qpts : num_nodes; in CeedBasisApplyGrad_Sycl()
341 const CeedInt P = is_transpose ? impl->num_qpts : impl->num_nodes; in CeedBasisApplyNonTensorInterp_Sycl()
342 const CeedInt Q = is_transpose ? impl->num_nodes : impl->num_qpts; in CeedBasisApplyNonTensorInterp_Sycl()
343 const CeedInt stride_0 = is_transpose ? 1 : impl->num_nodes; in CeedBasisApplyNonTensorInterp_Sycl()
[all …]
|
| H A D | ceed-sycl-restriction.sycl.cpp | 112 const CeedInt num_nodes = impl->num_nodes; in CeedElemRestrictionOffsetTranspose_Sycl() local
121 sycl::range<1> kernel_range(num_nodes * num_comp); in CeedElemRestrictionOffsetTranspose_Sycl()
127 const CeedInt node = id % num_nodes; in CeedElemRestrictionOffsetTranspose_Sycl()
128 const CeedInt comp = id / num_nodes; in CeedElemRestrictionOffsetTranspose_Sycl()
254 …CeedInt num_elem, elem_size, num_comp, num_nodes = 0, *ind_to_offset, *l_vec_ind… in CeedElemRestrictionOffset_Sycl() local
269 for (CeedInt i = 0; i < l_size; i++) num_nodes += is_node[i]; in CeedElemRestrictionOffset_Sycl()
270 impl->num_nodes = num_nodes; in CeedElemRestrictionOffset_Sycl()
274 CeedCallBackend(CeedCalloc(num_nodes, &l_vec_indices)); in CeedElemRestrictionOffset_Sycl()
284 const CeedInt size_offsets = num_nodes + 1; in CeedElemRestrictionOffset_Sycl()
314 …CeedCallSycl(ceed, impl->d_l_vec_indices = sycl::malloc_device<CeedInt>(num_nodes, data->sycl_devi… in CeedElemRestrictionOffset_Sycl()
[all …]
|
| H A D | ceed-sycl-ref-operator.sycl.cpp | 750 CeedInt num_nodes, num_qpts; in CeedOperatorAssembleDiagonalSetup_Sycl() local
751 CeedCallBackend(CeedBasisGetNumNodes(basis_in, &num_nodes)); in CeedOperatorAssembleDiagonalSetup_Sycl()
754 diag->num_nodes = num_nodes; 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()
834 const CeedSize num_nodes = diag->num_nodes; in CeedOperatorLinearDiagonal_Sycl() local
847 sycl::range<1> kernel_range(num_elem * num_nodes); in CeedOperatorLinearDiagonal_Sycl()
854 const CeedInt tid = idx % num_nodes; in CeedOperatorLinearDiagonal_Sycl()
[all …]
|
| /libCEED/examples/rust/mesh/src/ |
| H A D | lib.rs | 53 let num_nodes = p.pow(dim as u32); // number of nodes per element localVariable
69 let mut elem_nodes = vec![0; num_elem * num_nodes];
77 let loc_offset = e * num_nodes;
78 for loc_nodes in 0..num_nodes {
94 num_nodes,
|
| /libCEED/backends/magma/ |
| H A D | ceed-magma-basis.c | 33 CeedInt dim, num_comp, num_nodes, P_1d, Q_1d, P, Q; in CeedBasisApplyCore_Magma() local
43 CeedCallBackend(CeedBasisGetNumNodes(basis, &num_nodes)); in CeedBasisApplyCore_Magma()
281 CeedInt num_comp, num_nodes, num_qpts, P, Q, N; in CeedBasisApplyNonTensorCore_Magma() local
290 CeedCallBackend(CeedBasisGetNumNodes(basis, &num_nodes)); in CeedBasisApplyNonTensorCore_Magma()
292 P = num_nodes; in CeedBasisApplyNonTensorCore_Magma()
642 int CeedBasisCreateH1_Magma(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts… in CeedBasisCreateH1_Magma() argument
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 …]
|
| H A D | ceed-magma.h | 80 CEED_INTERN int CeedBasisCreateH1_Magma(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, Ceed…
82 CEED_INTERN int CeedBasisCreateHdiv_Magma(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, Ce…
84 CEED_INTERN int CeedBasisCreateHcurl_Magma(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, C…
|
| /libCEED/backends/hip-ref/ |
| H A D | ceed-hip-ref-basis.c | 257 CeedInt num_nodes, num_qpts; in CeedBasisApplyNonTensorCore_Hip() local
268 CeedCallBackend(CeedBasisGetNumNodes(basis, &num_nodes)); in CeedBasisApplyNonTensorCore_Hip()
285 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Hip()
295 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Hip()
305 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Hip()
315 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Hip()
445 int CeedBasisCreateH1_Hip(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts, … in CeedBasisCreateH1_Hip() argument
463 const CeedInt interp_bytes = q_bytes * num_nodes * q_comp_interp; in CeedBasisCreateH1_Hip()
469 const CeedInt grad_bytes = q_bytes * num_nodes * q_comp_grad; in CeedBasisCreateH1_Hip()
479 … basis_kernel_source, &data->module, 5, "BASIS_Q", num_qpts, "BASIS_P", num_nodes, "BASIS_Q_COMP_I… in CeedBasisCreateH1_Hip()
[all …]
|
| /libCEED/backends/cuda-ref/ |
| H A D | ceed-cuda-ref-basis.c | 258 CeedInt num_nodes, num_qpts; in CeedBasisApplyNonTensorCore_Cuda() local
269 CeedCallBackend(CeedBasisGetNumNodes(basis, &num_nodes)); in CeedBasisApplyNonTensorCore_Cuda()
286 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Cuda()
296 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Cuda()
306 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Cuda()
316 const int block_size_x = is_transpose ? num_nodes : num_qpts; in CeedBasisApplyNonTensorCore_Cuda()
446 int CeedBasisCreateH1_Cuda(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts,… in CeedBasisCreateH1_Cuda() argument
464 const CeedInt interp_bytes = q_bytes * num_nodes * q_comp_interp; in CeedBasisCreateH1_Cuda()
470 const CeedInt grad_bytes = q_bytes * num_nodes * q_comp_grad; in CeedBasisCreateH1_Cuda()
480 … basis_kernel_source, &data->module, 5, "BASIS_Q", num_qpts, "BASIS_P", num_nodes, "BASIS_Q_COMP_I… in CeedBasisCreateH1_Cuda()
[all …]
|
| /libCEED/backends/ref/ |
| H A D | ceed-ref.h | 72 CEED_INTERN int CeedBasisCreateH1_Ref(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedIn…
74 CEED_INTERN int CeedBasisCreateHdiv_Ref(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, Ceed…
76 CEED_INTERN int CeedBasisCreateHcurl_Ref(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, Cee…
|
| /libCEED/examples/rust-qfunctions/ |
| H A D | ex1-volume.c | 331 CeedInt num_nodes = CeedIntPow(p, dim); // number of scalar nodes per element in BuildCartesianRestriction() local
344 CeedInt *elem_nodes = malloc(sizeof(CeedInt) * num_elem * num_nodes); in BuildCartesianRestriction()
353 CeedInt *local_elem_nodes = elem_nodes + e * num_nodes; in BuildCartesianRestriction()
355 for (CeedInt l_nodes = 0; l_nodes < num_nodes; l_nodes++) { in BuildCartesianRestriction()
366 …CeedElemRestrictionCreate(ceed, num_elem, num_nodes, num_comp, scalar_size, num_comp * scalar_size… in BuildCartesianRestriction()
423 CeedInt num_nodes = mesh_size / dim; in TransformMeshCoords() local
425 for (CeedInt i = 0; i < num_nodes; i++) { in TransformMeshCoords()
428 CeedScalar u = coords[i], v = coords[i + num_nodes]; in TransformMeshCoords()
433 coords[i + num_nodes] = u * sin(v); in TransformMeshCoords()
|
| /libCEED/examples/ceed/ |
| H A D | ex3-volume.c | 310 CeedInt num_nodes = CeedIntPow(p, dim); // number of scalar nodes per element in BuildCartesianRestriction() local
323 CeedInt *elem_nodes = malloc(sizeof(CeedInt) * num_elem * num_nodes); in BuildCartesianRestriction()
332 CeedInt *local_elem_nodes = elem_nodes + e * num_nodes; in BuildCartesianRestriction()
334 for (CeedInt l_nodes = 0; l_nodes < num_nodes; l_nodes++) { in BuildCartesianRestriction()
346 …CeedElemRestrictionCreate(ceed, num_elem, num_nodes, num_comp, scalar_size, num_comp * scalar_size… in BuildCartesianRestriction()
403 CeedInt num_nodes = mesh_size / dim; in TransformMeshCoords() local
404 for (CeedInt i = 0; i < num_nodes; i++) { in TransformMeshCoords()
407 CeedScalar u = coords[i], v = coords[i + num_nodes]; in TransformMeshCoords()
412 coords[i + num_nodes] = u * sin(v); in TransformMeshCoords()
|
| H A D | ex1-volume.c | 321 CeedInt num_nodes = CeedIntPow(p, dim); // number of scalar nodes per element in BuildCartesianRestriction() local
334 CeedInt *elem_nodes = malloc(sizeof(CeedInt) * num_elem * num_nodes); in BuildCartesianRestriction()
343 CeedInt *local_elem_nodes = elem_nodes + e * num_nodes; in BuildCartesianRestriction()
345 for (CeedInt l_nodes = 0; l_nodes < num_nodes; l_nodes++) { in BuildCartesianRestriction()
356 …CeedElemRestrictionCreate(ceed, num_elem, num_nodes, num_comp, scalar_size, num_comp * scalar_size… in BuildCartesianRestriction()
413 CeedInt num_nodes = mesh_size / dim; in TransformMeshCoords() local
415 for (CeedInt i = 0; i < num_nodes; i++) { in TransformMeshCoords()
418 CeedScalar u = coords[i], v = coords[i + num_nodes]; in TransformMeshCoords()
423 coords[i + num_nodes] = u * sin(v); in TransformMeshCoords()
|