Changeset 72

Timestamp:

06/21/16 11:42:45 (8 years ago)

Author:

dumas

Message:

OpenMP parallelization in diagno-L2_mod.f90, eq_ellipt_sgbsv_mod-0.2.f90 and flottab2-0.7.f90

Location:

trunk/SOURCES

Files:

• Fichiers-parametres/Makefile.tof-lsce3130.inc (modified) (1 diff)
• New-remplimat/diagno-L2_mod.f90 (modified) (4 diffs)
• New-remplimat/eq_ellipt_sgbsv_mod-0.2.f90 (modified) (9 diffs)
• flottab2-0.7.f90 (modified) (23 diffs)
• velocities-polyn-0.3.f90 (modified) (5 diffs)

trunk/SOURCES/Fichiers-parametres/Makefile.tof-lsce3130.inc

@@ -39 +39 @@
         IFORT= ifort
 
-        ARITHMi    = -O2 -fp-model precise  -heap-arrays -traceback -mcmodel=medium #  -warn -traceback   -CB  -g  # options pour une meilleure arithmetique
+        ARITHMi    = -O2 -fp-model precise  -heap-arrays -traceback -mcmodel=medium  #-openmp #  -warn -traceback   -CB  -g  # options pour une meilleure arithmetique
                                                                       # (normalement reproductible)
     ifeq ($(debug), 1)

trunk/SOURCES/New-remplimat/diagno-L2_mod.f90

@@ -110 +110 @@
   call calc_pvi                    ! calcule les viscosites integrees
 
+!$OMP PARALLEL
+!$OMP WORKSHARE
    where (flot(:,:).and.(H.gt.1))  ! valeur analytique pour les shelfs
       pvi(:,:) = (4./coef_Sflot/rog)**2/btt(:,:,1,1)/H(:,:)
    end where
-
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
   ! avec couplage thermomecanique
 !  write(166,*) ' apres call calc_pvi',itour_pvi
@@ -178 +181 @@
 
   !      nouvelles vitesses         
-
+  !$OMP PARALLEL
+  !$OMP WORKSHARE
   uxbar(:,:)=uxb1(:,:)
   uybar(:,:)=uyb1(:,:)
-
+  !$OMP END WORKSHARE
 
 
   !    calcul de tobmx et tobmy (frottement basal) apres calcul des vitesses
   !    ---------------------------------------------------------------------
+  !$OMP DO
   do j=1,ny
      do i=1,nx
@@ -192 +197 @@
      enddo
   enddo
+  !$OMP END DO
 
   ! Mise ne réserve des vitesses flgzmx et flgzmy
+  !$OMP WORKSHARE
   where (flgzmx(:,:))
      uxflgz(:,:)=uxbar(:,:)
@@ -205 +212 @@
      uyflgz(:,:)=0.
   endwhere
-
+  !$OMP END WORKSHARE
+  !$OMP END PARALLEL
   !i=92
   !j=152

trunk/SOURCES/New-remplimat/eq_ellipt_sgbsv_mod-0.2.f90

@@ -8 +8 @@
 
 ! declaration des variables
-
+!$ USE OMP_LIB
 use remplimat_declar   ! toutes les routines de ce module connaissent les Tuij....
 
@@ -154 +154 @@
       !   write(6,*) 'sortie ok de sgbsv'
       !   write(6,*)'-----------------------'
-
+      !$OMP PARALLEL
+      !$OMP DO
       do j = ny1,ny2
          do i = nx1,nx2
@@ -171 +172 @@
          end do
       end do
-
-
+      !$OMP END DO
+      !$OMP END PARALLEL
 
    else
@@ -181 +182 @@
 
 else                 ! nblig=0, pas de passage par sgbsv
+   !$OMP PARALLEL
+   !$OMP WORKSHARE
    uxnew(nx1:nx2,ny1:ny2)=uxprec(nx1:nx2,ny1:ny2)
    uynew(nx1:nx2,ny1:ny2)=uyprec(nx1:nx2,ny1:ny2)
+   !$OMP END WORKSHARE
+   !$OMP END PARALLEL
 end if
 
@@ -209 +214 @@
 ! numerotation des lignes 
 !--------------------------
-
+!$OMP PARALLEL
+!$OMP WORKSHARE
 ligu_L2(:,:) = -9999
 ligv_L2(:,:) = -9999
 pos_ligu(:,:) = -9999
 pos_ligv(:,:) = -9999
