Mercurial > octave-antonio

diff scripts/linear-algebra/qzhess.m @ 11471:994e2a93a8e2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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	Sun Jan 09 13:44:15 2011 -0800
+++ b/scripts/linear-algebra/qzhess.m	Sun Jan 09 16:01:05 2011 -0800
@@ -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)