Skip to content

Commit

Permalink
Cleanups.
Browse files Browse the repository at this point in the history
  • Loading branch information
semi-h committed Feb 22, 2024
1 parent f5b24eb commit adbc5ab
Showing 1 changed file with 2 additions and 9 deletions.
11 changes: 2 additions & 9 deletions src/poisson_fft.f90
Original file line number Diff line number Diff line change
Expand Up @@ -55,8 +55,6 @@ subroutine base_init(self, xdirps, ydirps, zdirps, sz)
class(dirps_t), intent(in) :: xdirps, ydirps, zdirps
integer, intent(in) :: sz

integer :: nx, ny, nz

self%nx = xdirps%n; self%ny = ydirps%n; self%nz = zdirps%n

allocate (self%ax(self%nx), self%bx(self%nx))
Expand Down Expand Up @@ -84,7 +82,7 @@ subroutine waves_set(self, xdirps, ydirps, zdirps, sz)
real(dp) :: w, wp, rlexs, rleys, rlezs, xtt, ytt, ztt, xt1, yt1, zt1
complex(dp) :: xt2, yt2, zt2, xyzk

integer :: i, j, k, ka, kb, ix, iy, iz
integer :: i, j, ka, kb, ix, iy, iz


nx = xdirps%n; ny = ydirps%n; nz = zdirps%n
Expand Down Expand Up @@ -156,12 +154,7 @@ subroutine waves_set(self, xdirps, ydirps, zdirps, sz)
zk2(i) = cmplx(1._dp, 1._dp, kind=dp)*(nz*wp/zdirps%L)**2
end do

!do k = !sp%xst(3), sp%xen(3)
!do j = sp%xst(2), sp%xen(2)
! kxyz array is in z pencil orientation, with nz=nz/2+1 !!nooo
! kxyz array is in x pencil orientation.
!do k = 1, nx*ny/SZ
print*, 'set kxyz array'
print*, 'set waves array'
! TODO: do loop ranges below are valid only for single rank runs
do ka = 1, nz/2 + 1
do kb = 1, ny/sz
Expand Down

0 comments on commit adbc5ab

Please sign in to comment.