Re: Generalised Eigenvalue of complex matrices - Help!!

[Paul Kienzle]

I wrote some docs on this in the wiki:

        http://wiki.octave.org/wiki.pl?OctaveFortran

In the declaration you need two F77_CHAR_ARG_LEN_DECL at the end (no 
commas around them).

In the call you need two F77_CHAR_ARG_LEN(n) specs (no commas around 
them).

See also: octave/src/DLD-FUNCTIONS/qz.cc for an example.

Are you sure work_length is a Complex in your first call?

Use F77_XFCN if you are going to test for exceptions.  You don't need 
to test for exceptions if you don't use f77_XFCN.

The wiki page doesn't mention f77_exception_encountered, except in the 
example.  Note that qz doesn't always use f77_exception_encountered, 
and I'm not sure why.  From a brief glance at the code in 
octave/libcruft/misc/f77-fcn.h, it looks like every call to F77_XFCN 
should be be followed by an f77_exception_encountered test.  Please 
update the wiki once you find out.

- Paul

On Jul 19, 2005, at 5:24 PM, Stefan van der Walt wrote:

> I put together this early draft version of of 'gev' for calculating
> generalized eigenvalues.  I have never interfaced Fortran with Octave
> before, so I'd appreciate any tips on how to improve it.



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------