Het hulpprogramma BDEIG voor het oplossen van algemene eigenwaarde problemen in bandvorm

(1)

Het hulpprogramma BDEIG voor het oplossen van algemene

eigenwaarde problemen in bandvorm

Citation for published version (APA):

Kraker, de, A. (1978). Het hulpprogramma BDEIG voor het oplossen van algemene eigenwaarde problemen in bandvorm. (DCT rapporten; Vol. 1978.003). Technische Hogeschool Eindhoven.

Document status and date: Gepubliceerd: 01/01/1978

Document Version:

Uitgevers PDF, ook bekend als Version of Record

Please check the document version of this publication:

• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.

• The final author version and the galley proof are versions of the publication after peer review.

• The final published version features the final layout of the paper including the volume, issue and page numbers.

Link to publication

General rights

Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.

• Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain

• You may freely distribute the URL identifying the publication in the public portal.

If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:

www.tue.nl/taverne Take down policy

If you believe that this document breaches copyright please contact us at: openaccess@tue.nl

providing details and we will investigate your claim.

(2)

(3)

(4)

(5)

tu,.

L

a

(6)

b

(7)

(8)

0225 'OVD';

0230°PRûCEZURE' EIGBVVALUECNr DIAGr CODIAGJ Klr a r M D A > I !E35 'VALUE' N r K ~ D K2;

a240 'ARRAY' DIAGr CODIAG, L M D &

0245 'REGIN' 'INTEGER' I r lb KSI KLRD )<Ui3f

0250

0255 'ARRAY' CQt0 : rJ

..

13, FktKl : K 2 3 f

'INTEGER' NI KI, K2t

'REAL' U r Vr LKr DETr EPCJ E P S I r L% W r FLBr FU!3r LA, SM;

(9)

0205

0270 m 7 5 0280 0285 0290 0295

0300

0305

031

0 0315

0320

0325

0330

0335

0340

034s

0350

0355

'

0360 0365

'REAL,' 'PROCEDURE' STURM;

*BEEIN'

wm*

s;

U := O; DE3 :=

1;

KS := 0;

'FOR' I :=

1

'STEP'

1

'WTIL' ^J 'KI'

'REGIM' L := (DIAGEI3

-

LK)

*

DET

-

CQCI

-

11 *

Ut S ? = 1; 'IF' ARSCV) > LA *THEN' S := sc?

*ELSE* *IF* ARS(V> SM 'THEW'

'BEGIN' 'IF' ABSCDET) SM 'THEN' S := L A 'END'; U := S

*

DET; DET := S

*

V i

'IF' DET >= 0 'EQUIV' U >=

0

'TEEN' KS := KS

+

1

'END' f

'IF' DET =

0

'AND' U et 'THEN' KS := KS

-

1;

STUWUI := DET

'EJD' ;

'VAtOE' A, BJ F b FBJ KAJ In3i 'REcIiL' AJ BJ FA, F B I

*PROCEDIRE* ISOLATED

z m e b

B~ FA, FR, ICA, KB);

'INTEGm' 1(%L KR;

'REGIX' 'REAL' MJ CJ FCI '1,NTEGER' KCJ

'

IF' ARSCFT?) <= AWCFA)

'THEN' 'BEGIN' C := At FC := FAI KC := KA 'EVD'

'

ELSE' 0370 Al: 'BMiIW' C := Rf FC := F8; KC 2= KBt I3 := Af 037 5 FR := FA; KB X= W A := C i FA := FC; KA := KC 03881

'

BidB' ; 0385 A2: E4

:=

<A + B V S f 0390 0395

ma0

0405

@411

m15 0420 0425 043B 0435 B44EIr f m 5 m59

c45s

0460 0465 0470 m75 0480 0485 0490 0495 05cIJO

0505

051a

0515

0520

a525

(10)

0S30 'EIEGIN' LAMBQAEKJ := LE; FWKI := FLS; KWCKI := KLB 'END';

0535 'FOR' K := K l 'STEP'

1

'LhJTIL' K2 'Do'

054QJ 'BEGIN' 'IF' KLB K 'THEN'

0545 'BEGIN' UR := LB; FUR := FLE; KUB := KLBf

0550

LB := LA!YBDACKI; FLB := Fk'tKI; KLR := KWEK3

0555 'END' i

056~3 TEST: 'IF' UR

-

LR > EPS 'THEN'

0565 'BEGIN' 'IF' KLR

-

KUB >

1

0570 'THEN'

0575 BIS: 'BEGIN' LK := CUI3 + LF3)/2;

s m ;

0580 'IF' Ks K

0585 'THEN' 'BEGIN' U 9 := LK; FUR := DET; KUR := Ks 'END'

