# HG changeset patch # User jwe # Date 1068826541 0 # Node ID ac4e4807acc5f1586a7e65c7fa6c28adcfb768f4 # Parent 22ca4cc02525a5e5257616694b125c799080e0cb [project @ 2003-11-14 16:14:31 by jwe] diff -r 22ca4cc02525 -r ac4e4807acc5 scripts/ChangeLog --- a/scripts/ChangeLog Fri Nov 14 03:55:04 2003 +0000 +++ b/scripts/ChangeLog Fri Nov 14 16:15:41 2003 +0000 @@ -1,3 +1,9 @@ +2003-11-14 Gabriele Pannocchia + + * 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 * general/reshape: Delete. diff -r 22ca4cc02525 -r ac4e4807acc5 scripts/linear-algebra/krylov.m --- a/scripts/linear-algebra/krylov.m Fri Nov 14 03:55:04 2003 +0000 +++ b/scripts/linear-algebra/krylov.m Fri Nov 14 16:15:41 2003 +0000 @@ -43,7 +43,7 @@ ## Author: A. Scottedward Hodel -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 diff -r 22ca4cc02525 -r ac4e4807acc5 scripts/linear-algebra/krylovb.m --- a/scripts/linear-algebra/krylovb.m Fri Nov 14 03:55:04 2003 +0000 +++ b/scripts/linear-algebra/krylovb.m Fri Nov 14 16:15:41 2003 +0000 @@ -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