diff --git a/SRC/cgbrfsx.f b/SRC/cgbrfsx.f
index ab142be46..25d24f8cc 100644
--- a/SRC/cgbrfsx.f
+++ b/SRC/cgbrfsx.f
@@ -513,7 +513,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -527,7 +527,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/cgerfsx.f b/SRC/cgerfsx.f
index e707e9cf2..624bfb31c 100644
--- a/SRC/cgerfsx.f
+++ b/SRC/cgerfsx.f
@@ -490,7 +490,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -504,7 +504,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS(LA_LINRX_ITHRESH_I) = ITHRESH
+            PARAMS(LA_LINRX_ITHRESH_I) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/cherfsx.f b/SRC/cherfsx.f
index c29cd5303..da541616c 100644
--- a/SRC/cherfsx.f
+++ b/SRC/cherfsx.f
@@ -475,7 +475,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -489,7 +489,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/cla_gbamv.f b/SRC/cla_gbamv.f
index d57ba063a..3ae20616f 100644
--- a/SRC/cla_gbamv.f
+++ b/SRC/cla_gbamv.f
@@ -287,7 +287,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/cla_geamv.f b/SRC/cla_geamv.f
index 6e81a989e..da3c522bb 100644
--- a/SRC/cla_geamv.f
+++ b/SRC/cla_geamv.f
@@ -274,7 +274,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/cla_heamv.f b/SRC/cla_heamv.f
index dc7413427..83aae5017 100644
--- a/SRC/cla_heamv.f
+++ b/SRC/cla_heamv.f
@@ -261,7 +261,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/cla_lin_berr.f b/SRC/cla_lin_berr.f
index 6341abb7d..b9305a8b9 100644
--- a/SRC/cla_lin_berr.f
+++ b/SRC/cla_lin_berr.f
@@ -126,7 +126,7 @@
       REAL               SAFE1
 *     ..
 *     .. Statement Functions ..
-      COMPLEX            CABS1
+      REAL               CABS1
 *     ..
 *     .. Statement Function Definitions ..
       CABS1( CDUM ) = ABS( REAL( CDUM ) ) + ABS( AIMAG( CDUM ) )
@@ -138,7 +138,7 @@
 *     to compute AYB.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (NZ+1)*SAFE1
+      SAFE1 = REAL( NZ+1 )*SAFE1
 
       DO J = 1, NRHS
          BERR(J) = 0.0
