From 2326ea03981cefadcbf94d330e173e6d1da2e12a Mon Sep 17 00:00:00 2001 From: igor175 Date: Thu, 1 Nov 2012 19:37:46 +0000 Subject: [PATCH] fixesd inconsistencies between 'L' and 'U' parts in (z,c)lahef.f, added an IF check to 'U' part IF(KP.GT.1) where column interchange occurs --- SRC/clahef.f | 3 ++- SRC/zlahef.f | 3 ++- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/SRC/clahef.f b/SRC/clahef.f index f10649f..6deeaab 100644 --- a/SRC/clahef.f +++ b/SRC/clahef.f @@ -347,7 +347,8 @@ CALL CCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ), $ LDA ) CALL CLACGV( KK-1-KP, A( KP, KP+1 ), LDA ) - CALL CCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 ) + IF( KP.GT.1 ) + $ CALL CCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 ) * * Interchange rows KK and KP in last KK columns of A and W * diff --git a/SRC/zlahef.f b/SRC/zlahef.f index 705e8fc..cbe7633 100644 --- a/SRC/zlahef.f +++ b/SRC/zlahef.f @@ -347,7 +347,8 @@ CALL ZCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ), $ LDA ) CALL ZLACGV( KK-1-KP, A( KP, KP+1 ), LDA ) - CALL ZCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 ) + IF( KP.GT.1 ) + $ CALL ZCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 ) * * Interchange rows KK and KP in last KK columns of A and W * -- 2.7.4