EISPACK solution of generalized eigenvalue prob.

Thread View: sci.math.num-analysis

1 messages

[T] Back to Threads [V] Tree View [G] Group Articles

1 total messages Started by Matt Wette Fri, 21 Dec 1990 22:06

#578

Author: Matt Wette
Date: Fri, 21 Dec 1990 22:06

25 lines
1499 bytes

This article concerns the problem of finding two orthogonal
matrices Q and Z such that 
	Q'*A*Z is in upper real Schur form
	Q'*B*Z is upper triangular
where A and B are square matrices and ' denotes transpose.

I have encountered problems with the EISPACK routine QZIT
not converging.  The problem I am trying to solve is of
order 18 and a `relative Crawford number' for the problem,
given by evaluation of the Matlab expression
	cond(A/norm(A) + j*B/norm(B)),
is about 10^2.

Does anyone know of software which is a better implementation 
of (a variant of) the QZ algorithm and/or a better working
algorithm for solving the same problem?

Matt


-- 
 _________________________________________________________________
 Matthew R. Wette           | Jet Propulsion Laboratory, 198-326
 mwette@csi.jpl.nasa.gov    | 4800 Oak Grove Dr, Pasadena,CA 91109
 -----------------------------------------------------------------

Thread Navigation

This is a paginated view of messages in the thread with full content displayed inline.

Messages are displayed in chronological order, with the original post highlighted in green.

Use pagination controls to navigate through all messages in large threads.

Back to All Threads