diff --git a/SRC/cla_syamv.f b/SRC/cla_syamv.f
index 36ec24d81..b793f401f 100644
--- a/SRC/cla_syamv.f
+++ b/SRC/cla_syamv.f
@@ -263,7 +263,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/cporfsx.f b/SRC/cporfsx.f
index 4b7c09372..7658762b3 100644
--- a/SRC/cporfsx.f
+++ b/SRC/cporfsx.f
@@ -467,7 +467,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -481,7 +481,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS(LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/csyrfsx.f b/SRC/csyrfsx.f
index ed38a1ed8..978fd0718 100644
--- a/SRC/csyrfsx.f
+++ b/SRC/csyrfsx.f
@@ -477,7 +477,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -491,7 +491,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/dgbrfsx.f b/SRC/dgbrfsx.f
index c601f0f59..fa1eeec6c 100644
--- a/SRC/dgbrfsx.f
+++ b/SRC/dgbrfsx.f
@@ -516,7 +516,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/dgerfsx.f b/SRC/dgerfsx.f
index 0beeff12f..4ed6b8e76 100644
--- a/SRC/dgerfsx.f
+++ b/SRC/dgerfsx.f
@@ -489,7 +489,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/dporfsx.f b/SRC/dporfsx.f
index a2549e125..3404a8d9a 100644
--- a/SRC/dporfsx.f
+++ b/SRC/dporfsx.f
@@ -468,7 +468,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/dsyrfsx.f b/SRC/dsyrfsx.f
index a0e28affb..0af3672c9 100644
--- a/SRC/dsyrfsx.f
+++ b/SRC/dsyrfsx.f
@@ -475,7 +475,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/sgbrfsx.f b/SRC/sgbrfsx.f
index bb827d31f..1a7f5598a 100644
--- a/SRC/sgbrfsx.f
+++ b/SRC/sgbrfsx.f
@@ -516,7 +516,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -530,7 +530,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/sgerfsx.f b/SRC/sgerfsx.f
index c8d518cc7..78eebacba 100644
--- a/SRC/sgerfsx.f
+++ b/SRC/sgerfsx.f
@@ -489,7 +489,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -503,7 +503,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/sla_gbamv.f b/SRC/sla_gbamv.f
index c1714cd18..9260ca0ae 100644
--- a/SRC/sla_gbamv.f
+++ b/SRC/sla_gbamv.f
@@ -277,7 +277,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/sla_geamv.f b/SRC/sla_geamv.f
index 957fe3ff5..1de40cf85 100644
--- a/SRC/sla_geamv.f
+++ b/SRC/sla_geamv.f
@@ -264,7 +264,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/sla_lin_berr.f b/SRC/sla_lin_berr.f
index 3742eb656..4c6a32e91 100644
--- a/SRC/sla_lin_berr.f
+++ b/SRC/sla_lin_berr.f
@@ -131,7 +131,7 @@
 *     to compute AYB.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (NZ+1)*SAFE1
+      SAFE1 = REAL( NZ+1 )*SAFE1
 
       DO J = 1, NRHS
          BERR(J) = 0.0
diff --git a/SRC/sla_syamv.f b/SRC/sla_syamv.f
index 39fd3062d..ac4604a77 100644
--- a/SRC/sla_syamv.f
+++ b/SRC/sla_syamv.f
@@ -252,7 +252,7 @@
 *     number of additions in each row.
 *
       SAFE1 = SLAMCH( 'Safe minimum' )
-      SAFE1 = (N+1)*SAFE1
+      SAFE1 = REAL( N+1 )*SAFE1
 *
 *     Form  y := alpha*abs(A)*abs(x) + beta*abs(y).
 *
diff --git a/SRC/sporfsx.f b/SRC/sporfsx.f
index abdbe16a1..448c66a9e 100644
--- a/SRC/sporfsx.f
+++ b/SRC/sporfsx.f
@@ -468,7 +468,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -482,7 +482,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/ssyrfsx.f b/SRC/ssyrfsx.f
index 6a47f1168..7dca5dea2 100644
--- a/SRC/ssyrfsx.f
+++ b/SRC/ssyrfsx.f
@@ -475,7 +475,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
@@ -489,7 +489,7 @@
 *
       IF ( NPARAMS.GE.LA_LINRX_ITHRESH_I ) THEN
          IF ( PARAMS( LA_LINRX_ITHRESH_I ).LT.0.0 ) THEN
-            PARAMS( LA_LINRX_ITHRESH_I ) = ITHRESH
+            PARAMS( LA_LINRX_ITHRESH_I ) = REAL( ITHRESH )
          ELSE
             ITHRESH = INT( PARAMS( LA_LINRX_ITHRESH_I ) )
          END IF
diff --git a/SRC/zgbrfsx.f b/SRC/zgbrfsx.f
index 2c309e694..260c27088 100644
--- a/SRC/zgbrfsx.f
+++ b/SRC/zgbrfsx.f
@@ -513,7 +513,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/zgerfsx.f b/SRC/zgerfsx.f
index 4b7a2c439..4af17c5be 100644
--- a/SRC/zgerfsx.f
+++ b/SRC/zgerfsx.f
@@ -490,7 +490,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/zherfsx.f b/SRC/zherfsx.f
index 77a215fd5..e202b7971 100644
--- a/SRC/zherfsx.f
+++ b/SRC/zherfsx.f
@@ -475,7 +475,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/zla_lin_berr.f b/SRC/zla_lin_berr.f
index 7f96ed4fa..69b58f303 100644
--- a/SRC/zla_lin_berr.f
+++ b/SRC/zla_lin_berr.f
@@ -126,7 +126,7 @@
       DOUBLE PRECISION   SAFE1
 *     ..
 *     .. Statement Functions ..
-      COMPLEX*16         CABS1
+      DOUBLE PRECISION   CABS1
 *     ..
 *     .. Statement Function Definitions ..
       CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( DIMAG( CDUM ) )
diff --git a/SRC/zporfsx.f b/SRC/zporfsx.f
index 7f2a8405e..928c7b7cb 100644
--- a/SRC/zporfsx.f
+++ b/SRC/zporfsx.f
@@ -467,7 +467,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/SRC/zsyrfsx.f b/SRC/zsyrfsx.f
index 11e4fd1fb..196a09884 100644
--- a/SRC/zsyrfsx.f
+++ b/SRC/zsyrfsx.f
@@ -477,7 +477,7 @@
          IF ( PARAMS( LA_LINRX_ITREF_I ) .LT. 0.0D+0 ) THEN
             PARAMS( LA_LINRX_ITREF_I ) = ITREF_DEFAULT
          ELSE
-            REF_TYPE = PARAMS( LA_LINRX_ITREF_I )
+            REF_TYPE = INT( PARAMS( LA_LINRX_ITREF_I ) )
          END IF
       END IF
 *
diff --git a/TESTING/LIN/cebchvxx.f b/TESTING/LIN/cebchvxx.f
index a10a8f6f7..55206c645 100644
--- a/TESTING/LIN/cebchvxx.f
+++ b/TESTING/LIN/cebchvxx.f
@@ -120,7 +120,6 @@
 *     .. Local Arrays ..
       REAL               TSTRAT(NTESTS), RINV(NMAX), PARAMS(NPARAMS),
      $                   S(NMAX), R(NMAX),C(NMAX),RWORK(3*NMAX),
-     $                   DIFF(NMAX, NMAX),
      $                   ERRBND_N(NMAX*3), ERRBND_C(NMAX*3)
       INTEGER            IPIV(NMAX)
       COMPLEX            A(NMAX,NMAX),INVHILB(NMAX,NMAX),X(NMAX,NMAX),
@@ -254,13 +253,6 @@
             END IF
          END IF
 
-*        Calculating the difference between C**SVXX's X and the true X.
-         DO I = 1,N
-            DO J =1,NRHS
-               DIFF(I,J) = X(I,J) - INVHILB(I,J)
-            END DO
-         END DO
-
 *        Calculating the RCOND
          RNORM = 0
          RINORM = 0
diff --git a/TESTING/LIN/zebchvxx.f b/TESTING/LIN/zebchvxx.f
index b57ce88ec..b2f18bc29 100644
--- a/TESTING/LIN/zebchvxx.f
+++ b/TESTING/LIN/zebchvxx.f
@@ -120,7 +120,6 @@
 *     .. Local Arrays ..
       DOUBLE PRECISION   TSTRAT(NTESTS), RINV(NMAX), PARAMS(NPARAMS),
      $                   S(NMAX),R(NMAX),C(NMAX),RWORK(3*NMAX),
-     $                   DIFF(NMAX, NMAX),
      $                   ERRBND_N(NMAX*3), ERRBND_C(NMAX*3)
       INTEGER            IPIV(NMAX)
       COMPLEX*16         A(NMAX,NMAX),INVHILB(NMAX,NMAX),X(NMAX,NMAX),
@@ -254,13 +253,6 @@
             END IF
          END IF
 
-*        Calculating the difference between Z**SVXX's X and the true X.
-         DO I = 1,N
-            DO J =1,NRHS
-               DIFF(I,J) = X(I,J) - INVHILB(I,J)
-            END DO
-         END DO
-
 *        Calculating the RCOND
          RNORM = 0
          RINORM = 0