Mercurial > octave

comparison scripts/optimization/fsolve.m @ 21578:683a1beee538

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Use "FIXME:" for all code blocks needing further attention. * files-dock-widget.cc, file-editor-tab.cc, file-editor-tab.h, file-editor.cc, octave-qscintilla.cc, octave-qt-link.cc, __ilu__.cc, debug.h, oct-stream.cc, pt-jit.cc, __magick_read__.cc, ov-classdef.cc, pt-stmt.cc, oct-sort.cc, inputParser.m, validateattributes.m, fminunc.m, fsolve.m, fzero.m, __scatter__.m, print.m: Use "FIXME:" for all code blocks needing further attention.
author Rik <rik@octave.org>
date Fri, 01 Apr 2016 14:55:13 -0700
parents 3d60ed163b70
children b571fc85953f
comparison
equal deleted inserted replaced
21577:31823239207e 21578:683a1beee538
581 581
582 ## Solve the double dogleg trust-region least-squares problem: 582 ## Solve the double dogleg trust-region least-squares problem:
583 ## Minimize norm(r*x-b) subject to the constraint norm(d.*x) <= delta, 583 ## Minimize norm(r*x-b) subject to the constraint norm(d.*x) <= delta,
584 ## x being a convex combination of the gauss-newton and scaled gradient. 584 ## x being a convex combination of the gauss-newton and scaled gradient.
585 585
586 ## TODO: error checks 586 ## FIXME: error checks
587 ## TODO: handle singularity, or leave it up to mldivide? 587 ## FIXME: handle singularity, or leave it up to mldivide?
588 588
589 function x = __dogleg__ (r, b, d, delta) 589 function x = __dogleg__ (r, b, d, delta)
590 ## Get Gauss-Newton direction. 590 ## Get Gauss-Newton direction.
591 x = r \ b; 591 x = r \ b;
592 xn = norm (d .* x); 592 xn = norm (d .* x);