Mercurial > octave

comparison scripts/sparse/ichol.m @ 31254:5d6b058a22dc

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: use sequence "%!#" to disable BIST tests. * interp1.m, logspace.m, triplequad.m, textread.m, __w2mpth__.m, ode15i.m, __print_parse_opts__.m, ichol.m, discrete_rnd.m, isstring.m: Use sequence "%!#" to disable BIST tests.
author Rik <rik@octave.org>
date Tue, 04 Oct 2022 07:34:47 -0700
parents 796f54d4ddbf
children 597f3ee61a48
comparison
equal deleted inserted replaced
31253:a40c0b7aa376 31254:5d6b058a22dc
362 ## Complex entry in the diagonal 362 ## Complex entry in the diagonal
363 %!error <non-real pivot> ichol (A7) 363 %!error <non-real pivot> ichol (A7)
364 364
365 ## ICHOLT tests 365 ## ICHOLT tests
366 366
367 #%!test 367 %!#test
368 %! opts.type = "ict"; 368 %! opts.type = "ict";
369 %! opts.droptol = 1e-1; 369 %! opts.droptol = 1e-1;
370 %! opts.michol = "off"; 370 %! opts.michol = "off";
371 %! L = ichol (A1, opts); 371 %! L = ichol (A1, opts);
372 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.2065, 1e-4); 372 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.2065, 1e-4);
374 %! U = ichol (A1, opts); 374 %! U = ichol (A1, opts);
375 %! assert (norm (A1 - U' * U, "fro") / norm (A1, "fro"), 0.2065, 1e-4); 375 %! assert (norm (A1 - U' * U, "fro") / norm (A1, "fro"), 0.2065, 1e-4);
376 %! opts.shape = "lower"; 376 %! opts.shape = "lower";
377 %! L = ichol (A1, opts); 377 %! L = ichol (A1, opts);
378 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.2065, 1e-4); 378 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.2065, 1e-4);
379 %! 379
380 #%!test 380 %!#test
381 %! opts.type = "ict"; 381 %! opts.type = "ict";
382 %! opts.droptol = 1e-1; 382 %! opts.droptol = 1e-1;
383 %! opts.michol = "on"; 383 %! opts.michol = "on";
384 %! L = ichol (A1, opts); 384 %! L = ichol (A1, opts);
385 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.3266, 1e-4); 385 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.3266, 1e-4);
388 %! assert (norm (A1 - U' * U, "fro") / norm (A1, "fro"), 0.3266, 1e-4); 388 %! assert (norm (A1 - U' * U, "fro") / norm (A1, "fro"), 0.3266, 1e-4);
389 %! opts.shape = "lower"; 389 %! opts.shape = "lower";
390 %! opts.diagcomp = 3e-3; 390 %! opts.diagcomp = 3e-3;
391 %! L = ichol (A1, opts); 391 %! L = ichol (A1, opts);
392 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.3266, 1e-4); 392 %! assert (norm (A1 - L * L', "fro") / norm (A1, "fro"), 0.3266, 1e-4);
393 %! 393
394 %!test 394 %!test
395 %! opts.type = "ict"; 395 %! opts.type = "ict";
396 %! opts.droptol = 1e-1; 396 %! opts.droptol = 1e-1;
397 %! opts.michol = "off"; 397 %! opts.michol = "off";
398 %! L = ichol (A2, opts); 398 %! L = ichol (A2, opts);
399 %! assert (norm (A2 - L*L', "fro") / norm (A2, "fro"), 0.0893, 1e-4); 399 %! assert (norm (A2 - L*L', "fro") / norm (A2, "fro"), 0.0893, 1e-4);
400 %! opts.michol = "on"; 400 %! opts.michol = "on";
401 %! L = ichol (A2, opts); 401 %! L = ichol (A2, opts);
402 %! assert (norm (A2 - L*L', "fro") / norm (A2, "fro"), 0.2377, 1e-4); 402 %! assert (norm (A2 - L*L', "fro") / norm (A2, "fro"), 0.2377, 1e-4);
403 %! 403
404 %!test 404 %!test
405 %! opts.type = "ict"; 405 %! opts.type = "ict";
406 %! opts.droptol = 1e-1; 406 %! opts.droptol = 1e-1;
407 %! opts.michol = "off"; 407 %! opts.michol = "off";
408 %! L = ichol (A3, opts); 408 %! L = ichol (A3, opts);
409 %! assert (norm (A3 - L*L', "fro") / norm (A3, "fro"), eps, eps); 409 %! assert (norm (A3 - L*L', "fro") / norm (A3, "fro"), eps, eps);
410 %! opts.michol = "on"; 410 %! opts.michol = "on";
411 %! L = ichol (A3, opts); 411 %! L = ichol (A3, opts);
412 %! assert (norm (A3 - L*L', "fro") / norm (A3, "fro"), eps, eps); 412 %! assert (norm (A3 - L*L', "fro") / norm (A3, "fro"), eps, eps);
413 %! 413
414 %!test 414 %!test
415 %! opts.type = "ict"; 415 %! opts.type = "ict";
416 %! opts.droptol = 1e-1; 416 %! opts.droptol = 1e-1;
417 %! opts.michol = "off"; 417 %! opts.michol = "off";
418 %! L = ichol (A4, opts); 418 %! L = ichol (A4, opts);
419 %! assert (norm (A4 - L*L', "fro") / norm (A4, "fro"), 0.1224, 1e-4); 419 %! assert (norm (A4 - L*L', "fro") / norm (A4, "fro"), 0.1224, 1e-4);
420 %! opts.michol = "on"; 420 %! opts.michol = "on";
421 %! L = ichol (A4, opts); 421 %! L = ichol (A4, opts);
422 %! assert (norm (A4 - L*L', "fro") / norm (A4, "fro"), 0.2118, 1e-4); 422 %! assert (norm (A4 - L*L', "fro") / norm (A4, "fro"), 0.2118, 1e-4);
423 %! 423
424 %!test 424 %!test
425 %! opts.type = "ict"; 425 %! opts.type = "ict";
426 %! opts.droptol = 1e-1; 426 %! opts.droptol = 1e-1;
427 %! opts.michol = "off"; 427 %! opts.michol = "off";
428 %! L = ichol (A5, opts); 428 %! L = ichol (A5, opts);