Ignore:
Timestamp:
Oct 29, 2009, 2:55:23 PM (15 years ago)
Author:
yann meurdesoif
Message:

Optimisations SX9

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F

    r774 r1250  
    1313      real sx((n-1)*incx+1),sy((n-1)*incy+1)
    1414c
     15      if (incx.eq.1.and.incy.eq.1) then
     16      do 10 i=1,n
     17         sy(i)=sx(i)
     1810    continue
     19      else
    1520      iy=1
    1621      ix=1
    17       do 10 i=1,n
     22      do 11 i=1,n
    1823         sy(iy)=sx(ix)
    1924         ix=ix+incx
    2025         iy=iy+incy
    21 10    continue
     2611    continue
     27      endif
    2228c
    2329      return
    … …  
    3238c
    3339      ssum=0.
     40      if (incx.eq.1) then
     41      do 10 i=1,n
     42         ssum=ssum+sx(i)
     4310    continue
     44      else
    3445      ix=1
    35       do 10 i=1,n
     46      do 11 i=1,n
    3647         ssum=ssum+sx(ix)
    3748         ix=ix+incx
    38 10    continue
     4911    continue
     50      endif
    3951c
    4052      return
Note: See TracChangeset for help on using the changeset viewer.