Changeset 1601 for LMDZ4/branches/LMDZ4_AR5/libf/filtrez/mod_fft_mkl.F90

Timestamp:

Dec 7, 2011, 8:05:43 AM (12 years ago)

Author:

Ehouarn Millour

Message:

Updates and upgrades of filter:

• add minor corrections (r1591 of trunk) in filtreg_mod.F90: some arrays were oversized and the computation of the indexes from which the filter is applied could go wrong in extreme cases.
• adapt filtreg_p.F (r1597 of trunk) so that BLAS routine DGEMM is only used if 'BLAS' preprocessing flag is set.
• update FFT filter routines to match current 'trunk' versions (mostly cosmetic changes, exept for the implementation of use of FFTW in mod_fft_fftw.F90).
• update "arch-PW6_VARGAS.fcm" to enable use of FFTW

NB: implemented FFTs assume working precision to be double precision ; trying to use them when working precision is single precision will clearly end in tragedy.
EM

File:

• LMDZ4/branches/LMDZ4_AR5/libf/filtrez/mod_fft_mkl.F90 (modified) (3 diffs)

LMDZ4/branches/LMDZ4_AR5/libf/filtrez/mod_fft_mkl.F90

@@ -24 +24 @@
     INTEGER :: nb_vect_max
     REAL    :: rtmp=1.
-    COMPLEX*16 :: ctmp
+    COMPLEX :: ctmp
     INTEGER :: itmp=1
     INTEGER :: isign=0
@@ -60 +60 @@
     INTEGER,INTENT(IN)  :: nb_vect
     REAL,INTENT(IN)     :: vect((vsize+inc)*nb_vect)
-    COMPLEX*16,INTENT(OUT) :: TF_vect((vsize/2+1)*nb_vect)
+    COMPLEX,INTENT(OUT) :: TF_vect((vsize/2+1)*nb_vect)
     REAL                :: work(4*vsize*nb_vect)
     INTEGER             :: ierr
@@ -102 +102 @@
     INTEGER,INTENT(IN)  :: nb_vect
     REAL,INTENT(OUT)    :: vect((vsize+inc)*nb_vect)
-    COMPLEX*16,INTENT(IN ) :: TF_vect((vsize/2+1)*nb_vect)
+    COMPLEX,INTENT(IN ) :: TF_vect((vsize/2+1)*nb_vect)
     REAL                :: work(4*vsize*nb_vect)
     INTEGER             :: ierr