I'm currently trying to use the uBlas library with the units library and I'm running into some problems.  The matrix and vector classes I've instantiated with unit constraints don't seem to play well with uBlas methods.  Has anyone found ways of using the two libraries together?  I've included some lines of what I'm trying to do and would really appreciate some help.

 

Ryan

 

using namespace boost::numeric::ublas;

using namespace boost::units;

 

//Works fine.  This is expected though since dimensionless quantity can accept

//regular built in types.

bounded_matrix<quantity<si::dimensionless, 3, 3>  rotate = identity_matrix<double>(3);

 

//How do I initialize a vector of dimension lengths with something like the zero_matrix?

bounded_vector<quantity<si::length>, 3> diff = zero_matrix<double>(3) * si::meter; //???

 

//Lets pretend the diff variable has values in it.  This line doesn't compile.

//The compiler is saying it can't deduce the template argument.

//What can I do to help the compiler out so I can get the correct result?

bounded_vector<quantity<si::length>, 3> result = prod(rotate, diff);