57 real(kind=
rp),
allocatable :: xh_to_yh(:,:), xh_to_yht(:,:)
59 real(kind=
rp),
allocatable :: yh_to_xh(:,:), yh_to_xht(:,:)
61 type(c_ptr) :: xh_yh_d = c_null_ptr
63 type(c_ptr) :: xh_yht_d = c_null_ptr
65 type(c_ptr) :: yh_xh_d = c_null_ptr
67 type(c_ptr) :: yh_xht_d = c_null_ptr
88 type(
space_t),
intent(inout),
target :: Xh
89 type(
space_t),
intent(inout),
target :: Yh
90 integer :: deg_derivate
94 allocate(this%Xh_to_Yh(yh%lx,xh%lx))
95 allocate(this%Xh_to_YhT(xh%lx,yh%lx))
96 allocate(this%Yh_to_Xh(xh%lx,yh%lx))
97 allocate(this%Yh_to_XhT(yh%lx,xh%lx))
99 if (xh%t .eq.
gll .and. yh%t .eq.
gll)
then
100 else if ((xh%t .eq.
gl .and. yh%t .eq.
gll) .or. &
101 (yh%t .eq.
gl .and. xh%t .eq.
gll))
then
103 call neko_error(
'Unsupported interpolation')
106 call setup_intp(this%Xh_to_Yh, this%Xh_to_YhT, &
107 yh%zg, xh%zg, yh%lx, xh%lx, deg_derivate)
108 call setup_intp(this%Yh_to_Xh, this%Yh_to_XhT, &
109 xh%zg, yh%zg, xh%lx, yh%lx, deg_derivate)
114 call device_map(this%Xh_to_Yh, this%Xh_Yh_d, yh%lx*xh%lx)
115 call device_map(this%Xh_to_YhT, this%Xh_YhT_d, yh%lx*xh%lx)
116 call device_map(this%Yh_to_Xh, this%Yh_Xh_d, yh%lx*xh%lx)
117 call device_map(this%Yh_to_XhT, this%Yh_XhT_d, yh%lx*xh%lx)
118 call device_memcpy(this%Xh_to_Yh, this%Xh_Yh_d, yh%lx*xh%lx, &
120 call device_memcpy(this%Xh_to_YhT, this%Xh_YhT_d, yh%lx*xh%lx, &
122 call device_memcpy(this%Yh_to_Xh, this%Yh_Xh_d, yh%lx*xh%lx, &
124 call device_memcpy(this%Yh_to_XhT, this%Yh_XhT_d, yh%lx*xh%lx, &
133 if (
allocated(this%Xh_to_Yh))
then
134 deallocate(this%Xh_to_Yh)
136 if (
allocated(this%Xh_to_YhT))
then
137 deallocate(this%Xh_to_YhT)
139 if (
allocated(this%Yh_to_Xh))
then
140 deallocate(this%Yh_to_Xh)
142 if (
allocated(this%Yh_to_XhT))
then
143 deallocate(this%Yh_to_XhT)
145 if (c_associated(this%Yh_Xh_d))
then
148 if (c_associated(this%Yh_XhT_d))
then
151 if (c_associated(this%Xh_Yh_d))
then
154 if (c_associated(this%Xh_YhT_d))
then
169 real(kind=
rp),
intent(inout) :: x(this%Xh%lx, this%Xh%lx, this%Xh%lx, nel)
170 real(kind=
rp),
intent(inout) :: y(this%Yh%lx, this%Yh%lx, this%Yh%lx, nel)
171 if (to_space .eq. this%Yh)
then
172 call tnsr3d(y, this%Yh%lx, x, &
173 this%Xh%lx,this%Yh_to_XhT, &
174 this%Yh_to_Xh, this%Yh_to_Xh, nel)
175 else if (to_space .eq. this%Xh)
then
176 call tnsr3d(y, this%Xh%lx, x, &
177 this%Yh%lx,this%Yh_to_Xh, &
178 this%Yh_to_XhT, this%Yh_to_XhT, nel)
180 call neko_error(
'Invalid interpolation')
195 real(kind=
rp),
intent(inout) :: x(this%Xh%lx, this%Xh%lx, this%Xh%lx, nel)
196 real(kind=
rp),
intent(inout) :: y(this%Yh%lx, this%Yh%lx, this%Yh%lx, nel)
197 if (to_space .eq. this%Yh)
then
199 this%Xh%lx,this%Yh_to_XhT, &
200 this%Yh_to_Xh, this%Yh_to_Xh, nel)
201 else if (to_space .eq. this%Xh)
then
203 this%Yh%lx,this%Yh_to_Xh, &
204 this%Yh_to_XhT, this%Yh_to_XhT, nel)
206 call neko_error(
'Invalid interpolation')
Map a Fortran array to a device (allocate and associate)
Copy data between host and device (or device and device)
Device abstraction, common interface for various accelerators.
integer, parameter, public host_to_device
subroutine, public device_free(x_d)
Deallocate memory on the device.
Fast diagonalization methods from NEKTON.
subroutine, public setup_intp(jh, jht, z_to, z_from, n_to, n_from, derivative)
Compute interpolation weights for points z_to using values at points z_from.
Routines to interpolate between different spaces.
subroutine interpolator_map_host(this, y, x, nel, to_space)
Interpolates an array to one of Xh or Yh on host.
subroutine interpolator_map(this, y, x, nel, to_space)
Interpolates an array to one of Xh or Yh.
subroutine interpolator_init(this, Xh, Yh)
Constructor to initialize with two different spaces.
subroutine interpolator_free(this)
integer, parameter neko_bcknd_device
integer, parameter, public rp
Global precision used in computations.
Defines a function space.
integer, parameter, public gll
integer, parameter, public gl
subroutine, public tnsr3d_cpu(v, nv, u, nu, A, Bt, Ct, nelv)
subroutine, public tnsr3d(v, nv, u, nu, A, Bt, Ct, nelv)
Tensor product performed on nelv elements.
Interpolation between two space::space_t.
The function space for the SEM solution fields.