Mercurial > octave-nkf

comparison scripts/linear-algebra/krylov.m @ 10791:3140cb7a05a1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add spellchecker scripts for Octave and run spellcheck of documentation interpreter/doccheck: New directory for spelling/grammar scripts. interpreter/doccheck/README: Instructions for using scripts. interpreter/doccheck/spellcheck: Script to spellcheck a Texinfo file. interpreter/doccheck/aspell.conf: GNU Aspell configuration file for Octave documentation. interpreter/doccheck/aspell-octave.en.pws: Private Aspell dictionary. interpreter/doccheck/add_to_aspell_dict: Script to add new Octave-specific words to private Aspell dictionary. interpreter/octave.texi: New @nospell macro which forces Aspell to ignore the word marked by the macro. interpreter/mk_doc_cache.m: Skip new @nospell macro when building doc_cache.
author Rik <octave@nomad.inbox5.com>
date Sat, 17 Jul 2010 19:53:01 -0700
parents d1978e7364ad
children eb9e0b597d61
comparison
equal deleted inserted replaced
10790:01f1643dfbb1 10791:3140cb7a05a1
36 ## 36 ##
37 ## If @var{v} is a vector and @var{k} is greater than 37 ## If @var{v} is a vector and @var{k} is greater than
38 ## @code{length(A)-1}, then @var{h} contains the Hessenberg matrix such 38 ## @code{length(A)-1}, then @var{h} contains the Hessenberg matrix such
39 ## that @code{a*u == u*h}. 39 ## that @code{a*u == u*h}.
40 ## 40 ##
41 ## The value of @var{nu} is the dimension of the span of the krylov 41 ## The value of @var{nu} is the dimension of the span of the Krylov
42 ## subspace (based on @var{eps1}). 42 ## subspace (based on @var{eps1}).
43 ## 43 ##
44 ## If @var{b} is a vector and @var{k} is greater than @var{m-1}, then 44 ## If @var{b} is a vector and @var{k} is greater than @var{m-1}, then
45 ## @var{h} contains the Hessenberg decomposition of @var{a}. 45 ## @var{h} contains the Hessenberg decomposition of @var{a}.
46 ## 46 ##
48 ## default value is 1e-12. 48 ## default value is 1e-12.
49 ## 49 ##
50 ## If the optional parameter @var{pflg} is nonzero, row pivoting is used 50 ## If the optional parameter @var{pflg} is nonzero, row pivoting is used
51 ## to improve numerical behavior. The default value is 0. 51 ## to improve numerical behavior. The default value is 0.
52 ## 52 ##
53 ## Reference: Hodel and Misra, "Partial Pivoting in the Computation of 53 ## Reference: A. Hodel, P. Misra, @cite{Partial Pivoting in the Computation of
54 ## Krylov Subspaces", to be submitted to Linear Algebra and its 54 ## Krylov Subspaces of Large Sparse Systems}, Proceedings of the 42nd IEEE
55 ## Applications 55 ## Conference on Decision and Control, December 2003.
56 ## @end deftypefn 56 ## @end deftypefn
57 57
58 ## Author: A. Scottedward Hodel <a.s.hodel@eng.auburn.edu> 58 ## Author: A. Scottedward Hodel <a.s.hodel@eng.auburn.edu>
59 59
60 function [Uret, H, nu] = krylov (A, V, k, eps1, pflg); 60 function [Uret, H, nu] = krylov (A, V, k, eps1, pflg);