From: Dima Sorkin (dsorkin_at_[hidden])
Date: 2006-01-29 02:34:08
Quoting Karl Meerbergen :
> All I can say is that the bindings are written for column major matrices.
> Your matrix is row_major. I do not believe it is sufficient to transpose it,
> since for a row_major matrix the leading dimension corresponds to the number
> of columns. Lapack assumes it corresponds to the number of rows.
> You could pass trans(M) as an argument to lapack. trans(M) creates an object
> of type matrix_unary2<M,scalar_identity<typename M::value_type> >. There are
> no bindings for this kind of object. These should be written.
> This would be the clean way of doing it.
Alternatively, I think that a column-major temporary can be created.