Changeset 1372

Timestamp:

Apr 24, 2014 6:29:32 AM (11 years ago)

Author:

raasch

Message:

openMP-bugfix for fftw: some arrays defined as threadprivate

File:

• palm/trunk/SOURCE/fft_xy.f90 (modified) (3 diffs)

palm/trunk/SOURCE/fft_xy.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! openMP-bugfix for fftw: some arrays defined as threadprivate
 ! 
 ! Former revisions:
@@ -168 +168 @@
     INTEGER(KIND=C_INT) ::  ny_c  !:
     
+    !$OMP THREADPRIVATE( x_out, y_out, x_in, y_in )
     COMPLEX(KIND=C_DOUBLE_COMPLEX), DIMENSION(:), ALLOCATABLE, SAVE ::         &
        x_out  !:
@@ -320 +321 @@
           nx_c = nx+1
           ny_c = ny+1
+          !$OMP PARALLEL
           ALLOCATE( x_in(0:nx+2), y_in(0:ny+2), x_out(0:(nx+1)/2),             &
                     y_out(0:(ny+1)/2) )
+          !$OMP END PARALLEL
           plan_xf = FFTW_PLAN_DFT_R2C_1D( nx_c, x_in, x_out, FFTW_ESTIMATE )
           plan_xi = FFTW_PLAN_DFT_C2R_1D( nx_c, x_out, x_in, FFTW_ESTIMATE )