APPLYING INTEL PATCHES sent to Julie on Feb 19th 2016 by Dima from INTEL (dmitry...

[PATCH 19/42] Fix ?GBRFSX - do not proceed if INFO.NE.0 after ?GBCON

- and fix description of B and BERR in ZGBRFSX

diff --git a/SRC/cgbrfsx.f b/SRC/cgbrfsx.f
index fc73496910bab6818370bde7c8e7be4c5a1f7c16..1caee82ac78b67f60684d396247eebb09f91ec50 100644 (file)
--- a/SRC/cgbrfsx.f
+++ b/SRC/cgbrfsx.f
@@ -642,7 +642,7 @@
 *
 *     Perform refinement on each right-hand side
 *
-      IF ( REF_TYPE .NE. 0 ) THEN
+      IF ( REF_TYPE .NE. 0 .AND. INFO .EQ. 0 ) THEN
 
          PREC_TYPE = ILAPREC( 'D' )
 

diff --git a/SRC/dgbrfsx.f b/SRC/dgbrfsx.f
index c96c6233886ba98145d40b2df03df7980501a3db..92ed4ecf47388e039750d6491a451fed78dada56 100644 (file)
--- a/SRC/dgbrfsx.f
+++ b/SRC/dgbrfsx.f
@@ -646,7 +646,7 @@
 *
 *     Perform refinement on each right-hand side
 *
-      IF (REF_TYPE .NE. 0) THEN
+      IF ( REF_TYPE .NE. 0 .AND. INFO .EQ. 0 ) THEN
 
          PREC_TYPE = ILAPREC( 'E' )
 

diff --git a/SRC/sgbrfsx.f b/SRC/sgbrfsx.f
index 7245f867c87badb6f5f820097ca2d674a1beeafe..8c9efc5a51d33db5904ce1d8771e269690f1a3ee 100644 (file)
--- a/SRC/sgbrfsx.f
+++ b/SRC/sgbrfsx.f
@@ -646,7 +646,7 @@
 *
 *     Perform refinement on each right-hand side
 *
-      IF (REF_TYPE .NE. 0) THEN
+      IF ( REF_TYPE .NE. 0 .AND. INFO .EQ. 0 ) THEN
 
          PREC_TYPE = ILAPREC( 'D' )
 

diff --git a/SRC/zgbrfsx.f b/SRC/zgbrfsx.f
index 2b81d403b98a3428c90bfb02a04bea345e79ced7..0114d0dc4732fdf50b5e81c70de84f4d7e1f8d48 100644 (file)
--- a/SRC/zgbrfsx.f
+++ b/SRC/zgbrfsx.f
@@ -195,7 +195,7 @@
 *>
 *> \param[in] B
 *> \verbatim
-*>          B is DOUBLE PRECISION array, dimension (LDB,NRHS)
+*>          B is COMPLEX*16 array, dimension (LDB,NRHS)
 *>     The right hand side matrix B.
 *> \endverbatim
 *>
@@ -232,7 +232,7 @@
 *>
 *> \param[out] BERR
 *> \verbatim
-*>          BERR is COMPLEX*16 array, dimension (NRHS)
+*>          BERR is DOUBLE PRECISION array, dimension (NRHS)
 *>     Componentwise relative backward error.  This is the
 *>     componentwise relative backward error of each solution vector X(j)
 *>     (i.e., the smallest relative change in any element of A or B that
@@ -642,7 +642,7 @@
 *
 *     Perform refinement on each right-hand side
 *
-      IF ( REF_TYPE .NE. 0 ) THEN
+      IF ( REF_TYPE .NE. 0 .AND. INFO .EQ. 0 ) THEN
 
          PREC_TYPE = ILAPREC( 'E' )