commons-notifications mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From l..@apache.org
Subject svn commit: r948061 [3/22] - in /websites/production/commons/content/proper/commons-math: ./ xref/ xref/org/apache/commons/math3/ xref/org/apache/commons/math3/analysis/ xref/org/apache/commons/math3/analysis/differentiation/ xref/org/apache/commons/ma...
Date Fri, 17 Apr 2015 19:17:18 GMT
Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/BrentSolver.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/BrentSolver.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/BrentSolver.html Fri Apr 17 19:17:13 2015
@@ -37,208 +37,219 @@
 <a class="jxr_linenumber" name="L29" href="#L29">29</a>  <em class="jxr_javadoccomment"> * The function should be continuous but not necessarily smooth.</em>
 <a class="jxr_linenumber" name="L30" href="#L30">30</a>  <em class="jxr_javadoccomment"> * The {@code solve} method returns a zero {@code x} of the function {@code f}</em>
 <a class="jxr_linenumber" name="L31" href="#L31">31</a>  <em class="jxr_javadoccomment"> * in the given interval {@code [a, b]} to within a tolerance</em>
-<a class="jxr_linenumber" name="L32" href="#L32">32</a>  <em class="jxr_javadoccomment"> * {@code 6 eps abs(x) + t} where {@code eps} is the relative accuracy and</em>
+<a class="jxr_linenumber" name="L32" href="#L32">32</a>  <em class="jxr_javadoccomment"> * {@code 2 eps abs(x) + t} where {@code eps} is the relative accuracy and</em>
 <a class="jxr_linenumber" name="L33" href="#L33">33</a>  <em class="jxr_javadoccomment"> * {@code t} is the absolute accuracy.</em>
 <a class="jxr_linenumber" name="L34" href="#L34">34</a>  <em class="jxr_javadoccomment"> * The given interval must bracket the root.</em>
-<a class="jxr_linenumber" name="L35" href="#L35">35</a>  <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L36" href="#L36">36</a>  <em class="jxr_javadoccomment"> */</em>
-<a class="jxr_linenumber" name="L37" href="#L37">37</a>  <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a> <strong class="jxr_keyword">extends</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/AbstractUnivariateSolver.html">AbstractUnivariateSolver</a> {
-<a class="jxr_linenumber" name="L38" href="#L38">38</a>  
-<a class="jxr_linenumber" name="L39" href="#L39">39</a>      <em class="jxr_javadoccomment">/** Default absolute accuracy. */</em>
-<a class="jxr_linenumber" name="L40" href="#L40">40</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> DEFAULT_ABSOLUTE_ACCURACY = 1e-6;
-<a class="jxr_linenumber" name="L41" href="#L41">41</a>  
-<a class="jxr_linenumber" name="L42" href="#L42">42</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L43" href="#L43">43</a>  <em class="jxr_javadoccomment">     * Construct a solver with default accuracy (1e-6).</em>
-<a class="jxr_linenumber" name="L44" href="#L44">44</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L45" href="#L45">45</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>() {
-<a class="jxr_linenumber" name="L46" href="#L46">46</a>          <strong class="jxr_keyword">this</strong>(DEFAULT_ABSOLUTE_ACCURACY);
-<a class="jxr_linenumber" name="L47" href="#L47">47</a>      }
-<a class="jxr_linenumber" name="L48" href="#L48">48</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L49" href="#L49">49</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
-<a class="jxr_linenumber" name="L50" href="#L50">50</a>  <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L51" href="#L51">51</a>  <em class="jxr_javadoccomment">     * @param absoluteAccuracy Absolute accuracy.</em>
-<a class="jxr_linenumber" name="L52" href="#L52">52</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L53" href="#L53">53</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> absoluteAccuracy) {
-<a class="jxr_linenumber" name="L54" href="#L54">54</a>          <strong class="jxr_keyword">super</strong>(absoluteAccuracy);
-<a class="jxr_linenumber" name="L55" href="#L55">55</a>      }
-<a class="jxr_linenumber" name="L56" href="#L56">56</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L57" href="#L57">57</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
-<a class="jxr_linenumber" name="L58" href="#L58">58</a>  <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L59" href="#L59">59</a>  <em class="jxr_javadoccomment">     * @param relativeAccuracy Relative accuracy.</em>
+<a class="jxr_linenumber" name="L35" href="#L35">35</a>  <em class="jxr_javadoccomment"> * &lt;p&gt;</em>
+<a class="jxr_linenumber" name="L36" href="#L36">36</a>  <em class="jxr_javadoccomment"> *  The reference implementation is given in chapter 4 of</em>
+<a class="jxr_linenumber" name="L37" href="#L37">37</a>  <em class="jxr_javadoccomment"> *  &lt;blockquote&gt;</em>
+<a class="jxr_linenumber" name="L38" href="#L38">38</a>  <em class="jxr_javadoccomment"> *   &lt;b&gt;Algorithms for Minimization Without Derivatives&lt;/b&gt;&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L39" href="#L39">39</a>  <em class="jxr_javadoccomment"> *   &lt;em&gt;Richard P. Brent&lt;/em&gt;&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L40" href="#L40">40</a>  <em class="jxr_javadoccomment"> *   Dover, 2002&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L41" href="#L41">41</a>  <em class="jxr_javadoccomment"> *  &lt;/blockquote&gt;</em>
+<a class="jxr_linenumber" name="L42" href="#L42">42</a>  <em class="jxr_javadoccomment"> * &lt;/p&gt;</em>
+<a class="jxr_linenumber" name="L43" href="#L43">43</a>  <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L44" href="#L44">44</a>  <em class="jxr_javadoccomment"> * @see BaseAbstractUnivariateSolver</em>
+<a class="jxr_linenumber" name="L45" href="#L45">45</a>  <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L46" href="#L46">46</a>  <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a> <strong class="jxr_keyword">extends</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/AbstractUnivariateSolver.html">AbstractUnivariateSolver</a> {
+<a class="jxr_linenumber" name="L47" href="#L47">47</a>  
+<a class="jxr_linenumber" name="L48" href="#L48">48</a>      <em class="jxr_javadoccomment">/** Default absolute accuracy. */</em>
+<a class="jxr_linenumber" name="L49" href="#L49">49</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> DEFAULT_ABSOLUTE_ACCURACY = 1e-6;
+<a class="jxr_linenumber" name="L50" href="#L50">50</a>  
+<a class="jxr_linenumber" name="L51" href="#L51">51</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L52" href="#L52">52</a>  <em class="jxr_javadoccomment">     * Construct a solver with default absolute accuracy (1e-6).</em>
+<a class="jxr_linenumber" name="L53" href="#L53">53</a>  <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L54" href="#L54">54</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>() {
+<a class="jxr_linenumber" name="L55" href="#L55">55</a>          <strong class="jxr_keyword">this</strong>(DEFAULT_ABSOLUTE_ACCURACY);
+<a class="jxr_linenumber" name="L56" href="#L56">56</a>      }
+<a class="jxr_linenumber" name="L57" href="#L57">57</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L58" href="#L58">58</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
+<a class="jxr_linenumber" name="L59" href="#L59">59</a>  <em class="jxr_javadoccomment">     *</em>
 <a class="jxr_linenumber" name="L60" href="#L60">60</a>  <em class="jxr_javadoccomment">     * @param absoluteAccuracy Absolute accuracy.</em>
 <a class="jxr_linenumber" name="L61" href="#L61">61</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L62" href="#L62">62</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> relativeAccuracy,
