Home
last modified time | relevance | path

Searched refs:q (Results 1 – 25 of 209) sorted by relevance

123456789

/libCEED/tests/
H A Dt532-operator-f.h2 subroutine setup_mass(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,&
9 integer q,ierr
11 do i=1,q
12 v1(i)=u2(i)*(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
18 subroutine setup_diff(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,&
26 integer q,ierr
28 do i=1,q
29 w=u2(i)/(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
30 v1(i+q*0)=w*(u1(i+q*2)*u1(i+q*2)+u1(i+q*3)*u1(i+q*3))
31 v1(i+q*1)=w*(u1(i+q*0)*u1(i+q*0)+u1(i+q*1)*u1(i+q*1))
[all …]
H A Dt531-operator-f.h2 subroutine setup(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
9 integer q,ierr
11 do i=1,q
12 w=u2(i)/(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
13 v1(i+q*0)=w*(u1(i+q*2)*u1(i+q*2)+u1(i+q*3)*u1(i+q*3))
14 v1(i+q*1)=w*(u1(i+q*0)*u1(i+q*0)+u1(i+q*1)*u1(i+q*1))
15 v1(i+q*2)=-w*(u1(i+q*0)*u1(i+q*2)+u1(i+q*2)*u1(i+q*3))
21 subroutine diff(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
28 integer q,ierr
30 do i=1,q
[all …]
H A Dt522-operator-f.h2 subroutine setup(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
9 integer q,ierr
11 do i=1,q
12 w = u1(i)/(u2(i+0*q)*u2(i+3*q)-u2(i+1*q)*u2(i+2*q))
13 v1(i+0*q)=w*(u2(i+2*q)*u2(i+2*q)+u2(i+3*q)*u2(i+3*q))
14 v1(i+1*q)=w*(u2(i+0*q)*u2(i+0*q)+u2(i+1*q)*u2(i+1*q))
15 v1(i+2*q)=w*(u2(i+0*q)*u2(i+2*q)+u2(i+1*q)*u2(i+3*q))
21 subroutine diff(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
27 integer q,ierr
29 do i=1,q
[all …]
H A Dt534-operator-f.h2 subroutine setup(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
9 integer q,ierr
11 do i=1,q
12 w=u2(i)/(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
13 v1(i+q*0)=w*(u1(i+q*2)*u1(i+q*2)+u1(i+q*3)*u1(i+q*3))
14 v1(i+q*1)=-w*(u1(i+q*0)*u1(i+q*2)+u1(i+q*2)*u1(i+q*3))
15 v1(i+q*2)=w*(u1(i+q*0)*u1(i+q*0)+u1(i+q*1)*u1(i+q*1))
21 subroutine diff(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
28 integer q,ierr
30 do i=1,q
[all …]
H A Dt535-operator-f.h2 subroutine setup_mass(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,&
9 integer q,ierr
11 do i=1,q
12 v1(i)=u2(i)*(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
18 subroutine setup_diff(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,&
26 integer q,ierr
28 do i=1,q
29 w=u2(i)/(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
30 v1(i+q*0)=w*(u1(i+q*2)*u1(i+q*2)+u1(i+q*3)*u1(i+q*3))
31 v1(i+q*1)=w*(u1(i+q*0)*u1(i+q*0)+u1(i+q*1)*u1(i+q*1))
[all …]
H A Dt323-basis.c20 const CeedInt p = 6, q = 4, dim = 2; in main() local
22 CeedScalar q_ref[dim * q], q_weight[q]; in main()
23 CeedScalar interp[p * q], grad[dim * p * q]; in main()
37 CeedVectorCreate(ceed, q * dim, &v); in main()
41 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
50 for (int i = 0; i < q; i++) { in main()
51 CeedScalar dfx = EvalGrad(x_q[0 * q + i], x_q[1 * q + i]); in main()
52 …if (fabs(v_array[0 * q + i] - dfx) > 100. * CEED_EPSILON) printf("[%" CeedInt_FMT "] %f != %f\n", … in main()
53 dfx = EvalGrad(x_q[1 * q + i], x_q[0 * q + i]); in main()
54 …if (fabs(v_array[1 * q + i] - dfx) > 100. * CEED_EPSILON) printf("[%" CeedInt_FMT "] %f != %f\n", … in main()
H A Dt410-qfunction.c12 CeedInt q = 8; in main() local
13 CeedScalar v_true[q]; in main()
17 CeedVectorCreate(ceed, q, &dx); in main()
18 CeedVectorCreate(ceed, q, &w); in main()
19 CeedVectorCreate(ceed, q, &u); in main()
21 CeedScalar dx_array[q], w_array[q], u_array[q]; in main()
23 for (CeedInt i = 0; i < q; i++) { in main()
24 CeedScalar x = 2. * i / (q - 1) - 1; in main()
34 CeedVectorCreate(ceed, q, &v); in main()
36 CeedVectorCreate(ceed, q, &q_data); in main()
[all …]
H A Dt410-qfunction-f.f909 integer q,i variable
10 parameter(q=8)
11 real*8 jj(q)
12 real*8 ww(q)
13 real*8 uu(q)
14 real*8 vv(q)
15 real*8 vvv(q)
26 do i=0,q-1
28 x=2.0*i/(q-1)-1
34 call ceedvectorcreate(ceed,q,j,err)
[all …]
H A Dt405-qfunction.c14 CeedInt q = 8; in main() local
15 CeedScalar v_true[q]; in main()
19 CeedVectorCreate(ceed, q, &w); in main()
20 CeedVectorCreate(ceed, q, &u); in main()
22 CeedScalar w_array[q], u_array[q]; in main()
24 for (CeedInt i = 0; i < q; i++) { in main()
25 CeedScalar x = 2. * i / (q - 1) - 1; in main()
33 CeedVectorCreate(ceed, q, &v); in main()
35 CeedVectorCreate(ceed, q, &q_data); in main()
44 CeedQFunctionApply(qf_setup, q, in, out); in main()
[all …]
H A Dt400-qfunction.c14 CeedInt q = 8; in main() local
15 CeedScalar v_true[q]; in main()
19 CeedVectorCreate(ceed, q, &w); in main()
20 CeedVectorCreate(ceed, q, &u); in main()
22 CeedScalar w_array[q], u_array[q]; in main()
24 for (CeedInt i = 0; i < q; i++) { in main()
25 CeedScalar x = 2. * i / (q - 1) - 1; in main()
33 CeedVectorCreate(ceed, q, &v); in main()
35 CeedVectorCreate(ceed, q, &q_data); in main()
44 CeedQFunctionApply(qf_setup, q, in, out); in main()
[all …]
H A Dt323-basis-f.f9027 integer p,q,d variable
29 parameter(q=4)
32 real*8 qref(d*q)
33 real*8 qweight(q)
34 real*8 interp(p*q)
35 real*8 grad(d*p*q)
36 real*8 xq(d*q)
39 real*8 ooutput(d*q)
57 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
71 call ceedvectorcreate(ceed,q*d,output,err)
[all …]
H A Dt400-qfunction-f.f9014 integer q,i
15 parameter(q=8)
16 real*8 ww(q)
17 real*8 uu(q)
18 real*8 vv(q)
19 real*8 vvv(q)
42 do i=0,q-1
43 x=2.0*i/(q-1)-1
49 call ceedvectorcreate(ceed,q,w,err)
52 call ceedvectorcreate(ceed,q,u,err)
[all …]
H A Dt401-qfunction.c16 CeedInt q = 8; in main() local
17 CeedScalar v_true[q], ctx_data[5] = {1, 2, 3, 4, 5}; in main()
21 CeedVectorCreate(ceed, q, &w); in main()
22 CeedVectorCreate(ceed, q, &u); in main()
24 CeedScalar w_array[q], u_array[q]; in main()
26 for (CeedInt i = 0; i < q; i++) { in main()
27 CeedScalar x = 2. * i / (q - 1) - 1; in main()
36 CeedVectorCreate(ceed, q, &v); in main()
38 CeedVectorCreate(ceed, q, &q_data); in main()
47 CeedQFunctionApply(qf_setup, q, in, out); in main()
[all …]
H A Dt406-qfunction.c17 CeedInt q = 8; in main() local
18 CeedScalar v_true[q]; in main()
30 CeedVectorCreate(ceed, q, &w); in main()
31 CeedVectorCreate(ceed, q, &u); in main()
33 CeedScalar w_array[q], u_array[q]; in main()
35 for (CeedInt i = 0; i < q; i++) { in main()
36 CeedScalar x = 2. * i / (q - 1) - 1; in main()
44 CeedVectorCreate(ceed, q, &v); in main()
46 CeedVectorCreate(ceed, q, &q_data); in main()
55 CeedQFunctionApply(qf_setup, q, in, out); in main()
[all …]
H A Dt330-basis.h44 static void BuildHdivQuadrilateral(CeedInt q, CeedScalar *q_ref, CeedScalar *q_weights, CeedScalar … in BuildHdivQuadrilateral() argument
46 CeedScalar q_ref_1d[q], q_weight_1d[q]; in BuildHdivQuadrilateral()
49 CeedGaussQuadrature(q, q_ref_1d, q_weight_1d); in BuildHdivQuadrilateral()
53 CeedLobattoQuadrature(q, q_ref_1d, q_weight_1d); in BuildHdivQuadrilateral()
64 for (CeedInt i = 0; i < q; i++) { in BuildHdivQuadrilateral()
65 for (CeedInt j = 0; j < q; j++) { in BuildHdivQuadrilateral()
66 CeedInt k1 = q * i + j; in BuildHdivQuadrilateral()
68 q_ref[k1 + q * q] = q_ref_1d[i]; in BuildHdivQuadrilateral()
75 interp[k1 * 8 + k + 8 * q * q] = By[k]; in BuildHdivQuadrilateral()
H A Dt401-qfunction-f.f9015 integer q,i
16 parameter(q=8)
17 real*8 ww(q)
18 real*8 uu(q)
19 real*8 vv(q)
20 real*8 vvv(q)
54 do i=0,q-1
55 x=2.0*i/(q-1)-1
61 call ceedvectorcreate(ceed,q,w,err)
64 call ceedvectorcreate(ceed,q,u,err)
[all …]
H A Dt540-operator-f.h2 subroutine setup_mass(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,&
9 integer q,ierr
11 do i=1,q
12 v1(i)=u2(i)*(u1(i+q*0)*u1(i+q*3)-u1(i+q*1)*u1(i+q*2))
18 subroutine apply(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,&
24 integer q,ierr
26 do i=1,q
/libCEED/examples/rust/ex3-volume-vector/src/
H A Dmain.rs145 let q = qdata.len() / 2; in example_3_vector() localVariable
146 for i in 0..q { in example_3_vector()
148 qdata[i + q * 0] = weights[i] * jacobian[i]; in example_3_vector()
150 qdata[i + q * 1] = weights[i] / jacobian[i]; in example_3_vector()
154 let q = qdata.len() / 4; in example_3_vector() localVariable
155 for i in 0..q { in example_3_vector()
156 let j11 = jacobian[i + q * 0]; in example_3_vector()
157 let j21 = jacobian[i + q * 1]; in example_3_vector()
158 let j12 = jacobian[i + q * 2]; in example_3_vector()
159 let j22 = jacobian[i + q * 3]; in example_3_vector()
[all …]
/libCEED/examples/rust/ex3-volume/src/
H A Dmain.rs143 let q = qdata.len() / 2; in example_3() localVariable
144 for i in 0..q { in example_3()
146 qdata[i + q * 0] = weights[i] * jacobian[i]; in example_3()
148 qdata[i + q * 1] = weights[i] / jacobian[i]; in example_3()
152 let q = qdata.len() / 4; in example_3() localVariable
153 for i in 0..q { in example_3()
154 let j11 = jacobian[i + q * 0]; in example_3()
155 let j21 = jacobian[i + q * 1]; in example_3()
156 let j12 = jacobian[i + q * 2]; in example_3()
157 let j22 = jacobian[i + q * 3]; in example_3()
[all …]
/libCEED/examples/ceed/
H A Dex1-volume-f.h2 subroutine build_mass(ctx, q, j, w, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12, u13, u14, u15, u16,&
12 integer q, ierr
19 do i = 1, q
24 do i = 1, q
25 qdata(i) = (j(0*q + i)*j(3*q + i) - j(1*q + i)*j(2*q + i)) * w(i)
29 do i = 1, q
30 qdata(i) = (j(0*q + i) * (j(4*q + i)*j(8*q + i) - j(5*q + i)*j(7*q + i)) -&
31 &j(1*q + i) * (j(3*q + i)*j(8*q + i) - j(5*q + i)*j(6*q + i)) +&
32 &j(2*q + i) * (j(3*q + i)*j(7*q + i) - j(4*q + i)*j(6*q + i))) * w(i)
39 subroutine apply_mass(ctx, q, u, qdata, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12, u13, u14, u15, u…
[all …]
/libCEED/python/tests/
H A Dtest-4-qfunction.py58 q = 8
60 w_array = np.zeros(q, dtype=ceed.scalar_type())
61 u_array = np.zeros(q, dtype=ceed.scalar_type())
62 v_true = np.zeros(q, dtype=ceed.scalar_type())
63 for i in range(q):
64 x = 2. * i / (q - 1) - 1
69 dx = ceed.Vector(q)
71 w = ceed.Vector(q)
73 u = ceed.Vector(q)
75 v = ceed.Vector(q)
[all …]
/libCEED/examples/rust/ex2-surface-vector/src/
H A Dmain.rs155 let q = qdata.len() / 3; in example_2_vector() localVariable
156 for i in 0..q { in example_2_vector()
157 let j11 = jacobian[i + q * 0]; in example_2_vector()
158 let j21 = jacobian[i + q * 1]; in example_2_vector()
159 let j12 = jacobian[i + q * 2]; in example_2_vector()
160 let j22 = jacobian[i + q * 3]; in example_2_vector()
162 qdata[i + q * 0] = qw * (j12 * j12 + j22 * j22); in example_2_vector()
163 qdata[i + q * 1] = qw * (j11 * j11 + j21 * j21); in example_2_vector()
164 qdata[i + q * 2] = -qw * (j11 * j12 + j21 * j22); in example_2_vector()
168 let q = qdata.len() / 6; in example_2_vector() localVariable
[all …]
/libCEED/examples/rust/ex2-surface/src/
H A Dmain.rs152 let q = qdata.len() / 3; in example_2() localVariable
153 for i in 0..q { in example_2()
154 let j11 = jacobian[i + q * 0]; in example_2()
155 let j21 = jacobian[i + q * 1]; in example_2()
156 let j12 = jacobian[i + q * 2]; in example_2()
157 let j22 = jacobian[i + q * 3]; in example_2()
159 qdata[i + q * 0] = qw * (j12 * j12 + j22 * j22); in example_2()
160 qdata[i + q * 1] = qw * (j11 * j11 + j21 * j21); in example_2()
161 qdata[i + q * 2] = -qw * (j11 * j12 + j21 * j22); in example_2()
165 let q = qdata.len() / 6; in example_2() localVariable
[all …]
/libCEED/julia/LibCEED.jl/src/
H A DQuadrature.jl2 gauss_quadrature(q)
4 Return the Gauss-Legendre quadrature rule with `q` points (integrates polynomials of degree
5 $2q-1$ exactly).
9 function gauss_quadrature(q) argument
10 x = zeros(CeedScalar, q)
11 w = zeros(CeedScalar, q)
12 C.CeedGaussQuadrature(q, x, w)
21 lobatto_quadrature(q, mode::Mode=Abscissa)
23 Return the Gauss-Lobatto quadrature rule with `q` points (integrates polynomials of degree
24 $2q-3$ exactly).
[all …]
/libCEED/examples/rust/ex1-volume-vector/src/
H A Dmain.rs147 let q = qdata.len(); in example_1_vector() localVariable
150 *qdata = (jacobian[i + q * 0] * jacobian[i + q * 3] in example_1_vector()
151 - jacobian[i + q * 1] * jacobian[i + q * 2]) in example_1_vector()
157 let q = qdata.len(); in example_1_vector() localVariable
160 *qdata = (jacobian[i + q * 0] in example_1_vector()
161 * (jacobian[i + q * 4] * jacobian[i + q * 8] in example_1_vector()
162 - jacobian[i + q * 5] * jacobian[i + q * 7]) in example_1_vector()
163 - jacobian[i + q * 1] in example_1_vector()
164 * (jacobian[i + q * 3] * jacobian[i + q * 8] in example_1_vector()
165 - jacobian[i + q * 5] * jacobian[i + q * 6]) in example_1_vector()
[all …]

123456789