Mercurial > hg > octave-lyh

comparison scripts/linear-algebra/condest.m @ 7239:7d7f05863f27

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-12-03 20:10:06 by jwe]
author jwe
date Mon, 03 Dec 2007 20:10:06 +0000
parents b48a21816f2e
children 26f42a14d4f6
comparison
equal deleted inserted replaced
7238:22f210ca8cc1 7239:7d7f05863f27
194 %!test 194 %!test
195 %! N = 6; 195 %! N = 6;
196 %! A = hilb (N); 196 %! A = hilb (N);
197 %! cA = condest (A); 197 %! cA = condest (A);
198 %! cA_test = norm (inv (A), 1) * norm (A, 1); 198 %! cA_test = norm (inv (A), 1) * norm (A, 1);
199 %! assert (cA, cA_test, 2^-12); 199 %! assert (cA, cA_test, 2^-8);
200 200
201 %!test 201 %!test
202 %! N = 6; 202 %! N = 6;
203 %! A = hilb (N); 203 %! A = hilb (N);
204 %! solve = @(x) A\x; solve_t = @(x) A'\x; 204 %! solve = @(x) A\x; solve_t = @(x) A'\x;
205 %! cA = condest (A, solve, solve_t); 205 %! cA = condest (A, solve, solve_t);
206 %! cA_test = norm (inv (A), 1) * norm (A, 1); 206 %! cA_test = norm (inv (A), 1) * norm (A, 1);
207 %! assert (cA, cA_test, 2^-12); 207 %! assert (cA, cA_test, 2^-8);
208 208
209 %!test 209 %!test
210 %! N = 6; 210 %! N = 6;
211 %! A = hilb (N); 211 %! A = hilb (N);
212 %! apply = @(x) A*x; apply_t = @(x) A'*x; 212 %! apply = @(x) A*x; apply_t = @(x) A'*x;