使用LAPACK查找增强矩阵的特征值

Question

我无法使用boost LAPACK绑定找到增强矩阵的特征值。这是我最好的猜测：

template<class T>
boost::numeric::ublas::vector<double> diagonalize(const boost::numeric::ublas::matrix<T>& input)
{
    // create a working copy of the input
    boost::numeric::ublas::banded_matrix<T,boost::numeric::ublas::column_major> A(input);

    fortran_int_t n = input.size1();

    // Calculate tridiagonal form
    boost::numeric::ublas::vector<double> d(n);     //diagonal elements of tridiagonal matrix
    boost::numeric::ublas::vector<double> e(n-1);   //offdiagonal elements of tridiagonal matrix
    boost::numeric::ublas::vector<T> tau(n-1);  //other stuff
    boost::numeric::bindings::lapack::hetrd( A, d, e, tau );

    //Calculate eigenvalues
    boost::numeric::bindings::lapack::sterf(n,d,e);

    return d;
}

编译失败的地方：

boost/numeric/bindings/detail/property_map.hpp:30: error: no type named 'property_map' in 'struct 
 boost::numeric::bindings::detail::adaptor_access<boost::numeric::ublas::banded_matrix<std::complex<double>, 
boost::numeric::ublas::basic_column_major<long unsigned int, long int>,
boost::numeric::ublas::unbounded_array<std::complex<double>, 
std::allocator<std::complex<double> > > >, void> 

我知道我可能在这里搞乱矩阵类型，但编译器错误让我有点无助，我在网上找不到任何例子。