Changeset 683 for palm/trunk/SOURCE/poisfft.f90

Timestamp:

Feb 9, 2011 2:25:15 PM (13 years ago)

Author:

raasch

Message:

New:
---

optional exchange of ghost points in synchronous mode via MPI_SENDRCV,
steered by d3par parameter synchronous_exchange
(cpu_statistics, exchange_horiz, modules, parin)

openMP-parallelization of pressure solver (fft-method) for 2d-domain-decomposition
(poisfft, transpose)

Changed:

---

Errors:

---

mpt bugfix for netCDF4 usage (mrun)

File:

• palm/trunk/SOURCE/poisfft.f90 (modified) (12 diffs)

palm/trunk/SOURCE/poisfft.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! openMP parallelization for 2d-domain-decomposition
 !
 ! Former revisions:
@@ -287 +288 @@
 !
 !--    Define constant elements of the tridiagonal matrix.
+!$OMP  PARALLEL PRIVATE ( k, i )
+!$OMP  DO
        DO  k = 0, nz-1
           DO  i = nxl_z, nxr_z
@@ -293 +296 @@
           ENDDO
        ENDDO
+!$OMP  END PARALLEL
 
 #if defined( __parallel )
 !
 !--    Repeat for all y-levels.
+!$OMP  PARALLEL FIRSTPRIVATE( tri ) PRIVATE ( ar1, j )
+!$OMP  DO
        DO  j = nys_z, nyn_z
           IF ( j <= nnyh )  THEN
@@ -306 +312 @@
           CALL substi( ar, ar1, tri, j )
        ENDDO
+!$OMP  END PARALLEL
 #else
 !
@@ -527 +534 @@
 !
 !--    Performing the fft with one of the methods implemented
+!$OMP  PARALLEL PRIVATE ( j, k )
+!$OMP  DO
        DO  k = nzb_x, nzt_x
           DO  j = nys_x, nyn_x
@@ -532 +541 @@
           ENDDO
        ENDDO
+!$OMP  END PARALLEL
 
     END SUBROUTINE fftxp
@@ -550 +560 @@
 !
 !--    Performing the fft with one of the methods implemented
+!$OMP  PARALLEL PRIVATE ( j, k )
+!$OMP  DO
        DO  k = 1, nz
           DO  j = 0, ny
@@ -555 +567 @@
           ENDDO
        ENDDO
+!$OMP  END PARALLEL
 
     END SUBROUTINE fftx
@@ -575 +588 @@
 !
 !--    Performing the fft with one of the methods implemented
+!$OMP  PARALLEL PRIVATE ( i, k )
+!$OMP  DO
        DO  k = nzb_y, nzt_y
           DO  i = nxl_y, nxr_y
@@ -580 +595 @@
           ENDDO
        ENDDO
+!$OMP  END PARALLEL
 
     END SUBROUTINE fftyp
@@ -598 +614 @@
 !
 !--    Performing the fft with one of the methods implemented
+!$OMP  PARALLEL PRIVATE ( i, k )
+!$OMP  DO
        DO  k = 1, nz
           DO  i = 0, nx
@@ -603 +621 @@
           ENDDO
        ENDDO
+!$OMP  END PARALLEL
 
     END SUBROUTINE ffty