Ignore:
Timestamp:
Jan 8, 2013, 4:38:44 PM (12 years ago)
Author:
Ehouarn Millour
Message:

Added arch files for ADA (IDRIS IBMx3750) and made the following code modifications:
phylmd/printflag.F : removed "print" of unset variable (radpas0)
dyn3dmem/integrd_loc.F : removed unecessary "include mpif.h"
dyn3dmem/leapfrog_loc.F : removed unecessary "include mpif.h" and allocate saved variables at first call
dyn3dmem/mod_filtreg_p.F : added matmul() alternatives to call to BLAS routine SGEMM (which was incorectly set as DGEMM; which would fail if running with -r4)
filtrez/filtreg.F: changed calls to DGEMM into calls to SGEMM, so that code works with either -r4 or -r8 (the later being used in conjunction with "BLAS SGEMV=DGEMV SGEMM=DGEMM" preprocessing statements)
EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3dmem/mod_filtreg_p.F

    r1632 r1705  
    210210               IF( ifiltre.EQ.-2 )   THEN
    211211                  DO j = jdfil,jffil
    212                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     212#ifdef BLAS
     213                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    213214     &                    matrinvn(1,1,j), iim,
    214215     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    215216     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     217#else
     218                     champ_fft(:,j,:)=
     219     &                    matmul(matrinvn(:,:,j),champ_loc(:iim,j,:))
     220#endif
    216221                  ENDDO
    217222                 
    218223               ELSE IF ( griscal )     THEN
    219224                  DO j = jdfil,jffil
    220                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     225#ifdef BLAS
     226                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    221227     &                    matriceun(1,1,j), iim,
    222228     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    223229     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     230#else
     231                     champ_fft(:,j,:)=
     232     &                    matmul(matriceun(:,:,j),champ_loc(:iim,j,:))
     233#endif
    224234                  ENDDO
    225235                 
    226236               ELSE
    227237                  DO j = jdfil,jffil
    228                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     238#ifdef BLAS
     239                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    229240     &                    matricevn(1,1,j), iim,
    230241     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    231242     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     243#else
     244                     champ_fft(:,j,:)=
     245     &                    matmul(matricevn(:,:,j),champ_loc(:iim,j,:))
     246#endif
    232247                  ENDDO
    233248                 
    … …  
    238253               IF( ifiltre.EQ.-2 )   THEN
    239254                  DO j = jdfil,jffil
    240                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     255#ifdef BLAS
     256                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    241257     &                    matrinvs(1,1,j-jfiltsu+1), iim,
    242258     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    243259     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     260#else
     261                     champ_fft(:,j,:)=
     262     &                    matmul(matrinvs(:,:,j-jfiltsu+1),
     263     &                           champ_loc(:iim,j,:))
     264#endif
    244265                  ENDDO
    245266                 
    … …  
    247268                 
    248269                  DO j = jdfil,jffil
    249                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     270#ifdef BLAS
     271                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    250272     &                    matriceus(1,1,j-jfiltsu+1), iim,
    251273     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    252274     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     275#else
     276                     champ_fft(:,j,:)=
     277     &                    matmul(matriceus(:,:,j-jfiltsu+1),
     278     &                           champ_loc(:iim,j,:))
     279#endif
    253280                  ENDDO
    254281                 
    … …  
    256283                 
    257284                  DO j = jdfil,jffil
    258                      CALL DGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
     285#ifdef BLAS
     286                     CALL SGEMM("N", "N", iim, nbniv_loc, iim, 1.0,
    259287     &                    matricevs(1,1,j-jfiltsv+1), iim,
    260288     &                    champ_loc(1,j,1), iip1*(jje-jjb+1), 0.0,
    261289     &                    champ_fft(1,j,1), iip1*(jje-jjb+1))
     290#else
     291                     champ_fft(:,j,:)=
     292     &                    matmul(matricevs(:,:,j-jfiltsv+1),
     293     &                           champ_loc(:iim,j,:))
     294#endif
    262295                  ENDDO
    263296                 
    … …  
    269302               
    270303c     !-------------------------------------!
    271 c     ! Dés-agregation des niveau verticaux !
     304c     ! Dés-agregation des niveau verticaux !
    272305c     ! uniquement necessaire pour une      !
    273306c     ! execution OpenMP                    !
    … …  
    402435      END SUBROUTINE filtreg_p
    403436      END MODULE mod_filtreg_p
     437
Note: See TracChangeset for help on using the changeset viewer.