[Boost-bugs] [Boost C++ Libraries] #7549: uBLAS banded storage does not match common BLAS layout (from netlib)

Subject: [Boost-bugs] [Boost C++ Libraries] #7549: uBLAS banded storage does not match common BLAS layout (from netlib)
From: Boost C++ Libraries (noreply_at_[hidden])
Date: 2012-10-21 20:38:12

#7549: uBLAS banded storage does not match common BLAS layout (from netlib)
----------------------------------------+-----------------------------------
 Reporter:  guwi17                      |       Owner:  guwi17 
     Type:  Bugs                        |      Status:  new    
Milestone:  To Be Determined            |   Component:  uBLAS  
  Version:  Boost 1.52.0                |    Severity:  Problem
 Keywords:  ublas banded matrix layout  |  
----------------------------------------+-----------------------------------
 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
 }}}

-- 
Ticket URL: <https://svn.boost.org/trac/boost/ticket/7549>
Boost C++ Libraries <http://www.boost.org/>
Boost provides free peer-reviewed portable C++ source libraries.