You are viewing a plain text version of this content. The canonical link for it is here.
Posted to commits@commons.apache.org by lu...@apache.org on 2013/04/07 09:52:13 UTC
svn commit: r857590 [28/48] - in
/websites/production/commons/content/proper/commons-math/xref-test: ./
org/apache/commons/math3/ org/apache/commons/math3/analysis/
org/apache/commons/math3/analysis/differentiation/
org/apache/commons/math3/analysis/fu...
Modified: websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html (original)
+++ websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html Sun Apr 7 07:52:05 2013
@@ -29,213 +29,241 @@
<a class="jxr_linenumber" name="19" href="#19">19</a>
<a class="jxr_linenumber" name="20" href="#20">20</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.analysis.MultivariateFunction;
<a class="jxr_linenumber" name="21" href="#21">21</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.MaxEval;
-<a class="jxr_linenumber" name="22" href="#22">22</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.nonlinear.scalar.GoalType;
-<a class="jxr_linenumber" name="23" href="#23">23</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.InitialGuess;
-<a class="jxr_linenumber" name="24" href="#24">24</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.PointValuePair;
-<a class="jxr_linenumber" name="25" href="#25">25</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.SimpleValueChecker;
-<a class="jxr_linenumber" name="26" href="#26">26</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.nonlinear.scalar.ObjectiveFunction;
-<a class="jxr_linenumber" name="27" href="#27">27</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.util.FastMath;
-<a class="jxr_linenumber" name="28" href="#28">28</a> <strong class="jxr_keyword">import</strong> org.junit.Assert;
-<a class="jxr_linenumber" name="29" href="#29">29</a> <strong class="jxr_keyword">import</strong> org.junit.Test;
-<a class="jxr_linenumber" name="30" href="#30">30</a>
-<a class="jxr_linenumber" name="31" href="#31">31</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">SimplexOptimizerMultiDirectionalTest</a> {
-<a class="jxr_linenumber" name="32" href="#32">32</a> @Test
-<a class="jxr_linenumber" name="33" href="#33">33</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMinimize1() {
-<a class="jxr_linenumber" name="34" href="#34">34</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
-<a class="jxr_linenumber" name="35" href="#35">35</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a>();
-<a class="jxr_linenumber" name="36" href="#36">36</a>
-<a class="jxr_linenumber" name="37" href="#37">37</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
-<a class="jxr_linenumber" name="38" href="#38">38</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
-<a class="jxr_linenumber" name="39" href="#39">39</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
-<a class="jxr_linenumber" name="40" href="#40">40</a> GoalType.MINIMIZE,
-<a class="jxr_linenumber" name="41" href="#41">41</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -3, 0 }),
-<a class="jxr_linenumber" name="42" href="#42">42</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
-<a class="jxr_linenumber" name="43" href="#43">43</a> Assert.assertEquals(fourExtrema.xM, optimum.getPoint()[0], 4e-6);
-<a class="jxr_linenumber" name="44" href="#44">44</a> Assert.assertEquals(fourExtrema.yP, optimum.getPoint()[1], 3e-6);
-<a class="jxr_linenumber" name="45" href="#45">45</a> Assert.assertEquals(fourExtrema.valueXmYp, optimum.getValue(), 8e-13);
-<a class="jxr_linenumber" name="46" href="#46">46</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
-<a class="jxr_linenumber" name="47" href="#47">47</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
-<a class="jxr_linenumber" name="48" href="#48">48</a> }
-<a class="jxr_linenumber" name="49" href="#49">49</a>
-<a class="jxr_linenumber" name="50" href="#50">50</a> @Test
-<a class="jxr_linenumber" name="51" href="#51">51</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMinimize2() {
-<a class="jxr_linenumber" name="52" href="#52">52</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
-<a class="jxr_linenumber" name="53" href="#53">53</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a>();
-<a class="jxr_linenumber" name="54" href="#54">54</a>
-<a class="jxr_linenumber" name="55" href="#55">55</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
-<a class="jxr_linenumber" name="56" href="#56">56</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
-<a class="jxr_linenumber" name="57" href="#57">57</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
-<a class="jxr_linenumber" name="58" href="#58">58</a> GoalType.MINIMIZE,
-<a class="jxr_linenumber" name="59" href="#59">59</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 1, 0 }),
-<a class="jxr_linenumber" name="60" href="#60">60</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
-<a class="jxr_linenumber" name="61" href="#61">61</a> Assert.assertEquals(fourExtrema.xP, optimum.getPoint()[0], 2e-8);
-<a class="jxr_linenumber" name="62" href="#62">62</a> Assert.assertEquals(fourExtrema.yM, optimum.getPoint()[1], 3e-6);
-<a class="jxr_linenumber" name="63" href="#63">63</a> Assert.assertEquals(fourExtrema.valueXpYm, optimum.getValue(), 2e-12);
-<a class="jxr_linenumber" name="64" href="#64">64</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
-<a class="jxr_linenumber" name="65" href="#65">65</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
-<a class="jxr_linenumber" name="66" href="#66">66</a> }
-<a class="jxr_linenumber" name="67" href="#67">67</a>
-<a class="jxr_linenumber" name="68" href="#68">68</a> @Test
-<a class="jxr_linenumber" name="69" href="#69">69</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMaximize1() {
-<a class="jxr_linenumber" name="70" href="#70">70</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
-<a class="jxr_linenumber" name="71" href="#71">71</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a>();
-<a class="jxr_linenumber" name="72" href="#72">72</a>
-<a class="jxr_linenumber" name="73" href="#73">73</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
-<a class="jxr_linenumber" name="74" href="#74">74</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
-<a class="jxr_linenumber" name="75" href="#75">75</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
-<a class="jxr_linenumber" name="76" href="#76">76</a> GoalType.MAXIMIZE,
-<a class="jxr_linenumber" name="77" href="#77">77</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -3.0, 0.0 }),
-<a class="jxr_linenumber" name="78" href="#78">78</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
-<a class="jxr_linenumber" name="79" href="#79">79</a> Assert.assertEquals(fourExtrema.xM, optimum.getPoint()[0], 7e-7);
-<a class="jxr_linenumber" name="80" href="#80">80</a> Assert.assertEquals(fourExtrema.yM, optimum.getPoint()[1], 3e-7);
-<a class="jxr_linenumber" name="81" href="#81">81</a> Assert.assertEquals(fourExtrema.valueXmYm, optimum.getValue(), 2e-14);
-<a class="jxr_linenumber" name="82" href="#82">82</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
-<a class="jxr_linenumber" name="83" href="#83">83</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
-<a class="jxr_linenumber" name="84" href="#84">84</a> }
+<a class="jxr_linenumber" name="22" href="#22">22</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.InitialGuess;
+<a class="jxr_linenumber" name="23" href="#23">23</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.PointValuePair;
+<a class="jxr_linenumber" name="24" href="#24">24</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.SimpleValueChecker;
+<a class="jxr_linenumber" name="25" href="#25">25</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.SimpleBounds;
+<a class="jxr_linenumber" name="26" href="#26">26</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.nonlinear.scalar.GoalType;
+<a class="jxr_linenumber" name="27" href="#27">27</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.optim.nonlinear.scalar.ObjectiveFunction;
+<a class="jxr_linenumber" name="28" href="#28">28</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.MathUnsupportedOperationException;
+<a class="jxr_linenumber" name="29" href="#29">29</a> <strong class="jxr_keyword">import</strong> org.apache.commons.math3.util.FastMath;
+<a class="jxr_linenumber" name="30" href="#30">30</a> <strong class="jxr_keyword">import</strong> org.junit.Assert;
+<a class="jxr_linenumber" name="31" href="#31">31</a> <strong class="jxr_keyword">import</strong> org.junit.Test;
+<a class="jxr_linenumber" name="32" href="#32">32</a>
+<a class="jxr_linenumber" name="33" href="#33">33</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">SimplexOptimizerMultiDirectionalTest</a> {
+<a class="jxr_linenumber" name="34" href="#34">34</a> @Test(expected=MathUnsupportedOperationException.<strong class="jxr_keyword">class</strong>)
+<a class="jxr_linenumber" name="35" href="#35">35</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testBoundsUnsupported() {
+<a class="jxr_linenumber" name="36" href="#36">36</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-10, 1e-30);
+<a class="jxr_linenumber" name="37" href="#37">37</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a>();
+<a class="jxr_linenumber" name="38" href="#38">38</a>
+<a class="jxr_linenumber" name="39" href="#39">39</a> optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(100),
+<a class="jxr_linenumber" name="40" href="#40">40</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
+<a class="jxr_linenumber" name="41" href="#41">41</a> GoalType.MINIMIZE,
+<a class="jxr_linenumber" name="42" href="#42">42</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -3, 0 }),
+<a class="jxr_linenumber" name="43" href="#43">43</a> <strong class="jxr_keyword">new</strong> NelderMeadSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }),
+<a class="jxr_linenumber" name="44" href="#44">44</a> <strong class="jxr_keyword">new</strong> SimpleBounds(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -5, -1 },
+<a class="jxr_linenumber" name="45" href="#45">45</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 5, 1 }));
+<a class="jxr_linenumber" name="46" href="#46">46</a> }
+<a class="jxr_linenumber" name="47" href="#47">47</a>
+<a class="jxr_linenumber" name="48" href="#48">48</a> @Test
+<a class="jxr_linenumber" name="49" href="#49">49</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMinimize1() {
+<a class="jxr_linenumber" name="50" href="#50">50</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
+<a class="jxr_linenumber" name="51" href="#51">51</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a>();
+<a class="jxr_linenumber" name="52" href="#52">52</a>
+<a class="jxr_linenumber" name="53" href="#53">53</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
+<a class="jxr_linenumber" name="54" href="#54">54</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
+<a class="jxr_linenumber" name="55" href="#55">55</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
+<a class="jxr_linenumber" name="56" href="#56">56</a> GoalType.MINIMIZE,
+<a class="jxr_linenumber" name="57" href="#57">57</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -3, 0 }),
+<a class="jxr_linenumber" name="58" href="#58">58</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
+<a class="jxr_linenumber" name="59" href="#59">59</a> Assert.assertEquals(fourExtrema.xM, optimum.getPoint()[0], 4e-6);
+<a class="jxr_linenumber" name="60" href="#60">60</a> Assert.assertEquals(fourExtrema.yP, optimum.getPoint()[1], 3e-6);
+<a class="jxr_linenumber" name="61" href="#61">61</a> Assert.assertEquals(fourExtrema.valueXmYp, optimum.getValue(), 8e-13);
+<a class="jxr_linenumber" name="62" href="#62">62</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
+<a class="jxr_linenumber" name="63" href="#63">63</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
+<a class="jxr_linenumber" name="64" href="#64">64</a>
+<a class="jxr_linenumber" name="65" href="#65">65</a> <em class="jxr_comment">// Check that the number of iterations is updated (MATH-949).</em>
+<a class="jxr_linenumber" name="66" href="#66">66</a> Assert.assertTrue(optimizer.getIterations() > 0);
+<a class="jxr_linenumber" name="67" href="#67">67</a> }
+<a class="jxr_linenumber" name="68" href="#68">68</a>
+<a class="jxr_linenumber" name="69" href="#69">69</a> @Test
+<a class="jxr_linenumber" name="70" href="#70">70</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMinimize2() {
+<a class="jxr_linenumber" name="71" href="#71">71</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
+<a class="jxr_linenumber" name="72" href="#72">72</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a>();
+<a class="jxr_linenumber" name="73" href="#73">73</a>
+<a class="jxr_linenumber" name="74" href="#74">74</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
+<a class="jxr_linenumber" name="75" href="#75">75</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
+<a class="jxr_linenumber" name="76" href="#76">76</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
+<a class="jxr_linenumber" name="77" href="#77">77</a> GoalType.MINIMIZE,
+<a class="jxr_linenumber" name="78" href="#78">78</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 1, 0 }),
+<a class="jxr_linenumber" name="79" href="#79">79</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
+<a class="jxr_linenumber" name="80" href="#80">80</a> Assert.assertEquals(fourExtrema.xP, optimum.getPoint()[0], 2e-8);
+<a class="jxr_linenumber" name="81" href="#81">81</a> Assert.assertEquals(fourExtrema.yM, optimum.getPoint()[1], 3e-6);
+<a class="jxr_linenumber" name="82" href="#82">82</a> Assert.assertEquals(fourExtrema.valueXpYm, optimum.getValue(), 2e-12);
+<a class="jxr_linenumber" name="83" href="#83">83</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
+<a class="jxr_linenumber" name="84" href="#84">84</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
<a class="jxr_linenumber" name="85" href="#85">85</a>
-<a class="jxr_linenumber" name="86" href="#86">86</a> @Test
-<a class="jxr_linenumber" name="87" href="#87">87</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMaximize2() {
-<a class="jxr_linenumber" name="88" href="#88">88</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(<strong class="jxr_keyword">new</strong> SimpleValueChecker(1e-15, 1e-30));
-<a class="jxr_linenumber" name="89" href="#89">89</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a>();
-<a class="jxr_linenumber" name="90" href="#90">90</a>
-<a class="jxr_linenumber" name="91" href="#91">91</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
-<a class="jxr_linenumber" name="92" href="#92">92</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
-<a class="jxr_linenumber" name="93" href="#93">93</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
-<a class="jxr_linenumber" name="94" href="#94">94</a> GoalType.MAXIMIZE,
-<a class="jxr_linenumber" name="95" href="#95">95</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 1, 0 }),
-<a class="jxr_linenumber" name="96" href="#96">96</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
-<a class="jxr_linenumber" name="97" href="#97">97</a> Assert.assertEquals(fourExtrema.xP, optimum.getPoint()[0], 2e-8);
-<a class="jxr_linenumber" name="98" href="#98">98</a> Assert.assertEquals(fourExtrema.yP, optimum.getPoint()[1], 3e-6);
-<a class="jxr_linenumber" name="99" href="#99">99</a> Assert.assertEquals(fourExtrema.valueXpYp, optimum.getValue(), 2e-12);
-<a class="jxr_linenumber" name="100" href="#100">100</a> Assert.assertTrue(optimizer.getEvaluations() > 180);
-<a class="jxr_linenumber" name="101" href="#101">101</a> Assert.assertTrue(optimizer.getEvaluations() < 220);
-<a class="jxr_linenumber" name="102" href="#102">102</a> }
-<a class="jxr_linenumber" name="103" href="#103">103</a>
-<a class="jxr_linenumber" name="104" href="#104">104</a> @Test
-<a class="jxr_linenumber" name="105" href="#105">105</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testRosenbrock() {
-<a class="jxr_linenumber" name="106" href="#106">106</a> MultivariateFunction rosenbrock
-<a class="jxr_linenumber" name="107" href="#107">107</a> = <strong class="jxr_keyword">new</strong> MultivariateFunction() {
-<a class="jxr_linenumber" name="108" href="#108">108</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] x) {
-<a class="jxr_linenumber" name="109" href="#109">109</a> ++count;
-<a class="jxr_linenumber" name="110" href="#110">110</a> <strong class="jxr_keyword">double</strong> a = x[1] - x[0] * x[0];
-<a class="jxr_linenumber" name="111" href="#111">111</a> <strong class="jxr_keyword">double</strong> b = 1.0 - x[0];
-<a class="jxr_linenumber" name="112" href="#112">112</a> <strong class="jxr_keyword">return</strong> 100 * a * a + b * b;
-<a class="jxr_linenumber" name="113" href="#113">113</a> }
-<a class="jxr_linenumber" name="114" href="#114">114</a> };
+<a class="jxr_linenumber" name="86" href="#86">86</a> <em class="jxr_comment">// Check that the number of iterations is updated (MATH-949).</em>
+<a class="jxr_linenumber" name="87" href="#87">87</a> Assert.assertTrue(optimizer.getIterations() > 0);
+<a class="jxr_linenumber" name="88" href="#88">88</a> }
+<a class="jxr_linenumber" name="89" href="#89">89</a>
+<a class="jxr_linenumber" name="90" href="#90">90</a> @Test
+<a class="jxr_linenumber" name="91" href="#91">91</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMaximize1() {
+<a class="jxr_linenumber" name="92" href="#92">92</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-11, 1e-30);
+<a class="jxr_linenumber" name="93" href="#93">93</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a>();
+<a class="jxr_linenumber" name="94" href="#94">94</a>
+<a class="jxr_linenumber" name="95" href="#95">95</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
+<a class="jxr_linenumber" name="96" href="#96">96</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
+<a class="jxr_linenumber" name="97" href="#97">97</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
+<a class="jxr_linenumber" name="98" href="#98">98</a> GoalType.MAXIMIZE,
+<a class="jxr_linenumber" name="99" href="#99">99</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -3.0, 0.0 }),
+<a class="jxr_linenumber" name="100" href="#100">100</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
+<a class="jxr_linenumber" name="101" href="#101">101</a> Assert.assertEquals(fourExtrema.xM, optimum.getPoint()[0], 7e-7);
+<a class="jxr_linenumber" name="102" href="#102">102</a> Assert.assertEquals(fourExtrema.yM, optimum.getPoint()[1], 3e-7);
+<a class="jxr_linenumber" name="103" href="#103">103</a> Assert.assertEquals(fourExtrema.valueXmYm, optimum.getValue(), 2e-14);
+<a class="jxr_linenumber" name="104" href="#104">104</a> Assert.assertTrue(optimizer.getEvaluations() > 120);
+<a class="jxr_linenumber" name="105" href="#105">105</a> Assert.assertTrue(optimizer.getEvaluations() < 150);
+<a class="jxr_linenumber" name="106" href="#106">106</a>
+<a class="jxr_linenumber" name="107" href="#107">107</a> <em class="jxr_comment">// Check that the number of iterations is updated (MATH-949).</em>
+<a class="jxr_linenumber" name="108" href="#108">108</a> Assert.assertTrue(optimizer.getIterations() > 0);
+<a class="jxr_linenumber" name="109" href="#109">109</a> }
+<a class="jxr_linenumber" name="110" href="#110">110</a>
+<a class="jxr_linenumber" name="111" href="#111">111</a> @Test
+<a class="jxr_linenumber" name="112" href="#112">112</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMaximize2() {
+<a class="jxr_linenumber" name="113" href="#113">113</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(<strong class="jxr_keyword">new</strong> SimpleValueChecker(1e-15, 1e-30));
+<a class="jxr_linenumber" name="114" href="#114">114</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> fourExtrema = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a>();
<a class="jxr_linenumber" name="115" href="#115">115</a>
-<a class="jxr_linenumber" name="116" href="#116">116</a> count = 0;
-<a class="jxr_linenumber" name="117" href="#117">117</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(-1, 1e-3);
-<a class="jxr_linenumber" name="118" href="#118">118</a> PointValuePair optimum
-<a class="jxr_linenumber" name="119" href="#119">119</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(100),
-<a class="jxr_linenumber" name="120" href="#120">120</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(rosenbrock),
-<a class="jxr_linenumber" name="121" href="#121">121</a> GoalType.MINIMIZE,
-<a class="jxr_linenumber" name="122" href="#122">122</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -1.2, 1 }),
-<a class="jxr_linenumber" name="123" href="#123">123</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[][] {
-<a class="jxr_linenumber" name="124" href="#124">124</a> { -1.2, 1.0 },
-<a class="jxr_linenumber" name="125" href="#125">125</a> { 0.9, 1.2 },
-<a class="jxr_linenumber" name="126" href="#126">126</a> { 3.5, -2.3 } }));
+<a class="jxr_linenumber" name="116" href="#116">116</a> <strong class="jxr_keyword">final</strong> PointValuePair optimum
+<a class="jxr_linenumber" name="117" href="#117">117</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(200),
+<a class="jxr_linenumber" name="118" href="#118">118</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(fourExtrema),
+<a class="jxr_linenumber" name="119" href="#119">119</a> GoalType.MAXIMIZE,
+<a class="jxr_linenumber" name="120" href="#120">120</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 1, 0 }),
+<a class="jxr_linenumber" name="121" href="#121">121</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 0.2, 0.2 }));
+<a class="jxr_linenumber" name="122" href="#122">122</a> Assert.assertEquals(fourExtrema.xP, optimum.getPoint()[0], 2e-8);
+<a class="jxr_linenumber" name="123" href="#123">123</a> Assert.assertEquals(fourExtrema.yP, optimum.getPoint()[1], 3e-6);
+<a class="jxr_linenumber" name="124" href="#124">124</a> Assert.assertEquals(fourExtrema.valueXpYp, optimum.getValue(), 2e-12);
+<a class="jxr_linenumber" name="125" href="#125">125</a> Assert.assertTrue(optimizer.getEvaluations() > 180);
+<a class="jxr_linenumber" name="126" href="#126">126</a> Assert.assertTrue(optimizer.getEvaluations() < 220);
<a class="jxr_linenumber" name="127" href="#127">127</a>
-<a class="jxr_linenumber" name="128" href="#128">128</a> Assert.assertEquals(count, optimizer.getEvaluations());
-<a class="jxr_linenumber" name="129" href="#129">129</a> Assert.assertTrue(optimizer.getEvaluations() > 50);
-<a class="jxr_linenumber" name="130" href="#130">130</a> Assert.assertTrue(optimizer.getEvaluations() < 100);
-<a class="jxr_linenumber" name="131" href="#131">131</a> Assert.assertTrue(optimum.getValue() > 1e-2);
-<a class="jxr_linenumber" name="132" href="#132">132</a> }
-<a class="jxr_linenumber" name="133" href="#133">133</a>
-<a class="jxr_linenumber" name="134" href="#134">134</a> @Test
-<a class="jxr_linenumber" name="135" href="#135">135</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testPowell() {
-<a class="jxr_linenumber" name="136" href="#136">136</a> MultivariateFunction powell
-<a class="jxr_linenumber" name="137" href="#137">137</a> = <strong class="jxr_keyword">new</strong> MultivariateFunction() {
-<a class="jxr_linenumber" name="138" href="#138">138</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] x) {
-<a class="jxr_linenumber" name="139" href="#139">139</a> ++count;
-<a class="jxr_linenumber" name="140" href="#140">140</a> <strong class="jxr_keyword">double</strong> a = x[0] + 10 * x[1];
-<a class="jxr_linenumber" name="141" href="#141">141</a> <strong class="jxr_keyword">double</strong> b = x[2] - x[3];
-<a class="jxr_linenumber" name="142" href="#142">142</a> <strong class="jxr_keyword">double</strong> c = x[1] - 2 * x[2];
-<a class="jxr_linenumber" name="143" href="#143">143</a> <strong class="jxr_keyword">double</strong> d = x[0] - x[3];
-<a class="jxr_linenumber" name="144" href="#144">144</a> <strong class="jxr_keyword">return</strong> a * a + 5 * b * b + c * c * c * c + 10 * d * d * d * d;
-<a class="jxr_linenumber" name="145" href="#145">145</a> }
-<a class="jxr_linenumber" name="146" href="#146">146</a> };
-<a class="jxr_linenumber" name="147" href="#147">147</a>
-<a class="jxr_linenumber" name="148" href="#148">148</a> count = 0;
-<a class="jxr_linenumber" name="149" href="#149">149</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(-1, 1e-3);
-<a class="jxr_linenumber" name="150" href="#150">150</a> PointValuePair optimum
-<a class="jxr_linenumber" name="151" href="#151">151</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(1000),
-<a class="jxr_linenumber" name="152" href="#152">152</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(powell),
-<a class="jxr_linenumber" name="153" href="#153">153</a> GoalType.MINIMIZE,
-<a class="jxr_linenumber" name="154" href="#154">154</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 3, -1, 0, 1 }),
-<a class="jxr_linenumber" name="155" href="#155">155</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(4));
+<a class="jxr_linenumber" name="128" href="#128">128</a> <em class="jxr_comment">// Check that the number of iterations is updated (MATH-949).</em>
+<a class="jxr_linenumber" name="129" href="#129">129</a> Assert.assertTrue(optimizer.getIterations() > 0);
+<a class="jxr_linenumber" name="130" href="#130">130</a> }
+<a class="jxr_linenumber" name="131" href="#131">131</a>
+<a class="jxr_linenumber" name="132" href="#132">132</a> @Test
+<a class="jxr_linenumber" name="133" href="#133">133</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testRosenbrock() {
+<a class="jxr_linenumber" name="134" href="#134">134</a> MultivariateFunction rosenbrock
+<a class="jxr_linenumber" name="135" href="#135">135</a> = <strong class="jxr_keyword">new</strong> MultivariateFunction() {
+<a class="jxr_linenumber" name="136" href="#136">136</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] x) {
+<a class="jxr_linenumber" name="137" href="#137">137</a> ++count;
+<a class="jxr_linenumber" name="138" href="#138">138</a> <strong class="jxr_keyword">double</strong> a = x[1] - x[0] * x[0];
+<a class="jxr_linenumber" name="139" href="#139">139</a> <strong class="jxr_keyword">double</strong> b = 1.0 - x[0];
+<a class="jxr_linenumber" name="140" href="#140">140</a> <strong class="jxr_keyword">return</strong> 100 * a * a + b * b;
+<a class="jxr_linenumber" name="141" href="#141">141</a> }
+<a class="jxr_linenumber" name="142" href="#142">142</a> };
+<a class="jxr_linenumber" name="143" href="#143">143</a>
+<a class="jxr_linenumber" name="144" href="#144">144</a> count = 0;
+<a class="jxr_linenumber" name="145" href="#145">145</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(-1, 1e-3);
+<a class="jxr_linenumber" name="146" href="#146">146</a> PointValuePair optimum
+<a class="jxr_linenumber" name="147" href="#147">147</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(100),
+<a class="jxr_linenumber" name="148" href="#148">148</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(rosenbrock),
+<a class="jxr_linenumber" name="149" href="#149">149</a> GoalType.MINIMIZE,
+<a class="jxr_linenumber" name="150" href="#150">150</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { -1.2, 1 }),
+<a class="jxr_linenumber" name="151" href="#151">151</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[][] {
+<a class="jxr_linenumber" name="152" href="#152">152</a> { -1.2, 1.0 },
+<a class="jxr_linenumber" name="153" href="#153">153</a> { 0.9, 1.2 },
+<a class="jxr_linenumber" name="154" href="#154">154</a> { 3.5, -2.3 } }));
+<a class="jxr_linenumber" name="155" href="#155">155</a>
<a class="jxr_linenumber" name="156" href="#156">156</a> Assert.assertEquals(count, optimizer.getEvaluations());
-<a class="jxr_linenumber" name="157" href="#157">157</a> Assert.assertTrue(optimizer.getEvaluations() > 800);
-<a class="jxr_linenumber" name="158" href="#158">158</a> Assert.assertTrue(optimizer.getEvaluations() < 900);
+<a class="jxr_linenumber" name="157" href="#157">157</a> Assert.assertTrue(optimizer.getEvaluations() > 50);
+<a class="jxr_linenumber" name="158" href="#158">158</a> Assert.assertTrue(optimizer.getEvaluations() < 100);
<a class="jxr_linenumber" name="159" href="#159">159</a> Assert.assertTrue(optimum.getValue() > 1e-2);
<a class="jxr_linenumber" name="160" href="#160">160</a> }
<a class="jxr_linenumber" name="161" href="#161">161</a>
<a class="jxr_linenumber" name="162" href="#162">162</a> @Test
-<a class="jxr_linenumber" name="163" href="#163">163</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMath283() {
-<a class="jxr_linenumber" name="164" href="#164">164</a> <em class="jxr_comment">// fails because MultiDirectional.iterateSimplex is looping forever</em>
-<a class="jxr_linenumber" name="165" href="#165">165</a> <em class="jxr_comment">// the while(true) should be replaced with a convergence check</em>
-<a class="jxr_linenumber" name="166" href="#166">166</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-14, 1e-14);
-<a class="jxr_linenumber" name="167" href="#167">167</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a> function = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a>(0, 0, 1);
-<a class="jxr_linenumber" name="168" href="#168">168</a> PointValuePair estimate = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(1000),
-<a class="jxr_linenumber" name="169" href="#169">169</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(function),
-<a class="jxr_linenumber" name="170" href="#170">170</a> GoalType.MAXIMIZE,
-<a class="jxr_linenumber" name="171" href="#171">171</a> <strong class="jxr_keyword">new</strong> InitialGuess(function.getMaximumPosition()),
-<a class="jxr_linenumber" name="172" href="#172">172</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(2));
-<a class="jxr_linenumber" name="173" href="#173">173</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> EPSILON = 1e-5;
-<a class="jxr_linenumber" name="174" href="#174">174</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expectedMaximum = function.getMaximum();
-<a class="jxr_linenumber" name="175" href="#175">175</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actualMaximum = estimate.getValue();
-<a class="jxr_linenumber" name="176" href="#176">176</a> Assert.assertEquals(expectedMaximum, actualMaximum, EPSILON);
-<a class="jxr_linenumber" name="177" href="#177">177</a>
-<a class="jxr_linenumber" name="178" href="#178">178</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] expectedPosition = function.getMaximumPosition();
-<a class="jxr_linenumber" name="179" href="#179">179</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] actualPosition = estimate.getPoint();
-<a class="jxr_linenumber" name="180" href="#180">180</a> Assert.assertEquals(expectedPosition[0], actualPosition[0], EPSILON );
-<a class="jxr_linenumber" name="181" href="#181">181</a> Assert.assertEquals(expectedPosition[1], actualPosition[1], EPSILON );
-<a class="jxr_linenumber" name="182" href="#182">182</a> }
-<a class="jxr_linenumber" name="183" href="#183">183</a>
-<a class="jxr_linenumber" name="184" href="#184">184</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerNelderMeadTest.html">FourExtrema</a> <strong class="jxr_keyword">implements</strong> MultivariateFunction {
-<a class="jxr_linenumber" name="185" href="#185">185</a> <em class="jxr_comment">// The following function has 4 local extrema.</em>
-<a class="jxr_linenumber" name="186" href="#186">186</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xM = -3.841947088256863675365;
-<a class="jxr_linenumber" name="187" href="#187">187</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> yM = -1.391745200270734924416;
-<a class="jxr_linenumber" name="188" href="#188">188</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xP = 0.2286682237349059125691;
-<a class="jxr_linenumber" name="189" href="#189">189</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> yP = -yM;
-<a class="jxr_linenumber" name="190" href="#190">190</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXmYm = 0.2373295333134216789769; <em class="jxr_comment">// Local maximum.</em>
-<a class="jxr_linenumber" name="191" href="#191">191</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXmYp = -valueXmYm; <em class="jxr_comment">// Local minimum.</em>
-<a class="jxr_linenumber" name="192" href="#192">192</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXpYm = -0.7290400707055187115322; <em class="jxr_comment">// Global minimum.</em>
-<a class="jxr_linenumber" name="193" href="#193">193</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXpYp = -valueXpYm; <em class="jxr_comment">// Global maximum.</em>
-<a class="jxr_linenumber" name="194" href="#194">194</a>
-<a class="jxr_linenumber" name="195" href="#195">195</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] variables) {
-<a class="jxr_linenumber" name="196" href="#196">196</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = variables[0];
-<a class="jxr_linenumber" name="197" href="#197">197</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> y = variables[1];
-<a class="jxr_linenumber" name="198" href="#198">198</a> <strong class="jxr_keyword">return</strong> (x == 0 || y == 0) ? 0 :
-<a class="jxr_linenumber" name="199" href="#199">199</a> FastMath.atan(x) * FastMath.atan(x + 2) * FastMath.atan(y) * FastMath.atan(y) / (x * y);
-<a class="jxr_linenumber" name="200" href="#200">200</a> }
-<a class="jxr_linenumber" name="201" href="#201">201</a> }
-<a class="jxr_linenumber" name="202" href="#202">202</a>
-<a class="jxr_linenumber" name="203" href="#203">203</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a> <strong class="jxr_keyword">implements</strong> MultivariateFunction {
-<a class="jxr_linenumber" name="204" href="#204">204</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] maximumPosition;
-<a class="jxr_linenumber" name="205" href="#205">205</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> std;
-<a class="jxr_linenumber" name="206" href="#206">206</a>
-<a class="jxr_linenumber" name="207" href="#207">207</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a>(<strong class="jxr_keyword">double</strong> xOpt, <strong class="jxr_keyword">double</strong> yOpt, <strong class="jxr_keyword">double</strong> std) {
-<a class="jxr_linenumber" name="208" href="#208">208</a> maximumPosition = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { xOpt, yOpt };
-<a class="jxr_linenumber" name="209" href="#209">209</a> <strong class="jxr_keyword">this</strong>.std = std;
-<a class="jxr_linenumber" name="210" href="#210">210</a> }
+<a class="jxr_linenumber" name="163" href="#163">163</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testPowell() {
+<a class="jxr_linenumber" name="164" href="#164">164</a> MultivariateFunction powell
+<a class="jxr_linenumber" name="165" href="#165">165</a> = <strong class="jxr_keyword">new</strong> MultivariateFunction() {
+<a class="jxr_linenumber" name="166" href="#166">166</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] x) {
+<a class="jxr_linenumber" name="167" href="#167">167</a> ++count;
+<a class="jxr_linenumber" name="168" href="#168">168</a> <strong class="jxr_keyword">double</strong> a = x[0] + 10 * x[1];
+<a class="jxr_linenumber" name="169" href="#169">169</a> <strong class="jxr_keyword">double</strong> b = x[2] - x[3];
+<a class="jxr_linenumber" name="170" href="#170">170</a> <strong class="jxr_keyword">double</strong> c = x[1] - 2 * x[2];
+<a class="jxr_linenumber" name="171" href="#171">171</a> <strong class="jxr_keyword">double</strong> d = x[0] - x[3];
+<a class="jxr_linenumber" name="172" href="#172">172</a> <strong class="jxr_keyword">return</strong> a * a + 5 * b * b + c * c * c * c + 10 * d * d * d * d;
+<a class="jxr_linenumber" name="173" href="#173">173</a> }
+<a class="jxr_linenumber" name="174" href="#174">174</a> };
+<a class="jxr_linenumber" name="175" href="#175">175</a>
+<a class="jxr_linenumber" name="176" href="#176">176</a> count = 0;
+<a class="jxr_linenumber" name="177" href="#177">177</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(-1, 1e-3);
+<a class="jxr_linenumber" name="178" href="#178">178</a> PointValuePair optimum
+<a class="jxr_linenumber" name="179" href="#179">179</a> = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(1000),
+<a class="jxr_linenumber" name="180" href="#180">180</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(powell),
+<a class="jxr_linenumber" name="181" href="#181">181</a> GoalType.MINIMIZE,
+<a class="jxr_linenumber" name="182" href="#182">182</a> <strong class="jxr_keyword">new</strong> InitialGuess(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 3, -1, 0, 1 }),
+<a class="jxr_linenumber" name="183" href="#183">183</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(4));
+<a class="jxr_linenumber" name="184" href="#184">184</a> Assert.assertEquals(count, optimizer.getEvaluations());
+<a class="jxr_linenumber" name="185" href="#185">185</a> Assert.assertTrue(optimizer.getEvaluations() > 800);
+<a class="jxr_linenumber" name="186" href="#186">186</a> Assert.assertTrue(optimizer.getEvaluations() < 900);
+<a class="jxr_linenumber" name="187" href="#187">187</a> Assert.assertTrue(optimum.getValue() > 1e-2);
+<a class="jxr_linenumber" name="188" href="#188">188</a> }
+<a class="jxr_linenumber" name="189" href="#189">189</a>
+<a class="jxr_linenumber" name="190" href="#190">190</a> @Test
+<a class="jxr_linenumber" name="191" href="#191">191</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testMath283() {
+<a class="jxr_linenumber" name="192" href="#192">192</a> <em class="jxr_comment">// fails because MultiDirectional.iterateSimplex is looping forever</em>
+<a class="jxr_linenumber" name="193" href="#193">193</a> <em class="jxr_comment">// the while(true) should be replaced with a convergence check</em>
+<a class="jxr_linenumber" name="194" href="#194">194</a> SimplexOptimizer optimizer = <strong class="jxr_keyword">new</strong> SimplexOptimizer(1e-14, 1e-14);
+<a class="jxr_linenumber" name="195" href="#195">195</a> <strong class="jxr_keyword">final</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a> function = <strong class="jxr_keyword">new</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a>(0, 0, 1);
+<a class="jxr_linenumber" name="196" href="#196">196</a> PointValuePair estimate = optimizer.optimize(<strong class="jxr_keyword">new</strong> MaxEval(1000),
+<a class="jxr_linenumber" name="197" href="#197">197</a> <strong class="jxr_keyword">new</strong> ObjectiveFunction(function),
+<a class="jxr_linenumber" name="198" href="#198">198</a> GoalType.MAXIMIZE,
+<a class="jxr_linenumber" name="199" href="#199">199</a> <strong class="jxr_keyword">new</strong> InitialGuess(function.getMaximumPosition()),
+<a class="jxr_linenumber" name="200" href="#200">200</a> <strong class="jxr_keyword">new</strong> MultiDirectionalSimplex(2));
+<a class="jxr_linenumber" name="201" href="#201">201</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> EPSILON = 1e-5;
+<a class="jxr_linenumber" name="202" href="#202">202</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expectedMaximum = function.getMaximum();
+<a class="jxr_linenumber" name="203" href="#203">203</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actualMaximum = estimate.getValue();
+<a class="jxr_linenumber" name="204" href="#204">204</a> Assert.assertEquals(expectedMaximum, actualMaximum, EPSILON);
+<a class="jxr_linenumber" name="205" href="#205">205</a>
+<a class="jxr_linenumber" name="206" href="#206">206</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] expectedPosition = function.getMaximumPosition();
+<a class="jxr_linenumber" name="207" href="#207">207</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] actualPosition = estimate.getPoint();
+<a class="jxr_linenumber" name="208" href="#208">208</a> Assert.assertEquals(expectedPosition[0], actualPosition[0], EPSILON );
+<a class="jxr_linenumber" name="209" href="#209">209</a> Assert.assertEquals(expectedPosition[1], actualPosition[1], EPSILON );
+<a class="jxr_linenumber" name="210" href="#210">210</a> }
<a class="jxr_linenumber" name="211" href="#211">211</a>
-<a class="jxr_linenumber" name="212" href="#212">212</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getMaximum() {
-<a class="jxr_linenumber" name="213" href="#213">213</a> <strong class="jxr_keyword">return</strong> value(maximumPosition);
-<a class="jxr_linenumber" name="214" href="#214">214</a> }
-<a class="jxr_linenumber" name="215" href="#215">215</a>
-<a class="jxr_linenumber" name="216" href="#216">216</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong>[] getMaximumPosition() {
-<a class="jxr_linenumber" name="217" href="#217">217</a> <strong class="jxr_keyword">return</strong> maximumPosition.clone();
-<a class="jxr_linenumber" name="218" href="#218">218</a> }
-<a class="jxr_linenumber" name="219" href="#219">219</a>
-<a class="jxr_linenumber" name="220" href="#220">220</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] point) {
-<a class="jxr_linenumber" name="221" href="#221">221</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = point[0], y = point[1];
-<a class="jxr_linenumber" name="222" href="#222">222</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> twoS2 = 2.0 * std * std;
-<a class="jxr_linenumber" name="223" href="#223">223</a> <strong class="jxr_keyword">return</strong> 1.0 / (twoS2 * FastMath.PI) * FastMath.exp(-(x * x + y * y) / twoS2);
-<a class="jxr_linenumber" name="224" href="#224">224</a> }
-<a class="jxr_linenumber" name="225" href="#225">225</a> }
-<a class="jxr_linenumber" name="226" href="#226">226</a>
-<a class="jxr_linenumber" name="227" href="#227">227</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> count;
-<a class="jxr_linenumber" name="228" href="#228">228</a> }
+<a class="jxr_linenumber" name="212" href="#212">212</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">FourExtrema</a> <strong class="jxr_keyword">implements</strong> MultivariateFunction {
+<a class="jxr_linenumber" name="213" href="#213">213</a> <em class="jxr_comment">// The following function has 4 local extrema.</em>
+<a class="jxr_linenumber" name="214" href="#214">214</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xM = -3.841947088256863675365;
+<a class="jxr_linenumber" name="215" href="#215">215</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> yM = -1.391745200270734924416;
+<a class="jxr_linenumber" name="216" href="#216">216</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xP = 0.2286682237349059125691;
+<a class="jxr_linenumber" name="217" href="#217">217</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> yP = -yM;
+<a class="jxr_linenumber" name="218" href="#218">218</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXmYm = 0.2373295333134216789769; <em class="jxr_comment">// Local maximum.</em>
+<a class="jxr_linenumber" name="219" href="#219">219</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXmYp = -valueXmYm; <em class="jxr_comment">// Local minimum.</em>
+<a class="jxr_linenumber" name="220" href="#220">220</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXpYm = -0.7290400707055187115322; <em class="jxr_comment">// Global minimum.</em>
+<a class="jxr_linenumber" name="221" href="#221">221</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> valueXpYp = -valueXpYm; <em class="jxr_comment">// Global maximum.</em>
+<a class="jxr_linenumber" name="222" href="#222">222</a>
+<a class="jxr_linenumber" name="223" href="#223">223</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] variables) {
+<a class="jxr_linenumber" name="224" href="#224">224</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = variables[0];
+<a class="jxr_linenumber" name="225" href="#225">225</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> y = variables[1];
+<a class="jxr_linenumber" name="226" href="#226">226</a> <strong class="jxr_keyword">return</strong> (x == 0 || y == 0) ? 0 :
+<a class="jxr_linenumber" name="227" href="#227">227</a> FastMath.atan(x) * FastMath.atan(x + 2) * FastMath.atan(y) * FastMath.atan(y) / (x * y);
+<a class="jxr_linenumber" name="228" href="#228">228</a> }
+<a class="jxr_linenumber" name="229" href="#229">229</a> }
+<a class="jxr_linenumber" name="230" href="#230">230</a>
+<a class="jxr_linenumber" name="231" href="#231">231</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a> <strong class="jxr_keyword">implements</strong> MultivariateFunction {
+<a class="jxr_linenumber" name="232" href="#232">232</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] maximumPosition;
+<a class="jxr_linenumber" name="233" href="#233">233</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> std;
+<a class="jxr_linenumber" name="234" href="#234">234</a>
+<a class="jxr_linenumber" name="235" href="#235">235</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../../../org/apache/commons/math3/optim/nonlinear/scalar/noderiv/SimplexOptimizerMultiDirectionalTest.html">Gaussian2D</a>(<strong class="jxr_keyword">double</strong> xOpt, <strong class="jxr_keyword">double</strong> yOpt, <strong class="jxr_keyword">double</strong> std) {
+<a class="jxr_linenumber" name="236" href="#236">236</a> maximumPosition = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { xOpt, yOpt };
+<a class="jxr_linenumber" name="237" href="#237">237</a> <strong class="jxr_keyword">this</strong>.std = std;
+<a class="jxr_linenumber" name="238" href="#238">238</a> }
+<a class="jxr_linenumber" name="239" href="#239">239</a>
+<a class="jxr_linenumber" name="240" href="#240">240</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getMaximum() {
+<a class="jxr_linenumber" name="241" href="#241">241</a> <strong class="jxr_keyword">return</strong> value(maximumPosition);
+<a class="jxr_linenumber" name="242" href="#242">242</a> }
+<a class="jxr_linenumber" name="243" href="#243">243</a>
+<a class="jxr_linenumber" name="244" href="#244">244</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong>[] getMaximumPosition() {
+<a class="jxr_linenumber" name="245" href="#245">245</a> <strong class="jxr_keyword">return</strong> maximumPosition.clone();
+<a class="jxr_linenumber" name="246" href="#246">246</a> }
+<a class="jxr_linenumber" name="247" href="#247">247</a>
+<a class="jxr_linenumber" name="248" href="#248">248</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> value(<strong class="jxr_keyword">double</strong>[] point) {
+<a class="jxr_linenumber" name="249" href="#249">249</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = point[0], y = point[1];
+<a class="jxr_linenumber" name="250" href="#250">250</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> twoS2 = 2.0 * std * std;
+<a class="jxr_linenumber" name="251" href="#251">251</a> <strong class="jxr_keyword">return</strong> 1.0 / (twoS2 * FastMath.PI) * FastMath.exp(-(x * x + y * y) / twoS2);
+<a class="jxr_linenumber" name="252" href="#252">252</a> }
+<a class="jxr_linenumber" name="253" href="#253">253</a> }
+<a class="jxr_linenumber" name="254" href="#254">254</a>
+<a class="jxr_linenumber" name="255" href="#255">255</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> count;
+<a class="jxr_linenumber" name="256" href="#256">256</a> }
</pre>
<hr/><div id="footer">This page was automatically generated by <a href="http://maven.apache.org/">Maven</a></div></body>
</html>