Using boostbindings with lapack

[Preben]

Hi,

don't know if this is the correct group, but try anyway.

I would like to use boostbindings from the boost sandbox together with
my project where I use boost::numeric: :ublas vectors and matrices.

I would guess that I could use these directly with the bindings in the
boost sandbox like:

----
#include <boost/numeric/bindings/lapack/gesv.hpp>
#include "MechComputer.h pp"

#include <iostream>

using namespace boost::numeric;
namespace lapack = boost::numeric: :bindings::lapa ck;
using namespace std;
namespace rw {
namespace dynamics {
MechComputer::M echComputer(con st MechSystem& msystem,
const ublas::vector<d ouble>& q0,
const ublas::vector<d ouble>& v0)
: _system(msystem ), _q(q0), _v(v0) {
}
const ublas::vector<d ouble>& MechComputer::s tep(double h) {
ublas::matrix<d oubledG = _system.DConstr aints(_q);
std::vector<ubl as::matrix<doub le ddG =
_system.DDConst raints(_q);
ublas::vector<d oubledV = _system.DV(_q);
ublas::vector<d oublem = _system.GetMass ();

ublas::matrix<d oublemM = M(dG, m);
ublas::vector<d oublevFT = mFmT(dG, ddG, dV, _v, m);

ublas::matrix<d oubleB(vFT.size (), 1);

lapack::gesv(M, B); // Result in B

return m; // Just return something
}
}
}
----

Compiling with this command-line:
----
preben@thinkgen ~/working/master/src $ icc -I
/home/preben/working/boost-sandbox -c MechComputer.cp p
MechComputer.cp p(31): error: no instance of overloaded function
"boost::numeric ::bindings::lap ack::gesv" matches the argument list
argument types are: (boost::numeric ::ublas::matrix <double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
(const boost::numeric: :ublas::matrix< double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
&, const boost::numeric: :ublas::vector< double,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
&), boost::numeric: :ublas::matrix< double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>)
lapack::gesv(M, B); // Result in B
^

compilation aborted for MechComputer.cp p (code 2)
preben@thinkgen ~/working/master/src $
----

I don't understand the error and why it's wrong.
g++ complains even more and the error message is even more unclear for me.
/ Preben

[mlimber]

Preben wrote:

don't know if this is the correct group, but try anyway.

I would like to use boostbindings from the boost sandbox together with
my project where I use boost::numeric: :ublas vectors and matrices.

I would guess that I could use these directly with the bindings in the
boost sandbox like:

----
#include <boost/numeric/bindings/lapack/gesv.hpp>
#include "MechComputer.h pp"

#include <iostream>

using namespace boost::numeric;
namespace lapack = boost::numeric: :bindings::lapa ck;
using namespace std;
namespace rw {
namespace dynamics {
MechComputer::M echComputer(con st MechSystem& msystem,
const ublas::vector<d ouble>& q0,
const ublas::vector<d ouble>& v0)
: _system(msystem ), _q(q0), _v(v0) {
}
const ublas::vector<d ouble>& MechComputer::s tep(double h) {
ublas::matrix<d oubledG = _system.DConstr aints(_q);
std::vector<ubl as::matrix<doub le ddG =
_system.DDConst raints(_q);
ublas::vector<d oubledV = _system.DV(_q);
ublas::vector<d oublem = _system.GetMass ();

ublas::matrix<d oublemM = M(dG, m);
ublas::vector<d oublevFT = mFmT(dG, ddG, dV, _v, m);

ublas::matrix<d oubleB(vFT.size (), 1);

lapack::gesv(M, B); // Result in B

return m; // Just return something
}
}
}
----

Compiling with this command-line:
----
preben@thinkgen ~/working/master/src $ icc -I
/home/preben/working/boost-sandbox -c MechComputer.cp p
MechComputer.cp p(31): error: no instance of overloaded function
"boost::numeric ::bindings::lap ack::gesv" matches the argument list
argument types are: (boost::numeric ::ublas::matrix <double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
(const boost::numeric: :ublas::matrix< double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
&, const boost::numeric: :ublas::vector< double,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>
&), boost::numeric: :ublas::matrix< double,
boost::numeric: :ublas::row_maj or,
boost::numeric: :ublas::unbound ed_array<double , std::allocator< double>>>)
lapack::gesv(M, B); // Result in B
^

compilation aborted for MechComputer.cp p (code 2)
preben@thinkgen ~/working/master/src $
----

I don't understand the error and why it's wrong.
g++ complains even more and the error message is even more unclear for me.

Probably better to try the Boost User list. Also, you might consider
using STLFilt to make template error messages more readable.

Cheers! --M