Merge pull request #71 from ORNL-Fusion/asym_index

Remove redundant code in fourier.
ORNL-Fusion · Oct 31, 2024 · 4fbc02b · 4fbc02b
2 parents d059893 + c2696b2
commit 4fbc02b
Showing 1 changed file with 2 additions and 14 deletions.
diff --git 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)