27 lines
869 B
Plaintext
27 lines
869 B
Plaintext
$NetBSD: patch-ac,v 1.1.1.1 2008/03/09 17:52:57 tnn Exp $
|
|
|
|
Loosen types (real -> complex) to allow building with f2c.
|
|
|
|
--- cdiag.f.orig 1993-04-20 05:38:10.000000000 +0400
|
|
+++ cdiag.f 2008-03-05 21:26:56.000000000 +0300
|
|
@@ -195,8 +195,8 @@
|
|
C
|
|
C TO FIND THE EIGENVALUES AND VECTORS OF A TRI-DIAGONAL
|
|
C HERMITIAN MATRIX.
|
|
- REAL VALUE(*),W(*),PCK(2),ONE,ZERO,VEC(*)
|
|
- COMPLEX A(*),B(*),DN,UPCK
|
|
+ REAL VALUE(*),PCK(2),ONE,ZERO
|
|
+ COMPLEX A(*),B(*),DN,UPCK,W(*),VEC(*)
|
|
EQUIVALENCE (PCK(1),UPCK)
|
|
C WE TREAT VEC AS IF IT IS DEFINED AS COMPLEX VEC(IV,N)
|
|
C IN THE CALLING PROGRAM.
|
|
@@ -395,7 +395,7 @@
|
|
90 RETURN
|
|
END
|
|
SUBROUTINE ME08B (A,Q,B,N,IA)
|
|
- REAL A(IA,*),Q(2,*),B(IA,*)
|
|
+ COMPLEX A(IA,*),Q(2,*),B(IA,*)
|
|
DO 10 I=1,N
|
|
A(1,I)=A(1,I) -Q(1,1)*B(1,I)+Q(2,1)*B(2,I)
|
|
1 -Q(1,I)*B(1,1)+Q(2,I)*B(2,1)
|