| /libCEED/tests/ |
| H A D | t531-operator.c | 19 CeedInt nx = 3, ny = 2; in main() local 20 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q * q; in main() 30 for (CeedInt i = 0; i < nx * 2 + 1; i++) { in main() 32 x_array[i + j * (nx * 2 + 1) + 0 * num_dofs] = (CeedScalar)i / (2 * nx) + 0.5 * j; in main() 33 x_array[i + j * (nx * 2 + 1) + 1 * num_dofs] = (CeedScalar)j / (2 * ny) + 0.5 * i; in main() 43 for (CeedInt i = 0; i < nx * 2 + 1; i++) { in main() 45 u_array[i + j * (nx * 2 + 1)] = i * nx + j * ny; in main() 57 col = i % nx; in main() 58 row = i / nx; in main() 59 offset = col * (p - 1) + row * (nx * 2 + 1) * (p - 1); in main() [all …]
|
| H A D | t531-operator-f.f90 | 25 integer ndofs,nqpts,nx,ny 26 parameter(nx=3) 28 parameter(ndofs=(nx*2+1)*(ny*2+1)) 43 do i=0,nx*2 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 54 do i=0,nx*2 56 uu(i+j*(nx*2+1)+1)=i*nx+j*ny 70 col=mod(i,nx) 71 row=i/nx [all …]
|
| H A D | t533-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny 25 parameter(nx=3) 27 parameter(ndofs=(nx*2+1)*(ny*2+1)) 42 do i=0,nx*2 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) 57 col=mod(i,nx) 58 row=i/nx 59 offset=col*(p-1)+row*(nx*2+1)*(p-1) 62 indx(p*(p*i+k)+j+1)=offset+k*(nx*2+1)+j
|
| H A D | t533-operator.c | 19 CeedInt nx = 3, ny = 2; in main() local 20 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q * q; in main() 31 for (CeedInt i = 0; i < nx * 2 + 1; i++) { in main() 33 x_array[i + j * (nx * 2 + 1) + 0 * num_dofs] = (CeedScalar)i / (2 * nx); in main() 34 x_array[i + j * (nx * 2 + 1) + 1 * num_dofs] = (CeedScalar)j / (2 * ny); in main() 46 col = i % nx; in main() 47 row = i / nx; in main() 48 offset = col * (p - 1) + row * (nx * 2 + 1) * (p - 1); in main() 50 for (CeedInt k = 0; k < p; k++) ind_x[p * (p * i + k) + j] = offset + k * (nx * 2 + 1) + j; in main()
|
| H A D | t534-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny 25 parameter(nx=3) 27 parameter(ndofs=(nx*2+1)*(ny*2+1)) 42 do i=0,nx*2 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) 57 col=mod(i,nx) 58 row=i/nx 59 offset=col*(p-1)+row*(nx*2+1)*(p-1) 62 indx(p*(p*i+k)+j+1)=offset+k*(nx*2+1)+j
|
| H A D | t530-operator.c | 19 CeedInt nx = 3, ny = 2; in main() local 20 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q * q; in main() 29 for (CeedInt i = 0; i < nx * 2 + 1; i++) { in main() 31 x_array[i + j * (nx * 2 + 1) + 0 * num_dofs] = (CeedScalar)i / (2 * nx); in main() 32 x_array[i + j * (nx * 2 + 1) + 1 * num_dofs] = (CeedScalar)j / (2 * ny); in main() 44 col = i % nx; in main() 45 row = i / nx; in main() 46 offset = col * (p - 1) + row * (nx * 2 + 1) * (p - 1); in main() 48 for (CeedInt k = 0; k < p; k++) ind_x[p * (p * i + k) + j] = offset + k * (nx * 2 + 1) + j; in main()
|
| H A D | t530-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny 25 parameter(nx=3) 27 parameter(ndofs=(nx*2+1)*(ny*2+1)) 43 do i=0,nx*2 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) 58 col=mod(i,nx) 59 row=i/nx 60 offset=col*(p-1)+row*(nx*2+1)*(p-1) 63 indx(p*(p*i+k)+j+1)=offset+k*(nx*2+1)+j
|
| H A D | t510-operator.c | 21 CeedInt nx = 3, ny = 2; in main() local 23 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q; in main() 35 x_array[i] = (1. / (nx * 2)) * (CeedScalar)(i % (nx * 2 + 1)); in main() 36 x_array[i + num_dofs] = (1. / (ny * 2)) * (CeedScalar)(i / (nx * 2 + 1)); in main() 46 col = i % nx; in main() 47 row = i / nx; in main() 48 offset = col * 2 + row * (nx * 2 + 1) * 2; in main()
|
| H A D | t511-operator.c | 20 CeedInt nx = 3, ny = 2; in main() local 22 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q; in main() 34 x_array[i] = (1. / (nx * 2)) * (CeedScalar)(i % (nx * 2 + 1)); in main() 35 x_array[i + num_dofs] = (1. / (ny * 2)) * (CeedScalar)(i / (nx * 2 + 1)); in main() 45 col = i % nx; in main() 46 row = i / nx; in main() 47 offset = col * 2 + row * (nx * 2 + 1) * 2; in main()
|
| H A D | t511-operator-f.f90 | 28 integer ndofs,nqpts,nx,ny 29 parameter(nx=3) 31 parameter(ndofs=(nx*2+1)*(ny*2+1)) 54 arrx(i+1)=mod(i,(nx*2+1)) 55 arrx(i+1)=arrx(i+1)*(1.d0/(nx*2.d0)) 56 val=(i/(nx*2+1)) 60 col=mod(i,nx) 61 row=i/nx 62 offset=col*2+row*(nx*2+1)*2
|
| H A D | t510-operator-f.f90 | 28 integer ndofs,nqpts,nx,ny 29 parameter(nx=3) 31 parameter(ndofs=(nx*2+1)*(ny*2+1)) 53 arrx(i+1)=mod(i,(nx*2+1)) 54 arrx(i+1)=arrx(i+1)*(1.d0/(nx*2.d0)) 55 val=(i/(nx*2+1)) 59 col=mod(i,nx) 60 row=i/nx 61 offset=col*2+row*(nx*2+1)*2
|
| H A D | t535-operator-f.f90 | 24 integer ndofs,nqpts,nx,ny 25 parameter(nx=3) 27 parameter(ndofs=(nx*2+1)*(ny*2+1)) 42 do i=0,nx*2 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) 58 col=mod(i,nx) 59 row=i/nx 60 offset=col*(p-1)+row*(nx*2+1)*(p-1) 63 indx(p*(p*i+k)+j+1)=offset+k*(nx*2+1)+j
|
| H A D | t582-operator.c | 21 CeedInt nx = 3, ny = 2; in main() local 23 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q; in main() 36 x_array[i] = (1. / (nx * 2)) * (CeedScalar)(i % (nx * 2 + 1)); in main() 37 x_array[i + num_dofs] = (1. / (ny * 2)) * (CeedScalar)(i / (nx * 2 + 1)); in main() 47 col = i % nx; in main() 48 row = i / nx; in main() 49 offset = col * 2 + row * (nx * 2 + 1) * 2; in main()
|
| H A D | t581-operator.c | 21 CeedInt nx = 3, ny = 2; in main() local 23 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q; in main() 36 x_array[i] = (1. / (nx * 2)) * (CeedScalar)(i % (nx * 2 + 1)); in main() 37 x_array[i + num_dofs] = (1. / (ny * 2)) * (CeedScalar)(i / (nx * 2 + 1)); in main() 47 col = i % nx; in main() 48 row = i / nx; in main() 49 offset = col * 2 + row * (nx * 2 + 1) * 2; in main()
|
| H A D | t500-operator-f.f90 | 22 integer nx,nu 23 parameter(nx=nelem+1) 27 real*8 arrx(nx) 39 do i=0,nx-1 40 arrx(i+1)=i/(nx-1.d0) 47 call ceedelemrestrictioncreate(ceed,nelem,2,1,1,nx,ceed_mem_host,& 84 call ceedvectorcreate(ceed,nx,x,err)
|
| H A D | t503-operator-f.f90 | 22 integer nx,nu 23 parameter(nx=nelem+1) 27 real*8 arrx(nx) 45 do i=0,nx-1 46 arrx(i+1)=i/(nx-1.d0) 48 call ceedvectorcreate(ceed,nx,x,err) 58 call ceedelemrestrictioncreate(ceed,nelem,2,1,1,nx,ceed_mem_host,&
|
| H A D | t502-operator-f.f90 | 22 integer nx,nu 23 parameter(nx=nelem+1) 27 real*8 arrx(nx) 40 do i=0,nx-1 41 arrx(i+1)=i/(nx-1.d0) 48 call ceedelemrestrictioncreate(ceed,nelem,2,1,1,nx,ceed_mem_host,& 85 call ceedvectorcreate(ceed,nx,x,err)
|
| H A D | t501-operator-f.f90 | 22 integer nx,nu 23 parameter(nx=nelem+1) 27 real*8 arrx(nx) 40 do i=0,nx-1 41 arrx(i+1)=i/(nx-1.d0) 48 call ceedelemrestrictioncreate(ceed,nelem,2,1,1,nx,ceed_mem_host,& 85 call ceedvectorcreate(ceed,nx,x,err)
|
| H A D | t532-operator-f.f90 | 25 integer ndofs,nqpts,nx,ny 26 parameter(nx=3) 28 parameter(ndofs=(nx*2+1)*(ny*2+1)) 44 do i=0,nx*2 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) 60 col=mod(i,nx) 61 row=i/nx 62 offset=col*(p-1)+row*(nx*2+1)*(p-1) 65 indx(p*(p*i+k)+j+1)=offset+k*(nx*2+1)+j
|
| H A D | t532-operator.c | 19 CeedInt nx = 3, ny = 2; in main() local 20 CeedInt num_dofs = (nx * 2 + 1) * (ny * 2 + 1), num_qpts = num_elem * q * q; in main() 30 for (CeedInt i = 0; i < nx * 2 + 1; i++) { in main() 32 x_array[i + j * (nx * 2 + 1) + 0 * num_dofs] = (CeedScalar)i / (2 * nx); in main() 33 x_array[i + j * (nx * 2 + 1) + 1 * num_dofs] = (CeedScalar)j / (2 * ny); in main() 46 col = i % nx; in main() 47 row = i / nx; in main() 48 offset = col * (p - 1) + row * (nx * 2 + 1) * (p - 1); in main() 50 for (CeedInt k = 0; k < p; k++) ind_x[p * (p * i + k) + j] = offset + k * (nx * 2 + 1) + j; in main()
|
| H A D | t505-operator-f.f90 | 22 integer nx,nu 23 parameter(nx=nelem+1) 27 real*8 arrx(nx) 40 do i=0,nx-1 41 arrx(i+1)=i/(nx-1.d0) 48 call ceedelemrestrictioncreate(ceed,nelem,2,1,1,nx,ceed_mem_host,& 85 call ceedvectorcreate(ceed,nx,x,err)
|
| H A D | t524-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex 33 parameter(nx=3) 38 parameter(ndofs=(nx*2+1)*(ny*2+1)) 64 do j=0,nx*2 66 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx) 80 col=mod(i,nx) 81 row=i/nx 82 offset=col*2+row*(nx*2+1)*2 154 col=mod(i,nx) 155 row=i/nx
|
| H A D | t521-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex 33 parameter(nx=3) 38 parameter(ndofs=(nx*2+1)*(ny*2+1)) 64 do j=0,nx*2 66 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx) 80 col=mod(i,nx) 81 row=i/nx 82 offset=col*2+row*(nx*2+1)*2 154 col=mod(i,nx) 155 row=i/nx
|
| H A D | t520-operator-f.f90 | 32 integer ndofs,nqptstet,nqptshex,nqpts,nx,ny,nxtet,nytet,nxhex 33 parameter(nx=3) 38 parameter(ndofs=(nx*2+1)*(ny*2+1)) 63 do j=0,nx*2 65 arrx(i+j*(ny*2+1)+1*ndofs+1)=1.d0*j/(2*nx) 79 col=mod(i,nx) 80 row=i/nx 81 offset=col*2+row*(nx*2+1)*2 152 col=mod(i,nx) 153 row=i/nx
|
| /libCEED/python/tests/ |
| H A D | test-5-operator.py | 47 nx = nelem + 1 51 x = ceed.Vector(nx) 52 x_array = np.zeros(nx) 53 for i in range(nx): 54 x_array[i] = i / (nx - 1.0) 62 indx = np.zeros(nx * 2, dtype="int32") 63 for i in range(nx): 66 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx, 136 nx = nelem + 1 140 x = ceed.Vector(nx) [all …]
|