# HG changeset patch # User jwe # Date 1068826541 0 # Node ID ac4e4807acc5f1586a7e65c7fa6c28adcfb768f4 # Parent 22ca4cc02525a5e5257616694b125c799080e0cb [project @ 2003-11-14 16:14:31 by jwe] diff --git a/scripts/ChangeLog b/scripts/ChangeLog --- a/scripts/ChangeLog +++ b/scripts/ChangeLog @@ -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 --git a/scripts/linear-algebra/krylov.m b/scripts/linear-algebra/krylov.m --- a/scripts/linear-algebra/krylov.m +++ b/scripts/linear-algebra/krylov.m @@ -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 --git a/scripts/linear-algebra/krylovb.m b/scripts/linear-algebra/krylovb.m --- 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