Subject: Re: [ublas] [bindings][lapack] Initial high-level solve
From: Jesse Perla (jesseperla_at_[hidden])
Date: 2009-03-24 13:14:04
Looks great. One thing that would be nice is if the B could support vectors
as well as matrices. Quite often you want to solve AX=b with b a ublas
vector and don't want to convert between the two types. I am not sure if
there is the ability for the traits to treat a vector as an Nx1 matrix. I
assume this could use some kind of matrix adaptor or decoration to say
whether it is an Nx1 vs. 1xN for conformability.
I assume that it will be based on some sort of traits similar to
<http://lists.boost.org/MailArchives/ublas/2009/03/3384.php>It would be
nice if there was an equivalent decoration to these traits for calling this
function so the user isn't forced to use a bunch of adaptors with every
I look forward to testing!
On Sun, Mar 22, 2009 at 7:08 AM, Rutger ter Borg <rutger_at_[hidden]> wrote:
> I've created an initial prototype for a high-level interface to solve. It
> uses Boost.Parameter to allow a user to pass parameters by keyword;
> low-level routines create a temporary pivot vector if not provided and
> needed by that routine.
> See below for the syntax, all expressions are valid.
> #include <boost/numeric/bindings/lapack/solve.hpp>
> using namespace lapack::keywords;
> // basic usage
> lapack::solve( A, B );
> lapack::solve( A, B, _pivot = ipiv );
> // this is also valid
> lapack::solve( _B = B, _A = A );
> lapack::solve( A, _pivot = ipiv, _B = B );
> What we could use next is a good set of keywords/type decorators for
> addressing all solve routines and expert driver options. Suggestions are
> very welcome.
> ublas mailing list