Fix incorrect call to GEMM in ctrevc3.f and ztrevc3.f: alpha should have complex...

diff --git a/SRC/ctrevc3.f b/SRC/ctrevc3.f
index 279f147..1340116 100644 (file)
--- a/SRC/ctrevc3.f
+++ b/SRC/ctrevc3.f
@@ -593,7 +593,7 @@
 *              When the number of vectors stored reaches NB,
 *              or if this was last vector, do the GEMM
                IF( (IV.EQ.NB) .OR. (KI.EQ.N) ) THEN
-                  CALL CGEMM( 'N', 'N', N, IV, N-KI+IV, ONE,
+                  CALL CGEMM( 'N', 'N', N, IV, N-KI+IV, CONE,
      $                        VL( 1, KI-IV+1 ), LDVL,
      $                        WORK( KI-IV+1 + (1)*N ), N,
      $                        CZERO,

diff --git a/SRC/ztrevc3.f b/SRC/ztrevc3.f
index ad23ae7..6656852 100644 (file)
--- a/SRC/ztrevc3.f
+++ b/SRC/ztrevc3.f
@@ -593,7 +593,7 @@
 *              When the number of vectors stored reaches NB,
 *              or if this was last vector, do the GEMM
                IF( (IV.EQ.NB) .OR. (KI.EQ.N) ) THEN
-                  CALL ZGEMM( 'N', 'N', N, IV, N-KI+IV, ONE,
+                  CALL ZGEMM( 'N', 'N', N, IV, N-KI+IV, CONE,
      $                        VL( 1, KI-IV+1 ), LDVL,
      $                        WORK( KI-IV+1 + (1)*N ), N,
      $                        CZERO,