-
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 ! Boucle sur les noeuds avec la dimension la plus grande a l'interieur
 
@@ -248 +255 @@
 end do count
 
+
 nblig = count_line-1
 
@@ -253 +261 @@
 ! calcul de la largeur de bande. depend de la numerotation des lignes
 !--------------------------------------------------------------------
+!$OMP PARALLEL
+!$OMP WORKSHARE
 lu_band(:) = 0
 ll_band(:) = 0
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 
 larg_band: do j = ny1,ny2
@@ -367 +379 @@
 ! ecriture de la matrice proprement dite et du vecteur bdr
 !------------------------------------------------------------
+!$OMP PARALLEL
+!$OMP WORKSHARE
 mmat(:,:) = 0.
 bdr(:,:) = 0.
-
-
+!$OMP END WORKSHARE
+
+!$OMP DO PRIVATE(ligne,ilmin,ilmax,jlmin,jlmax,colonne)
 ij_loop : do j = ny1,ny2         ! attention, on rempli 2 lignes par 2 lignes (u,v)
    do i = nx1,nx2
@@ -481 +496 @@
  end do
 end do ij_loop
+!$OMP END DO
+!$OMP END PARALLEL
 
 

trunk/SOURCES/flottab2-0.7.f90

@@ -630 +630 @@
 
 !On compte comme englacé uniquement les calottes dont une partie est posée      
-!$OMP PARALLEL PRIVATE(smax_,smax_coord,smax_i,smax_j)
+!$OMP PARALLEL PRIVATE(smax_,smax_coord,smax_i,smax_j,mask_tache_ij)
 !$OMP DO
-do i=3,nx-2
-   do j=3,ny-2
+do j=3,ny-2
+   do i=3,nx-2
 test1:  if (.not.iceberg(table_out(i,j))) then ! on est pas sur un iceberg 
          if (nb_pts_tache(table_out(i,j)).ge.1) then
@@ -748 +748 @@
 subroutine determin_tache
 
-!!$ USE OMP_LIB
+!$ USE OMP_LIB
 
 implicit none
@@ -769 +769 @@
 enddo
 !      table_in  = .false.
-!!$OMP PARALLEL
-!!$OMP WORKSHARE
+!$OMP PARALLEL
+!$OMP WORKSHARE
 table_out(:,:) = 0
 iceberg(:)  = .true.
 icetrim (:) = .true.
 nb_pts_tache(:) = 0
-!!$OMP END WORKSHARE
-!!$OMP END PARALLEL
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 !    open(unit=100,file="tache.data",status='replace')
 
 ! 2-reperage des taches
 !----------------------
-!!$OMP PARALLEL PRIVATE(mask,label,indice)
-!!$OMP DO
-do i=2,nx-1
-  do j=2,ny-1
+
+do j=2,ny-1
+  do i=2,nx-1
 
 
@@ -852 +851 @@
   enddo
 enddo
-!!$OMP END DO
-!!$OMP END PARALLEL
+
 
 
@@ -869 +867 @@
 enddo
 
-!!$OMP PARALLEL
-!!$OMP DO REDUCTION(+:nb_pts_tache)
-do i=1,nx
-   do j=1,ny
+!$OMP PARALLEL
+!$OMP DO REDUCTION(+:nb_pts_tache)
+do j=1,ny
+   do i=1,nx
       if (table_out(i,j).ne.0) then
          table_out(i,j)=compt(table_out(i,j))  
@@ -879 +877 @@
    enddo
 enddo
-!!$OMP END DO
-!!$OMP END PARALLEL
+!$OMP END DO
+!$OMP END PARALLEL
 
 
@@ -915 +913 @@
 integer :: j_moins1,j_plus1,j_plus2
       
-      !!$OMP PARALLEL
-      !!$OMP DO
+      !$OMP PARALLEL
+      !$OMP DO
       do j=3,ny-2
         do i=3,nx-2
@@ -970 +968 @@
        end do
       end do
-      !!$OMP END DO
+      !$OMP END DO
+      !$OMP ENd PARALLEL
 
 !!$call detect_assym(nx,ny,0,41,1,0,1,0,H,itestf)
@@ -985 +984 @@
        
 baies: do k=1,2
-!!$OMP DO PRIVATE(i_moins1,j_moins1,i_plus1,j_plus1,i_plus2,j_plus2)
+!$OMP PARALLEL
+!$OMP DO PRIVATE(i_moins1,j_moins1,i_plus1,j_plus1,i_plus2,j_plus2)
 do j=1,ny
    do i=1,nx
