This is related to issue #42.

This is a fix for double precision. Need to fix S, C and Z now.

diff --git a/TESTING/dbal.in b/TESTING/dbal.in
index 94268a117b0d3b3d674dbff3ee0e5d62e869b958..c00fe3e83b823b1ddfd9c9842f2099c3c2603266 100644 (file)
--- a/TESTING/dbal.in
+++ b/TESTING/dbal.in
@@ -87,13 +87,13 @@ DBL:  Tests DGEBAL
   0.0000D+00  0.8192D+04  0.0000D+00  0.0000D+00  0.8000D+01
 
    1   5
-  1.0000D+00     0.0000D-03     0.0000D-03     0.0000D-03   250.0000D-03
+  1.0000D+00     0.0000D-03     0.0000D-03     0.0000D-03     2.0000D+00
   0.0000D-03     2.0000D+00     1.0240D+03    16.0000D+00    16.0000D+00
-  256.0000D-03     1.0000D-03     4.0000D+00     0.0000D-03     2.0480D+03
+  3.2000D-02     1.0000D-03     4.0000D+00     0.0000D-03     2.0480D+03
   0.0000D-03   250.0000D-03    16.0000D+00     4.0000D+00     4.0000D+00
   0.0000D-03     2.0480D+03     0.0000D-03     0.0000D-03     8.0000D+00
 
-  64.0000D+00  500.0000D-03  62.5000D-03  4.0000D+00  2.0000D+00
+  8.0000D+00  500.0000D-03  62.5000D-03  4.0000D+00  2.0000D+00
 
   4
   0.1000D+01  0.1000D+07  0.1000D+07  0.1000D+07
@@ -167,15 +167,15 @@ DBL:  Tests DGEBAL
   0.0000D+00  0.8000D+01  0.0000D+00  0.4000D-02  0.1250D+00 -0.2000D+00 0.3000D+01
 
   2   5
-  6.4000D+01   2.5000D-01   5.00000D-01   0.0000D+00   0.0000D+00   1.0000D+00  -2.0000D+00
-  0.0000D+00   4.0000D+00   2.00000D+00   4.0960D+00   1.6000D+00   0.0000D+00   1.0240D+01
-  0.0000D+00   5.0000D-01   3.00000D+00   4.0960D+00   1.0000D+00   0.0000D+00  -6.4000D+00
-  0.0000D+00   1.0000D+00  -3.90625D+00   1.0000D+00  -3.1250D+00   0.0000D+00   8.0000D+00
-  0.0000D+00  -2.0000D+00   4.00000D+00   1.6000D+00   2.0000D+00  -8.0000D+00   8.0000D+00
+  6.4000D+01   1.0000D+00   5.00000D-01   0.0000D+00   0.0000D+00   1.0000D+00  -2.0000D+00
+  0.0000D+00   4.0000D+00   5.00000D-01   1.0240D+00   8.0000D-01   0.0000D+00   2.5600D+00
+  0.0000D+00   2.0000D+00   3.00000D+00   4.0960D+00   2.0000D+00   0.0000D+00  -6.4000D+00
+  0.0000D+00   4.0000D+00  -3.90625D+00   1.0000D+00  -6.2500D+00   0.0000D+00   8.0000D+00
+  0.0000D+00  -4.0000D+00   2.00000D+00   8.0000D-01   2.0000D+00  -4.0000D+00   4.0000D+00
   0.0000D+00   0.0000D+00   0.00000D+00   0.0000D+00   0.0000D+00   6.0000D+00   1.0000D+00
   0.0000D+00   0.0000D+00   0.00000D+00   0.0000D+00   0.0000D+00   0.0000D+00   0.0000D+00
 
-  3.0000D+00  1.953125D-03  3.1250D-02  3.2000D+01  2.5000D-01  1.0000D+00 6.0000D+00
+  3.0000D+00  7.812500D-03  3.1250D-02  3.2000D+01  5.0000D-01  1.0000D+00 6.0000D+00
 
   5
   0.1000D+04  0.2000D+01  0.3000D+01  0.4000D+01  0.5000D+06
