[BACK]Return to idamax.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / blas

Diff for /rpl/lapack/blas/idamax.f between versions 1.6 and 1.7

version 1.6, 2010/12/21 13:51:26 version 1.7, 2011/07/22 07:38:02
Line 31 Line 31
       IF (N.LT.1 .OR. INCX.LE.0) RETURN        IF (N.LT.1 .OR. INCX.LE.0) RETURN
       IDAMAX = 1        IDAMAX = 1
       IF (N.EQ.1) RETURN        IF (N.EQ.1) RETURN
       IF (INCX.EQ.1) GO TO 20        IF (INCX.EQ.1) THEN
 *  *
 *        code for increment not equal to 1  *        code for increment equal to 1
 *  *
       IX = 1           DMAX = DABS(DX(1))
       DMAX = DABS(DX(1))           DO I = 2,N
       IX = IX + INCX              IF (DABS(DX(I)).GT.DMAX) THEN
       DO 10 I = 2,N                 IDAMAX = I
           IF (DABS(DX(IX)).LE.DMAX) GO TO 5                 DMAX = DABS(DX(I))
           IDAMAX = I              END IF
           DMAX = DABS(DX(IX))           END DO
     5     IX = IX + INCX        ELSE
    10 CONTINUE  
       RETURN  
 *  *
 *        code for increment equal to 1  *        code for increment not equal to 1
 *  *
    20 DMAX = DABS(DX(1))           IX = 1
       DO 30 I = 2,N           DMAX = DABS(DX(1))
           IF (DABS(DX(I)).LE.DMAX) GO TO 30           IX = IX + INCX
           IDAMAX = I           DO I = 2,N
           DMAX = DABS(DX(I))              IF (DABS(DX(IX)).GT.DMAX) THEN
    30 CONTINUE                 IDAMAX = I
                  DMAX = DABS(DX(IX))
               END IF
               IX = IX + INCX
            END DO
         END IF
       RETURN        RETURN
       END        END
Removed from v.1.6  
changed lines
  Added in v.1.7

CVSweb interface <joel.bertrand@systella.fr>