BUG: sparse/isolve: fix errors in Fortran code pointed to by gcc's wa…

…rnings
fish2000 · Jul 7, 2011 · 566fc2f · 566fc2f
1 parent e86854b
commit 566fc2f
Show file tree

Hide file tree

Showing 4 changed files with 7 additions and 4 deletions.
diff --git a/scipy/sparse/linalg/isolve/iterative/BiCGSTABREVCOM.f.src b/scipy/sparse/linalg/isolve/iterative/BiCGSTABREVCOM.f.src
@@ -107,7 +107,7 @@
 *  ==============================================================
 *
 *     .. Parameters ..
-      DOUBLE PRECISION   ZERO, ONE
+      <rt>             ZERO, ONE
       PARAMETER        ( ZERO = 0.0D+0, ONE = 1.0D+0 )
 *     ..
 *     .. Local Scalars ..

diff --git a/scipy/sparse/linalg/isolve/iterative/CGREVCOM.f.src b/scipy/sparse/linalg/isolve/iterative/CGREVCOM.f.src
@@ -94,7 +94,7 @@
 *  ============================================================
 *
 *     .. Parameters ..
-      DOUBLE PRECISION   ZERO, ONE
+      <rt>             ZERO, ONE
       PARAMETER        ( ZERO = 0.0D+0, ONE = 1.0D+0 )
 *     ..
 *     .. Local Scalars ..

diff --git a/scipy/sparse/linalg/isolve/iterative/GMRESREVCOM.f.src b/scipy/sparse/linalg/isolve/iterative/GMRESREVCOM.f.src
@@ -106,6 +106,7 @@
      $                    NEED1, NEED2
       <_t>  <xdot=sdot,ddot,cdotc,zdotc>
       <_t>  toz
+      <_t>    TMPVAL
       <rt>    BNRM2, RNORM, TOL,  
      $     <rc=s,d,sc,dz>NRM2,
      $     <rc>APPROXRES
@@ -267,7 +268,8 @@
          RNORM = <rc>NRM2( N, WORK( 1,V ), 1 )
          toz = ONE/RNORM
          CALL <_c>SCAL( N, toz, WORK( 1,V ), 1 )
-         CALL <_c>ELEMVEC( 1, N, RNORM, WORK( 1,S ) )
+         TMPVAL = RNORM
+         CALL <_c>ELEMVEC( 1, N, TMPVAL, WORK( 1,S ) )
 *
 *         DO 50 I = 1, RESTRT
          i = 1

diff --git a/scipy/sparse/linalg/isolve/iterative/QMRREVCOM.f.src b/scipy/sparse/linalg/isolve/iterative/QMRREVCOM.f.src
@@ -120,7 +120,8 @@
 *     .. Local Scalars ..
       INTEGER          R, D, P, PTLD, Q, S, V, VTLD, W, WTLD, Y, YTLD,
      $                 Z, ZTLD, MAXIT, NEED1, NEED2
-      <rt>   TOL, BNRM2, RHOTOL, BETATOL, GAMMATOL, DELTATOL, 
+      <rt>   TOL, BNRM2, RHOTOL, BETATOL,
+     $       GAMMATOL, DELTATOL,
      $       EPSTOL, XITOL,
      $       <sdsd=s,d,s,d>GETBREAK, 
      $       <rc=s,d,sc,dz>NRM2