LS,
I would really
appreciate it when someone could explain to me why the code snippet below does
not work; I would like to extract unit_lower or strict_lower triangular matrices
from a general matrix, but it seems I am doing something wrong. Any help would
be highly appreciated!
Kind
regards,
Pieter Ober
//
**********
matrix<double> M (3,
3);
triangular_matrix<double,lower>
L;
triangular_matrix<double,unit_lower> UL;
triangular_matrix<double,strict_lower> SL;
L =
triangular_adaptor<matrix<double>,
lower> (M);
UL =
triangular_adaptor<matrix<double>, unit_lower> (M); // throws a "boost::numeric::ublas::bad_index"
exception
SL =
triangular_adaptor<matrix<double>, strict_lower> (M); // throws a
"boost::numeric::ublas::bad_index"
exception