Mercurial > hg > octave-nkf
changeset 4609:ac4e4807acc5
[project @ 2003-11-14 16:14:31 by jwe]
author | jwe |
---|---|
date | Fri, 14 Nov 2003 16:15:41 +0000 |
parents | 22ca4cc02525 |
children | 02d2fcf835fc |
files | scripts/ChangeLog scripts/linear-algebra/krylov.m scripts/linear-algebra/krylovb.m |
diffstat | 3 files changed, 9 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/ChangeLog +++ b/scripts/ChangeLog @@ -1,3 +1,9 @@ +2003-11-14 Gabriele Pannocchia <g.pannocchia@ing.unipi.it> + + * linear-algebra/krylov.m: Return H = [] in Vnrm == 0 case. + + * linear-algebra/krylovb.m: Fix typo in usage message. + 2003-10-29 John W. Eaton <jwe@bevo.che.wisc.edu> * general/reshape: Delete.
--- a/scripts/linear-algebra/krylov.m +++ b/scripts/linear-algebra/krylov.m @@ -43,7 +43,7 @@ ## Author: A. Scottedward Hodel <a.s.hodel@eng.auburn.edu> -function [Uret,H,nu] = krylov(A,V,k,eps1,pflg); +function [Uret, H, nu] = krylov (A, V, k, eps1, pflg); defeps = 1e-12; @@ -81,6 +81,7 @@ ## check for trivial solution if (Vnrm == 0) Uret = []; + H = []; nu = 0; return; endif
--- a/scripts/linear-algebra/krylovb.m +++ b/scripts/linear-algebra/krylovb.m @@ -31,7 +31,7 @@ case (5) [Uret, H, Ucols] = krylov (A, V, k, eps1, pflg); otherwise - usage ("[Uret, Ucols] = krylovb (A, V, k ,eps1, pflg)); + usage ("[Uret, Ucols] = krylovb (A, V, k ,eps1, pflg)"); endswitch endfunction