comparison liboctave/CMatrix.cc @ 7532:493bb0de3199

avoid another xGELSD workspace query bug
author John W. Eaton <jwe@octave.org>
date Tue, 26 Feb 2008 02:47:56 -0500
parents 8c32f95c2639
children f9983d2761df
comparison
equal deleted inserted replaced
7531:c9a476b1e664 7532:493bb0de3199
2480 2480
2481 if (n-3*m > addend) 2481 if (n-3*m > addend)
2482 addend = n-3*m; 2482 addend = n-3*m;
2483 2483
2484 const octave_idx_type lworkaround = 4*m + m*m + addend; 2484 const octave_idx_type lworkaround = 4*m + m*m + addend;
2485
2486 if (std::real (work(0)) < lworkaround)
2487 work(0) = lworkaround;
2488 }
2489 else if (m >= n)
2490 {
2491 octave_idx_type lworkaround = 2*m + m*nrhs;
2485 2492
2486 if (std::real (work(0)) < lworkaround) 2493 if (std::real (work(0)) < lworkaround)
2487 work(0) = lworkaround; 2494 work(0) = lworkaround;
2488 } 2495 }
2489 2496