@@ -185,13 +185,13 @@ DBL:  Tests DGEBAL
   0.6000D+01  0.2000D+03  0.1000D+01  0.6000D+03  0.3000D+01
 
   1   5
-  1.0000D+03   3.1250D-02   3.7500D-01   6.2500D-02   3.90625D+03
-  5.7600D+02   0.0000D+00   1.6000D-03   1.0000D+00   1.5000D+00
-  0.0000D+00  -3.7500D+01   2.0000D+00   1.2500D-01   6.2500D-02
-  5.7600D+02   2.0000D-03   8.0000D+00   1.0000D+00  -5.0000D+02
-  7.6800D+02   4.0000D+02   1.6000D+01   1.2000D+03   3.0000D+00
+  1.0000D+03     3.1250D-02     3.7500D-01     3.1250D-02     1.95312500D+03
+  5.7600D+02     0.0000D+00     1.6000D-03     5.0000D-01     7.50000000D-01
+  0.0000D+00    -3.7500D+01     2.0000D+00     6.2500D-02     3.12500000D-02
+  1.1520D+03     4.0000D-03     1.6000D+01     1.0000D+00    -5.00000000D+02
+  1.5360D+03     8.0000D+02     3.2000D+01     1.2000D+03     3.00000000D+00
 
-  1.2800D+02  2.0000D+00  1.6000D+01  2.0000D+00  1.0000D+00
+  3.2000D+01  5.0000D-01   4.0000D+00   2.5000D-01  1.2500D-01
 
   6
   1.0000D+00  1.0000D+120 0.0000D+00  0.0000D+00  0.0000D+00 0.0000D+00
@@ -202,14 +202,13 @@ DBL:  Tests DGEBAL
   0.0000D+00  0.0000D+00  0.0000D+00  0.0000D+00  1.0000D-120 1.0000D+00
 
    1   6
-   1.000000000000000000D+00  6.344854593289122931D+03  0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00
-   1.576080247855779135D-04  1.000000000000000000D+00  6.344854593289122931D+03  0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00
-   0.000000000000000000D+00  1.576080247855779135D-04  1.000000000000000000D+00  3.172427296644561466D+03  0.000000000000000000D+00  0.000000000000000000D+00
-   0.000000000000000000D+00  0.000000000000000000D+00  3.152160495711558270D-04  1.000000000000000000D+00  1.586213648322280733D+03  0.000000000000000000D+00
-   0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00  6.304320991423116539D-04  1.000000000000000000D+00  1.586213648322280733D+03
-   0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00  0.000000000000000000D+00  6.304320991423116539D-04  1.000000000000000000D+00
-
-  2.494800386918399765D+291  1.582914569427869018D+175  1.004336277661868922D+59  3.186183822264904554D-58  5.053968264940243633D-175  8.016673440035891112D-292
-
+0.10000000000000000000D+01 0.63448545932891229313D+04 0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00
+0.15760802478557791348D-03 0.10000000000000000000D+01 0.63448545932891229313D+04 0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00
+0.00000000000000000000D+00 0.15760802478557791348D-03 0.10000000000000000000D+01 0.31724272966445614657D+04 0.00000000000000000000D+00 0.00000000000000000000D+00
+0.00000000000000000000D+00 0.00000000000000000000D+00 0.31521604957115582695D-03 0.10000000000000000000D+01 0.15862136483222807328D+04 0.00000000000000000000D+00
+0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00 0.63043209914231165391D-03 0.10000000000000000000D+01 0.79310682416114036641D+03
+0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00 0.00000000000000000000D+00 0.12608641982846233078D-02 0.10000000000000000000D+01
+
+  2.494800386918399765D+291  1.582914569427869018D+175  1.004336277661868922D+59  3.186183822264904554D-58  5.053968264940243633D-175  0.40083367200179455560D-291;
 
 0