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