Mercurial > hg > octave-nkf

comparison liboctave/dMatrix.cc @ 10813:2c2d4a2f1047

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix workspace bug workaround for xGELSD (since 7486:6a6d2abe51ff)
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 22 Jul 2010 09:54:03 +0200
parents 7c542263a92a
children 3622db30ff05
comparison
equal deleted inserted replaced
10812:5b68000faac1 10813:2c2d4a2f1047
2292 2292
2293 // The workspace query is broken in at least LAPACK 3.0.0 2293 // The workspace query is broken in at least LAPACK 3.0.0
2294 // through 3.1.1 when n >= mnthr. The obtuse formula below 2294 // through 3.1.1 when n >= mnthr. The obtuse formula below
2295 // should provide sufficient workspace for DGELSD to operate 2295 // should provide sufficient workspace for DGELSD to operate
2296 // efficiently. 2296 // efficiently.
2297 if (n >= mnthr) 2297 if (n > m && n >= mnthr)
2298 { 2298 {
2299 const octave_idx_type wlalsd 2299 const octave_idx_type wlalsd
2300 = 9*m + 2*m*smlsiz + 8*m*nlvl + m*nrhs + (smlsiz+1)*(smlsiz+1); 2300 = 9*m + 2*m*smlsiz + 8*m*nlvl + m*nrhs + (smlsiz+1)*(smlsiz+1);
2301 2301
2302 octave_idx_type addend = m; 2302 octave_idx_type addend = m;