[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Solving Ax=B where A is a complex, singular matrix

From: Marco Maggi
Subject: Re: [Help-gsl] Solving Ax=B where A is a complex, singular matrix
Date: Mon, 04 Jul 2011 10:29:26 +0200

Brian Gough wrote:

> Thanks for your email.  I  can confirm that the library is
> lacking the  complex SVD and QR  unfortunately.  If anyone
> can write either I will gladly include them.

In general, is  it maybe time to consider  to rely on LAPACK
now that they have made these changes to the C API?  I still
have  to try  the  latest  LAPACK, but  from  the README  of

  This implementation introduces:

  - row-major and  column-major matrix layout  controlled by
    the first function parameter;

  - an  implementation  with  working  arrays  (middle-level
    interface) as well as without working arrays (high-level

  - input scalars passed by value;

  - error  code  as  a  return  value instead  of  the  INFO

  I do  not think GSL will  ever be able to  provide as many
routines as  LAPACK does;  so if "officially"  GSL documents
how  to  use its  data  types  with  LAPACK's routines,  GSL
maintainers and users can stop worrying about missing things
(we have seen multiple  requests for complex matrix routines
on the GSL mailing list in recent years).

  Such  documentation work  should  be much  less work  than
implementing new routines for GSL.
Marco Maggi

reply via email to

[Prev in Thread] Current Thread [Next in Thread]