id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
7549,uBLAS banded storage does not match common BLAS layout (from netlib),Gunter,Gunter,"A banded_matrix stores the elements in an unusual way. This makes it impossible to call standard BLAS-libraries with these matrix type.

Example data taken from http://www.netlib.org/lapack/lug/node124.html
{{{
Running banded_matrix < column_major > test
Full matrix
[5,5]((11,12,0,0,0),(21,22,23,0,0),(31,32,33,34,0),(0,42,43,44,45),(0,0,53,54,55))
data() of matrix
0 12 23 34 45 11 22 33 44 55 21 32 43 54 0 31 42 53 0 0
Expected data() of matrix
0 11 21 31 12 22 32 42 23 33 43 53 34 44 54 0 45 55 0 0
}}}

{{{
Running banded_matrix < row_major > test
Full matrix
[5,5]((11,12,0,0,0),(21,22,23,0,0),(31,32,33,34,0),(0,42,43,44,45),(0,0,53,54,55))
data() of matrix
0 11 21 31 12 22 32 42 23 33 43 53 34 44 54 0 45 55 0 0
Expected data() of matrix
0 0 11 12 0 21 22 23 31 32 33 34 42 43 44 45 53 54 55 0
}}}
",Bugs,new,To Be Determined,uBLAS,Boost 1.52.0,Problem,,ublas banded matrix layout,Rutger ter Borg <rutger@…>