Mercurial > octave-nkf

comparison scripts/control/zgscal.m @ 3240:2e74d8aa1a20

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1999-04-07 18:33:23 by jwe]
author jwe
date Wed, 07 Apr 1999 18:34:20 +0000
parents 98e15955107e
children 6dd06d525de6
comparison
equal deleted inserted replaced
3239:d54ef0aa6e82 3240:2e74d8aa1a20
47 F = eye(nmp); 47 F = eye(nmp);
48 for kk=1:nmp 48 for kk=1:nmp
49 F(1:nmp,kk) = zgfmul(a,b,c,d,F(:,kk)); 49 F(1:nmp,kk) = zgfmul(a,b,c,d,F(:,kk));
50 endfor 50 endfor
51 51
52 [U,H,k1] = krylov(F,z,nmp,1e-12); 52 [U,H,k1] = krylov(F,z,nmp,1e-12,1);
53 if(!is_square(H)) 53 if(!is_square(H))
54 if(columns(H) != k1) 54 if(columns(H) != k1)
55 error("zgscal(tzero): k1=%d, columns(H)=%d",k1,columns(H)); 55 error("zgscal(tzero): k1=%d, columns(H)=%d",k1,columns(H));
56 elseif(rows(H) != k1+1) 56 elseif(rows(H) != k1+1)
57 error("zgscal: k1=%d, rows(H) = %d",k1,rows(H)); 57 error("zgscal: k1=%d, rows(H) = %d",k1,rows(H));