Commit f156fed2 authored by Gregor Michalicek's avatar Gregor Michalicek

Introduce oUnit to files in fft directory

parent 4f5701b8
......@@ -56,7 +56,8 @@ CONTAINS
! the following two constants should agree with the array dimensions
#ifdef CPP_Singleton
SUBROUTINE cfft(a, b, ntot, n, nspan, isn)
use m_juDFT
USE m_juDFT
USE m_constants
IMPLICIT NONE
! .. Scalar Arguments ..
INTEGER :: isn, n, nspan, ntot
......@@ -75,9 +76,6 @@ CONTAINS
REAL, ALLOCATABLE :: at(:), bt(:), ck(:), sk(:)
INTEGER, ALLOCATABLE :: nfac(:), np(:)
! ..
! .. Intrinsic Functions ..
INTRINSIC cos, real, mod, sin, sqrt
! ..
! .. Equivalences ..
EQUIVALENCE(i, ii)
! ..
......@@ -591,7 +589,7 @@ CONTAINS
! error finish, insufficient array storage
590 CONTINUE
! isn = 0
WRITE (6, FMT=8000)
WRITE (oUnit, FMT=8000)
CALL juDFT_error('array bounds exceeded', calledby='cfft')
8000 FORMAT('array bounds exceeded within subroutine cft')
667 CONTINUE
......
......@@ -32,6 +32,7 @@ CONTAINS
! G.Bihlmayer (UniWien)
! **********************************************************************
USE m_constants
USE m_cfft
IMPLICIT NONE
......@@ -54,15 +55,15 @@ CONTAINS
IF ((isn/=-1) .AND. (isn /= 1)) CALL juDFT_error("choose isn=+/- 1" &
,calledby ="rfft")
IF ((n1d < n1) .OR. (n2d < n2) .OR. (n3d < n3)) THEN
WRITE (6,*) 'n1d,n2d,n3d =',n1d,n2d,n3d
WRITE (6,*) 'n1 ,n2 ,n3 =',n1 ,n2 ,n3
WRITE (oUnit,*) 'n1d,n2d,n3d =',n1d,n2d,n3d
WRITE (oUnit,*) 'n1 ,n2 ,n3 =',n1 ,n2 ,n3
CALL juDFT_error("n(i) > n(i)d",calledby ="rfft")
ENDIF
IF ((n1 <= 2*nw1+1) .OR. &
(n2 <= 2*nw2+1) .OR. &
(n3 <= 2*nw3+1)) THEN
! WRITE (6,*) 'n1 ,n2 ,n3 =',n1 ,n2 ,n3
! WRITE (6,*) 'nw1,nw2,nw3 =',nw1,nw2,nw3
! WRITE (oUnit,*) 'n1 ,n2 ,n3 =',n1 ,n2 ,n3
! WRITE (oUnit,*) 'nw1,nw2,nw3 =',nw1,nw2,nw3
l_nopad= .TRUE.
ELSE
l_nopad= .FALSE.
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment