Boost logo

Ublas :

From: Hugo Duncan (hugoduncan_at_[hidden])
Date: 2007-03-17 14:41:57


Hello Karl,

>>> Any chance of getting the leading dimension fixed in the bindings
>>> traits?
>>>
>>> http://www.nabble.com/Lapack-bindings-t3237071.html
>>> http://www.nabble.com/ublas-traits-bug-----leading_dimension-p9469619
>>> .html

> In my opinion, there is no error in the bindings traits classes.
> To be sure, can someone tell me what the difference is between size_m
> and size_M? The leading dimension is the number of rows of ublas::matrix
> for column_major matrices.

The nearest I can find to a canonical description of leading dimension
(LDA):
   http://www.netlib.org/lapack/lug/node116.html

Although this is more informative.
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=217&sid=79b90021ff88b6cc1d18db411b2f7ad4

As I am sure you are all too aware, and you state above, these all show
the leading dimension as the number of rows. So far I think we agree.

I find the use of m and M non-intuitive, but believe:

size_m = minor size
size_M = major size

For a column major matrix, I believe the major size is the number of
columns, so the leading dimension should be the minor size. The debugger
leads me to confirm these beliefs. But I get thoroughly confused by major
and minor sizes, so I would not be too categoric in stating this.

> Note that for row_major matrices, lapack bindings do not work.

Agreed.

Thanks for looking at this,

Hugo