--- rpl/lapack/lapack/zheequb.f	2018/05/29 07:18:19	1.15
+++ rpl/lapack/lapack/zheequb.f	2023/08/07 08:39:23	1.16
@@ -117,8 +117,6 @@
 *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.
 *
-*> \date April 2012
-*
 *> \ingroup complex16HEcomputational
 *
 *> \par References:
@@ -132,10 +130,9 @@
 *  =====================================================================
       SUBROUTINE ZHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK, INFO )
 *
-*  -- LAPACK computational routine (version 3.8.0) --
+*  -- LAPACK computational routine --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     April 2012
 *
 *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, N
@@ -266,12 +263,12 @@
 *        avg = s^T beta / n
          AVG = 0.0D0
          DO I = 1, N
-            AVG = AVG + S( I )*WORK( I )
+            AVG = AVG + DBLE( S( I )*WORK( I ) )
          END DO
          AVG = AVG / N
 
          STD = 0.0D0
-         DO I = N+1, N
+         DO I = N+1, 2*N
             WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
          END DO
          CALL ZLASSQ( N, WORK( N+1 ), 1, SCALE, SUMSQ )
@@ -283,8 +280,8 @@
             T = CABS1( A( I, I ) )
             SI = S( I )
             C2 = ( N-1 ) * T
-            C1 = ( N-2 ) * ( WORK( I ) - T*SI )
-            C0 = -(T*SI)*SI + 2*WORK( I )*SI - N*AVG
+            C1 = ( N-2 ) * ( DBLE( WORK( I ) ) - T*SI )
+            C0 = -(T*SI)*SI + 2 * DBLE( WORK( I ) ) * SI - N*AVG
             D = C1*C1 - 4*C0*C2
 
             IF ( D .LE. 0 ) THEN
@@ -319,7 +316,7 @@
                END DO
             END IF
 
-            AVG = AVG + ( U + WORK( I ) ) * D / N
+            AVG = AVG + ( U + DBLE( WORK( I ) ) ) * D / N
             S( I ) = SI
          END DO
       END DO