| /libCEED/tests/ |
| H A D | t340-basis.h | 12 …uildHcurl2DSimplex(CeedScalar *q_ref, CeedScalar *q_weight, CeedScalar *interp, CeedScalar *curl) { in BuildHcurl2DSimplex() argument
49 curl[i * P + 0] = 24. * x1 - 6.; in BuildHcurl2DSimplex()
50 curl[i * P + 1] = 24. * x2 - 6.; in BuildHcurl2DSimplex()
51 curl[i * P + 2] = 24. * x1 + 24. * x2 - 18.; in BuildHcurl2DSimplex()
52 curl[i * P + 3] = -24. * x2 + 6.; in BuildHcurl2DSimplex()
53 curl[i * P + 4] = -24. * x1 - 24. * x2 + 18.; in BuildHcurl2DSimplex()
54 curl[i * P + 5] = 24. * x1 - 6.; in BuildHcurl2DSimplex()
55 curl[i * P + 6] = 24. * x1 + 48. * x2 - 24.; in BuildHcurl2DSimplex()
56 curl[i * P + 7] = -48. * x1 - 24. * x1 + 24.; in BuildHcurl2DSimplex()
|
| H A D | t342-basis.c | 15 CeedScalar interp[dim * p * q], curl[p * q]; in main() local
20 BuildHcurl2DSimplex(q_ref, q_weight, interp, curl); in main()
21 …CeedBasisCreateHcurl(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, curl, q_ref, q_weight, &basis); in main()
29 …if (fabs(curl[i] - curl_in_basis[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", curl[i], curl_in_… in main()
36 row_sum[i] += curl[j + i * p]; in main()
42 column_sum[i] += curl[i + j * p]; in main()
|
| H A D | t340-basis.c | 13 CeedScalar interp[dim * p * q], curl[p * q]; in main() local
20 BuildHcurl2DSimplex(q_ref, q_weight, interp, curl); in main()
21 …CeedBasisCreateHcurl(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, curl, q_ref, q_weight, &basis); in main()
|
| H A D | t341-basis.c | 15 CeedScalar interp[dim * p * q], curl[p * q]; in main() local
20 BuildHcurl2DSimplex(q_ref, q_weight, interp, curl); in main()
21 …CeedBasisCreateHcurl(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, curl, q_ref, q_weight, &basis); in main()
|
| H A D | README.md | 17 4. CeedBasis non-tensor H(curl) basis tests\
32 8. CeedOperator H(div) and H(curl) tests\
|
| /libCEED/tests/output/ |
| H A D | t340-basis.out | 1 CeedBasis in a H(curl) space on a triangle element
17 curl:
22 CeedBasis in a H(curl) space on a triangle element
38 curl:
|
| /libCEED/julia/LibCEED.jl/test/ |
| H A D | buildmats.jl | 45 curl = Array{T}(undef, 1, Q, P)
72 curl[1, i, 1] = 2.0
73 curl[1, i, 2] = -2.0
74 curl[1, i, 3] = -2.0
77 return interp, curl
|
| /libCEED/python/tests/ |
| H A D | buildmats.py | 95 curl = np.empty(P * Q, dtype=mat_dtype)
122 curl[i * P + 0] = 2.
123 curl[i * P + 1] = -2.
124 curl[i * P + 2] = -2.
126 return interp, curl
|
| H A D | test-3-basis.py | 352 interp, curl = bm.buildmatshcurl(qref, qweight, libceed.scalar_types[
355 b = ceed.BasisHcurl(libceed.TRIANGLE, 1, P, Q, interp, curl, qref, qweight)
|
| /libCEED/include/ceed/jit-source/cuda/ |
| H A D | cuda-ref-operator-assemble-diagonal.h | 22 … const CeedScalar *interp, const CeedScalar *grad, const CeedScalar *div, const CeedScalar *curl) { in GetBasisPointer() argument
37 *basis_ptr = curl; in GetBasisPointer()
|
| /libCEED/include/ceed/jit-source/hip/ |
| H A D | hip-ref-operator-assemble-diagonal.h | 22 … const CeedScalar *interp, const CeedScalar *grad, const CeedScalar *div, const CeedScalar *curl) { in GetBasisPointer() argument
37 *basis_ptr = curl; in GetBasisPointer()
|
| /libCEED/julia/LibCEED.jl/src/ |
| H A D | Basis.jl | 231 create_hcurl_basis(c::Ceed, topo::Topology, ncomp, nnodes, nqpts, interp, curl, qref, qweight)
233 Create a non tensor-product basis for H(curl) discretizations
243 - `curl`: Matrix of size `(curlcomp, nqpts, nnodes)`, `curlcomp = 1 if dim < 3 else dim`)
244 matrix expressing curl of basis functions at quadrature points.
257 curl::AbstractArray{CeedScalar},
264 @assert size(curl) == (curlcomp, nqpts, nnodes)
270 curl_rowmajor = permutedims(curl, [3, 2, 1])
466 getnumqpts(b), getnumnodes(b))` for a given vector $H(div)$ or $H(curl)$ basis.
543 Get the curl matrix of the given [`Basis`](@ref). Returns a tensor of size
|
| /libCEED/backends/ref/ |
| H A D | ceed-ref-basis.c | 225 const CeedScalar *curl; in CeedBasisApplyCore_Ref() local
227 CeedCallBackend(CeedBasisGetCurl(basis, &curl)); in CeedBasisApplyCore_Ref()
228 …eedTensorContractStridedApply(contract, num_comp, P, num_elem, q_comp, Q, curl, t_mode, add, u, v)… in CeedBasisApplyCore_Ref()
356 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis) { in CeedBasisCreateHcurl_Ref() argument
|
| H A D | ceed-ref.h | 77 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis);
|
| /libCEED/interface/ |
| H A D | ceed-basis.c | 1792 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis *basis) { in CeedBasisCreateHcurl() argument
1800 …CeedCall(CeedBasisCreateHcurl(delegate, topo, num_comp, num_nodes, num_qpts, interp, curl, q_ref, … in CeedBasisCreateHcurl()
1826 CeedCall(CeedMalloc(curl_comp * Q * P, &(*basis)->curl)); in CeedBasisCreateHcurl()
1828 if (curl) memcpy((*basis)->curl, curl, curl_comp * Q * P * sizeof(curl[0])); in CeedBasisCreateHcurl()
1829 CeedCall(ceed->BasisCreateHcurl(topo, dim, P, Q, interp, curl, q_ref, q_weight, *basis)); in CeedBasisCreateHcurl()
2006 const CeedScalar *q_ref, *q_weight, *interp, *grad, *div, *curl; in CeedBasisView() local
2016 CeedCall(CeedBasisGetCurl(basis, &curl)); in CeedBasisView()
2030 if (curl) { in CeedBasisView()
2032 CeedCall(CeedScalarView("curl", "\t% 12.8f", q_comp * Q, P, curl, tabs, stream)); in CeedBasisView()
2456 int CeedBasisGetCurl(CeedBasis basis, const CeedScalar **curl) { in CeedBasisGetCurl() argument
[all …]
|
| /libCEED/backends/magma/ |
| H A D | ceed-magma.h | 85 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis);
|
| H A D | ceed-magma-basis.c | 771 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis) { in CeedBasisCreateHcurl_Magma() argument
792 if (curl) { in CeedBasisCreateHcurl_Magma()
796 …ackend(magma_malloc((void **)&impl->d_curl, num_qpts * num_nodes * q_comp_curl * sizeof(curl[0]))); in CeedBasisCreateHcurl_Magma()
797 …magma_setvector(num_qpts * num_nodes * q_comp_curl, sizeof(curl[0]), curl, 1, impl->d_curl, 1, dat… in CeedBasisCreateHcurl_Magma()
|
| /libCEED/rust/libceed/src/ |
| H A D | basis.rs | 294 curl: &[crate::Scalar], in create_Hcurl()
313 curl.as_ptr(), in create_Hcurl()
|
| H A D | lib.rs | 1062 curl: &[crate::Scalar], in basis_Hcurl()
1067 self, topo, ncomp, nnodes, nqpts, interp, curl, qref, qweight, in basis_Hcurl()
|
| /libCEED/python/ |
| H A D | ceed.py | 541 interp, curl, qref, qweight):
565 interp, curl, qref, qweight)
|
| /libCEED/backends/cuda-ref/ |
| H A D | ceed-cuda-ref.h | 164 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis);
|
| H A D | ceed-cuda-ref-basis.c | 557 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis) { in CeedBasisCreateHcurl_Cuda() argument
579 if (curl) { in CeedBasisCreateHcurl_Cuda()
583 CeedCallCuda(ceed, cudaMemcpy(data->d_curl, curl, curl_bytes, cudaMemcpyHostToDevice)); in CeedBasisCreateHcurl_Cuda()
|
| /libCEED/backends/hip-ref/ |
| H A D | ceed-hip-ref.h | 169 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis);
|
| H A D | ceed-hip-ref-basis.c | 556 … const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis) { in CeedBasisCreateHcurl_Hip() argument
578 if (curl) { in CeedBasisCreateHcurl_Hip()
582 CeedCallHip(ceed, hipMemcpy(data->d_curl, curl, curl_bytes, hipMemcpyHostToDevice)); in CeedBasisCreateHcurl_Hip()
|
| /libCEED/include/ |
| H A D | ceed-impl.h | 229 …CeedScalar *curl; /* row-major matrix of shape [curl_dim * Q, P], curl_dim = 1 if dim < 3 else dim… member
|