Changeset 1372 for palm


Ignore:
Timestamp:
Apr 24, 2014 6:29:32 AM (8 years ago)
Author:
raasch
Message:

openMP-bugfix for fftw: some arrays defined as threadprivate

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/fft_xy.f90

    r1354 r1372  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! openMP-bugfix for fftw: some arrays defined as threadprivate
    2323!
    2424! Former revisions:
     
    168168    INTEGER(KIND=C_INT) ::  ny_c  !:
    169169   
     170    !$OMP THREADPRIVATE( x_out, y_out, x_in, y_in )
    170171    COMPLEX(KIND=C_DOUBLE_COMPLEX), DIMENSION(:), ALLOCATABLE, SAVE ::         &
    171172       x_out  !:
     
    320321          nx_c = nx+1
    321322          ny_c = ny+1
     323          !$OMP PARALLEL
    322324          ALLOCATE( x_in(0:nx+2), y_in(0:ny+2), x_out(0:(nx+1)/2),             &
    323325                    y_out(0:(ny+1)/2) )
     326          !$OMP END PARALLEL
    324327          plan_xf = FFTW_PLAN_DFT_R2C_1D( nx_c, x_in, x_out, FFTW_ESTIMATE )
    325328          plan_xi = FFTW_PLAN_DFT_C2R_1D( nx_c, x_out, x_in, FFTW_ESTIMATE )
Note: See TracChangeset for help on using the changeset viewer.