Re: LAPACK vs. GSL matrix algebra

[C J Kenneth Tan -- OptimaNumerics <cjtan at OptimaNumerics dot com>]

Linas,

You don't have to remember how to code in Fortran to use LAPACK. 

BLAS and LAPACK are available for most platforms these days.  If you
use BLAS and LAPACK, you can always use the free versions from Netlib.
When you need speed, you can switch to optimized versions of BLAS and
LAPACK.  High performance implementations of these libraries can
deliver very significant speedup.  We have seen improvements of
hundreds to thousands of percents depending on the routine and
architecture, compared to hardware vendors' code!


Kenneth Tan
-----------------------------------------------------------------------
C. J. Kenneth Tan, Ph.D.
OptimaNumerics Ltd.
E-mail: cjtan@OptimaNumerics.com      Telephone: +44 798 941 7838    
Web: http://www.OptimaNumerics.com    Facsimile: +44 289 066 3015 
-----------------------------------------------------------------------
This e-mail (and any attachments) is confidential and privileged.  It
is intended only for the addressee(s) stated above.  If you are not an
addressee, please accept my apologies and please do not use,
disseminate, disclose, copy, publish or distribute information in this
e-mail nor take any action through knowledge of its contents: to do so
is strictly prohibited and may be unlawful.  Please inform me that
this e-mail has gone astray, and delete this e-mail from your system.
Thank you for your co-operation.
-----------------------------------------------------------------------




On 2004-09-27 11:20 -0500 Linas Vepstas (linas@austin.ibm.com) wrote:

> Date: Mon, 27 Sep 2004 11:20:18 -0500
> From: Linas Vepstas <linas@austin.ibm.com>
> To: gsl-discuss@sources.redhat.com
> Subject: LAPACK vs. GSL matrix algebra
> 
> 
> 
> Hi,
> 
> What is the current thinking with adding to the matrix functionality 
> in GSL?   I need to find the eigenvectors of a non-symmetric real
> matrix, which is something that GSL doesn't currently support.
> 
> My choices seem to be
> -- try to remember how to code in fortran
> -- try to implement something in C, with an eye towards submitting 
>    it to GSL.
> 
> Any suggestions, words of encouragement to go one way or the other?
> 
> 
> --linas
> 
> 
>