Re: [Help-gsl] Re: gsl c++ wrapper ?

[Rhys Ulerich]

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

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

Quite true.  Underlying my statement was that your favorite C++
matrix/vector library can store data in a raw row-major format
compatible with GSL methods. Either ublas or eigen2 is capable of
row-major storage.

If you're after LAPACK, either ublas or eigen2 is capable of
column-major storage as well and probably could interop with LAPACK.

Hell, you might even use a mixture of storage types and get some
actual work done.

- Rhys