Mercurial > hg > octave-lyh
changeset 9759:97de6c916498
Fix test for valid threshold value in sparse LU factorization
author | David Bateman <dbateman@free.fr> |
---|---|
date | Tue, 27 Oct 2009 20:11:47 +0100 |
parents | 09da0bd91412 |
children | c6b97d218345 |
files | liboctave/ChangeLog liboctave/SparseCmplxLU.cc liboctave/SparsedbleLU.cc |
diffstat | 3 files changed, 16 insertions(+), 4 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/ChangeLog +++ b/liboctave/ChangeLog @@ -1,3 +1,15 @@ +2009-10-27 David Bateman <dbateman@free.fr> + + * SparsedbleLU.cc (SparseLU::SparseLU (const SparseMatrix&, + const Matrix&, bool), SparseLU::SparseLU (const SparseMatrix&, + const ColumnVector&, const Matrix&, boo, bool, double, bool, + bool)): Fix test for valid pivot tolerance. + * SparseCmplxLU.cc (SparseComplexLU::SparseComplexLU + (const SparseComplexMatrix&, const Matrix&, bool), + SparseComplexLU::SparseComplexLU (const SparseComplexMatrix&, + const ColumnVector&, const Matrix&, boo, bool, double, bool, + bool)): Ditto. + 2009-10-23 Jaroslav Hajek <highegg@gmail.com> * idx-vector.cc (idx_vector::idx_range_rep::idx_range_rep
--- a/liboctave/SparseCmplxLU.cc +++ b/liboctave/SparseCmplxLU.cc @@ -57,7 +57,7 @@ double tmp = octave_sparse_params::get_key ("spumoni"); if (!xisnan (tmp)) Control (UMFPACK_PRL) = tmp; - if (piv_thres.nelem() != 2) + if (piv_thres.nelem() == 2) { tmp = (piv_thres (0) > 1. ? 1. : piv_thres (0)); if (!xisnan (tmp)) @@ -266,7 +266,7 @@ double tmp = octave_sparse_params::get_key ("spumoni"); if (!xisnan (tmp)) Control (UMFPACK_PRL) = tmp; - if (piv_thres.nelem() != 2) + if (piv_thres.nelem() == 2) { tmp = (piv_thres (0) > 1. ? 1. : piv_thres (0)); if (!xisnan (tmp))
--- a/liboctave/SparsedbleLU.cc +++ b/liboctave/SparsedbleLU.cc @@ -57,7 +57,7 @@ if (!xisnan (tmp)) Control (UMFPACK_PRL) = tmp; - if (piv_thres.nelem() != 2) + if (piv_thres.nelem() == 2) { tmp = (piv_thres (0) > 1. ? 1. : piv_thres (0)); if (!xisnan (tmp)) @@ -253,7 +253,7 @@ if (!xisnan (tmp)) Control (UMFPACK_PRL) = tmp; - if (piv_thres.nelem() != 2) + if (piv_thres.nelem() == 2) { tmp = (piv_thres (0) > 1. ? 1. : piv_thres (0)); if (!xisnan (tmp))