Re: [ublas] lapack bindings: ppsv returns info=2

Subject: Re: [ublas] lapack bindings: ppsv returns info=2
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2011-04-23 07:58:39

On 2011-04-22 22:59, Rutger ter Borg wrote:

> You're right, the pp-prefix should be fine in this case. I sounds like
> something isn't passed right -- it has been a while since I've checked
> the details of packed storage, but have you tried column_major ordering
> for the symmetric matrix?
>
> ublas::symmetric_matrix<double, ..., column_major >
>
> If this happens to make things work, it is a bug in the bindings. In
> that case, it sounds like it is missing a static assertion / check for
> this case.
>

I can confirm that it is probably due to insufficient checking of the 
packed data structure (adding column_major does the trick). If it is 
possible to do stuff like swapping lower/row_major with 
upper/column_major, we should add that trickery to uplo_tag.

Cheers,

Rutger