--- rpl/lapack/lapack/zhetri2x.f 2018/05/29 07:18:21 1.12 +++ rpl/lapack/lapack/zhetri2x.f 2023/08/07 08:39:25 1.13 @@ -113,17 +113,14 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date June 2017 -* *> \ingroup complex16HEcomputational * * ===================================================================== SUBROUTINE ZHETRI2X( UPLO, N, A, LDA, IPIV, WORK, NB, INFO ) * -* -- LAPACK computational routine (version 3.7.1) -- +* -- LAPACK computational routine -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* June 2017 * * .. Scalar Arguments .. CHARACTER UPLO @@ -242,8 +239,8 @@ ELSE * 2 x 2 diagonal NNB T = ABS ( WORK(K+1,1) ) - AK = REAL ( A( K, K ) ) / T - AKP1 = REAL ( A( K+1, K+1 ) ) / T + AK = DBLE ( A( K, K ) ) / T + AKP1 = DBLE ( A( K+1, K+1 ) ) / T AKKP1 = WORK(K+1,1) / T D = T*( AK*AKP1-ONE ) WORK(K,INVD) = AKP1 / D @@ -419,8 +416,8 @@ ELSE * 2 x 2 diagonal NNB T = ABS ( WORK(K-1,1) ) - AK = REAL ( A( K-1, K-1 ) ) / T - AKP1 = REAL ( A( K, K ) ) / T + AK = DBLE ( A( K-1, K-1 ) ) / T + AKP1 = DBLE ( A( K, K ) ) / T AKKP1 = WORK(K-1,1) / T D = T*( AK*AKP1-ONE ) WORK(K-1,INVD) = AKP1 / D