From c2696b24e526021df51a83b74afe27fbc40cc477 Mon Sep 17 00:00:00 2001 From: cianciosa Date: Thu, 31 Oct 2024 16:39:30 -0400 Subject: [PATCH] Remove redundant code in fourier. --- Sources/fourier.f90 | 16 ++-------------- 1 file changed, 2 insertions(+), 14 deletions(-) diff --git a/Sources/fourier.f90 b/Sources/fourier.f90 index ace2fc9..0c17b63 100644 --- a/Sources/fourier.f90 +++ b/Sources/fourier.f90 @@ -751,13 +751,7 @@ PURE SUBROUTINE fourier_tomnsp_2d_v(this, xmn, parity) xmn(:,:) = 0 IF (parity .eq. f_cos) THEN DO j = 1, SIZE(this%cosnv, 1) ! nzeta - workcos = this%workmj1(m0,j)*this%cosnv(j,n0) - worksin = this%workmj2(m0,j)*this%sinnv(j,n0) - - xmn(m0 + moff,n0 + noff) = xmn(m0 + moff,n0 + noff) & - + workcos - worksin - - DO m = m1, UBOUND(this%cosmu, 2) ! mpol + DO m = m0, UBOUND(this%cosmu, 2) ! mpol workcos = this%workmj1(m,j)*this%cosnv(j,n0) worksin = this%workmj2(m,j)*this%sinnv(j,n0) @@ -788,13 +782,7 @@ PURE SUBROUTINE fourier_tomnsp_2d_v(this, xmn, parity) xmn(m0 + moff,:-n1 + noff) = zero ELSE DO j = 1, SIZE(this%cosnv, 1) ! nzeta - workcos = this%workmj2(m0,j)*this%cosnv(j,n0) - worksin = this%workmj1(m0,j)*this%sinnv(j,n0) - - xmn(m0 + moff,n0 + noff) = xmn(m0 + moff,n0 + noff) & - + workcos + worksin - - DO m = m1, UBOUND(this%cosmu, 2) ! mpol + DO m = m0, UBOUND(this%cosmu, 2) ! mpol workcos = this%workmj2(m,j)*this%cosnv(j,n0) worksin = this%workmj1(m,j)*this%sinnv(j,n0)