Mercurial > hg > octave-lyh
diff scripts/linear-algebra/qzhess.m @ 11471:994e2a93a8e2
Use uppercase 'A' to refer to matrix inputs in m-files.
author | Rik <octave@nomad.inbox5.com> |
---|---|
date | Sun, 09 Jan 2011 16:01:05 -0800 |
parents | 693e22af08ae |
children | fd0a3ac60b0e |
line wrap: on
line diff
--- a/scripts/linear-algebra/qzhess.m +++ b/scripts/linear-algebra/qzhess.m @@ -18,11 +18,11 @@ ## <http://www.gnu.org/licenses/>. ## -*- texinfo -*- -## @deftypefn {Function File} {[@var{aa}, @var{bb}, @var{q}, @var{z}] =} qzhess (@var{a}, @var{b}) +## @deftypefn {Function File} {[@var{aa}, @var{bb}, @var{q}, @var{z}] =} qzhess (@var{A}, @var{B}) ## Compute the Hessenberg-triangular decomposition of the matrix pencil -## @code{(@var{a}, @var{b})}, returning -## @code{@var{aa} = @var{q} * @var{a} * @var{z}}, -## @code{@var{bb} = @var{q} * @var{b} * @var{z}}, with @var{q} and @var{z} +## @code{(@var{A}, @var{B})}, returning +## @code{@var{aa} = @var{q} * @var{A} * @var{z}}, +## @code{@var{bb} = @var{q} * @var{B} * @var{z}}, with @var{q} and @var{z} ## orthogonal. For example: ## ## @example @@ -46,22 +46,22 @@ ## Created: August 1993 ## Adapted-By: jwe -function [aa, bb, q, z] = qzhess (a, b) +function [aa, bb, q, z] = qzhess (A, B) if (nargin != 2) print_usage (); endif - [na, ma] = size (a); - [nb, mb] = size (b); + [na, ma] = size (A); + [nb, mb] = size (B); if (na != ma || na != nb || nb != mb) error ("qzhess: incompatible dimensions"); endif ## Reduce to hessenberg-triangular form. - [q, bb] = qr (b); - aa = q' * a; + [q, bb] = qr (B); + aa = q' * A; q = q'; z = eye (na); for j = 1:(na-2)