[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Generalised Eigenvalue of complex matrices - Help!!
From: |
Paul Kienzle |
Subject: |
Re: Generalised Eigenvalue of complex matrices - Help!! |
Date: |
Tue, 19 Jul 2005 21:26:06 -0400 |
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
-------------------------------------------------------------