| /libCEED/tests/ |
| H A D | t533-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny
27 parameter(ndofs=(nx*2+1)*(ny*2+1))
30 real*8 arrx(d*ndofs),aa(nqpts),uu(ndofs),vv(ndofs),atrue(ndofs)
44 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)
45 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)
48 call ceedvectorcreate(ceed,d*ndofs,x,err)
68 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
71 call ceedelemrestrictioncreate(ceed,nelem,p*p,1,1,ndofs,&
123 call ceedvectorcreate(ceed,ndofs,a,err)
128 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t534-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny
27 parameter(ndofs=(nx*2+1)*(ny*2+1))
30 real*8 arrx(d*ndofs),aa(nqpts),uu(ndofs),vv(ndofs),atrue(ndofs)
44 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)
45 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)
48 call ceedvectorcreate(ceed,d*ndofs,x,err)
68 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
71 call ceedelemrestrictioncreate(ceed,nelem,p*p,1,1,ndofs,&
127 call ceedvectorcreate(ceed,ndofs,a,err)
132 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t535-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny
27 parameter(ndofs=(nx*2+1)*(ny*2+1))
30 real*8 arrx(d*ndofs),aa(nqpts),uu(ndofs),vv(ndofs),atrue(ndofs)
44 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)
45 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)
48 call ceedvectorcreate(ceed,d*ndofs,x,err)
69 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
72 call ceedelemrestrictioncreate(ceed,nelem,p*p,1,1,ndofs,&
158 call ceedvectorcreate(ceed,ndofs,a,err)
163 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t531-operator-f.f90 | 25 integer ndofs,nqpts,nx,ny
28 parameter(ndofs=(nx*2+1)*(ny*2+1))
31 real*8 arrx(d*ndofs),uu(ndofs),vv(ndofs),vvlin(ndofs)
45 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)+j*0.5
46 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)+i*0.5
49 call ceedvectorcreate(ceed,d*ndofs,x,err)
59 call ceedvectorcreate(ceed,ndofs,u,err)
62 call ceedvectorcreate(ceed,ndofs,v,err)
63 call ceedvectorcreate(ceed,ndofs,v_lin,err)
81 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
[all …]
|
| H A D | t511-operator-f.f90 | 28 integer ndofs,nqpts,nx,ny
31 parameter(ndofs=(nx*2+1)*(ny*2+1))
34 real*8 arrx(d*ndofs)
43 real*8 hv(ndofs)
53 do i=0,ndofs-1
57 arrx(i+1+ndofs)=val*(1.d0/(ny*2.d0))
79 call ceedelemrestrictioncreate(ceed,nelem,p,d,ndofs,d*ndofs,&
82 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,ndofs,ceed_mem_host,&
114 call ceedvectorcreate(ceed,d*ndofs,x,err)
134 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t510-operator-f.f90 | 28 integer ndofs,nqpts,nx,ny
31 parameter(ndofs=(nx*2+1)*(ny*2+1))
34 real*8 arrx(d*ndofs)
42 real*8 hv(ndofs)
52 do i=0,ndofs-1
56 arrx(i+1+ndofs)=val*(1.d0/(ny*2.d0))
78 call ceedelemrestrictioncreate(ceed,nelem,p,d,ndofs,d*ndofs,&
81 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,ndofs,ceed_mem_host,&
113 call ceedvectorcreate(ceed,d*ndofs,x,err)
133 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t536-operator-f.f90 | 28 integer ndofs,nqpts,nx,ny
31 parameter(ndofs=(nx*2+1)*(ny*2+1))
34 real*8 arrx(d*ndofs),aa(nqpts),uu(ndofs),vv(ndofs),atrue(ndofs)
51 do i=0,ndofs-1
55 arrx(i+1+ndofs)=val*(1.d0/(ny*2.d0))
57 call ceedvectorcreate(ceed,d*ndofs,x,err)
87 call ceedelemrestrictioncreate(ceed,nelem,p,d,ndofs,d*ndofs,&
90 call ceedelemrestrictioncreate(ceed,nelem,p,1,1,ndofs,&
180 call ceedvectorcreate(ceed,ndofs,a,err)
185 call ceedvectorcreate(ceed,ndofs,u,err)
[all …]
|
| H A D | t530-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny
27 parameter(ndofs=(nx*2+1)*(ny*2+1))
30 real*8 arrx(d*ndofs),aa(nqpts),qq(nqpts),vv(ndofs)
45 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)
46 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)
49 call ceedvectorcreate(ceed,d*ndofs,x,err)
69 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
72 call ceedelemrestrictioncreate(ceed,nelem,p*p,1,1,ndofs,&
142 call ceedvectorcreate(ceed,ndofs,u,err)
144 call ceedvectorcreate(ceed,ndofs,v,err)
[all …]
|
| H A D | t524-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex
38 parameter(ndofs=(nx*2+1)*(ny*2+1))
43 real*8 arrx(d*ndofs)
51 real*8 hv(ndofs)
65 arrx(i+j*(ny*2+1)+0*ndofs+1)=1.d0*i/(2*ny)
66 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx)
70 call ceedvectorcreate(ceed,d*ndofs,x,err)
100 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,d,ndofs,d*ndofs,&
103 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,1,1,ndofs,&
165 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,d,ndofs,d*ndofs,&
[all …]
|
| H A D | t521-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex
38 parameter(ndofs=(nx*2+1)*(ny*2+1))
43 real*8 arrx(d*ndofs)
51 real*8 hv(ndofs)
65 arrx(i+j*(ny*2+1)+0*ndofs+1)=1.d0*i/(2*ny)
66 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx)
70 call ceedvectorcreate(ceed,d*ndofs,x,err)
100 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,d,ndofs,d*ndofs,&
103 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,1,1,ndofs,&
165 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,d,ndofs,&
[all …]
|
| H A D | t520-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex
38 parameter(ndofs=(nx*2+1)*(ny*2+1))
43 real*8 arrx(d*ndofs)
51 real*8 hv(ndofs)
64 arrx(i+j*(ny*2+1)+0*ndofs+1)=1.d0*i/(2*ny)
65 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx)
69 call ceedvectorcreate(ceed,d*ndofs,x,err)
99 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,d,ndofs,d*ndofs,&
101 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,1,1,ndofs,&
163 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,d,ndofs,&
[all …]
|
| H A D | t522-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex
38 parameter(ndofs=(nx*2+1)*(ny*2+1))
43 real*8 arrx(d*ndofs)
51 real*8 hv(ndofs)
65 arrx(i+j*(ny*2+1)+0*ndofs+1)=1.d0*i/(2*ny)
66 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx)
70 call ceedvectorcreate(ceed,d*ndofs,x,err)
100 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,d,ndofs,d*ndofs,&
103 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,1,1,ndofs,&
166 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,d,ndofs,&
[all …]
|
| H A D | t532-operator-f.f90 | 25 integer ndofs,nqpts,nx,ny
28 parameter(ndofs=(nx*2+1)*(ny*2+1))
31 real*8 arrx(d*ndofs),vv(ndofs)
46 arrx(i+j*(nx*2+1)+0*ndofs+1)=1.d0*i/(2*nx)
47 arrx(i+j*(nx*2+1)+1*ndofs+1)=1.d0*j/(2*ny)
50 call ceedvectorcreate(ceed,d*ndofs,x,err)
71 call ceedelemrestrictioncreate(ceed,nelem,p*p,d,ndofs,d*ndofs,&
74 call ceedelemrestrictioncreate(ceed,nelem,p*p,1,1,ndofs,&
160 call ceedvectorcreate(ceed,ndofs,u,err)
162 call ceedvectorcreate(ceed,ndofs,v,err)
[all …]
|
| H A D | t540-operator-f.f90 | 23 integer ndofs,nqpts
24 parameter(ndofs=p*p)
26 real*8 arrx(d*nelem*2*2),uu(ndofs)
57 call ceedelemrestrictioncreatestrided(ceed,nelem,p*p,1,ndofs,&
109 call ceedvectorcreate(ceed,ndofs,u,err)
111 call ceedvectorcreate(ceed,ndofs,v,err)
124 do i=1,ndofs
|
| H A D | t523-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex
38 parameter(ndofs=(nx*2+1)*(ny*2+1))
60 call ceedvectorcreate(ceed,d*ndofs,x,err)
88 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,d,ndofs,d*ndofs,&
91 call ceedelemrestrictioncreate(ceed,nelemtet,ptet,1,1,ndofs,&
155 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,d,ndofs,d*ndofs,&
158 call ceedelemrestrictioncreate(ceed,nelemhex,phex*phex,1,1,ndofs,&
|
| /libCEED/python/tests/ |
| H A D | test-5-operator.py | 590 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
594 x = ceed.Vector(dim * ndofs)
595 x_array = np.zeros(dim * ndofs)
596 for i in range(ndofs):
598 x_array[i + ndofs] = (1. / (ny * 2)) * (i / (nx * 2 + 1))
602 u = ceed.Vector(ndofs)
603 v = ceed.Vector(ndofs)
626 rx = ceed.ElemRestriction(nelem, p, dim, ndofs, dim * ndofs, indx,
629 ru = ceed.ElemRestriction(nelem, p, 1, 1, ndofs, indx,
681 for i in range(ndofs):
[all …]
|
| /libCEED/rust/libceed/src/ |
| H A D | lib.rs | 1191 let ndofs = p * nelem - nelem + 1; in ceed_t501() localVariable
1197 let mut u = ceed.vector(ndofs)?; in ceed_t501()
1199 let mut v = ceed.vector(ndofs)?; in ceed_t501()
1216 let ru = ceed.elem_restriction(nelem, p, 1, 1, ndofs, MemType::Host, &indu)?; in ceed_t501()
|