[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Re: gsl c++ wrapper ?
From: |
Leo Razoumov |
Subject: |
Re: [Help-gsl] Re: gsl c++ wrapper ? |
Date: |
Mon, 31 Aug 2009 18:44:08 -0400 |
On 8/31/09, Rhys Ulerich <address@hidden> wrote:
> > Designing a good C++ wrapper to any matrix
> > library is difficult provided one tries to keep C++ semantics.
>
>
> It may be easier, if you have a favorite C++ matrix/vector library, to
> design a GSL-wrapper for the matrix/vector library itself.
>
> For example, write wrappers that accept ublas [1] or eigen2 [2]
> matrices and vectors, wrap the same storage with gsl_matrix and
> gsl_vector views, and then call the GSL methods directly. Such
> GSL-ublas or GSL-eigen2 bindings would be independent of both
> projects, mainly boilerplate, and easy to test for correctness.
>
> - Rhys
Many good and tested libraries like LAPACK have come from Fortran and
use column-major matrix storage incompatible with GSL. One cannot
simply wrap them up in some sort of GSL objects that can be passed to
GSL functions.
--Leo--