Re: [ublas] [bindings] How do I call the new gees? Why is it different from the old one?

Subject: Re: [ublas] [bindings] How do I call the new gees? Why is it different from the old one?
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2010-04-19 11:11:11

Paul Leopardi wrote:

> #elif defined(_GLUCAT_USE_BINDINGS)
>     lapack::gees ('N', 'N', (logical_t*)0, T, fortran_int_t(0), lambda, V
>     );
> #else

[snip]

> In particular, I was expecting the call to match
> 
> gees( const char jobvs, const char sort, logical_t* select, MatrixA& a,
>         fortran_int_t& sdim, VectorW& w, MatrixVS& vs )
> 

Isn't the non-match due to sdim being an LValue? I.e., please try 

fortran_int_t sdim;
lapack::gees( 'N', 'N', 0, sdim, lambda, V );

Cheers,

Rutger