Return to isamax.f CVS log

Up to [Research Unix] / researchv10no / cmd / view2d / num

Annotation of researchv10no/cmd/view2d/num/isamax.f, revision 1.1

1.1     ! root        1:       integer function isamax(n,sx,incx)
        !             2: c
        !             3: c     finds the index of element having max. absolute value.
        !             4: c     jack dongarra, linpack, 3/11/78.
        !             5: c
        !             6:       real sx(1),smax
        !             7:       integer i,incx,ix,n
        !             8: c
        !             9:       isamax = 0
        !            10:       if( n .lt. 1 ) return
        !            11:       isamax = 1
        !            12:       if(n.eq.1)return
        !            13:       if(incx.eq.1)go to 20
        !            14: c
        !            15: c        code for increment not equal to 1
        !            16: c
        !            17:       ix = 1
        !            18:       smax = abs(sx(1))
        !            19:       ix = ix + incx
        !            20:       do 10 i = 2,n
        !            21:          if(abs(sx(ix)).le.smax) go to 5
        !            22:          isamax = i
        !            23:          smax = abs(sx(ix))
        !            24:     5    ix = ix + incx
        !            25:    10 continue
        !            26:       return
        !            27: c
        !            28: c        code for increment equal to 1
        !            29: c
        !            30:    20 smax = abs(sx(1))
        !            31:       do 30 i = 2,n
        !            32:          if(abs(sx(i)).le.smax) go to 30
        !            33:          isamax = i
        !            34:          smax = abs(sx(i))
        !            35:    30 continue
        !            36:       return
        !            37:       end

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.