0590 'ELSE'

B95 'BEGIN' LB := LK; FLB :=

DET;

PCLB := KS; I := K:

0600 'FOR' I := I 9 I 'WHILE' I e= K2 'MD' I e= KS 'Do'

0S05 'REXIN' LAMBDACII := LK;

0610

F W C X J := DET; %t13 := KS 0615 'END'

0620

'END' ; W25 'GOTO' TEST 0630 'END' 0635 'ELSE' 0640 'IF' Ft@ = B 18645 'THF3iJ' LK := UB

0650 'ELSE' 'IF' FLB = @ 'THEXI' 'GOTO' RIS

0655 'ELSE' ISOLATED ZERO<Lásr UBI FLBr FUBr. f d R r IGJR) 0660 'END'

€3665 'ELSE' LK := <UB + L)7)/2;

0670 ASSIGN: LAMBDACKI := LK 0ó75 'mD' 'END';

Q680 'PROCEDURE' EIGE?Jt.ECTOR(piJr DIAGD CODfAGr Kl0 K20 k D A t

2);

0685 'VALUE' Ns K l s K2;

0690 'INTEGER' N r Klr )c2; 'ARRAY' DIMI CODIAGD LMBEIFlr

22

0695 'BEEIN' 'INTEGER@ Ir KI Ft

e -,.a.- ...--e11 --.-.-II--.-

klim H W L r L m i I I ~ L I W ~

@7I0 'ARRAY' MI Pr B r Rfl t NJr ZKtI : N + 21,'

0715 'BOOLEt%J' 'ARRAY' ICCl : h3

-

13;

j 0720 'BOOLEAN' 'PROCEDIBE' Q R T ~ O ~ ~ ~ ~ E ~

i3725 'BMjIN'

'

IWTEEER' I 0 KK; 'REAL' NZr VZ0r NZ€MZ

07318 NZB0 := N% t = SQRT(IhJPROD(1r

1~

ND ZKCIJ, ZKEII

0735 ORTHONOAMALIZE := eTRüE':

@?YB AF K KI 'I"M33'

m4s

'BEGIN' 'ARRp;Y' GtKl 2 K

-

1 3 3

@75@ 'FOR' IM := KI 'STEP'

1

'UNTIL' K

-

1

"W'

0755 GCKKJ := -INPROD(I l r M D ZKEIIr ZCIr K K I s 0):

0760 'FOR' I := i 'STEP'

1

'IEVTIL' N 'Co'

e 9 e

0765 ZKEI3 := INPRODCKKr Kir K I r ZEI, KKJr G C K K J r & K t I J ) ;

0770 NZO := %Zi NZ := SQRT(INmEOD<Ir I r Nr ZKtIIr ZKcf3r G#>>t

0775

*rv

NZ <=

0.001 *

y200 ~TIIEIV*

0780 'BEGIN' ORTXBIORMALIZE := 'FALSE'; 'GOTO'

READY

*E?JD';

0785 'IF' NZ C= 6.5

*

NZO 'THEN' 'GOTO' AGAIN

0790 'END';

(11)

0800 READY: e805 'END'; 0810 'REAL' 'PROCEWRE' R4VlXBíf 0815 *BEEIN' FFGXVD := FRPND

*

1021;

0820 := FRAVD := FR14'JD

-

ENTIER(FR&VDI 0825 'END'; 0830 MOW :=

2

t (-332): FRAND := @ e 5 + 1/18445763 @335 'FOR' I :=

1

'STEP'

1

'UIVTIL' N 'W'

@840 'BEGIN' U := ABSCDIAGCI.1) +

0845 ('IF' I >

1 'W'

F1RSCCODIAGCI

-

1 3 ) 'ELSE'

a )

+

0850 ('IF' I < N 'THEN' ARS<CODIAGtIJ) 'ELSE'

0);

0855 'IF' U r NORM 'THEN* NOW4 := U

9860 *END';

0865 EPS I= SORM

*

MACHEPC; LK := NORPI + EPSS

0870 QCNI := REN3 := ZKCN 4 1 3 := ZKCN + 21 := 0:

0975 'FOR' K := K l 'STEP'

1

*WTIL' K2 'Do'

@?80 'BEZE!' F := 0; LK := LK

-

