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

Diff for /rpl/lapack/lapack/dlaebz.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:28 version 1.8, 2011/07/22 07:38:06
Line 2 Line 2
      $                   RELTOL, PIVMIN, D, E, E2, NVAL, AB, C, MOUT,       $                   RELTOL, PIVMIN, D, E, E2, NVAL, AB, C, MOUT,
      $                   NAB, WORK, IWORK, INFO )       $                   NAB, WORK, IWORK, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary routine (version 3.3.1) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     November 2006  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            IJOB, INFO, MINP, MMAX, MOUT, N, NBMIN, NITMAX        INTEGER            IJOB, INFO, MINP, MMAX, MOUT, N, NBMIN, NITMAX
Line 251 Line 251
 *        Compute the number of eigenvalues in the initial intervals.  *        Compute the number of eigenvalues in the initial intervals.
 *  *
          MOUT = 0           MOUT = 0
 *DIR$ NOVECTOR  
          DO 30 JI = 1, MINP           DO 30 JI = 1, MINP
             DO 20 JP = 1, 2              DO 20 JP = 1, 2
                TMP1 = D( 1 ) - AB( JI, JP )                 TMP1 = D( 1 ) - AB( JI, JP )
Line 407 Line 406
                   TMP2 = MIN( TMP2, -PIVMIN )                    TMP2 = MIN( TMP2, -PIVMIN )
                END IF                 END IF
 *  *
 *              A series of compiler directives to defeat vectorization  
 *              for the next loop  
 *  
 *$PL$ CMCHAR=' '  
 CDIR$          NEXTSCALAR  
 C$DIR          SCALAR  
 CDIR$          NEXT SCALAR  
 CVD$L          NOVECTOR  
 CDEC$          NOVECTOR  
 CVD$           NOVECTOR  
 *VDIR          NOVECTOR  
 *VOCL          LOOP,SCALAR  
 CIBM           PREFER SCALAR  
 *$PL$ CMCHAR='*'  
 *  
                DO 90 J = 2, N                 DO 90 J = 2, N
                   TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1                    TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1
                   IF( TMP2.LE.PIVMIN ) THEN                    IF( TMP2.LE.PIVMIN ) THEN
Line 487  CIBM           PREFER SCALAR Line 471  CIBM           PREFER SCALAR
   100       CONTINUE    100       CONTINUE
             KL = KLNEW              KL = KLNEW
 *  *
 *           End of Serial Version of the loop  
 *  
          END IF           END IF
 *  *
 *        Check for convergence  *        Check for convergence
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>