|
Ublas : |
Subject: Re: [ublas] [bindings] New traits system
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2010-01-09 07:52:16
Thomas Klimpel wrote:
> I don't think that this has anything to do with MKL, and I think the crash
> happens earlier before you are in the MKL on an assert. This is again an
> inconsistency in the documentation, because the optimal workspace
> determined by lapack for gelsd is smaller than the minimal required
> workspace mentioned in the documentation. So I cheated when updating the
> regression tests, and did some hand-editing to gelsd.hpp. The updated
> bindings have of course reverted the hand editing, so it's no surprise
> that it no longer works.
There was an error in the computation of nlvl. The correct computation was
in the Fortran code.Fixed and checked in. Gelsd lapack regression test
passes.
Cheers,
Rutger