Re: use LAPACK routine for triangular systems?

[Evan Monroig]

On 11/25/05, Fredrik Lingvall <address@hidden> wrote:
> >In my algorithm I end with a triangular system to solve, and I know
> >that LAPACK has some specific routines for triangular systems (like
> >this one http://www.netlib.org/lapack/double/dtrtrs.f). Is there a way
> >to take advantage of these with octave?
> >
> I have made I few mex-files for some of the BLAS and LAPACK routines
> for that purpose. I have attached  an example  using the LAPACK POTRI
> routine. You can use the mex-tools in octave-forge to build the
> corresponding
> oct-files for octave (you need to link to your BLAS and LAPACK libs).

Thanks, I will try your method. If I understand well, the process to do that is:

1) get the LAPACK f-file
2) wrap up the file in the way you did to have a mex-file
3) use mex-tools to convert that to an oct-file

Am I right?

Evan



-------------------------------------------------------------
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
-------------------------------------------------------------