Fix scaling in the transpose case of extra-precise refinement.

The row-scaling factors need passed in the transpose case.  Apparently,
none of our release tests include ill-scaled matrices.

Signed-off-by: Jason Riedy <ejr@cs.berkeley.edu>

diff --git a/SRC/cgbrfsx.f b/SRC/cgbrfsx.f
index 63abc5d7be4bf625f2bc9bfd426036f73488e17e..e357b8c743456996f8e9b9bf82732952837c6586 100644 (file)
--- a/SRC/cgbrfsx.f
+++ b/SRC/cgbrfsx.f
@@ -515,7 +515,7 @@
      $           INFO )
          ELSE
             CALL CLA_GBRFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N, KL, KU,
-     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, C, B,
+     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), RWORK, WORK(1), RWORK, RCOND,
      $           ITHRESH, RTHRESH, UNSTABLE_THRESH, IGNORE_CWISE,

diff --git a/SRC/cgerfsx.f b/SRC/cgerfsx.f
index 3e1288f10e802a1038798ba621fb1d3a2149bc92..4979a907e2a8bbf0b79661d29b9b43a4d1470483 100644 (file)
--- a/SRC/cgerfsx.f
+++ b/SRC/cgerfsx.f
@@ -500,7 +500,7 @@
      $           INFO )
          ELSE
             CALL CLA_GERFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N,
-     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, C, B,
+     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), RWORK, WORK(1), RWORK, RCOND,
      $           ITHRESH, RTHRESH, UNSTABLE_THRESH, IGNORE_CWISE,

diff --git a/SRC/dgbrfsx.f b/SRC/dgbrfsx.f
index fa0f8db59fe80a3c7be1019493175eaaf3e378da..d6dbd184b785d5e1342577ca7262e26070203b53 100644 (file)
--- a/SRC/dgbrfsx.f
+++ b/SRC/dgbrfsx.f
@@ -519,7 +519,7 @@
      $           IGNORE_CWISE, INFO )
          ELSE
             CALL DLA_GBRFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N, KL, KU,
-     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, C, B,
+     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK( N+1 ), WORK( 1 ), WORK( 2*N+1 ),
      $           WORK( 1 ), RCOND, ITHRESH, RTHRESH, UNSTABLE_THRESH,

diff --git a/SRC/dgerfsx.f b/SRC/dgerfsx.f
index 6ecfd642951cd08bd715568730620a67d8cb87b2..bb1878cbf498041dc51c7fc085ac9d6e0b24a7bd 100644 (file)
--- a/SRC/dgerfsx.f
+++ b/SRC/dgerfsx.f
@@ -500,7 +500,7 @@
      $           IGNORE_CWISE, INFO )
          ELSE
             CALL DLA_GERFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N,
-     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, C, B,
+     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), WORK(1), WORK(2*N+1),
      $           WORK(1), RCOND, ITHRESH, RTHRESH, UNSTABLE_THRESH,

diff --git a/SRC/sgbrfsx.f b/SRC/sgbrfsx.f
index 828297f0b8e1fcfdf52201c022d24bdbb969c984..b43b63d53acb7781f66b79a6a1df619473c5cfd6 100644 (file)
--- a/SRC/sgbrfsx.f
+++ b/SRC/sgbrfsx.f
@@ -519,7 +519,7 @@
      $           IGNORE_CWISE, INFO )
          ELSE
             CALL SLA_GBRFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N, KL, KU,
-     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, C, B,
+     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK( N+1 ), WORK( 1 ), WORK( 2*N+1 ),
      $           WORK( 1 ), RCOND, ITHRESH, RTHRESH, UNSTABLE_THRESH,

diff --git a/SRC/sgerfsx.f b/SRC/sgerfsx.f
index 7933c7e5d95e6832804690c336c9a3ac289c828c..67625446ec45eee392519ea1a7803fe5c02bbd88 100644 (file)
--- a/SRC/sgerfsx.f
+++ b/SRC/sgerfsx.f
@@ -500,7 +500,7 @@
      $           IGNORE_CWISE, INFO )
          ELSE
             CALL SLA_GERFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N,
-     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, C, B,
+     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), WORK(1), WORK(2*N+1),
      $           WORK(1), RCOND, ITHRESH, RTHRESH, UNSTABLE_THRESH,

diff --git a/SRC/zgbrfsx.f b/SRC/zgbrfsx.f
index 8c4bdd89fe72917940407cd78ea405f30842a4b6..3251138b589640f1e86eb02e092ffb3ddd3d72ee 100644 (file)
--- a/SRC/zgbrfsx.f
+++ b/SRC/zgbrfsx.f
@@ -515,7 +515,7 @@
      $           INFO )
          ELSE
             CALL ZLA_GBRFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N, KL, KU,
-     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, C, B,
+     $           NRHS, AB, LDAB, AFB, LDAFB, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), RWORK, WORK(1), RWORK, RCOND,
      $           ITHRESH, RTHRESH, UNSTABLE_THRESH, IGNORE_CWISE,

diff --git a/SRC/zgerfsx.f b/SRC/zgerfsx.f
index edb1bcaef5e017c68e2d099c6ce1079fa28f48c7..f5e1de436b8d17c899950500dde395b343e98dd7 100644 (file)
--- a/SRC/zgerfsx.f
+++ b/SRC/zgerfsx.f
@@ -500,7 +500,7 @@
      $           INFO )
          ELSE
             CALL ZLA_GERFSX_EXTENDED( PREC_TYPE, TRANS_TYPE,  N,
-     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, C, B,
+     $           NRHS, A, LDA, AF, LDAF, IPIV, ROWEQU, R, B,
      $           LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM,
      $           ERR_BNDS_COMP, WORK(N+1), RWORK, WORK(1), RWORK, RCOND,
      $           ITHRESH, RTHRESH, UNSTABLE_THRESH, IGNORE_CWISE,