| /libCEED/julia/LibCEED.jl/test/ |
| H A D | buildmats.jl | 3 interp = Array{T}(undef, dim, Q, P) 24 interp[1, i, 1] = 0.0 25 interp[2, i, 1] = 1.0 - x2 26 interp[1, i, 2] = x1 - 1.0 27 interp[2, i, 2] = 0.0 28 interp[1, i, 3] = -x1 29 interp[2, i, 3] = 0.0 30 interp[1, i, 4] = 0.0 31 interp[2, i, 4] = x2 39 return interp, div [all …]
|
| /libCEED/python/tests/ |
| H A D | buildmats.py | 6 interp = np.empty(P * Q, dtype=mat_dtype) 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.) 31 interp[i * P + 4] = 4. * x1 * x2 32 interp[i * P + 5] = 2. * x2 * (x2 - 1. / 2.) 49 return interp, grad 54 interp = np.empty(dim * P * Q, dtype=mat_dtype) 75 interp[(i + 0) * P + 0] = 0. [all …]
|
| /libCEED/tests/ |
| H A D | t340-basis.h | 12 static void BuildHcurl2DSimplex(CeedScalar *q_ref, CeedScalar *q_weight, CeedScalar *interp, CeedSc… in BuildHcurl2DSimplex() argument 32 interp[(i + 0) * P + 0] = 2. * x2 * (1. - 4. * x1); in BuildHcurl2DSimplex() 33 interp[(i + Q) * P + 0] = 4. * x1 * (2. * x1 - 1.); in BuildHcurl2DSimplex() 34 interp[(i + 0) * P + 1] = 4. * x2 * (1. - 2. * x2); in BuildHcurl2DSimplex() 35 interp[(i + Q) * P + 1] = 2. * x1 * (4. * x2 - 1.); in BuildHcurl2DSimplex() 36 interp[(i + 0) * P + 2] = 2. * x2 * (-4. * x1 - 4. * x2 + 3.); in BuildHcurl2DSimplex() 37 interp[(i + Q) * P + 2] = 8. * x1 * x1 + 8. * x1 * x2 - 12. * x1 - 6. * x2 + 4.; in BuildHcurl2DSimplex() 38 interp[(i + 0) * P + 3] = 4. * x2 * (2. * x2 - 1.); in BuildHcurl2DSimplex() 39 interp[(i + Q) * P + 3] = -8. * x1 * x2 + 2. * x1 + 6. * x2 - 2.; in BuildHcurl2DSimplex() 40 interp[(i + 0) * P + 4] = 8. * x1 * x2 - 6. * x1 + 8. * x2 * x2 - 12. * x2 + 4.; in BuildHcurl2DSimplex() [all …]
|
| H A D | t320-basis.h | 10 static void Build2DSimplex(CeedScalar *q_ref, CeedScalar *q_weight, CeedScalar *interp, CeedScalar … in Build2DSimplex() argument 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() 34 interp[i * P + 4] = 4. * x1 * x2; in Build2DSimplex() 35 interp[i * P + 5] = 2. * x2 * (x2 - 1. / 2.); in Build2DSimplex()
|
| H A D | t320-basis-f.h | 9 subroutine buildmats(qref,qweight,interp,grad) 17 real*8 interp(p*q) 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.); 32 interp(i*P+5)=4.*x1*x2; 33 interp(i*P+6)=2.*x2*(x2-1./2.);
|
| 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() 29 …if (fabs(interp[i] - interp_in_basis[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", interp[i], in… in main() 36 row_sum[i] += interp[j + i * p]; in main() 42 column_sum[i] += interp[i + j * p]; in main()
|
| H A D | t331-basis.c | 16 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main() local 20 BuildHdivQuadrilateral(q, q_ref, q_weights, interp, div, CEED_GAUSS); in main() 21 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main() 29 …if (fabs(interp[i] - interp_in_basis[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", interp[i], in… in main()
|
| H A D | t319-basis.h | 60 …ation(CeedQuadMode quad_mode, CeedInt *P, CeedInt *Q, const CeedScalar **interp, const CeedScalar … in GetCellToFaceTabulation() argument 65 *interp = (const CeedScalar *)Q1_interp_gauss; in GetCellToFaceTabulation() 69 *interp = (const CeedScalar *)Q1_interp_gauss_lobatto; in GetCellToFaceTabulation()
|
| H A D | t320-basis-f.f90 | 19 real*8 interp(p*q) 30 call buildmats(qref,qweight,interp,grad) 32 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
|
| 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 | t320-basis.c | 13 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 20 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 21 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t330-basis.c | 14 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main() local 21 BuildHdivQuadrilateral(q, q_ref, q_weights, interp, div, CEED_GAUSS); in main() 22 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main()
|
| H A D | t324-basis.c | 16 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 26 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 27 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t363-basis.c | 16 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 26 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 27 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t319-basis.c | 159 const CeedScalar *interp, *grad; in main() local 165 CeedBasisGetInterp(basis_from, &interp); in main() 167 …CeedBasisCreateH1(ceed, topo, num_comp, num_nodes, num_qpts, interp, grad, NULL, NULL, &basis_from… in main() 173 CeedBasisGetInterp(basis_to, &interp); in main() 175 …CeedBasisCreateH1(ceed, topo, num_comp, num_nodes, num_qpts, interp, grad, NULL, NULL, &basis_to_n… in main() 212 const CeedScalar *grad, *interp; in main() local 214 GetCellToFaceTabulation(CEED_GAUSS, &P, &Q, &interp, &grad); in main() 217 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_HEX, 1, P, Q, (CeedScalar *)interp, (CeedScalar *)grad, q_re… in main()
|
| H A D | t362-basis.c | 17 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 29 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 30 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t321-basis.c | 20 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 37 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 38 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t510-operator.c | 26 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 71 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 72 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p, q, interp, grad, q_ref, q_weight, &basis_x… in main() 74 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 75 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis_u); in main()
|
| H A D | t511-operator.c | 25 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 70 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 71 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p, q, interp, grad, q_ref, q_weight, &basis_x… in main() 73 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 74 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis_u); in main()
|
| H A D | t332-basis.c | 16 CeedScalar interp[dim * p * num_qpts], div[p * num_qpts]; in main() local 20 BuildHdivQuadrilateral(q, q_ref, q_weights, interp, div, CEED_GAUSS); in main() 21 …CeedBasisCreateHdiv(ceed, CEED_TOPOLOGY_QUAD, 1, p, num_qpts, interp, div, q_ref, q_weights, &basi… in main()
|
| H A D | t322-basis.c | 20 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 38 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 39 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| /libCEED/tests/output/ |
| H A D | t320-basis.out | 8 interp: 29 interp:
|
| H A D | t320-basis-f.out | 8 interp: 29 interp:
|
| H A D | t330-basis.out | 8 interp: 44 interp:
|
| H A D | t340-basis.out | 8 interp: 29 interp:
|