-<a class="jxr_linenumber" name="L63" href="#L63">63</a>                         <strong class="jxr_keyword">double</strong> absoluteAccuracy) {
-<a class="jxr_linenumber" name="L64" href="#L64">64</a>          <strong class="jxr_keyword">super</strong>(relativeAccuracy, absoluteAccuracy);
-<a class="jxr_linenumber" name="L65" href="#L65">65</a>      }
-<a class="jxr_linenumber" name="L66" href="#L66">66</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L67" href="#L67">67</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
-<a class="jxr_linenumber" name="L68" href="#L68">68</a>  <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L69" href="#L69">69</a>  <em class="jxr_javadoccomment">     * @param relativeAccuracy Relative accuracy.</em>
-<a class="jxr_linenumber" name="L70" href="#L70">70</a>  <em class="jxr_javadoccomment">     * @param absoluteAccuracy Absolute accuracy.</em>
-<a class="jxr_linenumber" name="L71" href="#L71">71</a>  <em class="jxr_javadoccomment">     * @param functionValueAccuracy Function value accuracy.</em>
-<a class="jxr_linenumber" name="L72" href="#L72">72</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L73" href="#L73">73</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> relativeAccuracy,
-<a class="jxr_linenumber" name="L74" href="#L74">74</a>                         <strong class="jxr_keyword">double</strong> absoluteAccuracy,
-<a class="jxr_linenumber" name="L75" href="#L75">75</a>                         <strong class="jxr_keyword">double</strong> functionValueAccuracy) {
-<a class="jxr_linenumber" name="L76" href="#L76">76</a>          <strong class="jxr_keyword">super</strong>(relativeAccuracy, absoluteAccuracy, functionValueAccuracy);
-<a class="jxr_linenumber" name="L77" href="#L77">77</a>      }
-<a class="jxr_linenumber" name="L78" href="#L78">78</a>  
-<a class="jxr_linenumber" name="L79" href="#L79">79</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L80" href="#L80">80</a>  <em class="jxr_javadoccomment">     * {@inheritDoc}</em>
-<a class="jxr_linenumber" name="L81" href="#L81">81</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L82" href="#L82">82</a>      @Override
-<a class="jxr_linenumber" name="L83" href="#L83">83</a>      <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> doSolve()
-<a class="jxr_linenumber" name="L84" href="#L84">84</a>          <strong class="jxr_keyword">throws</strong> NoBracketingException,
-<a class="jxr_linenumber" name="L85" href="#L85">85</a>                 TooManyEvaluationsException,
-<a class="jxr_linenumber" name="L86" href="#L86">86</a>                 <a href="../../../../../../org/apache/commons/math3/exception/NumberIsTooLargeException.html">NumberIsTooLargeException</a> {
-<a class="jxr_linenumber" name="L87" href="#L87">87</a>          <strong class="jxr_keyword">double</strong> min = getMin();
-<a class="jxr_linenumber" name="L88" href="#L88">88</a>          <strong class="jxr_keyword">double</strong> max = getMax();
-<a class="jxr_linenumber" name="L89" href="#L89">89</a>          <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> initial = getStartValue();
-<a class="jxr_linenumber" name="L90" href="#L90">90</a>          <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> functionValueAccuracy = getFunctionValueAccuracy();
-<a class="jxr_linenumber" name="L91" href="#L91">91</a>  
-<a class="jxr_linenumber" name="L92" href="#L92">92</a>          verifySequence(min, initial, max);
-<a class="jxr_linenumber" name="L93" href="#L93">93</a>  
-<a class="jxr_linenumber" name="L94" href="#L94">94</a>          <em class="jxr_comment">// Return the initial guess if it is good enough.</em>
-<a class="jxr_linenumber" name="L95" href="#L95">95</a>          <strong class="jxr_keyword">double</strong> yInitial = computeObjectiveValue(initial);
-<a class="jxr_linenumber" name="L96" href="#L96">96</a>          <strong class="jxr_keyword">if</strong> (FastMath.abs(yInitial) &lt;= functionValueAccuracy) {
-<a class="jxr_linenumber" name="L97" href="#L97">97</a>              <strong class="jxr_keyword">return</strong> initial;
-<a class="jxr_linenumber" name="L98" href="#L98">98</a>          }
-<a class="jxr_linenumber" name="L99" href="#L99">99</a>  
-<a class="jxr_linenumber" name="L100" href="#L100">100</a>         <em class="jxr_comment">// Return the first endpoint if it is good enough.</em>
-<a class="jxr_linenumber" name="L101" href="#L101">101</a>         <strong class="jxr_keyword">double</strong> yMin = computeObjectiveValue(min);
-<a class="jxr_linenumber" name="L102" href="#L102">102</a>         <strong class="jxr_keyword">if</strong> (FastMath.abs(yMin) &lt;= functionValueAccuracy) {
-<a class="jxr_linenumber" name="L103" href="#L103">103</a>             <strong class="jxr_keyword">return</strong> min;
-<a class="jxr_linenumber" name="L104" href="#L104">104</a>         }
-<a class="jxr_linenumber" name="L105" href="#L105">105</a> 
-<a class="jxr_linenumber" name="L106" href="#L106">106</a>         <em class="jxr_comment">// Reduce interval if min and initial bracket the root.</em>
-<a class="jxr_linenumber" name="L107" href="#L107">107</a>         <strong class="jxr_keyword">if</strong> (yInitial * yMin &lt; 0) {
-<a class="jxr_linenumber" name="L108" href="#L108">108</a>             <strong class="jxr_keyword">return</strong> brent(min, initial, yMin, yInitial);
+<a class="jxr_linenumber" name="L62" href="#L62">62</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> absoluteAccuracy) {
+<a class="jxr_linenumber" name="L63" href="#L63">63</a>          <strong class="jxr_keyword">super</strong>(absoluteAccuracy);
+<a class="jxr_linenumber" name="L64" href="#L64">64</a>      }
+<a class="jxr_linenumber" name="L65" href="#L65">65</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L66" href="#L66">66</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
+<a class="jxr_linenumber" name="L67" href="#L67">67</a>  <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L68" href="#L68">68</a>  <em class="jxr_javadoccomment">     * @param relativeAccuracy Relative accuracy.</em>
+<a class="jxr_linenumber" name="L69" href="#L69">69</a>  <em class="jxr_javadoccomment">     * @param absoluteAccuracy Absolute accuracy.</em>
+<a class="jxr_linenumber" name="L70" href="#L70">70</a>  <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L71" href="#L71">71</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> relativeAccuracy,
+<a class="jxr_linenumber" name="L72" href="#L72">72</a>                         <strong class="jxr_keyword">double</strong> absoluteAccuracy) {
+<a class="jxr_linenumber" name="L73" href="#L73">73</a>          <strong class="jxr_keyword">super</strong>(relativeAccuracy, absoluteAccuracy);
+<a class="jxr_linenumber" name="L74" href="#L74">74</a>      }
+<a class="jxr_linenumber" name="L75" href="#L75">75</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L76" href="#L76">76</a>  <em class="jxr_javadoccomment">     * Construct a solver.</em>
+<a class="jxr_linenumber" name="L77" href="#L77">77</a>  <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L78" href="#L78">78</a>  <em class="jxr_javadoccomment">     * @param relativeAccuracy Relative accuracy.</em>
+<a class="jxr_linenumber" name="L79" href="#L79">79</a>  <em class="jxr_javadoccomment">     * @param absoluteAccuracy Absolute accuracy.</em>
+<a class="jxr_linenumber" name="L80" href="#L80">80</a>  <em class="jxr_javadoccomment">     * @param functionValueAccuracy Function value accuracy.</em>
+<a class="jxr_linenumber" name="L81" href="#L81">81</a>  <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L82" href="#L82">82</a>  <em class="jxr_javadoccomment">     * @see BaseAbstractUnivariateSolver#BaseAbstractUnivariateSolver(double,double,double)</em>
+<a class="jxr_linenumber" name="L83" href="#L83">83</a>  <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L84" href="#L84">84</a>      <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/BrentSolver.html">BrentSolver</a>(<strong class="jxr_keyword">double</strong> relativeAccuracy,
+<a class="jxr_linenumber" name="L85" href="#L85">85</a>                         <strong class="jxr_keyword">double</strong> absoluteAccuracy,
+<a class="jxr_linenumber" name="L86" href="#L86">86</a>                         <strong class="jxr_keyword">double</strong> functionValueAccuracy) {
+<a class="jxr_linenumber" name="L87" href="#L87">87</a>          <strong class="jxr_keyword">super</strong>(relativeAccuracy, absoluteAccuracy, functionValueAccuracy);
+<a class="jxr_linenumber" name="L88" href="#L88">88</a>      }
+<a class="jxr_linenumber" name="L89" href="#L89">89</a>  
+<a class="jxr_linenumber" name="L90" href="#L90">90</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L91" href="#L91">91</a>  <em class="jxr_javadoccomment">     * {@inheritDoc}</em>
+<a class="jxr_linenumber" name="L92" href="#L92">92</a>  <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L93" href="#L93">93</a>      @Override
+<a class="jxr_linenumber" name="L94" href="#L94">94</a>      <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> doSolve()
+<a class="jxr_linenumber" name="L95" href="#L95">95</a>          <strong class="jxr_keyword">throws</strong> NoBracketingException,
+<a class="jxr_linenumber" name="L96" href="#L96">96</a>                 TooManyEvaluationsException,
+<a class="jxr_linenumber" name="L97" href="#L97">97</a>                 <a href="../../../../../../org/apache/commons/math3/exception/NumberIsTooLargeException.html">NumberIsTooLargeException</a> {
+<a class="jxr_linenumber" name="L98" href="#L98">98</a>          <strong class="jxr_keyword">double</strong> min = getMin();
+<a class="jxr_linenumber" name="L99" href="#L99">99</a>          <strong class="jxr_keyword">double</strong> max = getMax();
+<a class="jxr_linenumber" name="L100" href="#L100">100</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> initial = getStartValue();
+<a class="jxr_linenumber" name="L101" href="#L101">101</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> functionValueAccuracy = getFunctionValueAccuracy();
+<a class="jxr_linenumber" name="L102" href="#L102">102</a> 
+<a class="jxr_linenumber" name="L103" href="#L103">103</a>         verifySequence(min, initial, max);
+<a class="jxr_linenumber" name="L104" href="#L104">104</a> 
+<a class="jxr_linenumber" name="L105" href="#L105">105</a>         <em class="jxr_comment">// Return the initial guess if it is good enough.</em>
+<a class="jxr_linenumber" name="L106" href="#L106">106</a>         <strong class="jxr_keyword">double</strong> yInitial = computeObjectiveValue(initial);
+<a class="jxr_linenumber" name="L107" href="#L107">107</a>         <strong class="jxr_keyword">if</strong> (FastMath.abs(yInitial) &lt;= functionValueAccuracy) {
+<a class="jxr_linenumber" name="L108" href="#L108">108</a>             <strong class="jxr_keyword">return</strong> initial;
 <a class="jxr_linenumber" name="L109" href="#L109">109</a>         }
 <a class="jxr_linenumber" name="L110" href="#L110">110</a> 
-<a class="jxr_linenumber" name="L111" href="#L111">111</a>         <em class="jxr_comment">// Return the second endpoint if it is good enough.</em>
-<a class="jxr_linenumber" name="L112" href="#L112">112</a>         <strong class="jxr_keyword">double</strong> yMax = computeObjectiveValue(max);
-<a class="jxr_linenumber" name="L113" href="#L113">113</a>         <strong class="jxr_keyword">if</strong> (FastMath.abs(yMax) &lt;= functionValueAccuracy) {
-<a class="jxr_linenumber" name="L114" href="#L114">114</a>             <strong class="jxr_keyword">return</strong> max;
+<a class="jxr_linenumber" name="L111" href="#L111">111</a>         <em class="jxr_comment">// Return the first endpoint if it is good enough.</em>
+<a class="jxr_linenumber" name="L112" href="#L112">112</a>         <strong class="jxr_keyword">double</strong> yMin = computeObjectiveValue(min);
+<a class="jxr_linenumber" name="L113" href="#L113">113</a>         <strong class="jxr_keyword">if</strong> (FastMath.abs(yMin) &lt;= functionValueAccuracy) {
+<a class="jxr_linenumber" name="L114" href="#L114">114</a>             <strong class="jxr_keyword">return</strong> min;
 <a class="jxr_linenumber" name="L115" href="#L115">115</a>         }
 <a class="jxr_linenumber" name="L116" href="#L116">116</a> 
-<a class="jxr_linenumber" name="L117" href="#L117">117</a>         <em class="jxr_comment">// Reduce interval if initial and max bracket the root.</em>
-<a class="jxr_linenumber" name="L118" href="#L118">118</a>         <strong class="jxr_keyword">if</strong> (yInitial * yMax &lt; 0) {
-<a class="jxr_linenumber" name="L119" href="#L119">119</a>             <strong class="jxr_keyword">return</strong> brent(initial, max, yInitial, yMax);
+<a class="jxr_linenumber" name="L117" href="#L117">117</a>         <em class="jxr_comment">// Reduce interval if min and initial bracket the root.</em>
+<a class="jxr_linenumber" name="L118" href="#L118">118</a>         <strong class="jxr_keyword">if</strong> (yInitial * yMin &lt; 0) {
+<a class="jxr_linenumber" name="L119" href="#L119">119</a>             <strong class="jxr_keyword">return</strong> brent(min, initial, yMin, yInitial);
 <a class="jxr_linenumber" name="L120" href="#L120">120</a>         }
 <a class="jxr_linenumber" name="L121" href="#L121">121</a> 
-<a class="jxr_linenumber" name="L122" href="#L122">122</a>         <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NoBracketingException.html">NoBracketingException</a>(min, max, yMin, yMax);
-<a class="jxr_linenumber" name="L123" href="#L123">123</a>     }
-<a class="jxr_linenumber" name="L124" href="#L124">124</a> 
-<a class="jxr_linenumber" name="L125" href="#L125">125</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L126" href="#L126">126</a> <em class="jxr_javadoccomment">     * Search for a zero inside the provided interval.</em>
-<a class="jxr_linenumber" name="L127" href="#L127">127</a> <em class="jxr_javadoccomment">     * This implementation is based on the algorithm described at page 58 of</em>
-<a class="jxr_linenumber" name="L128" href="#L128">128</a> <em class="jxr_javadoccomment">     * the book</em>
-<a class="jxr_linenumber" name="L129" href="#L129">129</a> <em class="jxr_javadoccomment">     * &lt;quote&gt;</em>
-<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_javadoccomment">     *  &lt;b&gt;Algorithms for Minimization Without Derivatives&lt;/b&gt;</em>
-<a class="jxr_linenumber" name="L131" href="#L131">131</a> <em class="jxr_javadoccomment">     *  &lt;it&gt;Richard P. Brent&lt;/it&gt;</em>
-<a class="jxr_linenumber" name="L132" href="#L132">132</a> <em class="jxr_javadoccomment">     *  Dover 0-486-41998-3</em>
-<a class="jxr_linenumber" name="L133" href="#L133">133</a> <em class="jxr_javadoccomment">     * &lt;/quote&gt;</em>
-<a class="jxr_linenumber" name="L134" href="#L134">134</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_javadoccomment">     * @param lo Lower bound of the search interval.</em>
-<a class="jxr_linenumber" name="L136" href="#L136">136</a> <em class="jxr_javadoccomment">     * @param hi Higher bound of the search interval.</em>
-<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment">     * @param fLo Function value at the lower bound of the search interval.</em>
-<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment">     * @param fHi Function value at the higher bound of the search interval.</em>
-<a class="jxr_linenumber" name="L139" href="#L139">139</a> <em class="jxr_javadoccomment">     * @return the value where the function is zero.</em>
-<a class="jxr_linenumber" name="L140" href="#L140">140</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L141" href="#L141">141</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> brent(<strong class="jxr_keyword">double</strong> lo, <strong class="jxr_keyword">double</strong> hi,
-<a class="jxr_linenumber" name="L142" href="#L142">142</a>                          <strong class="jxr_keyword">double</strong> fLo, <strong class="jxr_keyword">double</strong> fHi) {
-<a class="jxr_linenumber" name="L143" href="#L143">143</a>         <strong class="jxr_keyword">double</strong> a = lo;
-<a class="jxr_linenumber" name="L144" href="#L144">144</a>         <strong class="jxr_keyword">double</strong> fa = fLo;
-<a class="jxr_linenumber" name="L145" href="#L145">145</a>         <strong class="jxr_keyword">double</strong> b = hi;
-<a class="jxr_linenumber" name="L146" href="#L146">146</a>         <strong class="jxr_keyword">double</strong> fb = fHi;
-<a class="jxr_linenumber" name="L147" href="#L147">147</a>         <strong class="jxr_keyword">double</strong> c = a;
-<a class="jxr_linenumber" name="L148" href="#L148">148</a>         <strong class="jxr_keyword">double</strong> fc = fa;
-<a class="jxr_linenumber" name="L149" href="#L149">149</a>         <strong class="jxr_keyword">double</strong> d = b - a;
-<a class="jxr_linenumber" name="L150" href="#L150">150</a>         <strong class="jxr_keyword">double</strong> e = d;
-<a class="jxr_linenumber" name="L151" href="#L151">151</a> 
-<a class="jxr_linenumber" name="L152" href="#L152">152</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> t = getAbsoluteAccuracy();
-<a class="jxr_linenumber" name="L153" href="#L153">153</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> eps = getRelativeAccuracy();
-<a class="jxr_linenumber" name="L154" href="#L154">154</a> 
-<a class="jxr_linenumber" name="L155" href="#L155">155</a>         <strong class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>) {
-<a class="jxr_linenumber" name="L156" href="#L156">156</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(fc) &lt; FastMath.abs(fb)) {
-<a class="jxr_linenumber" name="L157" href="#L157">157</a>                 a = b;
-<a class="jxr_linenumber" name="L158" href="#L158">158</a>                 b = c;
-<a class="jxr_linenumber" name="L159" href="#L159">159</a>                 c = a;
-<a class="jxr_linenumber" name="L160" href="#L160">160</a>                 fa = fb;
-<a class="jxr_linenumber" name="L161" href="#L161">161</a>                 fb = fc;
-<a class="jxr_linenumber" name="L162" href="#L162">162</a>                 fc = fa;
-<a class="jxr_linenumber" name="L163" href="#L163">163</a>             }
-<a class="jxr_linenumber" name="L164" href="#L164">164</a> 
-<a class="jxr_linenumber" name="L165" href="#L165">165</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol = 2 * eps * FastMath.abs(b) + t;
-<a class="jxr_linenumber" name="L166" href="#L166">166</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> m = 0.5 * (c - b);
-<a class="jxr_linenumber" name="L167" href="#L167">167</a> 
-<a class="jxr_linenumber" name="L168" href="#L168">168</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(m) &lt;= tol ||
-<a class="jxr_linenumber" name="L169" href="#L169">169</a>                 Precision.equals(fb, 0))  {
-<a class="jxr_linenumber" name="L170" href="#L170">170</a>                 <strong class="jxr_keyword">return</strong> b;
-<a class="jxr_linenumber" name="L171" href="#L171">171</a>             }
-<a class="jxr_linenumber" name="L172" href="#L172">172</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(e) &lt; tol ||
-<a class="jxr_linenumber" name="L173" href="#L173">173</a>                 FastMath.abs(fa) &lt;= FastMath.abs(fb)) {
-<a class="jxr_linenumber" name="L174" href="#L174">174</a>                 <em class="jxr_comment">// Force bisection.</em>
-<a class="jxr_linenumber" name="L175" href="#L175">175</a>                 d = m;
-<a class="jxr_linenumber" name="L176" href="#L176">176</a>                 e = d;
-<a class="jxr_linenumber" name="L177" href="#L177">177</a>             } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L178" href="#L178">178</a>                 <strong class="jxr_keyword">double</strong> s = fb / fa;
-<a class="jxr_linenumber" name="L179" href="#L179">179</a>                 <strong class="jxr_keyword">double</strong> p;
-<a class="jxr_linenumber" name="L180" href="#L180">180</a>                 <strong class="jxr_keyword">double</strong> q;
-<a class="jxr_linenumber" name="L181" href="#L181">181</a>                 <em class="jxr_comment">// The equality test (a == c) is intentional,</em>
-<a class="jxr_linenumber" name="L182" href="#L182">182</a>                 <em class="jxr_comment">// it is part of the original Brent's method and</em>
-<a class="jxr_linenumber" name="L183" href="#L183">183</a>                 <em class="jxr_comment">// it should NOT be replaced by proximity test.</em>
-<a class="jxr_linenumber" name="L184" href="#L184">184</a>                 <strong class="jxr_keyword">if</strong> (a == c) {
-<a class="jxr_linenumber" name="L185" href="#L185">185</a>                     <em class="jxr_comment">// Linear interpolation.</em>
-<a class="jxr_linenumber" name="L186" href="#L186">186</a>                     p = 2 * m * s;
-<a class="jxr_linenumber" name="L187" href="#L187">187</a>                     q = 1 - s;
-<a class="jxr_linenumber" name="L188" href="#L188">188</a>                 } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L189" href="#L189">189</a>                     <em class="jxr_comment">// Inverse quadratic interpolation.</em>
-<a class="jxr_linenumber" name="L190" href="#L190">190</a>                     q = fa / fc;
-<a class="jxr_linenumber" name="L191" href="#L191">191</a>                     <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> r = fb / fc;
-<a class="jxr_linenumber" name="L192" href="#L192">192</a>                     p = s * (2 * m * q * (q - r) - (b - a) * (r - 1));
-<a class="jxr_linenumber" name="L193" href="#L193">193</a>                     q = (q - 1) * (r - 1) * (s - 1);
-<a class="jxr_linenumber" name="L194" href="#L194">194</a>                 }
-<a class="jxr_linenumber" name="L195" href="#L195">195</a>                 <strong class="jxr_keyword">if</strong> (p &gt; 0) {
-<a class="jxr_linenumber" name="L196" href="#L196">196</a>                     q = -q;
-<a class="jxr_linenumber" name="L197" href="#L197">197</a>                 } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L198" href="#L198">198</a>                     p = -p;
-<a class="jxr_linenumber" name="L199" href="#L199">199</a>                 }
-<a class="jxr_linenumber" name="L200" href="#L200">200</a>                 s = e;
-<a class="jxr_linenumber" name="L201" href="#L201">201</a>                 e = d;
-<a class="jxr_linenumber" name="L202" href="#L202">202</a>                 <strong class="jxr_keyword">if</strong> (p &gt;= 1.5 * m * q - FastMath.abs(tol * q) ||
-<a class="jxr_linenumber" name="L203" href="#L203">203</a>                     p &gt;= FastMath.abs(0.5 * s * q)) {
-<a class="jxr_linenumber" name="L204" href="#L204">204</a>                     <em class="jxr_comment">// Inverse quadratic interpolation gives a value</em>
-<a class="jxr_linenumber" name="L205" href="#L205">205</a>                     <em class="jxr_comment">// in the wrong direction, or progress is slow.</em>
-<a class="jxr_linenumber" name="L206" href="#L206">206</a>                     <em class="jxr_comment">// Fall back to bisection.</em>
-<a class="jxr_linenumber" name="L207" href="#L207">207</a>                     d = m;
-<a class="jxr_linenumber" name="L208" href="#L208">208</a>                     e = d;
-<a class="jxr_linenumber" name="L209" href="#L209">209</a>                 } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L210" href="#L210">210</a>                     d = p / q;
-<a class="jxr_linenumber" name="L211" href="#L211">211</a>                 }
-<a class="jxr_linenumber" name="L212" href="#L212">212</a>             }
-<a class="jxr_linenumber" name="L213" href="#L213">213</a>             a = b;
-<a class="jxr_linenumber" name="L214" href="#L214">214</a>             fa = fb;
-<a class="jxr_linenumber" name="L215" href="#L215">215</a> 
-<a class="jxr_linenumber" name="L216" href="#L216">216</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(d) &gt; tol) {
-<a class="jxr_linenumber" name="L217" href="#L217">217</a>                 b += d;
-<a class="jxr_linenumber" name="L218" href="#L218">218</a>             } <strong class="jxr_keyword">else</strong> <strong class="jxr_keyword">if</strong> (m &gt; 0) {
-<a class="jxr_linenumber" name="L219" href="#L219">219</a>                 b += tol;
-<a class="jxr_linenumber" name="L220" href="#L220">220</a>             } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L221" href="#L221">221</a>                 b -= tol;
-<a class="jxr_linenumber" name="L222" href="#L222">222</a>             }
-<a class="jxr_linenumber" name="L223" href="#L223">223</a>             fb = computeObjectiveValue(b);
-<a class="jxr_linenumber" name="L224" href="#L224">224</a>             <strong class="jxr_keyword">if</strong> ((fb &gt; 0 &amp;&amp; fc &gt; 0) ||
-<a class="jxr_linenumber" name="L225" href="#L225">225</a>                 (fb &lt;= 0 &amp;&amp; fc &lt;= 0)) {
-<a class="jxr_linenumber" name="L226" href="#L226">226</a>                 c = a;
-<a class="jxr_linenumber" name="L227" href="#L227">227</a>                 fc = fa;
-<a class="jxr_linenumber" name="L228" href="#L228">228</a>                 d = b - a;
-<a class="jxr_linenumber" name="L229" href="#L229">229</a>                 e = d;
-<a class="jxr_linenumber" name="L230" href="#L230">230</a>             }
-<a class="jxr_linenumber" name="L231" href="#L231">231</a>         }
-<a class="jxr_linenumber" name="L232" href="#L232">232</a>     }
-<a class="jxr_linenumber" name="L233" href="#L233">233</a> }
+<a class="jxr_linenumber" name="L122" href="#L122">122</a>         <em class="jxr_comment">// Return the second endpoint if it is good enough.</em>
+<a class="jxr_linenumber" name="L123" href="#L123">123</a>         <strong class="jxr_keyword">double</strong> yMax = computeObjectiveValue(max);
+<a class="jxr_linenumber" name="L124" href="#L124">124</a>         <strong class="jxr_keyword">if</strong> (FastMath.abs(yMax) &lt;= functionValueAccuracy) {
+<a class="jxr_linenumber" name="L125" href="#L125">125</a>             <strong class="jxr_keyword">return</strong> max;
+<a class="jxr_linenumber" name="L126" href="#L126">126</a>         }
+<a class="jxr_linenumber" name="L127" href="#L127">127</a> 
+<a class="jxr_linenumber" name="L128" href="#L128">128</a>         <em class="jxr_comment">// Reduce interval if initial and max bracket the root.</em>
+<a class="jxr_linenumber" name="L129" href="#L129">129</a>         <strong class="jxr_keyword">if</strong> (yInitial * yMax &lt; 0) {
+<a class="jxr_linenumber" name="L130" href="#L130">130</a>             <strong class="jxr_keyword">return</strong> brent(initial, max, yInitial, yMax);
+<a class="jxr_linenumber" name="L131" href="#L131">131</a>         }
+<a class="jxr_linenumber" name="L132" href="#L132">132</a> 
+<a class="jxr_linenumber" name="L133" href="#L133">133</a>         <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NoBracketingException.html">NoBracketingException</a>(min, max, yMin, yMax);
+<a class="jxr_linenumber" name="L134" href="#L134">134</a>     }
+<a class="jxr_linenumber" name="L135" href="#L135">135</a> 
+<a class="jxr_linenumber" name="L136" href="#L136">136</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment">     * Search for a zero inside the provided interval.</em>
+<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment">     * This implementation is based on the algorithm described at page 58 of</em>
+<a class="jxr_linenumber" name="L139" href="#L139">139</a> <em class="jxr_javadoccomment">     * the book</em>
+<a class="jxr_linenumber" name="L140" href="#L140">140</a> <em class="jxr_javadoccomment">     * &lt;blockquote&gt;</em>
+<a class="jxr_linenumber" name="L141" href="#L141">141</a> <em class="jxr_javadoccomment">     *  &lt;b&gt;Algorithms for Minimization Without Derivatives&lt;/b&gt;</em>
+<a class="jxr_linenumber" name="L142" href="#L142">142</a> <em class="jxr_javadoccomment">     *  &lt;it&gt;Richard P. Brent&lt;/it&gt;</em>
+<a class="jxr_linenumber" name="L143" href="#L143">143</a> <em class="jxr_javadoccomment">     *  Dover 0-486-41998-3</em>
+<a class="jxr_linenumber" name="L144" href="#L144">144</a> <em class="jxr_javadoccomment">     * &lt;/blockquote&gt;</em>
+<a class="jxr_linenumber" name="L145" href="#L145">145</a> <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L146" href="#L146">146</a> <em class="jxr_javadoccomment">     * @param lo Lower bound of the search interval.</em>
+<a class="jxr_linenumber" name="L147" href="#L147">147</a> <em class="jxr_javadoccomment">     * @param hi Higher bound of the search interval.</em>
+<a class="jxr_linenumber" name="L148" href="#L148">148</a> <em class="jxr_javadoccomment">     * @param fLo Function value at the lower bound of the search interval.</em>
+<a class="jxr_linenumber" name="L149" href="#L149">149</a> <em class="jxr_javadoccomment">     * @param fHi Function value at the higher bound of the search interval.</em>
+<a class="jxr_linenumber" name="L150" href="#L150">150</a> <em class="jxr_javadoccomment">     * @return the value where the function is zero.</em>
+<a class="jxr_linenumber" name="L151" href="#L151">151</a> <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L152" href="#L152">152</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> brent(<strong class="jxr_keyword">double</strong> lo, <strong class="jxr_keyword">double</strong> hi,
+<a class="jxr_linenumber" name="L153" href="#L153">153</a>                          <strong class="jxr_keyword">double</strong> fLo, <strong class="jxr_keyword">double</strong> fHi) {
+<a class="jxr_linenumber" name="L154" href="#L154">154</a>         <strong class="jxr_keyword">double</strong> a = lo;
+<a class="jxr_linenumber" name="L155" href="#L155">155</a>         <strong class="jxr_keyword">double</strong> fa = fLo;
+<a class="jxr_linenumber" name="L156" href="#L156">156</a>         <strong class="jxr_keyword">double</strong> b = hi;
+<a class="jxr_linenumber" name="L157" href="#L157">157</a>         <strong class="jxr_keyword">double</strong> fb = fHi;
+<a class="jxr_linenumber" name="L158" href="#L158">158</a>         <strong class="jxr_keyword">double</strong> c = a;
+<a class="jxr_linenumber" name="L159" href="#L159">159</a>         <strong class="jxr_keyword">double</strong> fc = fa;
+<a class="jxr_linenumber" name="L160" href="#L160">160</a>         <strong class="jxr_keyword">double</strong> d = b - a;
+<a class="jxr_linenumber" name="L161" href="#L161">161</a>         <strong class="jxr_keyword">double</strong> e = d;
+<a class="jxr_linenumber" name="L162" href="#L162">162</a> 
+<a class="jxr_linenumber" name="L163" href="#L163">163</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> t = getAbsoluteAccuracy();
+<a class="jxr_linenumber" name="L164" href="#L164">164</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> eps = getRelativeAccuracy();
+<a class="jxr_linenumber" name="L165" href="#L165">165</a> 
+<a class="jxr_linenumber" name="L166" href="#L166">166</a>         <strong class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>) {
+<a class="jxr_linenumber" name="L167" href="#L167">167</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(fc) &lt; FastMath.abs(fb)) {
+<a class="jxr_linenumber" name="L168" href="#L168">168</a>                 a = b;
+<a class="jxr_linenumber" name="L169" href="#L169">169</a>                 b = c;
+<a class="jxr_linenumber" name="L170" href="#L170">170</a>                 c = a;
+<a class="jxr_linenumber" name="L171" href="#L171">171</a>                 fa = fb;
+<a class="jxr_linenumber" name="L172" href="#L172">172</a>                 fb = fc;
+<a class="jxr_linenumber" name="L173" href="#L173">173</a>                 fc = fa;
+<a class="jxr_linenumber" name="L174" href="#L174">174</a>             }
+<a class="jxr_linenumber" name="L175" href="#L175">175</a> 
+<a class="jxr_linenumber" name="L176" href="#L176">176</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol = 2 * eps * FastMath.abs(b) + t;
+<a class="jxr_linenumber" name="L177" href="#L177">177</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> m = 0.5 * (c - b);
+<a class="jxr_linenumber" name="L178" href="#L178">178</a> 
+<a class="jxr_linenumber" name="L179" href="#L179">179</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(m) &lt;= tol ||
+<a class="jxr_linenumber" name="L180" href="#L180">180</a>                 Precision.equals(fb, 0))  {
+<a class="jxr_linenumber" name="L181" href="#L181">181</a>                 <strong class="jxr_keyword">return</strong> b;
+<a class="jxr_linenumber" name="L182" href="#L182">182</a>             }
+<a class="jxr_linenumber" name="L183" href="#L183">183</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(e) &lt; tol ||
+<a class="jxr_linenumber" name="L184" href="#L184">184</a>                 FastMath.abs(fa) &lt;= FastMath.abs(fb)) {
+<a class="jxr_linenumber" name="L185" href="#L185">185</a>                 <em class="jxr_comment">// Force bisection.</em>
+<a class="jxr_linenumber" name="L186" href="#L186">186</a>                 d = m;
+<a class="jxr_linenumber" name="L187" href="#L187">187</a>                 e = d;
+<a class="jxr_linenumber" name="L188" href="#L188">188</a>             } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L189" href="#L189">189</a>                 <strong class="jxr_keyword">double</strong> s = fb / fa;
+<a class="jxr_linenumber" name="L190" href="#L190">190</a>                 <strong class="jxr_keyword">double</strong> p;
+<a class="jxr_linenumber" name="L191" href="#L191">191</a>                 <strong class="jxr_keyword">double</strong> q;
+<a class="jxr_linenumber" name="L192" href="#L192">192</a>                 <em class="jxr_comment">// The equality test (a == c) is intentional,</em>
+<a class="jxr_linenumber" name="L193" href="#L193">193</a>                 <em class="jxr_comment">// it is part of the original Brent's method and</em>
+<a class="jxr_linenumber" name="L194" href="#L194">194</a>                 <em class="jxr_comment">// it should NOT be replaced by proximity test.</em>
+<a class="jxr_linenumber" name="L195" href="#L195">195</a>                 <strong class="jxr_keyword">if</strong> (a == c) {
+<a class="jxr_linenumber" name="L196" href="#L196">196</a>                     <em class="jxr_comment">// Linear interpolation.</em>
+<a class="jxr_linenumber" name="L197" href="#L197">197</a>                     p = 2 * m * s;
+<a class="jxr_linenumber" name="L198" href="#L198">198</a>                     q = 1 - s;
+<a class="jxr_linenumber" name="L199" href="#L199">199</a>                 } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L200" href="#L200">200</a>                     <em class="jxr_comment">// Inverse quadratic interpolation.</em>
+<a class="jxr_linenumber" name="L201" href="#L201">201</a>                     q = fa / fc;
+<a class="jxr_linenumber" name="L202" href="#L202">202</a>                     <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> r = fb / fc;
+<a class="jxr_linenumber" name="L203" href="#L203">203</a>                     p = s * (2 * m * q * (q - r) - (b - a) * (r - 1));
+<a class="jxr_linenumber" name="L204" href="#L204">204</a>                     q = (q - 1) * (r - 1) * (s - 1);
+<a class="jxr_linenumber" name="L205" href="#L205">205</a>                 }
+<a class="jxr_linenumber" name="L206" href="#L206">206</a>                 <strong class="jxr_keyword">if</strong> (p &gt; 0) {
+<a class="jxr_linenumber" name="L207" href="#L207">207</a>                     q = -q;
+<a class="jxr_linenumber" name="L208" href="#L208">208</a>                 } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L209" href="#L209">209</a>                     p = -p;
+<a class="jxr_linenumber" name="L210" href="#L210">210</a>                 }
+<a class="jxr_linenumber" name="L211" href="#L211">211</a>                 s = e;
+<a class="jxr_linenumber" name="L212" href="#L212">212</a>                 e = d;
+<a class="jxr_linenumber" name="L213" href="#L213">213</a>                 <strong class="jxr_keyword">if</strong> (p &gt;= 1.5 * m * q - FastMath.abs(tol * q) ||
+<a class="jxr_linenumber" name="L214" href="#L214">214</a>                     p &gt;= FastMath.abs(0.5 * s * q)) {
+<a class="jxr_linenumber" name="L215" href="#L215">215</a>                     <em class="jxr_comment">// Inverse quadratic interpolation gives a value</em>
+<a class="jxr_linenumber" name="L216" href="#L216">216</a>                     <em class="jxr_comment">// in the wrong direction, or progress is slow.</em>
+<a class="jxr_linenumber" name="L217" href="#L217">217</a>                     <em class="jxr_comment">// Fall back to bisection.</em>
+<a class="jxr_linenumber" name="L218" href="#L218">218</a>                     d = m;
+<a class="jxr_linenumber" name="L219" href="#L219">219</a>                     e = d;
+<a class="jxr_linenumber" name="L220" href="#L220">220</a>                 } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L221" href="#L221">221</a>                     d = p / q;
+<a class="jxr_linenumber" name="L222" href="#L222">222</a>                 }
+<a class="jxr_linenumber" name="L223" href="#L223">223</a>             }
+<a class="jxr_linenumber" name="L224" href="#L224">224</a>             a = b;
+<a class="jxr_linenumber" name="L225" href="#L225">225</a>             fa = fb;
+<a class="jxr_linenumber" name="L226" href="#L226">226</a> 
+<a class="jxr_linenumber" name="L227" href="#L227">227</a>             <strong class="jxr_keyword">if</strong> (FastMath.abs(d) &gt; tol) {
+<a class="jxr_linenumber" name="L228" href="#L228">228</a>                 b += d;
+<a class="jxr_linenumber" name="L229" href="#L229">229</a>             } <strong class="jxr_keyword">else</strong> <strong class="jxr_keyword">if</strong> (m &gt; 0) {
+<a class="jxr_linenumber" name="L230" href="#L230">230</a>                 b += tol;
+<a class="jxr_linenumber" name="L231" href="#L231">231</a>             } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L232" href="#L232">232</a>                 b -= tol;
+<a class="jxr_linenumber" name="L233" href="#L233">233</a>             }
+<a class="jxr_linenumber" name="L234" href="#L234">234</a>             fb = computeObjectiveValue(b);
+<a class="jxr_linenumber" name="L235" href="#L235">235</a>             <strong class="jxr_keyword">if</strong> ((fb &gt; 0 &amp;&amp; fc &gt; 0) ||
+<a class="jxr_linenumber" name="L236" href="#L236">236</a>                 (fb &lt;= 0 &amp;&amp; fc &lt;= 0)) {
+<a class="jxr_linenumber" name="L237" href="#L237">237</a>                 c = a;
+<a class="jxr_linenumber" name="L238" href="#L238">238</a>                 fc = fa;
+<a class="jxr_linenumber" name="L239" href="#L239">239</a>                 d = b - a;
+<a class="jxr_linenumber" name="L240" href="#L240">240</a>                 e = d;
+<a class="jxr_linenumber" name="L241" href="#L241">241</a>             }
+<a class="jxr_linenumber" name="L242" href="#L242">242</a>         }
+<a class="jxr_linenumber" name="L243" href="#L243">243</a>     }
+<a class="jxr_linenumber" name="L244" href="#L244">244</a> }
 </pre>
 <hr/>
 <div id="footer">Copyright &#169; 2003&#x2013;2015 <a href="http://www.apache.org/">The Apache Software Foundation</a>. All rights reserved.</div>



Mime
View raw message