| /libCEED/tests/ |
| H A D | t340-basis.h | 30 CeedScalar x1 = q_ref[0 * Q + i], x2 = q_ref[1 * Q + i]; in BuildHcurl2DSimplex() local 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() 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() 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() 41 interp[(i + Q) * P + 4] = 2. * x1 * (-4. * x1 - 4. * x2 + 3.); in BuildHcurl2DSimplex() 42 interp[(i + 0) * P + 5] = -8. * x1 * x2 + 6. * x1 + 2. * x2 - 2.; in BuildHcurl2DSimplex() [all …]
|
| H A D | t320-basis.h | 28 CeedScalar x1 = q_ref[0 * Q + i], x2 = q_ref[1 * Q + i]; in Build2DSimplex() local 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() 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() [all …]
|
| H A D | t320-basis-f.h | 19 real*8 x1,x2 25 x1 = qref(0*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.); 32 interp(i*P+5)=4.*x1*x2; 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.); [all …]
|
| H A D | t323-basis-f.f90 | 7 subroutine feval(x1,x2,val) argument 8 real*8 x1,x2,val local 10 val=x1*x1+x2*x2+x1*x2+1 14 subroutine dfeval(x1,x2,val) argument 15 real*8 x1,x2,val local 17 val=2*x1+x2 41 real*8 x1,x2 variable 61 x1=xr(0*p+i) 63 call feval(x1,x2,val) 78 x1=xq(0*q+i) [all …]
|
| H A D | t322-basis-f.f90 | 7 subroutine feval(x1,x2,val) argument 8 real*8 x1,x2,val local 10 val=x1*x1+x2*x2+x1*x2+1 34 real*8 x1,x2 variable 54 x1=xr(0*p+i) 56 call feval(x1,x2,val)
|
| H A D | t323-basis.c | 11 static CeedScalar Eval(CeedScalar x1, CeedScalar x2) { return x1 * x1 + x2 * x2 + x1 * x2 + 1; } in Eval() argument 14 CeedScalar EvalGrad(CeedScalar x1, CeedScalar x2) { return 2 * x1 + x2; } in EvalGrad() argument
|
| H A D | t321-basis.c | 11 static CeedScalar Eval(CeedScalar x1, CeedScalar x2) { return x1 * x1 + x2 * x2 + x1 * x2 + 1; } in Eval() argument
|
| H A D | t322-basis.c | 11 static CeedScalar Eval(CeedScalar x1, CeedScalar x2) { return x1 * x1 + x2 * x2 + x1 * x2 + 1; } in Eval() argument
|
| /libCEED/python/tests/ |
| H A D | buildmats.py | 24 x1 = qref[0 * 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.) 31 interp[i * P + 4] = 4. * x1 * x2 35 (1. * (x1 + x2 - 1. / 2.) + (x1 + x2 - 1.) * 1.) 37 (1. * (x1 + x2 - 1. / 2.) + (x1 + x2 - 1.) * 1.) 38 grad[(i + 0) * P + 1] = -4. * (1. * (x1 + x2 - 1.) + x1 * 1.) 39 grad[(i + Q) * P + 1] = -4. * (x1 * 1.) [all …]
|
| H A D | test-3-basis.py | 33 def feval(x1, x2): argument 34 return x1 * x1 + x2 * x2 + x1 * x2 + 1 37 def dfeval(x1, x2): argument 38 return 2 * x1 + x2
|
| /libCEED/julia/LibCEED.jl/test/ |
| H A D | buildmats.jl | 21 x1 = qref[1, i] 26 interp[1, i, 2] = x1 - 1.0 28 interp[1, i, 3] = -x1 62 x1 = qref[1, i] 66 interp[2, i, 1] = x1 68 interp[2, i, 2] = 1.0 - x1 70 interp[2, i, 3] = x1
|
| /libCEED/julia/LibCEED.jl/docs/src/ |
| H A D | Misc.md | 3 User Q-functions often perform small (1x1, 2x2, or 3x3) linear algebra
|
| /libCEED/examples/nek/bps/ |
| H A D | bps.usr | 360 x1 = 1 361 call rescale_x(xm1,x0,x1) 362 call rescale_x(ym1,x0,x1) 363 call rescale_x(zm1,x0,x1)
|