|
Ublas : |
Subject: Re: [ublas] [bindings] New traits system
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2010-01-09 09:56:44
Jesse wrote:
> Yup, Thomas was right and I am an idiot thinking it was a problem inside
> of MKL. I commented out the assertions on the workspace (line 142 of
> gelsd.hpp) and it worked great for the lapack::optimal_workspace().
>
> If these kinds of bugs are common in automatic generation, may I make a
> proposal? Add on a tag to detail::workspace2<> to say whether an optimal
> workspace was used or not.
>
> Then for the assertion checks on the size of the workspace, only check if
> NOT an optimal workspace?
>
Please check the latest gelsd in svn -- it's fixed. The generator is full of
templates which override stuff found in the lapack documentation (see
libs/numeric/bindings/tools/templates). Let's say this is not the first
inconsistency found :-)
Cheers,
Rutger