Home
last modified time | relevance | path

Searched full:1 (Results 1 – 25 of 859) sorted by relevance

12345678910>>...35

/libCEED/backends/magma/tuning/
H A Dv100_rtc.h6 {3, 1, 1024, 1, 1 },
7 {3, 1, 1024, 2, 1 },
8 {3, 1, 5120, 1, 1 },
9 {3, 1, 5120, 2, 1 },
10 {3, 1, 10240, 1, 1 },
11 {3, 1, 10240, 2, 1 },
12 {3, 1, 51200, 1, 1 },
13 {3, 1, 51200, 2, 1 },
14 {3, 1, 102400, 1, 1 },
15 {3, 1, 102400, 2, 1 },
[all …]
H A Dmi100_rtc.h6 {3, 1, 1024, 1, 2 },
7 {3, 1, 1024, 2, 1 },
8 {3, 1, 5120, 1, 1 },
9 {3, 1, 5120, 2, 1 },
10 {3, 1, 10240, 1, 1 },
11 {3, 1, 10240, 2, 1 },
12 {3, 1, 51200, 1, 1 },
13 {3, 1, 51200, 2, 1 },
14 {3, 1, 102400, 1, 1 },
15 {3, 1, 102400, 2, 1 },
[all …]
H A Dh100_rtc.h6 {3, 1, 1024, 1, 1 },
7 {3, 1, 1024, 2, 1 },
8 {3, 1, 5120, 1, 1 },
9 {3, 1, 5120, 2, 1 },
10 {3, 1, 10240, 1, 1 },
11 {3, 1, 10240, 2, 1 },
12 {3, 1, 51200, 1, 1 },
13 {3, 1, 51200, 2, 1 },
14 {3, 1, 102400, 1, 1 },
15 {3, 1, 102400, 2, 1 },
[all …]
H A Dmi250x_rtc.h6 {3, 1, 1024, 1, 1},
7 {3, 1, 1024, 2, 1},
8 {3, 1, 5120, 1, 2},
9 {3, 1, 5120, 2, 2},
10 {3, 1, 10240, 1, 4},
11 {3, 1, 10240, 2, 4},
12 {3, 1, 51200, 1, 1},
13 {3, 1, 51200, 2, 1},
14 {3, 1, 102400, 1, 1},
15 {3, 1, 102400, 2, 1},
[all …]
H A Da100_rtc.h6 {3, 1, 1024, 1, 1 },
7 {3, 1, 1024, 2, 1 },
8 {3, 1, 5120, 1, 1 },
9 {3, 1, 5120, 2, 1 },
10 {3, 1, 10240, 1, 1 },
11 {3, 1, 10240, 2, 1 },
12 {3, 1, 51200, 1, 1 },
13 {3, 1, 51200, 2, 2 },
14 {3, 1, 102400, 1, 1 },
15 {3, 1, 102400, 2, 1 },
[all …]
H A Dv100.h6 {3, 512, 1, 64, 0},
7 {3, 1024, 1, 32, 0},
8 {3, 1536, 1, 32, 0},
9 {3, 2048, 1, 64, 0},
10 {3, 2560, 1, 128, 0},
11 {3, 3072, 1, 64, 0},
12 {3, 3584, 1, 64, 0},
13 {3, 4096, 1, 64, 0},
14 {3, 4608, 1, 32, 0},
15 {3, 5120, 1, 32, 0},
[all …]
H A Da100.h6 {3, 512, 1, 256, 0},
7 {3, 1024, 1, 32, 0},
8 {3, 1536, 1, 32, 0},
9 {3, 2048, 1, 2048, 0},
10 {3, 2560, 1, 256, 0},
11 {3, 3072, 1, 32, 0},
12 {3, 3584, 1, 32, 0},
13 {3, 4096, 1, 512, 0},
14 {3, 4608, 1, 512, 0},
15 {3, 5120, 1, 64, 0},
[all …]
H A Dmi100.h6 {3, 512, 1, 32, 0},
7 {3, 1024, 1, 512, 0},
8 {3, 1536, 1, 256, 0},
9 {3, 2048, 1, 512, 0},
10 {3, 2560, 1, 512, 0},
11 {3, 3072, 1, 512, 0},
12 {3, 3584, 1, 32, 1},
13 {3, 4096, 1, 512, 0},
14 {3, 4608, 1, 512, 0},
15 {3, 5120, 1, 32, 1},
[all …]
H A Dmi250x.h6 {3, 512, 1, 128, 0},
7 {3, 1024, 1, 64, 1},
8 {3, 1536, 1, 64, 1},
9 {3, 2048, 1, 1024, 1},
10 {3, 2560, 1, 256, 1},
11 {3, 3072, 1, 1024, 1},
12 {3, 3584, 1, 512, 1},
13 {3, 4096, 1, 2048, 0},
14 {3, 4608, 1, 512, 1},
15 {3, 5120, 1, 256, 1},
[all …]
/libCEED/julia/LibCEED.jl/test/
H A Dbuildmats.jl6 qref[1, 1] = -1.0/sqrt(3.0)
7 qref[1, 2] = qref[1, 1]
8 qref[1, 3] = qref[1, 1]
9 qref[1, 4] = -qref[1, 1]
10 qref[2, 1] = -qref[1, 1]
11 qref[2, 2] = -qref[1, 1]
12 qref[2, 3] = qref[1, 1]
13 qref[2, 4] = qref[1, 1]
14 qweight[1] = 1.0
20 for i = 1:Q
[all …]
H A Druntests.jl69 @test axes(v) == (1:n,)
70 @test ndims(v) == 1
71 @test ndims(CeedVector) == 1
81 for p ∈ [1, 2, Inf]
140 ncomp = 1
155 @test w1d ≈ CeedScalar[1/3, 4/3, 1/3]
164 w1d = CeedScalar[1/3, 4/3, 1/3]
167 d2d[1, :, :] = kron(b1d, d1d)
171 b2 = create_tensor_h1_basis(c, dim2, 1, p, q, b1d, d1d, q1d, w1d)
180 1,
[all …]
/libCEED/python/tests/
H A Dbuildmats.py10 qref[1] = 0.6
11 qref[2] = 1. / 3.
15 qref[6] = 1. / 3.
18 qweight[1] = 25. / 96.
25 x2 = qref[1 * Q + i]
27 interp[i * P + 0] = 2. * (x1 + x2 - 1.) * (x1 + x2 - 1. / 2.)
28 interp[i * P + 1] = -4. * x1 * (x1 + x2 - 1.)
29 interp[i * P + 2] = 2. * x1 * (x1 - 1. / 2.)
30 interp[i * P + 3] = -4. * x2 * (x1 + x2 - 1.)
32 interp[i * P + 5] = 2. * x2 * (x2 - 1. / 2.)
[all …]
H A Dtest-5-operator.py47 nx = nelem + 1
48 nu = nelem * (p - 1) + 1
65 indx[2 * i + 1] = i + 1
66 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
72 indu[p * i + j] = i * (p - 1) + j
73 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
75 strides = np.array([1, q, q], dtype="int32")
76 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
79 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
80 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
[all …]
/libCEED/tests/
H A Dt320-basis-f.h21 qref=(/2.d-1,6.d-1,1.d0/3.d0,2.d-1,2.d-1,2.d-1,1.d0/3.d0,6.d-1/)
24 do i=0,q-1
25 x1 = qref(0*q+i+1)
26 x2 = qref(1*q+i+1);
28 interp(i*P+1)=2.*(x1+x2-1.)*(x1+x2-1./2.);
29 interp(i*P+2)=-4.*x1*(x1+x2-1.);
30 interp(i*P+3)=2.*x1*(x1-1./2.);
31 interp(i*P+4)=-4.*x2*(x1+x2-1.);
33 interp(i*P+6)=2.*x2*(x2-1./2.);
35 grad((i+0)*P+1)=2.*(1.*(x1+x2-1./2.)+(x1+x2-1.)*1.);
[all …]
H A Dt320-basis.h14 q_ref[1] = 0.6; in Build2DSimplex()
15 q_ref[2] = 1. / 3.; in Build2DSimplex()
19 q_ref[6] = 1. / 3.; in Build2DSimplex()
22 q_weight[1] = 25. / 96.; in Build2DSimplex()
28 CeedScalar x1 = q_ref[0 * Q + i], x2 = q_ref[1 * Q + i]; in Build2DSimplex()
30 interp[i * P + 0] = 2. * (x1 + x2 - 1.) * (x1 + x2 - 1. / 2.); in Build2DSimplex()
31 interp[i * P + 1] = -4. * x1 * (x1 + x2 - 1.); in Build2DSimplex()
32 interp[i * P + 2] = 2. * x1 * (x1 - 1. / 2.); in Build2DSimplex()
33 interp[i * P + 3] = -4. * x2 * (x1 + x2 - 1.); in Build2DSimplex()
35 interp[i * P + 5] = 2. * x2 * (x2 - 1. / 2.); in Build2DSimplex()
[all …]
H A Dt593-operator.c2 /// Test 1D mass matrix operator at points with heterogeneous points per element
3 /// \test Test 1D mass matrix operator at points with heterogeneous points per element
13 CeedInt num_elem = 3, dim = 1, p = 3, q = 5; in main()
14 …CeedInt num_nodes_x = num_elem + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem =… in main()
15 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points]; in main()
24 CeedInit(argv[1], &ceed); in main()
27 for (CeedInt i = 0; i < num_nodes_x; i++) x_array_mesh[i] = (CeedScalar)i / (num_nodes_x - 1); in main()
30 ind_x[2 * i + 1] = i + 1; in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_nodes_x, CEED_MEM_HOST, CEED_USE_POINTER, i… in main()
39 ind_u[p * i + j] = i * (p - 1) + j; in main()
[all …]
/libCEED/examples/solids/qfunctions/
H A Dlinear.h30 …CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasResidual_Linear()
35 // (*grad_u)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[1]; in ElasResidual_Linear()
46 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasResidual_Linear()
47 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasResidual_Linear()
48 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasResidual_Linear()
53 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasResidual_Linear()
59 // dXdx = (dx/dX)^(-1) in ElasResidual_Linear()
70 // e = 1/2 (grad u + (grad u)^T) in ElasResidual_Linear()
73 …{(grad_u[0][0] + grad_u[0][0]) / 2., (grad_u[0][1] + grad_u[1][0]) / 2., (grad_u[0][2] + grad_u[2]… in ElasResidual_Linear()
74 …{(grad_u[1][0] + grad_u[0][1]) / 2., (grad_u[1][1] + grad_u[1][1]) / 2., (grad_u[1][2] + grad_u[2]… in ElasResidual_Linear()
[all …]
/libCEED/examples/ceed/
H A Dex3-volume.h20 // in[1] is quadrature weights, size (Q) in build_mass_diff()
21 const CeedScalar *w = in[1]; in build_mass_diff()
28 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_mass_diff()
35 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
43 // 1 3 2 1 -J10 J00 in build_mass_diff()
45 const CeedScalar J10 = J[0][1][i]; in build_mass_diff()
46 const CeedScalar J01 = J[1][0][i]; in build_mass_diff()
47 const CeedScalar J11 = J[1][1][i]; in build_mass_diff()
54 q_data[1][i] = qw * (J01 * J01 + J11 * J11); in build_mass_diff()
69 // A[k][j] = J[j+1][k+1]*J[j+2][k+2] - J[j+1][k+2]*J[j+2][k+1] in build_mass_diff()
[all …]
H A Dex1-volume-f.f908 ! libCEED Example 1
11 ! Arbitrary mesh and solution degrees in 1D, 2D and 3D are supported from the same code.
28 ! Test in 1D-3D
29 ! TESTARGS(name = "1D User QFunction") -ceed {ceed_resource} -d 1 -t
32 ! TESTARGS(name = "1D Gallery QFunction") -ceed {ceed_resource} -d 1 -t -g
58 do while (num_elem > 1)
60 s = s + 1
64 do d = 1, fe_dim
67 sd = sd + 1
68 r = r - 1
[all …]
H A Dex2-surface.h20 // in[1] is quadrature weights, size (Q) in build_diff()
21 const CeedScalar *w = in[1]; in build_diff()
28 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_diff()
37 // 1 3 2 1 -J10 J00 in build_diff()
39 const CeedScalar J10 = J[0][1][i]; in build_diff()
40 const CeedScalar J01 = J[1][0][i]; in build_diff()
41 const CeedScalar J11 = J[1][1][i]; in build_diff()
45 q_data[1][i] = qw * (J00 * J00 + J10 * J10); in build_diff()
59 // A[k][j] = J[j+1][k+1]*J[j+2][k+2] - J[j+1][k+2]*J[j+2][k+1] in build_diff()
61 …J[(k + 1) % 3][(j + 1) % 3][i] * J[(k + 2) % 3][(j + 2) % 3][i] - J[(k + 2) % 3][(j + 1) % 3][i] *… in build_diff()
[all …]
/libCEED/examples/python/qfunctions/
H A Dex3-volume.h15 // in[1] is quadrature weights, size (Q) in build_mass_diff()
16 const CeedScalar *w = in[1]; in build_mass_diff()
23 case 11: { // dim = 1, space_dim = 1 in build_mass_diff()
24 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_mass_diff()
31 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
39 // 1 3 2 1 -J10 J00 in build_mass_diff()
41 const CeedScalar J10 = J[0][1][i]; in build_mass_diff()
42 const CeedScalar J01 = J[1][0][i]; in build_mass_diff()
43 const CeedScalar J11 = J[1][1][i]; in build_mass_diff()
50 q_data[1][i] = qw * (J01 * J01 + J11 * J11); in build_mass_diff()
[all …]
H A Dex2-surface.h15 // in[1] is quadrature weights, size (Q) in build_diff()
16 const CeedScalar *w = in[1]; in build_diff()
24 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_diff()
33 // 1 3 2 1 -J10 J00 in build_diff()
35 const CeedScalar J10 = J[0][1][i]; in build_diff()
36 const CeedScalar J01 = J[1][0][i]; in build_diff()
37 const CeedScalar J11 = J[1][1][i]; in build_diff()
41 q_data[1][i] = qw * (J00 * J00 + J10 * J10); in build_diff()
55 // A[k][j] = J[j+1][k+1]*J[j+2][k+2] - J[j+1][k+2]*J[j+2][k+1] in build_diff()
57 …J[(k + 1) % 3][(j + 1) % 3][i] * J[(k + 2) % 3][(j + 2) % 3][i] - J[(k + 2) % 3][(j + 1) % 3][i] *… in build_diff()
[all …]
/libCEED/include/ceed/jit-source/gallery/
H A Dceed-poisson3dbuild.h16 // in[1] is quadrature weights, size (Q) in Poisson3DBuild()
17 const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], *w = in[1]; in Poisson3DBuild()
30 // A[k][j] = J[k+1][j+1]*J[k+2][j+2] - J[k+2][j+1]*J[k+1][j+2] in Poisson3DBuild()
31 A[k][j] = J[(k + 1) % dim][(j + 1) % dim][i] * J[(k + 2) % dim][(j + 2) % dim][i] - in Poisson3DBuild()
32 J[(k + 2) % dim][(j + 1) % dim][i] * J[(k + 1) % dim][(j + 2) % dim][i]; in Poisson3DBuild()
35 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in Poisson3DBuild()
40 // 5 1 3 in Poisson3DBuild()
42 q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in Poisson3DBuild()
43 q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in Poisson3DBuild()
44 q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in Poisson3DBuild()
[all …]
/libCEED/tests/output/
H A Dt523-operator.out3 SubOperator 1 - quadrilateral elements
11 Size: 1
14 Input field 1:
19 1 output field:
22 Size: 1
25 SubOperator 1 - quadrilateral elements:
31 Size: 1
34 Input field 1:
39 1 output field:
42 Size: 1
[all …]
/libCEED/examples/nek/bps/
H A Dbps.usr10 C_TESTARGS(name="BP1") -c {ceed_resource} -e bp1 -n 1 -b 4 -test
11 C_TESTARGS(name="BP3") -c {ceed_resource} -e bp3 -n 1 -b 4 -test
20 real*8 ctx(1)
31 do i=1,q
36 a21=u2(i+q*1)
64 $ +u1(i+q*1)*u1(i+q*1)
77 real*8 ctx(1)
83 do i=1,q
96 real*8 ctx(1)
108 do i=1,q
[all …]

12345678910>>...35