Mercurial > octave-nkf

comparison src/DLD-FUNCTIONS/chol.cc @ 8562:a6edd5c23cb5

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
use replacement methods if qrupdate is not available
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 22 Jan 2009 11:10:47 +0100
parents d66c9b6e506a
children fb1c929dbbb7
comparison
equal deleted inserted replaced
8561:66165de2cc42 8562:a6edd5c23cb5
574 print_usage (); 574 print_usage ();
575 575
576 return retval; 576 return retval;
577 } 577 }
578 578
579 #ifdef HAVE_QRUPDATE
580
581 DEFUN_DLD (cholupdate, args, nargout, 579 DEFUN_DLD (cholupdate, args, nargout,
582 "-*- texinfo -*-\n\ 580 "-*- texinfo -*-\n\
583 @deftypefn {Loadable Function} {[@var{R1}, @var{info}] =} cholupdate (@var{R}, @var{u}, @var{op})\n\ 581 @deftypefn {Loadable Function} {[@var{R1}, @var{info}] =} cholupdate (@var{R}, @var{u}, @var{op})\n\
584 Update or downdate a Cholesky factorization. Given an upper triangular\n\ 582 Update or downdate a Cholesky factorization. Given an upper triangular\n\
585 matrix @var{R} and a column vector @var{u}, attempt to determine another\n\ 583 matrix @var{R} and a column vector @var{u}, attempt to determine another\n\
1274 %! 1272 %!
1275 %! assert(norm(triu(R1)-R1,Inf) == 0) 1273 %! assert(norm(triu(R1)-R1,Inf) == 0)
1276 %! assert(norm(R1'*R1 - single(Ac(p,p)),Inf) < 1e1*eps('single')) 1274 %! assert(norm(R1'*R1 - single(Ac(p,p)),Inf) < 1e1*eps('single'))
1277 */ 1275 */
1278 1276
1279 #endif
1280
1281 /* 1277 /*
1282 ;;; Local Variables: *** 1278 ;;; Local Variables: ***
1283 ;;; mode: C++ *** 1279 ;;; mode: C++ ***
1284 ;;; End: *** 1280 ;;; End: ***
1285 */ 1281 */