# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1236242279 -3600
# Node ID 57c3155754d69c8a01d2ec90fe37b5914a6f5de9
# Parent  eeed8d6498118146dfc43e9d1f2515ff0c0ff758
fix condest

diff -r eeed8d649811 -r 57c3155754d6 scripts/ChangeLog
--- a/scripts/ChangeLog	Thu Mar 05 02:56:37 2009 -0500
+++ b/scripts/ChangeLog	Thu Mar 05 09:37:59 2009 +0100
@@ -1,3 +1,8 @@
+2009-03-05  Jaroslav Hajek  <highegg@gmail.com>
+
+	* linear-algebra/condest.m: Use lu rather than splu. Use relative tols
+	for tests.
+
 2009-03-05  John W. Eaton  <jwe@octave.org>
 
 	* deprecated/spdiag.m: Ensure sparse result.
diff -r eeed8d649811 -r 57c3155754d6 scripts/linear-algebra/condest.m
--- a/scripts/linear-algebra/condest.m	Thu Mar 05 02:56:37 2009 -0500
+++ b/scripts/linear-algebra/condest.m	Thu Mar 05 09:37:59 2009 +0100
@@ -163,7 +163,7 @@
 
   if (! have_solve)
     if (issparse (A))
-      [L, U, P, Pc] = splu (A);
+      [L, U, P, Pc] = lu (A);
       solve = @(x) Pc' * (U \ (L \ (P * x)));
       solve_t = @(x) P' * (L' \ (U' \ (Pc * x)));
     else
@@ -203,7 +203,7 @@
 %!  A = hilb (N);
 %!  cA = condest (A);
 %!  cA_test = norm (inv (A), 1) * norm (A, 1);
-%!  assert (cA, cA_test, 2^-8);
+%!  assert (cA, cA_test, -2^-8);
 
 %!test
 %!  N = 6;
@@ -211,7 +211,7 @@
 %!  solve = @(x) A\x; solve_t = @(x) A'\x;
 %!  cA = condest (A, solve, solve_t);
 %!  cA_test = norm (inv (A), 1) * norm (A, 1);
-%!  assert (cA, cA_test, 2^-8);
+%!  assert (cA, cA_test, -2^-8);
 
 %!test
 %!  N = 6;
@@ -220,7 +220,7 @@
 %!  solve = @(x) A\x; solve_t = @(x) A'\x;
 %!  cA = condest (apply, apply_t, solve, solve_t, N);
 %!  cA_test = norm (inv (A), 1) * norm (A, 1);
-%!  assert (cA, cA_test, 2^-6);
+%!  assert (cA, cA_test, -2^-6);
 
 %!test
 %!  N = 12;