1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
$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)
|
