diff scripts/linear-algebra/condest.m @ 8912:57c3155754d6

fix condest
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 05 Mar 2009 09:37:59 +0100
parents e07e93c04080
children eb63fbe60fab
line wrap: on
line diff
--- 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;