[Scipy-svn] r4929 - trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK

[scipy-svn at scipy.org]

Author: cdavid
Date: 2008-11-02 19:25:53 -0600 (Sun, 02 Nov 2008)
New Revision: 4929

Modified:
   trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/clahqr.f
   trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/zlahqr.f
Log:
Use wrapped wcladiv/wzladiv instead of direct LAPACK calls.

Modified: trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/clahqr.f
===================================================================
--- trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/clahqr.f	2008-11-03 01:25:00 UTC (rev 4928)
+++ trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/clahqr.f	2008-11-03 01:25:53 UTC (rev 4929)
@@ -104,8 +104,8 @@
 *     ..
 *     .. External Functions ..
       REAL               CLANHS, SLAMCH
-      COMPLEX            CLADIV
-      EXTERNAL           CLANHS, SLAMCH, CLADIV
+      COMPLEX            WCLADIV
+      EXTERNAL           CLANHS, SLAMCH, WCLADIV
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           CCOPY, CLARFG, CSCAL
@@ -221,7 +221,7 @@
                Y = SQRT( X*X+U )
                IF( REAL( X )*REAL( Y )+AIMAG( X )*AIMAG( Y ).LT.RZERO )
      $            Y = -Y
-               T = T - CLADIV( U, ( X+Y ) )
+               T = T - WCLADIV( U, ( X+Y ) )
             END IF
          END IF
 *

Modified: trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/zlahqr.f
===================================================================
--- trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/zlahqr.f	2008-11-03 01:25:00 UTC (rev 4928)
+++ trunk/scipy/sparse/linalg/eigen/arpack/ARPACK/LAPACK/zlahqr.f	2008-11-03 01:25:53 UTC (rev 4929)
@@ -104,8 +104,8 @@
 *     ..
 *     .. External Functions ..
       DOUBLE PRECISION   ZLANHS, DLAMCH
-      COMPLEX*16         ZLADIV
-      EXTERNAL           ZLANHS, DLAMCH, ZLADIV
+      COMPLEX*16         WZLADIV
+      EXTERNAL           ZLANHS, DLAMCH, WZLADIV
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           ZCOPY, ZLARFG, ZSCAL
@@ -221,7 +221,7 @@
                Y = SQRT( X*X+U )
                IF( DBLE( X )*DBLE( Y )+DIMAG( X )*DIMAG( Y ).LT.RZERO )
      $            Y = -Y
-               T = T - ZLADIV( U, ( X+Y ) )
+               T = T - WZLADIV( U, ( X+Y ) )
             END IF
          END IF
 *