# HG changeset patch
# User David Bateman <dbateman@free.fr>
# Date 1259046032 -3600
# Node ID ae70ce64beb7968ea0417dab36a62dddd0941ca3
# Parent  38f3c198ba630c26880db38bc978703e5fdbd426
Fix test for valid threshold value in sparse LU factorization

diff -r 38f3c198ba63 -r ae70ce64beb7 liboctave/ChangeLog
--- a/liboctave/ChangeLog	Tue Nov 24 07:58:44 2009 +0100
+++ b/liboctave/ChangeLog	Tue Nov 24 08:00:32 2009 +0100
@@ -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-09-18  Jaroslav Hajek  <highegg@gmail.com>
 
 	Version 3.2.3 released.
diff -r 38f3c198ba63 -r ae70ce64beb7 liboctave/SparseCmplxLU.cc
--- a/liboctave/SparseCmplxLU.cc	Tue Nov 24 07:58:44 2009 +0100
+++ b/liboctave/SparseCmplxLU.cc	Tue Nov 24 08:00:32 2009 +0100
@@ -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))
diff -r 38f3c198ba63 -r ae70ce64beb7 liboctave/SparsedbleLU.cc
--- a/liboctave/SparsedbleLU.cc	Tue Nov 24 07:58:44 2009 +0100
+++ b/liboctave/SparsedbleLU.cc	Tue Nov 24 08:00:32 2009 +0100
@@ -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))