[boost] ublas matrix range problem

From: boost (boost_at_[hidden])
Date: 2002-04-05 09:52:58

Salut,

currently I'm fighting with the followinf ublas problem:

typedef numerics::matrix<double, numerics::row_major,
    numerics::unbounded_array<double> > TpMatrix;
typedef numerics::matrix_range< TpMatrix > TpMatrixRange;
    
 TpMatrix U( 1, 2 );
 TpMatrixRange r1  = TpMatrixRange( U, numerics::range( 0,1),
    numerics::range( 0,1) ); // o.k.
       
TpMatrixRange r2 ;
r2 = TpMatrixRange( U, numerics::range( 0,1), numerics::range( 0,1) ); // ??? 

Works for r1, but not for r2 on my Debian/Linux g++2.95.4:

t: /home/peter/DMRG/boost/numeric/ublas/storage.h:41: const size_t & 
numerics::common_impl_ex<size_t>(const size_t &, const size_t &, const char 
*, int): Assertion `size1 == size2' failed.

Is this intended ?

Best wishes,
Peter