[ublas] LAPACK/BLAS bindings

Subject: [ublas] LAPACK/BLAS bindings
From: Peter Man (plwm2_at_[hidden])
Date: 2010-06-06 10:05:19

Dear all,

Would anyone know what the problem with the following code is?

///////////////////////////////////////////////////////////////////////////////////////////////////////////
Code 1 - inside a function called "void UTILS::LinAlgBasicTest()"
///////////////////////////////////////////////////////////////////////////////////////////////////////////
namespace ublas = boost::numeric::ublas;
namespace blas = boost::numeric::bindings::blas;
namespace lapack = boost::numeric::bindings::lapack;

ublas::vector<double> uv (18);
for(unsigned int i=0; i<uv.size(); ++i)
{
uv(i) = i;
}
std::cout << blas::nrm2(uv) << std::endl;
ublas::vector_range<ublas::vector<double> > ur (uv, ublas::range (3, 15));
std::cout << blas::nrm2(ur) << std::endl;

This produces errors (in cygwin's g++):
/cygdrive/c/Program 
Files/boost/boost_1_42/boost/numeric/bindings/detail/property_map.hpp: 
In instantiation of 
‘boost::numeric::bindings::detail::property_at<boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double, 
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > 
 > >, boost::numeric::bindings::tag::value_type>’:
/cygdrive/c/Program 
Files/boost/boost_1_42/boost/numeric/bindings/value_type.hpp:20: 
instantiated from 
‘boost::numeric::bindings::value_type<boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double, 
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > 
 > > >’
./source/utils_linalgtest.cpp:224: instantiated from here
/cygdrive/c/Program 
Files/boost/boost_1_42/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::vector_range<boost::numeric::ublas::vector<double, 
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > 
 > >, void>’
./source/utils_linalgtest.cpp: In function ‘void UTILS::LinAlgBasicTest()’:
./source/utils_linalgtest.cpp:224: error: no matching function for call 
to 
‘nrm2(boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double, 
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > 
 > >&)’

or (in Microsoft Visual Studio 2010):
c:\program 
files\boost\boost_1_42\boost\numeric\bindings\detail\property_map.hpp(30): 
error C2039: 'property_map' : is not a member of 
'boost::numeric::bindings::detail::adaptor_access<T>'
1> with
1> [
1> 
T=boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double,boost::numeric::ublas::unbounded_array<double>>>
1> ]
1> c:\program 
files\boost\boost_1_42\boost\numeric\bindings\value_type.hpp(20) : see 
reference to class template instantiation 
'boost::numeric::bindings::detail::property_at<T,Key>' being compiled
1> with
1> [
1> 
T=boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double,boost::numeric::ublas::unbounded_array<double>>>,
1> Key=boost::numeric::bindings::tag::value_type
1> ]
1> c:\users\george\desktop\wholeproject 
new\utils\source\utils_linalgtest.cpp(224) : see reference to class 
template instantiation 'boost::numeric::bindings::value_type<T>' being 
compiled
1> with
1> [
1> 
T=boost::numeric::ublas::vector_range<boost::numeric::ublas::vector<double,boost::numeric::ublas::unbounded_array<double>>>
1> ]
1>c:\program 
files\boost\boost_1_42\boost\numeric\bindings\detail\property_map.hpp(30): 
error C2146: syntax error : missing ',' before identifier 'property_map'
1>c:\program 
files\boost\boost_1_42\boost\numeric\bindings\detail\property_map.hpp(30): 
error C2065: 'property_map' : undeclared identifier
1>c:\program 
files\boost\boost_1_42\boost\numeric\bindings\detail\property_map.hpp(30): 
error C2955: 'boost::mpl::at' : use of class template requires template 
argument list
1> c:\program files\boost\boost_1_42\boost\mpl\at.hpp(32) : see 
declaration of 'boost::mpl::at'
1>c:\program 
files\boost\boost_1_42\boost\numeric\bindings\detail\property_map.hpp(30): 
error C2039: 'type' : is not a member of 'boost::mpl::at'
1> c:\program files\boost\boost_1_42\boost\mpl\at.hpp(32) : see 
declaration of 'boost::mpl::at'
1>c:\users\george\desktop\wholeproject 
new\utils\source\utils_linalgtest.cpp(224): error C2893: Failed to 
specialize function template 
'nrm2_impl<bindings::value_type<T>::type>::result_type 
boost::numeric::bindings::blas::nrm2(const VectorX &)'
1> With the following template arguments:
1> 'boost::numeric::ublas::vector_range<V>'
1> with
1> [
1> 
V=boost::numeric::ublas::vector<double,boost::numeric::ublas::unbounded_array<double>>
1> ]

Also, when trying the ublas_getri.cc file in libs/numeric/bindings/atlas 
folder, it runs, but produces the following runtime error (in cygwin g++):

assertion "bindings::size(work.select(real_type())) >= min_size_work( 
bindings::size_column(a) )" failed: file "/cygdrive/c/Program 
Files/boost/boost_1_42/boost/numeric/bindings/lapack/computational/getri.hpp", 
line 202, function: static ptrdiff_t 
boost::numeric::bindings::lapack::getri_impl<Value, typename 
boost::enable_if<boost::numeric::bindings::is_real<T>, 
void>::type>::invoke(MatrixA&, const VectorIPIV&, 
boost::numeric::bindings::lapack::detail::workspace1<WORK>) [with 
MatrixA = boost::numeric::ublas::matrix<double, 
boost::numeric::ublas::basic_column_major<unsigned int, int>, 
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > 
 >, VectorIPIV = std::vector<int, std::allocator<int> >, WORK = 
boost::numeric::bindings::detail::array<double>, Value = double]
Aborted (core dumped)

Thank you very much in advance, and I apologise for the long email!!

Peter