[ublas] Plans for bindings/lapack/computational/gbsv.hpp?

Subject: [ublas] Plans for bindings/lapack/computational/gbsv.hpp?
From: Larry Evans (cppljevans_at_[hidden])
Date: 2011-11-18 13:29:59

Currently, there's no gbsv.hpp in that directory.
Are there any plans to add one?

There is a test driver with a name that suggests the use of gbsv;

http://svn.boost.org/svn/boost/sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_gbsv.cpp

however, it only uses:
  lapack::gbtrf
  lapack::gbtrs
and when MatrB is a ublas::matrix instead of ublas::vector, the runtime
assertion:

boost/numeric/bindings/lapack/computational/gbtrs.hpp:158: static
std::ptrdiff_t
boost::numeric::bindings::lapack::gbtrs_impl<Value>::invoke(const
MatrixAB&, const VectorIPIV&, MatrixB&) [with MatrixAB =
boost::numeric::ublas::banded_matrix<double,
boost::numeric::ublas::basic_row_major<>,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> >
>; VectorIPIV = std::vector<int>; MatrixB =
boost::numeric::ublas::matrix<double,
boost::numeric::ublas::basic_row_major<>,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> >
>; Value = double; std::ptrdiff_t = long int]: Assertion
`(bindings::bandwidth_upper_op(ab, trans())-
bindings::bandwidth_lower_op(ab, trans())) >= 0' failed.

occurs.

TIA.

-Larry