Boost logo

Ublas :

From: Pieter (PieterB_at_[hidden])
Date: 2007-09-25 03:20:17


Using boost version 1.34.1 and Visual C++ 2005, the following code
throws a bad_index exception when used with unit_lower triangular
matrices, while the same works fine for regular triangular matrices:

matrix<double> A(3,3);
triangular_matrix<double,unit_lower> UL(3,3);
// fill A
UL = getUnitLower(A);


template<class T>
getUnitLower(matrix<double> M) const
        triangular_adaptor<const matrix<T>, unit_lower> UL(M);
        return UL;

>From the documentation, I concluded that this should be legal to do, but
it seems the code attempts to assign the diagonal of UL. Is anyone aware
of this problem?

Many thanks,