@@ -1020 +1020 @@
    end do
 end do
-!!$OMP END DO
+!$OMP END DO
+!$OMP END PARALLEL
 end do baies
 
@@ -1032 +1033 @@
 !!$end if
 
-!!$OMP DO
+!$OMP PARALLEL
+!$OMP DO
 do j=2,ny-1
    do i=2,nx-1
@@ -1046 +1048 @@
    end do
 end do
-!!$OMP END DO
+!$OMP END DO
 
 ! traitement des bords. On considere que l'exterieur n'a pas de glace
 ! attention ce n'est vrai que sur la grande grille
 
-!!$OMP DO PRIVATE(i)
+!$OMP DO PRIVATE(i)
 do j=2,ny-1
    i=1
@@ -1058 +1060 @@
    front(i,j)=(ice(i-1,j)+ice(i,j+1)+ice(i,j-1))
 end do
-!!$OMP END DO 
-
-!!$OMP DO PRIVATE(j)
+!$OMP END DO 
+
+!$OMP DO PRIVATE(j)
 do i=2,nx-1
    j=1 
@@ -1067 +1069 @@
    front(i,j)=(ice(i-1,j)+ice(i+1,j)+ice(i,j-1))
 end do
-!!$OMP END DO 
+!$OMP END DO 
+
 ! traitement des coins
 
@@ -1089 +1092 @@
 !   si ces deux taches sont posées (ou une des deux), il n'y a pas assez de conditions aux limites
 
-!!$OMP DO 
+!$OMP DO 
 do j=1,ny
    do i=1,nx-1
@@ -1102 +1105 @@
    end do
 end do
-!!$OMP END DO
-
-!!$OMP DO 
+!$OMP END DO
+
+!$OMP DO 
 do j=1,ny-1
    do i=1,nx
@@ -1117 +1120 @@
    end do
 end do
-!!$OMP END DO
+!$OMP END DO
 
 !isolx signifie pas de voisins en x
@@ -1125 +1128 @@
 
 ! calcul de frontfacex et isolx
-!!$OMP DO
+!$OMP DO
 do j=1,ny
    do i=2,nx-1
@@ -1145 +1148 @@
    end do
 end do
-!!$OMP END DO
+!$OMP END DO
 
 ! calcul de frontfacey et isoly
-!!$OMP DO
+!$OMP DO
 do j=2,ny-1
    do i=1,nx
@@ -1168 +1171 @@
    end do
 end do
-!!$OMP END DO
+!$OMP END DO
 
 
@@ -1174 +1177 @@
 ! attention ce n'est vrai que sur la grande grille
 
-!!$OMP DO PRIVATE(i)
+!$OMP DO PRIVATE(i)
 do j=2,ny-1
    i=1
@@ -1193 +1196 @@
    end if
 end do
-!!$OMP END DO
-
-!!$OMP DO PRIVATE(j)
+!$OMP END DO
+
+!$OMP DO PRIVATE(j)
 do i=2,nx-1
    j=1 
@@ -1214 +1217 @@
    end if
 end do
-!!$OMP END DO
-!!OMP END PARALLEL
+!$OMP END DO
+!$OMP END PARALLEL
 
 return

trunk/SOURCES/velocities-polyn-0.3.f90

@@ -23 +23 @@
 !     Attention la vitesse verticale est "super reduite"
 ! ==========================================================================
-
+!$ USE OMP_LIB
 use module3d_phy
 !use deform_declar 
@@ -32 +32 @@
 if (itracebug.eq.1)  call tracebug(' Entree dans routine SIA_velocity')
 
-    
+!$OMP PARALLEL
+!$OMP DO COLLAPSE(2)    
 do k=1,nz
    do j=2,ny-1
@@ -101 +102 @@
    end do
 end do
-
+!$OMP END DO
 
 !       *************************** Z VELOCITIES ******************
 
-
+!$OMP DO
 do j=2,ny-1
    do i=2,nx-1
@@ -144 +145 @@
      end do
  end do
+!$OMP END DO
 
+!$OMP DO
   do j=2,ny-1
      do i=2,nx-1  
@@ -151 +154 @@
      end do
   end do
-
+!$OMP END DO
+!$OMP END PARALLEL
 
 end subroutine SIA_velocities