From 633d65dc6afe9b12342f3d6add4bfc1a561a0b7f Mon Sep 17 00:00:00 2001 From: Eugene Chereshnev Date: Mon, 29 Aug 2016 13:31:18 +0700 Subject: [PATCH] Fix incorrect call to GEMM in ctrevc3.f and ztrevc3.f: alpha should have complex type --- SRC/ctrevc3.f | 2 +- SRC/ztrevc3.f | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/SRC/ctrevc3.f b/SRC/ctrevc3.f index 279f147..1340116 100644 --- 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 --- 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, -- 2.7.4