FPSf i 08R5 @895 @90@ 09@5 891 0 091 5 0920 0925 093@ @935 0948 0945 095@ 0955 a960 Qi965 09763 , (21975 @98(21 0985 a980 SET: 0995

100m

mm:

1005

l a l 0 l a 1 5

1

B20

1025

1030

1035

1040 1045

1050

2 055 lB6D

1

@65 i

'IF' LK > LAMRIXtKl 'THDJ' LK := L74MI3MCYIi

U := D I A G t l l

-

LK;

v

:= COD11 :=

'IF' ARS(CODIAGt.13) .c EP'S 'THEN' EPS 'ELSE' CODIAGtlJs 'FOR' I :=

1

'STEP'

1

'WTIL' N

-

1 'MI'

N' COD1 := COD115 COD11 := IF* I = N

-

1

'THEN' 0

'ELSE' IF' ABSGODIAGCI +

11 >

e EPS 'THEN' WS

'IF' ARS(CODI> <= ARS(U3 'TFEX'

'REGf3' ICE13 := 'FALSE'i X I 3 := U; QCIJ := V; REI3 := Sl

'ELSE' CODIAGC I

+

1

l i

M t I l := CODI/U; U := DIAGCI +

11 -

LK

-

MEI3

*

t J

V := COD11

'Ewe 'ELSE'

'REGIg' ICE13 := 'TRUE'; PCIJ := COD11

QLII := DI4GCI + 1 3

-

LK; RCIJ := CODXI; PIC11 :=

'IF' U =

0

'AVD' AQSCCODI) e= EPS ' T W '

1

U := L'

-

M t I l

*

QICIIi V t= -Méf3

*

RF-I3

*ELSE*

u/com

i

'END' 'END';

PtVl := 'IF' U =

0

'TFEN' EPS eELSE' Ui

'FOR' I :=

1

'STEP' i 'LNTIL' 3 'Do' ZKCII := RAVDOM; FIRST ITEXATIOX := 'TRUE';

'FOR' I := N 'STEP'

-1

'UC\ITIL'

1

'DO' ZKEII :=

(ZKEII

-

Q L I I

*

ZKCI +

11 -

RCI3

*

ZKCI + 23>/PCfl;

'

IF' 'NOT' ORmûRMALIZE

'

T%EhJ'

'BEGI3J' 'IF' F =c 5 'THEN' 'BEGIV' F := F

+

I ; 'GOTQ' SET 'END'

'ELSE' 'IF' 'NOT' FIRST ITERATION ' T W '

*FOR' I :=

1

'STEP'

1

'WTIL' N 'D3' ZWEI3 := 0 EiJD' i

'IF' FIRST ITERATION 'THEN'

'REGIN' 'FOR' I :=

1

'STEP'

1

'UNTIL' N

-

1

'Do'

'REGIN' 'IF' ICCI3 'THEN'

'BEGIN' 0 := ZKEIJ; ZKCII := ZKCI + 13:

'

mD' t

ZKCI + 1 3 := Z K f I +

11 -

MCIJ

*

ZKEfl ZKCI + 1 3 := U

(12)

1070 'END' f

1075 FIRST ITE24TION := 'FALSE'; 'GOTO' BACK

1080

'

EVD'

:

1085 'FOR' I 1 'STEP' 1 'LNTIL' EV 'El' ZCI, KI := ZF(Ci3

1090 'FSD' 'END';

1095 'PROCECiURE* BACKT&XVSFORMATION@YN, Hs Kl, KCfr E); 1100 'VALUE' NI Klr K21 'INTEGER' NI Kl, K2; 'ARRAY' Ft S i

1105 'REGIN' 'INTEGEX' 11 KJ KKf 'REK' S3

1 1 1@ 'FOR' KK := K1 'STEP' 1 'DJTIL' K2 'Do' 1115 'FOR' K := N

-

2 'STEP' -1 'UVTIL' 1 'El'

1120 'IF' WCKJ KI 'NB' 0 'TIIEN'

1 I25 'BEGIN' S := FNPRODCIr K + 1, N a HtI, KI, ZCI, K K l r 6331

1130

alt

KI1

1135 'FOR' I := K + 1 TEP' 1 'WTIL' N 'W'

1140 Z t f ,

ml

:= ZCI, 3

-

s

*

m b

KI

1145 'END'

I, 11, r 2 ~ WDW AI, RI,

c;

i 'IJIJTILQ 12

*mg

c

:=

c

+ AI

*

Rr;

1180 'PROCEWRE' CHOLBD DECOMPOSITI(SN WM,& FAIL);

1185 'VAtuE' NJM; 'lNVTEGER' NIM; ' ~ 3 " Af 'LABEL' FAfLf

ll9# '3EGI1v' 'INTEGER' áClM1,IsJf 'REAt' Xis

1 i95 M i :=MI

1200 'FüR' K:=l 'STEP' 1 'WTIL' N

'Do'

1205 eBBEGIN"IF' AtK,QI<@ 'OR' ACK,al=0 '"WN"GOT0' FAILS

1210 ACKsG3 :~SQRT(ACKJ@~

>;

1215 OIF9 WN-K *"pw\fw MirM-KS

122@ @FORe f:=l 'STEP' 1 'LEVTIL' M1 'Do'

i225

123@

1240 'FOR' J:=l 'STEF)' 1 'UNTIL' M1 'IXI'

AC% Ilt=ACKo IIbACKoBI;

Arw *-.-it

i235 rlLn# $ 4 * - a s

1845 'BEGDJ' X:=AEKnJ3i 'IF' X *NQv @ 'THEX'

1255 Af K+JJ 1 3 Z SACK+ JI 1 I-X*At Ka

'STZF'

1

'"TIL@

'ma

'-n. _run-i i = f g i i I

6, 1850 'FOR' I:=@ 'STEP' 1 'WTIL' M1-3 'XI*

1 3

1268 'END'

1265

'

END'

i27Q

'mm'

DEcmmJERaJ VAf Dl2

C w m o I l x

WGS. C?-!QLSKf*

1275 GEBRUIK MRF[END VAN BAYCMATRICE% 1260

12135

129B CFBIL, FAIL, VERDER 5;

1295 'VALUE' N,MaKl,KSP

'

INTmEX' N,Mif(l,K21 'ARRAY' A J Q I L M D A J X ; 130@ 'F#)OLEANs EIVEGCROL;

'LAREL'

FAIL,VMDERP

1305 'BEEIN' 'INTEGER' I,J,MlrMFbPaRf 'REAL' St 1310 'ARRAY' DIACL CODIAGs XlCltNI;

1315

'

IF' CML 'THEN' CHDLRD DECOMFOSITION (NIM~QIFAIL);

1320 MR:=1;

1325 'FOR' R:=l 'STEP' 1 'INTIL' N 'DD' 1330 'BEGIN' 'IF' f B M 'THEN' MR:=R-M;

1335 'FOR' &=MR OSSTEP' 1 'C1NTIL' N 'Do' 'PROCEDURE' BDEIG <N,MJ AJ QI LAMRDAJ Xs K10 K2, E I W *

(13)

1340 1345 1350 1355 1360 1365 137@ 1375 1380 1385 1390 1395 1400 14@5 1410 1415 1420 1425 143@ 1435 f44a 1445 l45@ 1455 146@ 1405 1470 I475 1480 1485 1490 1495

1500

1505 1510 1515

1520

1525

1535 1540 1545 1558 1555 156@ 1530 X~CJJ:=A~JIRJ;

'FOR' I:=MR 'STEP' I 'UNTIL' N 'Do'

'I3EGI.J' Ml:=*IF' I>M 'THEN' I-M 'ELSE'

1;

P:='IF' Ml>MR 'THEN' M l 'ELSE' MR;

S:=XlCI1;

'FOR' J:=P 'LiTFP'

1

'UVTIL'

1-1

*XI'

X1 t IJ t _SACI, RI :=S/QC Ia03f

S:=S-QE 3, I-J3

*x1

r:

J J ;

'END'

'END';

'FORe R:=l 'STEP' 1 'WTIL' N 'Do'

'BEGXN' 'FOR'

&=1

'STEP'

1

'WTIL' R 'D3'

'FOR' I:=l 'STEP' i 'WTIL' i? ' W e

'BEX5IN' Ml:=eIIF' P M 'THEN' I-M 'ELSE'

1;

XlCS3:sACR,JI; s : = X l t I l ;

'FOR' J:=MI 'STEP'

1

'WTIL'

1-1

'IM'

St=S-GC 4, P-JI*Xl C 53;

XI E I I 3 =At R a I 1 8 =At Ir R3 :=S/QL

11

03 f

(14)

17

(15)

(16)

i

\

(17)

(18)

(19)

i , +2 +m8@57Pl9*+ B -*533402'+ # -*754344'+ 0 -*533402'+ 0 +* 833B€30'* 8 +m5334@2'+ 0 +.754344*+ $3

(20)

+3 +014@313'+ 1 +o458043'+ 0 +o237100'+ 0 -0335311'* # -o410670'+ 0 +0122732'+ 0 +o474281 * + 0 +4 +oS08274*+ 1 -0351469'+ El +0181933'+ 0 +0257293'+ @ -0315118'+ 0 -0941757'- 1 +0363867*+ 0 +5 +.281465'+ 1 +0220942'+ @ -031246&1'+ @ +022#942'+ @ -019992D'- 7 -e220942'+ 5 +o312460'+ 5 i +6 +033774@'4 1 +.753593'- 1 -c145583'+ 0 +0205885'* @ -0252157'4 0 *0281245'+ 0 -0291166'+ &I

WILT U NOG MEER C E V U N DOORRMENEEJ TIK IXN 1s ZO NIET TIK M U 0