|
|
Ublas : |
From: Nico Galoppo (ngaloppo_at_[hidden])
Date: 2008-03-13 22:05:44
Hi,
Upon calling the lapack gesvd in numeric/bindings/lapack/gesvd.cpp,
with JobU = 'O', I'm getting an assertion error on line 418:
417 assert ((jobu == 'N' && traits::leading_dimension (u) >= 1)
418 || (jobu == 'O' && traits::leading_dimension (u) >= 1)
419 || (jobu == 'A' && traits::leading_dimension (u) >= m)
420 || (jobu == 'S' && traits::leading_dimension (u) >= m));
I don't understand why this assertion is correct. jobu == 'O' means
that the left-hand vectors should be written to input matrix A, not to
u, hence, the u input to gesvd should essentially be able to be a null
pointer (or zero-size matrix).
Am I missing something?
Thanks!
--nico
-- Nico Galoppo :: http://www.ngaloppo.org