Home
last modified time | relevance | path

Searched refs:b (Results 1 – 25 of 76) sorted by relevance

1234

/libCEED/backends/avx/
H A Dceed-avx-tensor.c21 #define fmadd(c, a, b) (c) = _mm256_fmadd_pd((a), (b), (c)) argument
23 #define fmadd(c, a, b) (c) += _mm256_mul_pd((a), (b)) argument
33 #define fmadd(c, a, b) (c) = _mm_fmadd_ps((a), (b), (c)) argument
35 #define fmadd(c, a, b) (c) += _mm_mul_ps((a), (b)) argument
60 for (CeedInt b = 0; b < B; b++) { in CeedTensorContract_Avx_Blocked() local
62 rtype tqv = set1(t[(j + jj) * t_stride_0 + b * t_stride_1]); in CeedTensorContract_Avx_Blocked()
64 fmadd(vv[jj][cc], tqv, loadu(&u[(a * B + b) * C + c + cc * 4])); in CeedTensorContract_Avx_Blocked()
83 for (CeedInt b = 0; b < B; b++) { in CeedTensorContract_Avx_Blocked() local
85 rtype tqv = set1(t[(j + jj) * t_stride_0 + b * t_stride_1]); in CeedTensorContract_Avx_Blocked()
88 fmadd(vv[jj][cc], tqv, loadu(&u[(a * B + b) * C + c + cc * 4])); in CeedTensorContract_Avx_Blocked()
[all …]
/libCEED/include/ceed/
H A Ddeprecated.h22 …int CeedCompositeOperatorCreate(Ceed a, CeedOperator *b) { return CeedOperatorCreateComposite(a, b… in CeedCompositeOperatorCreate() argument
24 …dCompositeOperatorAddSub(CeedOperator a, CeedOperator b) { return CeedOperatorCompositeAddSub(a, b… in CeedCompositeOperatorAddSub() argument
26 …dCompositeOperatorGetNumSub(CeedOperator a, CeedInt *b) { return CeedOperatorCompositeGetNumSub(a,… in CeedCompositeOperatorGetNumSub() argument
28 …teOperatorGetSubList(CeedOperator a, CeedOperator **b) { return CeedOperatorCompositeGetSubList(a,… in CeedCompositeOperatorGetSubList() argument
30 static inline int CeedCompositeOperatorGetSubByName(CeedOperator a, const char *b, CeedOperator *c)… in CeedCompositeOperatorGetSubByName() argument
31 return CeedOperatorCompositeGetSubByName(a, b, c); in CeedCompositeOperatorGetSubByName()
34 static inline int CeedCompositeOperatorGetMultiplicity(CeedOperator a, CeedInt b, CeedInt *c, CeedV… in CeedCompositeOperatorGetMultiplicity() argument
35 return CeedOperatorCompositeGetMultiplicity(a, b, c, d); in CeedCompositeOperatorGetMultiplicity()
/libCEED/julia/LibCEED.jl/src/
H A DBasis.jl34 destroy(b::Basis) = C.CeedBasisDestroy(b.ref) # COV_EXCL_LINE
35 Base.getindex(b::Basis) = b.ref[]
36 Base.show(io::IO, ::MIME"text/plain", b::Basis) = ceed_show(io, b, C.CeedBasisView)
290 …apply!(b::Basis, nelem, tmode::TransposeMode, emode::EvalMode, u::AbstractCeedVector, v::AbstractC…
308 b::Basis,
315 C.CeedBasisApply(b[], nelem, tmode, emode, u[], v[])
319 apply(b::Basis, u::AbstractVector; nelem=1, tmode=NOTRANSPOSE, emode=EVAL_INTERP)
321 Performs the same function as the above-defined [`apply!`](@ref apply!(b::Basis, nelem,
327 function apply(b::Basis, u::AbstractVector; nelem=1, tmode=NOTRANSPOSE, emode=EVAL_INTERP)
329 ccall((:CeedBasisGetCeed, C.libceed), Cint, (C.CeedBasis, Ptr{C.Ceed}), b[], ceed_ref)
[all …]
/libCEED/tests/
H A Dt300-basis-f.f907 integer b variable
14 call ceedbasiscreatetensorh1lagrange(ceed,1,1,4,4,ceed_gauss_lobatto,b,&
16 call ceedbasisview(b,err)
17 call ceedbasisdestroy(b,err)
19 call ceedbasiscreatetensorh1lagrange(ceed,1,1,4,4,ceed_gauss,b,err)
20 call ceedbasisview(b,err)
21 call ceedbasissetnumviewtabs(b,1,err)
22 call ceedbasisview(b,err)
23 call ceedbasisdestroy(b,err)
H A Dt302-basis-f.f907 integer b,p variable
21 call ceedbasiscreatetensorh1lagrange(ceed,1,1,p,p,ceed_gauss_lobatto,b,&
23 call ceedbasisgetcollocatedgrad(b,collograd1d,err)
24 call ceedbasisgetgrad1d(b,grad1d,gradoffset,err)
35 call ceedbasisdestroy(b,err)
38 call ceedbasiscreatetensorh1lagrange(ceed,1,1,p,p,ceed_gauss,b,err)
39 call ceedbasisgetcollocatedgrad(b,collograd1d,err)
41 call ceedbasisgetqref(b,qref,qoffset,err)
58 call ceedbasisdestroy(b,err)
61 call ceedbasiscreatetensorh1lagrange(ceed,1,1,p,p+2,ceed_gauss,b,err)
[all …]
H A Dt101-vector-f.f909 real*8 b(10) variable
30 call ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
33 diff=b(boffset+i)-10-i
36 write(*,*) 'Error reading array b(',i,')=',b(boffset+i)
41 call ceedvectorrestorearrayread(x,b,boffset,err)
44 call ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
47 diff=b(boffset+i)-val
50 write(*,*) 'Error reading array b(',i,')=',b(boffset+i)
55 call ceedvectorrestorearrayread(x,b,boffset,err)
H A Dt320-basis-f.f9022 integer b
33 & b,err)
34 call ceedbasisview(b,err)
35 call ceedbasissetnumviewtabs(b,1,err)
36 call ceedbasisview(b,err)
38 call ceedbasisdestroy(b,err)
H A Dt306-basis-f.f907 integer b variable
15 call ceedbasiscreatetensorh1lagrange(ceed,3,1,4,5,ceed_gauss_lobatto,b,&
18 call ceedbasisgetnumnodes(b,p,err)
19 call ceedbasisgetnumquadraturepoints(b,q,err)
32 call ceedbasisdestroy(b,err)
H A Dt100-vector-f.f9010 real*8 b(10) variable
28 call ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
31 diff=b(i+boffset)-10-i
34 write(*,*) 'Error reading array b(',i,')=',b(i+boffset)
39 call ceedvectorrestorearrayread(x,b,boffset,err)
H A Dt119-vector-f.f9010 real*8 b(10) variable
30 call ceedvectorgetarrayread(x,ceed_mem_host,b,boffset,err)
32 diff=1./(real(10+i,8))-b(i+boffset)
35 write(*,*) 'Error reading array b(',i,')=',b(i+boffset),diff
40 call ceedvectorrestorearrayread(x,b,boffset,err)
H A Dt106-vector-f.f909 real*8 b(10) variable
31 b(i)=0
34 call ceedvectorsetarray(y,ceed_mem_host,ceed_use_pointer,b,boffset,err)
42 diff = b(i+boffset)-(9+i)
45 write(*,*) 'Error reading array b(',i,') = ',b(i+boffset),' != ',9+i
H A Dt104-vector-f.f909 real*8 b(10) variable
28 call ceedvectorgetarray(x,ceed_mem_host,b,boffset,err)
29 b(boffset+3) = -3.14
30 call ceedvectorrestorearray(x,b,boffset,err)
H A Dt109-vector-f.f909 real*8 b(10) variable
41 call ceedvectorgetarraywrite(x,ceed_mem_host,b,boffset,err)
42 b(boffset+5) = -3.14
43 call ceedvectorrestorearray(x,b,boffset,err)
H A Dt105-vector-f.f909 real*8 b(10) variable
30 call ceedvectorgetarrayread(x,ceed_mem_device,b,boffset,err)
31 call ceedvectorsetarray(y,ceed_mem_device,ceed_copy_values,b,boffset,err)
32 call ceedvectorrestorearrayread(x,b,boffset,err)
H A Dt110-vector.c12 CeedScalar *a, *b; in main() local
21 CeedVectorGetArray(x, CEED_MEM_HOST, &b); in main()
25 CeedVectorRestoreArray(x, &b); in main()
/libCEED/examples/nek/
H A Dnek-examples.sh272 if [ ! -f b$i/b$i.map ]; then
279 mkdir -p b$i
280 sed "5s/.*/-$nex -$ney -$nez/" ${CEED_DIR}/examples/nek/boxes/b.box > b$i/b$i.box
281 cp ${CEED_DIR}/examples/nek/boxes/b1e.rea b$i/
283 cd b$i
284 genbb b$i
285 genbb b$i &> log || {
308 if [ ! -f ${nek_box_dir}/b${nek_box}/b${nek_box}.rea ] || \
309 [ ! -f ${nek_box_dir}/b${nek_box}/b${nek_box}.map ]; then
320 echo b${nek_box} > SESSION.NAME
[all …]
/libCEED/python/tests/
H A Dtest-1-vector.py24 with x.array_read() as b:
25 for i in range(len(b)):
26 assert b[i] == value
42 with x.array_read() as b:
44 assert b[i] == 10 + i
59 with x.array() as b:
60 for i in range(len(b)):
61 assert b[i] == 10 + i
86 b = x.get_array_read()
128 with x.array() as b:
[all …]
H A Dtest-3-basis.py48 b = ceed.BasisTensorH1Lagrange(1, 1, 4, 4, libceed.GAUSS_LOBATTO)
49 print(b)
50 del b
52 b = ceed.BasisTensorH1Lagrange(1, 1, 4, 4, libceed.GAUSS)
53 print(b)
54 del b
70 b = ceed.BasisTensorH1Lagrange(3, 1, 4, 5, libceed.GAUSS_LOBATTO)
72 p = b.get_num_nodes()
73 q = b.get_num_quadrature_points()
215 b = ceed.BasisH1(libceed.TRIANGLE, 1, P, Q, interp, grad, qref, qweight)
[all …]
/libCEED/include/ceed/jit-source/cuda/
H A Dcuda-ref-basis-tensor-at-points.h90 …for (CeedInt b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) … in InterpAtPoints() local
115 for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c]; in InterpAtPoints() local
209 …for (CeedInt b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) … in InterpTransposeAtPoints() local
270 …for (CeedInt b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) … in GradAtPoints() local
299 for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c]; in GradAtPoints() local
400 …for (CeedInt b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) … in GradTransposeAtPoints() local
/libCEED/include/ceed/jit-source/hip/
H A Dhip-ref-basis-tensor-at-points.h90 …for (CeedInt b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) … in InterpAtPoints() local
115 for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c]; in InterpAtPoints() local
209 …for (CeedInt b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) … in InterpTransposeAtPoints() local
270 …for (CeedInt b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) … in GradAtPoints() local
299 for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c]; in GradAtPoints() local
400 …for (CeedInt b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) … in GradTransposeAtPoints() local
/libCEED/backends/ref/
H A Dceed-ref-tensor.c30 for (CeedInt b = 0; b < B; b++) { in CeedTensorContractApply_Ref() local
32 CeedScalar tq = t[j * t_stride_0 + b * t_stride_1]; in CeedTensorContractApply_Ref()
33 for (CeedInt c = 0; c < C; c++) v[(a * J + j) * C + c] += tq * u[(a * B + b) * C + c]; in CeedTensorContractApply_Ref()
/libCEED/rust/libceed/src/
H A Doperator.rs542 pub fn create<'b>( in create()
544 qf: impl Into<QFunctionOpt<'b>>, in create() argument
545 dqf: impl Into<QFunctionOpt<'b>>, in create() argument
546 dqfT: impl Into<QFunctionOpt<'b>>, in create() argument
803 pub fn field<'b>( in field()
806 r: impl Into<ElemRestrictionOpt<'b>>, in field() argument
807 b: impl Into<BasisOpt<'b>>, in field() argument
808 v: impl Into<VectorOpt<'b>>, in field() argument
817 b.into().to_raw(), in field()
1708 pub fn create_multigrid_level<'b>( in create_multigrid_level()
[all …]
/libCEED/include/
H A Dceed-fortran-name.h21 #define TOKEN_PASTE_(a, b) a##b argument
22 #define TOKEN_PASTE(a, b) TOKEN_PASTE_(a, b) argument
/libCEED/backends/opt/
H A Dceed-opt-tensor.c27 for (CeedInt b = 0; b < B; b++) { in CeedTensorContractApply_Core_Opt() local
29 CeedScalar tq = t[j * t_stride_0 + b * t_stride_1]; in CeedTensorContractApply_Core_Opt()
30 for (CeedInt c = 0; c < C; c++) v[(a * J + j) * C + c] += tq * u[(a * B + b) * C + c]; in CeedTensorContractApply_Core_Opt()
/libCEED/examples/fluids/qfunctions/
H A Dutils.h18 CEED_QFUNCTION_HELPER CeedScalar Max(CeedScalar a, CeedScalar b) { return a < b ? b : a; } in Max() argument
19 CEED_QFUNCTION_HELPER CeedScalar Min(CeedScalar a, CeedScalar b) { return a < b ? a : b; } in Min() argument
21 CEED_QFUNCTION_HELPER void SwapScalar(CeedScalar *a, CeedScalar *b) { in SwapScalar() argument
23 *a = *b; in SwapScalar()
24 *b = temp; in SwapScalar()
72 CeedScalar *b) { in MatVecNM() argument
75 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M); in MatVecNM()
78 …CeedInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[i] += A[j * M + i] *… in MatVecNM()
84 … CeedScalar A[3][3], const CeedScalar x[3], const CeedTransposeMode transpose_A, CeedScalar b[3]) { in MatVec3()
85 MatVecNM((const CeedScalar *)A, (const CeedScalar *)x, 3, 3, transpose_A, (CeedScalar *)b); in MatVec3()

1234