| /libCEED/tests/ |
| H A D | t320-basis.h | 10 …oid Build2DSimplex(CeedScalar *q_ref, CeedScalar *q_weight, CeedScalar *interp, CeedScalar *grad) { in Build2DSimplex() argument 37 grad[(i + 0) * P + 0] = 2. * (1. * (x1 + x2 - 1. / 2.) + (x1 + x2 - 1.) * 1.); in Build2DSimplex() 38 grad[(i + Q) * P + 0] = 2. * (1. * (x1 + x2 - 1. / 2.) + (x1 + x2 - 1.) * 1.); in Build2DSimplex() 39 grad[(i + 0) * P + 1] = -4. * (1. * (x1 + x2 - 1.) + x1 * 1.); in Build2DSimplex() 40 grad[(i + Q) * P + 1] = -4. * (x1 * 1.); in Build2DSimplex() 41 grad[(i + 0) * P + 2] = 2. * (1. * (x1 - 1. / 2.) + x1 * 1.); in Build2DSimplex() 42 grad[(i + Q) * P + 2] = 2. * 0.; in Build2DSimplex() 43 grad[(i + 0) * P + 3] = -4. * (x2 * 1.); in Build2DSimplex() 44 grad[(i + Q) * P + 3] = -4. * (1. * (x1 + x2 - 1.) + x2 * 1.); in Build2DSimplex() 45 grad[(i + 0) * P + 4] = 4. * (1. * x2); in Build2DSimplex() [all …]
|
| H A D | t320-basis-f.h | 9 subroutine buildmats(qref,qweight,interp,grad) 18 real*8 grad(d*p*q) 35 grad((i+0)*P+1)=2.*(1.*(x1+x2-1./2.)+(x1+x2-1.)*1.); 36 grad((i+Q)*P+1)=2.*(1.*(x1+x2-1./2.)+(x1+x2-1.)*1.); 37 grad((i+0)*P+2)=-4.*(1.*(x1+x2-1.)+x1*1.); 38 grad((i+Q)*P+2)=-4.*(x1*1.); 39 grad((i+0)*P+3)=2.*(1.*(x1-1./2.)+x1*1.); 40 grad((i+Q)*P+3)=2.*0.; 41 grad((i+0)*P+4)=-4.*(x2*1.); 42 grad((i+Q)*P+4)=-4.*(1.*(x1+x2-1.)+x2*1.); [all …]
|
| H A D | t319-basis.h | 60 …edQuadMode quad_mode, CeedInt *P, CeedInt *Q, const CeedScalar **interp, const CeedScalar **grad) { in GetCellToFaceTabulation() argument 66 *grad = (const CeedScalar *)Q1_grad_gauss; in GetCellToFaceTabulation() 70 *grad = (const CeedScalar *)Q1_grad_gauss_lobatto; in GetCellToFaceTabulation()
|
| 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() 35 column_sum[i] += grad[i + j * p]; 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() 36 column_sum[i] += grad[i + j * p]; in main()
|
| H A D | t320-basis-f.f90 | 20 real*8 grad(d*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 | 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 | t325-basis.c | 16 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 21 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 22 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, num_comp, p, q, interp, grad, q_ref, q_weight, &ba… in main() 45 column_sum[i] += grad[i + j * p]; in main()
|
| H A D | t319-basis.c | 159 const CeedScalar *interp, *grad; in main() local 166 CeedBasisGetGrad(basis_from, &grad); in main() 167 …CeedBasisCreateH1(ceed, topo, num_comp, num_nodes, num_qpts, interp, grad, NULL, NULL, &basis_from… in main() 174 CeedBasisGetGrad(basis_to, &grad); 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 | 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()
|
| H A D | t526-operator.c | 35 CeedScalar interp[p_tet * q_tet], grad[dim * p_tet * q_tet]; in main() local 72 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 73 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p_tet, q_tet, interp, grad, q_ref, q_weight, … in main() 75 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 76 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p_tet, q_tet, interp, grad, q_ref, q_weight, &b… in main()
|
| H A D | t582-operator.c | 27 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 79 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 80 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p, q, interp, grad, q_ref, q_weight, &basis_x… in main() 82 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 83 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis_u); in main()
|
| H A D | t581-operator.c | 27 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 76 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 77 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p, q, interp, grad, q_ref, q_weight, &basis_x… in main() 79 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 80 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis_u); in main()
|
| H A D | t323-basis.c | 23 CeedScalar interp[p * q], grad[dim * p * q]; in main() local 40 Build2DSimplex(q_ref, q_weight, interp, grad); in main() 41 CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, 1, p, q, interp, grad, q_ref, q_weight, &basis); in main()
|
| H A D | t322-basis-f.f90 | 28 real*8 grad(d*p*q) variable 46 call buildmats(qref,qweight,interp,grad) 50 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
|
| H A D | t511-operator-f.f90 | 40 real*8 grad(d*p*q) 88 call buildmats(qref,qweight,interp,grad) 89 call ceedbasiscreateh1(ceed,ceed_triangle,d,p,q,interp,grad,qref,qweight,& 91 call buildmats(qref,qweight,interp,grad) 92 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
|
| H A D | t510-operator-f.f90 | 40 real*8 grad(d*p*q) 87 call buildmats(qref,qweight,interp,grad) 88 call ceedbasiscreateh1(ceed,ceed_triangle,d,p,q,interp,grad,qref,qweight,& 90 call buildmats(qref,qweight,interp,grad) 91 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
|
| H A D | t323-basis-f.f90 | 35 real*8 grad(d*p*q) variable 53 call buildmats(qref,qweight,interp,grad) 57 call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,&
|
| /libCEED/python/tests/ |
| H A D | buildmats.py | 7 grad = np.empty(dim * P * Q, dtype=mat_dtype) 34 grad[(i + 0) * P + 0] = 2. * \ 36 grad[(i + Q) * P + 0] = 2. * \ 38 grad[(i + 0) * P + 1] = -4. * (1. * (x1 + x2 - 1.) + x1 * 1.) 39 grad[(i + Q) * P + 1] = -4. * (x1 * 1.) 40 grad[(i + 0) * P + 2] = 2. * (1. * (x1 - 1. / 2.) + x1 * 1.) 41 grad[(i + Q) * P + 2] = 2. * 0. 42 grad[(i + 0) * P + 3] = -4. * (x2 * 1.) 43 grad[(i + Q) * P + 3] = -4. * (1. * (x1 + x2 - 1.) + x2 * 1.) 44 grad[(i + 0) * P + 4] = 4. * (1. * x2) [all …]
|
| /libCEED/tests/output/ |
| H A D | t320-basis.out | 13 grad: 34 grad:
|
| H A D | t320-basis-f.out | 13 grad: 34 grad:
|