[ublas] bindings for gbsv

Subject: [ublas] bindings for gbsv
From: Karl Meerbergen (karl.meerbergen_at_[hidden])
Date: 2009-04-09 07:10:52

Hi,

The bindings for gbtrf state that the banded matrix should be row_major. 
This is surprising since all LAPACK routines should take column major 
matrices (with implicitly assumed transposition row_major can be allowed).

I suppose the reason for this is that the ublas banded matrices are not 
compatible with the lapack storage format. Ublas matrices should be 
passed in row_major format to be interpreted as column major by lapack. 
The traits specialization for ublas should thus take a row_major matrix, 
but interpret it as a column_major matrix. This would solve this strange 
situation.

Any thoughts?

Karl