Rutger/Karl:
Again, sorry for taking so long to test. Travel and a computer rebuild.
> This is right, trtri is meant for normal dense matrices (=storage format) of
> which you say they are triangular. I've added support for uplo detection in
> trtri, to use it now please use ublas' triangular adaptor to indicate
>whether it is upper/lower.
Thanks for reminding me about triangular storage formats. I think I get it now, how to use the triangular adaptor if I have dense storage.
I tried the following code:
ublas::triangular_matrix<double, ublas::lower, ublas::column_major> A(2,2);
lapack::tptri('N', A);
And got a bunch of missing traits:
1>C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/tptri.hpp(75): error: class "boost::numeric::bindings::traits::matrix_traits<boost::numeric::ublas::triangular_matrix<double, boost::numeric::ublas::lower, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<QuantLib::Real={double}>>>>" has no member "value_type"
1> typedef typename traits::matrix_traits< MatrixAP >::value_type value_type;
1>C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/tptri.hpp(56): error: class "boost::numeric::bindings::traits::type_traits<<error-type>>" has no member "real_type"
1> typedef typename traits::type_traits<ValueType>::real_type real_type;
1>C:\working\libraries\boost\boost/numeric/bindings/traits/matrix_traits.hpp(139): error: class "boost::numeric::bindings::traits::matrix_traits<boost::numeric::ublas::triangular_matrix<double, boost::numeric::ublas::lower, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<QuantLib::Real={double}>>>>" has no member "uplo_type"
1> typedef typename matrix_traits<SymmM>::uplo_type uplo_t;
It sure would be nice if when the higher level interfaces are developed they would dispatch to the correct LAPACK routine based on the packed vs. dense storage.
Thanks,
Jesse