| |
| | [No title] (Site not responding. Last check: 2007-10-08) |
 | | > # Rayleigh Quotient Iteration w/ deflation > # of complex mxm C, > # to get _ALL_ eigenvalues > # followed by inverse iteration > # to get all eigenvectors. |
 | | > M:=evalf(evalm(C)): > printf(`\n RAYLEIGH QUOTIENT ITERATION w/ DEFLATION\n`); RAYLEIGH QUOTIENT ITERATION w/ DEFLATION > n:=m; n := 5 > for k from 1 to m-1 do > Idn:=diag('1.0' $ ('i'=1..n)); > xv:=randvector(n): > yv:=randvector(n): > zv:=matadd(xv,yv,1,I): > dis:=infinity; > test:=0.; > for it from 1 to 25 > while (dis>1.0 or test |
 | | > X:=evalf(augment(uv,Idn)); > QRdecomp(X,Q='U'); > U:=evalf(evalm(U)); > M:=multiply(htranspose(U),M,U); > M:=submatrix(M,2..n,2..n); > n:=n-1; > od: lam=-43.5766450696602522140604+106.358451291029228444776*I, it=12, dis=3.54e-01, test=1.22e+00 lam=-108.067839327715693469085+23.8134396333860157517779*I, it=11, dis=2.88e-01, test=2.01e+00 lam=33.8351778994824143244665+108.487127124126530217014*I, it=9, dis=4.27e-01, test=3.89e+00 lam=134.899622987860321272466-91.0738602668532689202781*I, it=7, dis=2.00e-01, test=1.58e+01 > printf(`\n`); > lambda[m]:=M[1,1]; lambda[5] := 48.9096835100332100862128 - 188.585157781688505493290 I > printf(`\n INVERSE ITERATION. |
| www.math.okstate.edu /~burchar/ralqdfl.txt (162 words) |
|