commons-notifications mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From l..@apache.org
Subject svn commit: r948061 [4/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/LaguerreSolver.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/LaguerreSolver.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/LaguerreSolver.html Fri Apr 17 19:17:13 2015
@@ -39,10 +39,10 @@
 <a class="jxr_linenumber" name="L31" href="#L31">31</a>  <em class="jxr_javadoccomment"> * Implements the &lt;a href="<a href="http://mathworld.wolfram.com/LaguerresMethod.html" target="alexandria_uri">http://mathworld.wolfram.com/LaguerresMethod.html</a>"&gt;</em>
 <a class="jxr_linenumber" name="L32" href="#L32">32</a>  <em class="jxr_javadoccomment"> * Laguerre's Method&lt;/a&gt; for root finding of real coefficient polynomials.</em>
 <a class="jxr_linenumber" name="L33" href="#L33">33</a>  <em class="jxr_javadoccomment"> * For reference, see</em>
-<a class="jxr_linenumber" name="L34" href="#L34">34</a>  <em class="jxr_javadoccomment"> * &lt;quote&gt;</em>
-<a class="jxr_linenumber" name="L35" href="#L35">35</a>  <em class="jxr_javadoccomment"> *  &lt;b&gt;A First Course in Numerical Analysis&lt;/b&gt;</em>
-<a class="jxr_linenumber" name="L36" href="#L36">36</a>  <em class="jxr_javadoccomment"> *  ISBN 048641454X, chapter 8.</em>
-<a class="jxr_linenumber" name="L37" href="#L37">37</a>  <em class="jxr_javadoccomment"> * &lt;/quote&gt;</em>
+<a class="jxr_linenumber" name="L34" href="#L34">34</a>  <em class="jxr_javadoccomment"> * &lt;blockquote&gt;</em>
+<a class="jxr_linenumber" name="L35" href="#L35">35</a>  <em class="jxr_javadoccomment"> *  &lt;b&gt;A First Course in Numerical Analysis&lt;/b&gt;&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L36" href="#L36">36</a>  <em class="jxr_javadoccomment"> *  ISBN 048641454X, chapter 8.&lt;br&gt;</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"> * Laguerre's method is global in the sense that it can start with any initial</em>
 <a class="jxr_linenumber" name="L39" href="#L39">39</a>  <em class="jxr_javadoccomment"> * approximation and be able to solve all roots from that point.</em>
 <a class="jxr_linenumber" name="L40" href="#L40">40</a>  <em class="jxr_javadoccomment"> * The algorithm requires a bracketing condition.</em>
@@ -189,9 +189,9 @@
 <a class="jxr_linenumber" name="L181" href="#L181">181</a> <em class="jxr_javadoccomment">     *</em>
 <a class="jxr_linenumber" name="L182" href="#L182">182</a> <em class="jxr_javadoccomment">     * @param coefficients Polynomial coefficients.</em>
 <a class="jxr_linenumber" name="L183" href="#L183">183</a> <em class="jxr_javadoccomment">     * @param initial Start value.</em>
-<a class="jxr_linenumber" name="L184" href="#L184">184</a> <em class="jxr_javadoccomment">     * @return the point at which the function value is zero.</em>
+<a class="jxr_linenumber" name="L184" href="#L184">184</a> <em class="jxr_javadoccomment">     * @return the full set of complex roots of the polynomial</em>
 <a class="jxr_linenumber" name="L185" href="#L185">185</a> <em class="jxr_javadoccomment">     * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
-<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded.</em>
+<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded when solving for one of the roots</em>
 <a class="jxr_linenumber" name="L187" href="#L187">187</a> <em class="jxr_javadoccomment">     * @throws NullArgumentException if the {@code coefficients} is</em>
 <a class="jxr_linenumber" name="L188" href="#L188">188</a> <em class="jxr_javadoccomment">     * {@code null}.</em>
 <a class="jxr_linenumber" name="L189" href="#L189">189</a> <em class="jxr_javadoccomment">     * @throws NoDataException if the {@code coefficients} array is empty.</em>
@@ -202,200 +202,250 @@
 <a class="jxr_linenumber" name="L194" href="#L194">194</a>         <strong class="jxr_keyword">throws</strong> NullArgumentException,
 <a class="jxr_linenumber" name="L195" href="#L195">195</a>                NoDataException,
 <a class="jxr_linenumber" name="L196" href="#L196">196</a>                <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
-<a class="jxr_linenumber" name="L197" href="#L197">197</a>         setup(Integer.MAX_VALUE,
-<a class="jxr_linenumber" name="L198" href="#L198">198</a>               <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/analysis/polynomials/PolynomialFunction.html">PolynomialFunction</a>(coefficients),
-<a class="jxr_linenumber" name="L199" href="#L199">199</a>               Double.NEGATIVE_INFINITY,
-<a class="jxr_linenumber" name="L200" href="#L200">200</a>               Double.POSITIVE_INFINITY,
-<a class="jxr_linenumber" name="L201" href="#L201">201</a>               initial);
-<a class="jxr_linenumber" name="L202" href="#L202">202</a>         <strong class="jxr_keyword">return</strong> complexSolver.solveAll(ComplexUtils.convertToComplex(coefficients),
-<a class="jxr_linenumber" name="L203" href="#L203">203</a>                                       <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(initial, 0d));
-<a class="jxr_linenumber" name="L204" href="#L204">204</a>     }
-<a class="jxr_linenumber" name="L205" href="#L205">205</a> 
-<a class="jxr_linenumber" name="L206" href="#L206">206</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L207" href="#L207">207</a> <em class="jxr_javadoccomment">     * Find a complex root for the polynomial with the given coefficients,</em>
-<a class="jxr_linenumber" name="L208" href="#L208">208</a> <em class="jxr_javadoccomment">     * starting from the given initial value.</em>
-<a class="jxr_linenumber" name="L209" href="#L209">209</a> <em class="jxr_javadoccomment">     * &lt;br/&gt;</em>
-<a class="jxr_linenumber" name="L210" href="#L210">210</a> <em class="jxr_javadoccomment">     * Note: This method is not part of the API of {@link BaseUnivariateSolver}.</em>
-<a class="jxr_linenumber" name="L211" href="#L211">211</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L212" href="#L212">212</a> <em class="jxr_javadoccomment">     * @param coefficients Polynomial coefficients.</em>
-<a class="jxr_linenumber" name="L213" href="#L213">213</a> <em class="jxr_javadoccomment">     * @param initial Start value.</em>
-<a class="jxr_linenumber" name="L214" href="#L214">214</a> <em class="jxr_javadoccomment">     * @return the point at which the function value is zero.</em>
-<a class="jxr_linenumber" name="L215" href="#L215">215</a> <em class="jxr_javadoccomment">     * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
-<a class="jxr_linenumber" name="L216" href="#L216">216</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded.</em>
-<a class="jxr_linenumber" name="L217" href="#L217">217</a> <em class="jxr_javadoccomment">     * @throws NullArgumentException if the {@code coefficients} is</em>
-<a class="jxr_linenumber" name="L218" href="#L218">218</a> <em class="jxr_javadoccomment">     * {@code null}.</em>
-<a class="jxr_linenumber" name="L219" href="#L219">219</a> <em class="jxr_javadoccomment">     * @throws NoDataException if the {@code coefficients} array is empty.</em>
-<a class="jxr_linenumber" name="L220" href="#L220">220</a> <em class="jxr_javadoccomment">     * @since 3.1</em>
-<a class="jxr_linenumber" name="L221" href="#L221">221</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L222" href="#L222">222</a>     <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> solveComplex(<strong class="jxr_keyword">double</strong>[] coefficients,
-<a class="jxr_linenumber" name="L223" href="#L223">223</a>                                 <strong class="jxr_keyword">double</strong> initial)
-<a class="jxr_linenumber" name="L224" href="#L224">224</a>         <strong class="jxr_keyword">throws</strong> NullArgumentException,
-<a class="jxr_linenumber" name="L225" href="#L225">225</a>                NoDataException,
-<a class="jxr_linenumber" name="L226" href="#L226">226</a>                <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
-<a class="jxr_linenumber" name="L227" href="#L227">227</a>         setup(Integer.MAX_VALUE,
-<a class="jxr_linenumber" name="L228" href="#L228">228</a>               <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/analysis/polynomials/PolynomialFunction.html">PolynomialFunction</a>(coefficients),
-<a class="jxr_linenumber" name="L229" href="#L229">229</a>               Double.NEGATIVE_INFINITY,
-<a class="jxr_linenumber" name="L230" href="#L230">230</a>               Double.POSITIVE_INFINITY,
-<a class="jxr_linenumber" name="L231" href="#L231">231</a>               initial);
-<a class="jxr_linenumber" name="L232" href="#L232">232</a>         <strong class="jxr_keyword">return</strong> complexSolver.solve(ComplexUtils.convertToComplex(coefficients),
-<a class="jxr_linenumber" name="L233" href="#L233">233</a>                                    <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(initial, 0d));
-<a class="jxr_linenumber" name="L234" href="#L234">234</a>     }
-<a class="jxr_linenumber" name="L235" href="#L235">235</a> 
-<a class="jxr_linenumber" name="L236" href="#L236">236</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L237" href="#L237">237</a> <em class="jxr_javadoccomment">     * Class for searching all (complex) roots.</em>
-<a class="jxr_linenumber" name="L238" href="#L238">238</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L239" href="#L239">239</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/LaguerreSolver.html">ComplexSolver</a> {
-<a class="jxr_linenumber" name="L240" href="#L240">240</a>         <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L241" href="#L241">241</a> <em class="jxr_javadoccomment">         * Check whether the given complex root is actually a real zero</em>
-<a class="jxr_linenumber" name="L242" href="#L242">242</a> <em class="jxr_javadoccomment">         * in the given interval, within the solver tolerance level.</em>
-<a class="jxr_linenumber" name="L243" href="#L243">243</a> <em class="jxr_javadoccomment">         *</em>
-<a class="jxr_linenumber" name="L244" href="#L244">244</a> <em class="jxr_javadoccomment">         * @param min Lower bound for the interval.</em>
-<a class="jxr_linenumber" name="L245" href="#L245">245</a> <em class="jxr_javadoccomment">         * @param max Upper bound for the interval.</em>
-<a class="jxr_linenumber" name="L246" href="#L246">246</a> <em class="jxr_javadoccomment">         * @param z Complex root.</em>
-<a class="jxr_linenumber" name="L247" href="#L247">247</a> <em class="jxr_javadoccomment">         * @return {@code true} if z is a real zero.</em>
-<a class="jxr_linenumber" name="L248" href="#L248">248</a> <em class="jxr_javadoccomment">         */</em>
-<a class="jxr_linenumber" name="L249" href="#L249">249</a>         <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">boolean</strong> isRoot(<strong class="jxr_keyword">double</strong> min, <strong class="jxr_keyword">double</strong> max, <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> z) {
-<a class="jxr_linenumber" name="L250" href="#L250">250</a>             <strong class="jxr_keyword">if</strong> (isSequence(min, z.getReal(), max)) {
-<a class="jxr_linenumber" name="L251" href="#L251">251</a>                 <strong class="jxr_keyword">double</strong> tolerance = FastMath.max(getRelativeAccuracy() * z.abs(), getAbsoluteAccuracy());
-<a class="jxr_linenumber" name="L252" href="#L252">252</a>                 <strong class="jxr_keyword">return</strong> (FastMath.abs(z.getImaginary()) &lt;= tolerance) ||
-<a class="jxr_linenumber" name="L253" href="#L253">253</a>                      (z.abs() &lt;= getFunctionValueAccuracy());
-<a class="jxr_linenumber" name="L254" href="#L254">254</a>             }
-<a class="jxr_linenumber" name="L255" href="#L255">255</a>             <strong class="jxr_keyword">return</strong> false;
-<a class="jxr_linenumber" name="L256" href="#L256">256</a>         }
-<a class="jxr_linenumber" name="L257" href="#L257">257</a> 
-<a class="jxr_linenumber" name="L258" href="#L258">258</a>         <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L259" href="#L259">259</a> <em class="jxr_javadoccomment">         * Find all complex roots for the polynomial with the given</em>
-<a class="jxr_linenumber" name="L260" href="#L260">260</a> <em class="jxr_javadoccomment">         * coefficients, starting from the given initial value.</em>
-<a class="jxr_linenumber" name="L261" href="#L261">261</a> <em class="jxr_javadoccomment">         *</em>
-<a class="jxr_linenumber" name="L262" href="#L262">262</a> <em class="jxr_javadoccomment">         * @param coefficients Polynomial coefficients.</em>
-<a class="jxr_linenumber" name="L263" href="#L263">263</a> <em class="jxr_javadoccomment">         * @param initial Start value.</em>
-<a class="jxr_linenumber" name="L264" href="#L264">264</a> <em class="jxr_javadoccomment">         * @return the point at which the function value is zero.</em>
-<a class="jxr_linenumber" name="L265" href="#L265">265</a> <em class="jxr_javadoccomment">         * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
-<a class="jxr_linenumber" name="L266" href="#L266">266</a> <em class="jxr_javadoccomment">         * if the maximum number of evaluations is exceeded.</em>
-<a class="jxr_linenumber" name="L267" href="#L267">267</a> <em class="jxr_javadoccomment">         * @throws NullArgumentException if the {@code coefficients} is</em>
-<a class="jxr_linenumber" name="L268" href="#L268">268</a> <em class="jxr_javadoccomment">         * {@code null}.</em>
-<a class="jxr_linenumber" name="L269" href="#L269">269</a> <em class="jxr_javadoccomment">         * @throws NoDataException if the {@code coefficients} array is empty.</em>
-<a class="jxr_linenumber" name="L270" href="#L270">270</a> <em class="jxr_javadoccomment">         */</em>
-<a class="jxr_linenumber" name="L271" href="#L271">271</a>         <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[] solveAll(<a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> coefficients[], <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> initial)
-<a class="jxr_linenumber" name="L272" href="#L272">272</a>             <strong class="jxr_keyword">throws</strong> NullArgumentException,
-<a class="jxr_linenumber" name="L273" href="#L273">273</a>                    NoDataException,
-<a class="jxr_linenumber" name="L274" href="#L274">274</a>                    <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
-<a class="jxr_linenumber" name="L275" href="#L275">275</a>             <strong class="jxr_keyword">if</strong> (coefficients == <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L276" href="#L276">276</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NullArgumentException.html">NullArgumentException</a>();
-<a class="jxr_linenumber" name="L277" href="#L277">277</a>             }
-<a class="jxr_linenumber" name="L278" href="#L278">278</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = coefficients.length - 1;
-<a class="jxr_linenumber" name="L279" href="#L279">279</a>             <strong class="jxr_keyword">if</strong> (n == 0) {
-<a class="jxr_linenumber" name="L280" href="#L280">280</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NoDataException.html">NoDataException</a>(LocalizedFormats.POLYNOMIAL);
-<a class="jxr_linenumber" name="L281" href="#L281">281</a>             }
-<a class="jxr_linenumber" name="L282" href="#L282">282</a>             <em class="jxr_comment">// Coefficients for deflated polynomial.</em>
-<a class="jxr_linenumber" name="L283" href="#L283">283</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> c[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n + 1];
-<a class="jxr_linenumber" name="L284" href="#L284">284</a>             <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt;= n; i++) {
-<a class="jxr_linenumber" name="L285" href="#L285">285</a>                 c[i] = coefficients[i];
-<a class="jxr_linenumber" name="L286" href="#L286">286</a>             }
-<a class="jxr_linenumber" name="L287" href="#L287">287</a> 
-<a class="jxr_linenumber" name="L288" href="#L288">288</a>             <em class="jxr_comment">// Solve individual roots successively.</em>
-<a class="jxr_linenumber" name="L289" href="#L289">289</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> root[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n];
-<a class="jxr_linenumber" name="L290" href="#L290">290</a>             <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; n; i++) {
-<a class="jxr_linenumber" name="L291" href="#L291">291</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> subarray[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n - i + 1];
-<a class="jxr_linenumber" name="L292" href="#L292">292</a>                 System.arraycopy(c, 0, subarray, 0, subarray.length);
-<a class="jxr_linenumber" name="L293" href="#L293">293</a>                 root[i] = solve(subarray, initial);
-<a class="jxr_linenumber" name="L294" href="#L294">294</a>                 <em class="jxr_comment">// Polynomial deflation using synthetic division.</em>
-<a class="jxr_linenumber" name="L295" href="#L295">295</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> newc = c[n - i];
-<a class="jxr_linenumber" name="L296" href="#L296">296</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> oldc = <strong class="jxr_keyword">null</strong>;
-<a class="jxr_linenumber" name="L297" href="#L297">297</a>                 <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = n - i - 1; j &gt;= 0; j--) {
-<a class="jxr_linenumber" name="L298" href="#L298">298</a>                     oldc = c[j];
-<a class="jxr_linenumber" name="L299" href="#L299">299</a>                     c[j] = newc;
-<a class="jxr_linenumber" name="L300" href="#L300">300</a>                     newc = oldc.add(newc.multiply(root[i]));
-<a class="jxr_linenumber" name="L301" href="#L301">301</a>                 }
-<a class="jxr_linenumber" name="L302" href="#L302">302</a>             }
-<a class="jxr_linenumber" name="L303" href="#L303">303</a> 
-<a class="jxr_linenumber" name="L304" href="#L304">304</a>             <strong class="jxr_keyword">return</strong> root;
-<a class="jxr_linenumber" name="L305" href="#L305">305</a>         }
-<a class="jxr_linenumber" name="L306" href="#L306">306</a> 
-<a class="jxr_linenumber" name="L307" href="#L307">307</a>         <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L308" href="#L308">308</a> <em class="jxr_javadoccomment">         * Find a complex root for the polynomial with the given coefficients,</em>
-<a class="jxr_linenumber" name="L309" href="#L309">309</a> <em class="jxr_javadoccomment">         * starting from the given initial value.</em>
-<a class="jxr_linenumber" name="L310" href="#L310">310</a> <em class="jxr_javadoccomment">         *</em>
-<a class="jxr_linenumber" name="L311" href="#L311">311</a> <em class="jxr_javadoccomment">         * @param coefficients Polynomial coefficients.</em>
-<a class="jxr_linenumber" name="L312" href="#L312">312</a> <em class="jxr_javadoccomment">         * @param initial Start value.</em>
-<a class="jxr_linenumber" name="L313" href="#L313">313</a> <em class="jxr_javadoccomment">         * @return the point at which the function value is zero.</em>
-<a class="jxr_linenumber" name="L314" href="#L314">314</a> <em class="jxr_javadoccomment">         * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
-<a class="jxr_linenumber" name="L315" href="#L315">315</a> <em class="jxr_javadoccomment">         * if the maximum number of evaluations is exceeded.</em>
-<a class="jxr_linenumber" name="L316" href="#L316">316</a> <em class="jxr_javadoccomment">         * @throws NullArgumentException if the {@code coefficients} is</em>
-<a class="jxr_linenumber" name="L317" href="#L317">317</a> <em class="jxr_javadoccomment">         * {@code null}.</em>
-<a class="jxr_linenumber" name="L318" href="#L318">318</a> <em class="jxr_javadoccomment">         * @throws NoDataException if the {@code coefficients} array is empty.</em>
-<a class="jxr_linenumber" name="L319" href="#L319">319</a> <em class="jxr_javadoccomment">         */</em>
-<a class="jxr_linenumber" name="L320" href="#L320">320</a>         <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> solve(<a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> coefficients[], <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> initial)
-<a class="jxr_linenumber" name="L321" href="#L321">321</a>             <strong class="jxr_keyword">throws</strong> NullArgumentException,
-<a class="jxr_linenumber" name="L322" href="#L322">322</a>                    NoDataException,
-<a class="jxr_linenumber" name="L323" href="#L323">323</a>                    <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
-<a class="jxr_linenumber" name="L324" href="#L324">324</a>             <strong class="jxr_keyword">if</strong> (coefficients == <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L325" href="#L325">325</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NullArgumentException.html">NullArgumentException</a>();
-<a class="jxr_linenumber" name="L326" href="#L326">326</a>             }
-<a class="jxr_linenumber" name="L327" href="#L327">327</a> 
+<a class="jxr_linenumber" name="L197" href="#L197">197</a>        <strong class="jxr_keyword">return</strong> solveAllComplex(coefficients, initial, Integer.MAX_VALUE);
+<a class="jxr_linenumber" name="L198" href="#L198">198</a>     }
+<a class="jxr_linenumber" name="L199" href="#L199">199</a> 
+<a class="jxr_linenumber" name="L200" href="#L200">200</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L201" href="#L201">201</a> <em class="jxr_javadoccomment">     * Find all complex roots for the polynomial with the given</em>
+<a class="jxr_linenumber" name="L202" href="#L202">202</a> <em class="jxr_javadoccomment">     * coefficients, starting from the given initial value.</em>
+<a class="jxr_linenumber" name="L203" href="#L203">203</a> <em class="jxr_javadoccomment">     * &lt;br/&gt;</em>
+<a class="jxr_linenumber" name="L204" href="#L204">204</a> <em class="jxr_javadoccomment">     * Note: This method is not part of the API of {@link BaseUnivariateSolver}.</em>
+<a class="jxr_linenumber" name="L205" href="#L205">205</a> <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L206" href="#L206">206</a> <em class="jxr_javadoccomment">     * @param coefficients polynomial coefficients</em>
+<a class="jxr_linenumber" name="L207" href="#L207">207</a> <em class="jxr_javadoccomment">     * @param initial start value</em>
+<a class="jxr_linenumber" name="L208" href="#L208">208</a> <em class="jxr_javadoccomment">     * @param maxEval maximum number of evaluations</em>
+<a class="jxr_linenumber" name="L209" href="#L209">209</a> <em class="jxr_javadoccomment">     * @return the full set of complex roots of the polynomial</em>
+<a class="jxr_linenumber" name="L210" href="#L210">210</a> <em class="jxr_javadoccomment">     * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
+<a class="jxr_linenumber" name="L211" href="#L211">211</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded when solving for one of the roots</em>
+<a class="jxr_linenumber" name="L212" href="#L212">212</a> <em class="jxr_javadoccomment">     * @throws NullArgumentException if the {@code coefficients} is</em>
+<a class="jxr_linenumber" name="L213" href="#L213">213</a> <em class="jxr_javadoccomment">     * {@code null}</em>
+<a class="jxr_linenumber" name="L214" href="#L214">214</a> <em class="jxr_javadoccomment">     * @throws NoDataException if the {@code coefficients} array is empty</em>
+<a class="jxr_linenumber" name="L215" href="#L215">215</a> <em class="jxr_javadoccomment">     * @since 3.5</em>
+<a class="jxr_linenumber" name="L216" href="#L216">216</a> <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L217" href="#L217">217</a>     <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[] solveAllComplex(<strong class="jxr_keyword">double</strong>[] coefficients,
+<a class="jxr_linenumber" name="L218" href="#L218">218</a>                                      <strong class="jxr_keyword">double</strong> initial, <strong class="jxr_keyword">int</strong> maxEval)
+<a class="jxr_linenumber" name="L219" href="#L219">219</a>         <strong class="jxr_keyword">throws</strong> NullArgumentException,
+<a class="jxr_linenumber" name="L220" href="#L220">220</a>                NoDataException,
+<a class="jxr_linenumber" name="L221" href="#L221">221</a>                <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
+<a class="jxr_linenumber" name="L222" href="#L222">222</a>         setup(maxEval,
+<a class="jxr_linenumber" name="L223" href="#L223">223</a>               <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/analysis/polynomials/PolynomialFunction.html">PolynomialFunction</a>(coefficients),
+<a class="jxr_linenumber" name="L224" href="#L224">224</a>               Double.NEGATIVE_INFINITY,
+<a class="jxr_linenumber" name="L225" href="#L225">225</a>               Double.POSITIVE_INFINITY,
+<a class="jxr_linenumber" name="L226" href="#L226">226</a>               initial);
+<a class="jxr_linenumber" name="L227" href="#L227">227</a>         <strong class="jxr_keyword">return</strong> complexSolver.solveAll(ComplexUtils.convertToComplex(coefficients),
+<a class="jxr_linenumber" name="L228" href="#L228">228</a>                                       <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(initial, 0d));
+<a class="jxr_linenumber" name="L229" href="#L229">229</a>     }
+<a class="jxr_linenumber" name="L230" href="#L230">230</a> 
+<a class="jxr_linenumber" name="L231" href="#L231">231</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L232" href="#L232">232</a> <em class="jxr_javadoccomment">     * Find a complex root for the polynomial with the given coefficients,</em>
+<a class="jxr_linenumber" name="L233" href="#L233">233</a> <em class="jxr_javadoccomment">     * starting from the given initial value.</em>
+<a class="jxr_linenumber" name="L234" href="#L234">234</a> <em class="jxr_javadoccomment">     * &lt;br/&gt;</em>
+<a class="jxr_linenumber" name="L235" href="#L235">235</a> <em class="jxr_javadoccomment">     * Note: This method is not part of the API of {@link BaseUnivariateSolver}.</em>
+<a class="jxr_linenumber" name="L236" href="#L236">236</a> <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L237" href="#L237">237</a> <em class="jxr_javadoccomment">     * @param coefficients Polynomial coefficients.</em>
+<a class="jxr_linenumber" name="L238" href="#L238">238</a> <em class="jxr_javadoccomment">     * @param initial Start value.</em>
+<a class="jxr_linenumber" name="L239" href="#L239">239</a> <em class="jxr_javadoccomment">     * @return a complex root of the polynomial</em>
+<a class="jxr_linenumber" name="L240" href="#L240">240</a> <em class="jxr_javadoccomment">     * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
+<a class="jxr_linenumber" name="L241" href="#L241">241</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded.</em>
+<a class="jxr_linenumber" name="L242" href="#L242">242</a> <em class="jxr_javadoccomment">     * @throws NullArgumentException if the {@code coefficients} is</em>
+<a class="jxr_linenumber" name="L243" href="#L243">243</a> <em class="jxr_javadoccomment">     * {@code null}.</em>
+<a class="jxr_linenumber" name="L244" href="#L244">244</a> <em class="jxr_javadoccomment">     * @throws NoDataException if the {@code coefficients} array is empty.</em>
+<a class="jxr_linenumber" name="L245" href="#L245">245</a> <em class="jxr_javadoccomment">     * @since 3.1</em>
+<a class="jxr_linenumber" name="L246" href="#L246">246</a> <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L247" href="#L247">247</a>     <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> solveComplex(<strong class="jxr_keyword">double</strong>[] coefficients,
+<a class="jxr_linenumber" name="L248" href="#L248">248</a>                                 <strong class="jxr_keyword">double</strong> initial)
+<a class="jxr_linenumber" name="L249" href="#L249">249</a>         <strong class="jxr_keyword">throws</strong> NullArgumentException,
+<a class="jxr_linenumber" name="L250" href="#L250">250</a>                NoDataException,
+<a class="jxr_linenumber" name="L251" href="#L251">251</a>                <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
+<a class="jxr_linenumber" name="L252" href="#L252">252</a>        <strong class="jxr_keyword">return</strong> solveComplex(coefficients, initial, Integer.MAX_VALUE);
+<a class="jxr_linenumber" name="L253" href="#L253">253</a>     }
+<a class="jxr_linenumber" name="L254" href="#L254">254</a> 
+<a class="jxr_linenumber" name="L255" href="#L255">255</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L256" href="#L256">256</a> <em class="jxr_javadoccomment">     * Find a complex root for the polynomial with the given coefficients,</em>
+<a class="jxr_linenumber" name="L257" href="#L257">257</a> <em class="jxr_javadoccomment">     * starting from the given initial value.</em>
+<a class="jxr_linenumber" name="L258" href="#L258">258</a> <em class="jxr_javadoccomment">     * &lt;br/&gt;</em>
+<a class="jxr_linenumber" name="L259" href="#L259">259</a> <em class="jxr_javadoccomment">     * Note: This method is not part of the API of {@link BaseUnivariateSolver}.</em>
+<a class="jxr_linenumber" name="L260" href="#L260">260</a> <em class="jxr_javadoccomment">     *</em>
+<a class="jxr_linenumber" name="L261" href="#L261">261</a> <em class="jxr_javadoccomment">     * @param coefficients polynomial coefficients</em>
+<a class="jxr_linenumber" name="L262" href="#L262">262</a> <em class="jxr_javadoccomment">     * @param initial start value</em>
+<a class="jxr_linenumber" name="L263" href="#L263">263</a> <em class="jxr_javadoccomment">     * @param maxEval maximum number of evaluations</em>
+<a class="jxr_linenumber" name="L264" href="#L264">264</a> <em class="jxr_javadoccomment">     * @return a complex root of the polynomial</em>
+<a class="jxr_linenumber" name="L265" href="#L265">265</a> <em class="jxr_javadoccomment">     * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
+<a class="jxr_linenumber" name="L266" href="#L266">266</a> <em class="jxr_javadoccomment">     * if the maximum number of evaluations is exceeded</em>
+<a class="jxr_linenumber" name="L267" href="#L267">267</a> <em class="jxr_javadoccomment">     * @throws NullArgumentException if the {@code coefficients} is</em>
+<a class="jxr_linenumber" name="L268" href="#L268">268</a> <em class="jxr_javadoccomment">     * {@code null}</em>
+<a class="jxr_linenumber" name="L269" href="#L269">269</a> <em class="jxr_javadoccomment">     * @throws NoDataException if the {@code coefficients} array is empty</em>
+<a class="jxr_linenumber" name="L270" href="#L270">270</a> <em class="jxr_javadoccomment">     * @since 3.1</em>
+<a class="jxr_linenumber" name="L271" href="#L271">271</a> <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L272" href="#L272">272</a>     <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> solveComplex(<strong class="jxr_keyword">double</strong>[] coefficients,
+<a class="jxr_linenumber" name="L273" href="#L273">273</a>                                 <strong class="jxr_keyword">double</strong> initial, <strong class="jxr_keyword">int</strong> maxEval)
+<a class="jxr_linenumber" name="L274" href="#L274">274</a>         <strong class="jxr_keyword">throws</strong> NullArgumentException,
+<a class="jxr_linenumber" name="L275" href="#L275">275</a>                NoDataException,
+<a class="jxr_linenumber" name="L276" href="#L276">276</a>                <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
+<a class="jxr_linenumber" name="L277" href="#L277">277</a>         setup(maxEval,
+<a class="jxr_linenumber" name="L278" href="#L278">278</a>               <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/analysis/polynomials/PolynomialFunction.html">PolynomialFunction</a>(coefficients),
+<a class="jxr_linenumber" name="L279" href="#L279">279</a>               Double.NEGATIVE_INFINITY,
+<a class="jxr_linenumber" name="L280" href="#L280">280</a>               Double.POSITIVE_INFINITY,
+<a class="jxr_linenumber" name="L281" href="#L281">281</a>               initial);
+<a class="jxr_linenumber" name="L282" href="#L282">282</a>         <strong class="jxr_keyword">return</strong> complexSolver.solve(ComplexUtils.convertToComplex(coefficients),
+<a class="jxr_linenumber" name="L283" href="#L283">283</a>                                    <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(initial, 0d));
+<a class="jxr_linenumber" name="L284" href="#L284">284</a>     }
+<a class="jxr_linenumber" name="L285" href="#L285">285</a> 
+<a class="jxr_linenumber" name="L286" href="#L286">286</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L287" href="#L287">287</a> <em class="jxr_javadoccomment">     * Class for searching all (complex) roots.</em>
+<a class="jxr_linenumber" name="L288" href="#L288">288</a> <em class="jxr_javadoccomment">     */</em>
+<a class="jxr_linenumber" name="L289" href="#L289">289</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/math3/analysis/solvers/LaguerreSolver.html">ComplexSolver</a> {
+<a class="jxr_linenumber" name="L290" href="#L290">290</a>         <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L291" href="#L291">291</a> <em class="jxr_javadoccomment">         * Check whether the given complex root is actually a real zero</em>
+<a class="jxr_linenumber" name="L292" href="#L292">292</a> <em class="jxr_javadoccomment">         * in the given interval, within the solver tolerance level.</em>
+<a class="jxr_linenumber" name="L293" href="#L293">293</a> <em class="jxr_javadoccomment">         *</em>
+<a class="jxr_linenumber" name="L294" href="#L294">294</a> <em class="jxr_javadoccomment">         * @param min Lower bound for the interval.</em>
+<a class="jxr_linenumber" name="L295" href="#L295">295</a> <em class="jxr_javadoccomment">         * @param max Upper bound for the interval.</em>
+<a class="jxr_linenumber" name="L296" href="#L296">296</a> <em class="jxr_javadoccomment">         * @param z Complex root.</em>
+<a class="jxr_linenumber" name="L297" href="#L297">297</a> <em class="jxr_javadoccomment">         * @return {@code true} if z is a real zero.</em>
+<a class="jxr_linenumber" name="L298" href="#L298">298</a> <em class="jxr_javadoccomment">         */</em>
+<a class="jxr_linenumber" name="L299" href="#L299">299</a>         <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">boolean</strong> isRoot(<strong class="jxr_keyword">double</strong> min, <strong class="jxr_keyword">double</strong> max, <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> z) {
+<a class="jxr_linenumber" name="L300" href="#L300">300</a>             <strong class="jxr_keyword">if</strong> (isSequence(min, z.getReal(), max)) {
+<a class="jxr_linenumber" name="L301" href="#L301">301</a>                 <strong class="jxr_keyword">double</strong> tolerance = FastMath.max(getRelativeAccuracy() * z.abs(), getAbsoluteAccuracy());
+<a class="jxr_linenumber" name="L302" href="#L302">302</a>                 <strong class="jxr_keyword">return</strong> (FastMath.abs(z.getImaginary()) &lt;= tolerance) ||
+<a class="jxr_linenumber" name="L303" href="#L303">303</a>                      (z.abs() &lt;= getFunctionValueAccuracy());
+<a class="jxr_linenumber" name="L304" href="#L304">304</a>             }
+<a class="jxr_linenumber" name="L305" href="#L305">305</a>             <strong class="jxr_keyword">return</strong> false;
+<a class="jxr_linenumber" name="L306" href="#L306">306</a>         }
+<a class="jxr_linenumber" name="L307" href="#L307">307</a> 
+<a class="jxr_linenumber" name="L308" href="#L308">308</a>         <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L309" href="#L309">309</a> <em class="jxr_javadoccomment">         * Find all complex roots for the polynomial with the given</em>
+<a class="jxr_linenumber" name="L310" href="#L310">310</a> <em class="jxr_javadoccomment">         * coefficients, starting from the given initial value.</em>
+<a class="jxr_linenumber" name="L311" href="#L311">311</a> <em class="jxr_javadoccomment">         *</em>
+<a class="jxr_linenumber" name="L312" href="#L312">312</a> <em class="jxr_javadoccomment">         * @param coefficients Polynomial coefficients.</em>
+<a class="jxr_linenumber" name="L313" href="#L313">313</a> <em class="jxr_javadoccomment">         * @param initial Start value.</em>
+<a class="jxr_linenumber" name="L314" href="#L314">314</a> <em class="jxr_javadoccomment">         * @return the point at which the function value is zero.</em>
+<a class="jxr_linenumber" name="L315" href="#L315">315</a> <em class="jxr_javadoccomment">         * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
+<a class="jxr_linenumber" name="L316" href="#L316">316</a> <em class="jxr_javadoccomment">         * if the maximum number of evaluations is exceeded.</em>
+<a class="jxr_linenumber" name="L317" href="#L317">317</a> <em class="jxr_javadoccomment">         * @throws NullArgumentException if the {@code coefficients} is</em>
+<a class="jxr_linenumber" name="L318" href="#L318">318</a> <em class="jxr_javadoccomment">         * {@code null}.</em>
+<a class="jxr_linenumber" name="L319" href="#L319">319</a> <em class="jxr_javadoccomment">         * @throws NoDataException if the {@code coefficients} array is empty.</em>
+<a class="jxr_linenumber" name="L320" href="#L320">320</a> <em class="jxr_javadoccomment">         */</em>
+<a class="jxr_linenumber" name="L321" href="#L321">321</a>         <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[] solveAll(<a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> coefficients[], <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> initial)
+<a class="jxr_linenumber" name="L322" href="#L322">322</a>             <strong class="jxr_keyword">throws</strong> NullArgumentException,
+<a class="jxr_linenumber" name="L323" href="#L323">323</a>                    NoDataException,
+<a class="jxr_linenumber" name="L324" href="#L324">324</a>                    <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
+<a class="jxr_linenumber" name="L325" href="#L325">325</a>             <strong class="jxr_keyword">if</strong> (coefficients == <strong class="jxr_keyword">null</strong>) {
+<a class="jxr_linenumber" name="L326" href="#L326">326</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NullArgumentException.html">NullArgumentException</a>();
+<a class="jxr_linenumber" name="L327" href="#L327">327</a>             }
 <a class="jxr_linenumber" name="L328" href="#L328">328</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = coefficients.length - 1;
 <a class="jxr_linenumber" name="L329" href="#L329">329</a>             <strong class="jxr_keyword">if</strong> (n == 0) {
 <a class="jxr_linenumber" name="L330" href="#L330">330</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NoDataException.html">NoDataException</a>(LocalizedFormats.POLYNOMIAL);
 <a class="jxr_linenumber" name="L331" href="#L331">331</a>             }
-<a class="jxr_linenumber" name="L332" href="#L332">332</a> 
-<a class="jxr_linenumber" name="L333" href="#L333">333</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> absoluteAccuracy = getAbsoluteAccuracy();
-<a class="jxr_linenumber" name="L334" href="#L334">334</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> relativeAccuracy = getRelativeAccuracy();
-<a class="jxr_linenumber" name="L335" href="#L335">335</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> functionValueAccuracy = getFunctionValueAccuracy();
-<a class="jxr_linenumber" name="L336" href="#L336">336</a> 
-<a class="jxr_linenumber" name="L337" href="#L337">337</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> nC  = <strong class="jxr_keyword">new</strong> Complex(n, 0);
-<a class="jxr_linenumber" name="L338" href="#L338">338</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> n1C = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(n - 1, 0);
-<a class="jxr_linenumber" name="L339" href="#L339">339</a> 
-<a class="jxr_linenumber" name="L340" href="#L340">340</a>             <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> z = initial;
-<a class="jxr_linenumber" name="L341" href="#L341">341</a>             <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> oldz = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(Double.POSITIVE_INFINITY,
-<a class="jxr_linenumber" name="L342" href="#L342">342</a>                                        Double.POSITIVE_INFINITY);
-<a class="jxr_linenumber" name="L343" href="#L343">343</a>             <strong class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>) {
-<a class="jxr_linenumber" name="L344" href="#L344">344</a>                 <em class="jxr_comment">// Compute pv (polynomial value), dv (derivative value), and</em>
-<a class="jxr_linenumber" name="L345" href="#L345">345</a>                 <em class="jxr_comment">// d2v (second derivative value) simultaneously.</em>
-<a class="jxr_linenumber" name="L346" href="#L346">346</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> pv = coefficients[n];
-<a class="jxr_linenumber" name="L347" href="#L347">347</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dv = Complex.ZERO;
-<a class="jxr_linenumber" name="L348" href="#L348">348</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> d2v = Complex.ZERO;
-<a class="jxr_linenumber" name="L349" href="#L349">349</a>                 <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = n-1; j &gt;= 0; j--) {
-<a class="jxr_linenumber" name="L350" href="#L350">350</a>                     d2v = dv.add(z.multiply(d2v));
-<a class="jxr_linenumber" name="L351" href="#L351">351</a>                     dv = pv.add(z.multiply(dv));
-<a class="jxr_linenumber" name="L352" href="#L352">352</a>                     pv = coefficients[j].add(z.multiply(pv));
-<a class="jxr_linenumber" name="L353" href="#L353">353</a>                 }
-<a class="jxr_linenumber" name="L354" href="#L354">354</a>                 d2v = d2v.multiply(<strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(2.0, 0.0));
-<a class="jxr_linenumber" name="L355" href="#L355">355</a> 
-<a class="jxr_linenumber" name="L356" href="#L356">356</a>                 <em class="jxr_comment">// Check for convergence.</em>
-<a class="jxr_linenumber" name="L357" href="#L357">357</a>                 <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tolerance = FastMath.max(relativeAccuracy * z.abs(),
-<a class="jxr_linenumber" name="L358" href="#L358">358</a>                                                       absoluteAccuracy);
-<a class="jxr_linenumber" name="L359" href="#L359">359</a>                 <strong class="jxr_keyword">if</strong> ((z.subtract(oldz)).abs() &lt;= tolerance) {
-<a class="jxr_linenumber" name="L360" href="#L360">360</a>                     <strong class="jxr_keyword">return</strong> z;
-<a class="jxr_linenumber" name="L361" href="#L361">361</a>                 }
-<a class="jxr_linenumber" name="L362" href="#L362">362</a>                 <strong class="jxr_keyword">if</strong> (pv.abs() &lt;= functionValueAccuracy) {
-<a class="jxr_linenumber" name="L363" href="#L363">363</a>                     <strong class="jxr_keyword">return</strong> z;
-<a class="jxr_linenumber" name="L364" href="#L364">364</a>                 }
-<a class="jxr_linenumber" name="L365" href="#L365">365</a> 
-<a class="jxr_linenumber" name="L366" href="#L366">366</a>                 <em class="jxr_comment">// Now pv != 0, calculate the new approximation.</em>
-<a class="jxr_linenumber" name="L367" href="#L367">367</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> G = dv.divide(pv);
-<a class="jxr_linenumber" name="L368" href="#L368">368</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> G2 = G.multiply(G);
-<a class="jxr_linenumber" name="L369" href="#L369">369</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> H = G2.subtract(d2v.divide(pv));
-<a class="jxr_linenumber" name="L370" href="#L370">370</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> delta = n1C.multiply((nC.multiply(H)).subtract(G2));
-<a class="jxr_linenumber" name="L371" href="#L371">371</a>                 <em class="jxr_comment">// Choose a denominator larger in magnitude.</em>
-<a class="jxr_linenumber" name="L372" href="#L372">372</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> deltaSqrt = delta.sqrt();
-<a class="jxr_linenumber" name="L373" href="#L373">373</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dplus = G.add(deltaSqrt);
-<a class="jxr_linenumber" name="L374" href="#L374">374</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dminus = G.subtract(deltaSqrt);
-<a class="jxr_linenumber" name="L375" href="#L375">375</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> denominator = dplus.abs() &gt; dminus.abs() ? dplus : dminus;
-<a class="jxr_linenumber" name="L376" href="#L376">376</a>                 <em class="jxr_comment">// Perturb z if denominator is zero, for instance,</em>
-<a class="jxr_linenumber" name="L377" href="#L377">377</a>                 <em class="jxr_comment">// p(x) = x^3 + 1, z = 0.</em>
-<a class="jxr_linenumber" name="L378" href="#L378">378</a>                 <strong class="jxr_keyword">if</strong> (denominator.equals(<strong class="jxr_keyword">new</strong> Complex(0.0, 0.0))) {
-<a class="jxr_linenumber" name="L379" href="#L379">379</a>                     z = z.add(<strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(absoluteAccuracy, absoluteAccuracy));
-<a class="jxr_linenumber" name="L380" href="#L380">380</a>                     oldz = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(Double.POSITIVE_INFINITY,
-<a class="jxr_linenumber" name="L381" href="#L381">381</a>                                        Double.POSITIVE_INFINITY);
-<a class="jxr_linenumber" name="L382" href="#L382">382</a>                 } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L383" href="#L383">383</a>                     oldz = z;
-<a class="jxr_linenumber" name="L384" href="#L384">384</a>                     z = z.subtract(nC.divide(denominator));
-<a class="jxr_linenumber" name="L385" href="#L385">385</a>                 }
-<a class="jxr_linenumber" name="L386" href="#L386">386</a>                 incrementEvaluationCount();
-<a class="jxr_linenumber" name="L387" href="#L387">387</a>             }
-<a class="jxr_linenumber" name="L388" href="#L388">388</a>         }
-<a class="jxr_linenumber" name="L389" href="#L389">389</a>     }
-<a class="jxr_linenumber" name="L390" href="#L390">390</a> }
+<a class="jxr_linenumber" name="L332" href="#L332">332</a>             <em class="jxr_comment">// Coefficients for deflated polynomial.</em>
+<a class="jxr_linenumber" name="L333" href="#L333">333</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> c[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n + 1];
+<a class="jxr_linenumber" name="L334" href="#L334">334</a>             <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt;= n; i++) {
+<a class="jxr_linenumber" name="L335" href="#L335">335</a>                 c[i] = coefficients[i];
+<a class="jxr_linenumber" name="L336" href="#L336">336</a>             }
+<a class="jxr_linenumber" name="L337" href="#L337">337</a> 
+<a class="jxr_linenumber" name="L338" href="#L338">338</a>             <em class="jxr_comment">// Solve individual roots successively.</em>
+<a class="jxr_linenumber" name="L339" href="#L339">339</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> root[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n];
+<a class="jxr_linenumber" name="L340" href="#L340">340</a>             <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; n; i++) {
+<a class="jxr_linenumber" name="L341" href="#L341">341</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> subarray[] = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>[n - i + 1];
+<a class="jxr_linenumber" name="L342" href="#L342">342</a>                 System.arraycopy(c, 0, subarray, 0, subarray.length);
+<a class="jxr_linenumber" name="L343" href="#L343">343</a>                 root[i] = solve(subarray, initial);
+<a class="jxr_linenumber" name="L344" href="#L344">344</a>                 <em class="jxr_comment">// Polynomial deflation using synthetic division.</em>
+<a class="jxr_linenumber" name="L345" href="#L345">345</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> newc = c[n - i];
+<a class="jxr_linenumber" name="L346" href="#L346">346</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> oldc = <strong class="jxr_keyword">null</strong>;
+<a class="jxr_linenumber" name="L347" href="#L347">347</a>                 <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = n - i - 1; j &gt;= 0; j--) {
+<a class="jxr_linenumber" name="L348" href="#L348">348</a>                     oldc = c[j];
+<a class="jxr_linenumber" name="L349" href="#L349">349</a>                     c[j] = newc;
+<a class="jxr_linenumber" name="L350" href="#L350">350</a>                     newc = oldc.add(newc.multiply(root[i]));
+<a class="jxr_linenumber" name="L351" href="#L351">351</a>                 }
+<a class="jxr_linenumber" name="L352" href="#L352">352</a>             }
+<a class="jxr_linenumber" name="L353" href="#L353">353</a> 
+<a class="jxr_linenumber" name="L354" href="#L354">354</a>             <strong class="jxr_keyword">return</strong> root;
+<a class="jxr_linenumber" name="L355" href="#L355">355</a>         }
+<a class="jxr_linenumber" name="L356" href="#L356">356</a> 
+<a class="jxr_linenumber" name="L357" href="#L357">357</a>         <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L358" href="#L358">358</a> <em class="jxr_javadoccomment">         * Find a complex root for the polynomial with the given coefficients,</em>
+<a class="jxr_linenumber" name="L359" href="#L359">359</a> <em class="jxr_javadoccomment">         * starting from the given initial value.</em>
+<a class="jxr_linenumber" name="L360" href="#L360">360</a> <em class="jxr_javadoccomment">         *</em>
+<a class="jxr_linenumber" name="L361" href="#L361">361</a> <em class="jxr_javadoccomment">         * @param coefficients Polynomial coefficients.</em>
+<a class="jxr_linenumber" name="L362" href="#L362">362</a> <em class="jxr_javadoccomment">         * @param initial Start value.</em>
+<a class="jxr_linenumber" name="L363" href="#L363">363</a> <em class="jxr_javadoccomment">         * @return the point at which the function value is zero.</em>
+<a class="jxr_linenumber" name="L364" href="#L364">364</a> <em class="jxr_javadoccomment">         * @throws org.apache.commons.math3.exception.TooManyEvaluationsException</em>
+<a class="jxr_linenumber" name="L365" href="#L365">365</a> <em class="jxr_javadoccomment">         * if the maximum number of evaluations is exceeded.</em>
+<a class="jxr_linenumber" name="L366" href="#L366">366</a> <em class="jxr_javadoccomment">         * @throws NullArgumentException if the {@code coefficients} is</em>
+<a class="jxr_linenumber" name="L367" href="#L367">367</a> <em class="jxr_javadoccomment">         * {@code null}.</em>
+<a class="jxr_linenumber" name="L368" href="#L368">368</a> <em class="jxr_javadoccomment">         * @throws NoDataException if the {@code coefficients} array is empty.</em>
+<a class="jxr_linenumber" name="L369" href="#L369">369</a> <em class="jxr_javadoccomment">         */</em>
+<a class="jxr_linenumber" name="L370" href="#L370">370</a>         <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> solve(<a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> coefficients[], <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> initial)
+<a class="jxr_linenumber" name="L371" href="#L371">371</a>             <strong class="jxr_keyword">throws</strong> NullArgumentException,
+<a class="jxr_linenumber" name="L372" href="#L372">372</a>                    NoDataException,
+<a class="jxr_linenumber" name="L373" href="#L373">373</a>                    <a href="../../../../../../org/apache/commons/math3/exception/TooManyEvaluationsException.html">TooManyEvaluationsException</a> {
+<a class="jxr_linenumber" name="L374" href="#L374">374</a>             <strong class="jxr_keyword">if</strong> (coefficients == <strong class="jxr_keyword">null</strong>) {
+<a class="jxr_linenumber" name="L375" href="#L375">375</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NullArgumentException.html">NullArgumentException</a>();
+<a class="jxr_linenumber" name="L376" href="#L376">376</a>             }
+<a class="jxr_linenumber" name="L377" href="#L377">377</a> 
+<a class="jxr_linenumber" name="L378" href="#L378">378</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = coefficients.length - 1;
+<a class="jxr_linenumber" name="L379" href="#L379">379</a>             <strong class="jxr_keyword">if</strong> (n == 0) {
+<a class="jxr_linenumber" name="L380" href="#L380">380</a>                 <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/exception/NoDataException.html">NoDataException</a>(LocalizedFormats.POLYNOMIAL);
+<a class="jxr_linenumber" name="L381" href="#L381">381</a>             }
+<a class="jxr_linenumber" name="L382" href="#L382">382</a> 
+<a class="jxr_linenumber" name="L383" href="#L383">383</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> absoluteAccuracy = getAbsoluteAccuracy();
+<a class="jxr_linenumber" name="L384" href="#L384">384</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> relativeAccuracy = getRelativeAccuracy();
+<a class="jxr_linenumber" name="L385" href="#L385">385</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> functionValueAccuracy = getFunctionValueAccuracy();
+<a class="jxr_linenumber" name="L386" href="#L386">386</a> 
+<a class="jxr_linenumber" name="L387" href="#L387">387</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> nC  = <strong class="jxr_keyword">new</strong> Complex(n, 0);
+<a class="jxr_linenumber" name="L388" href="#L388">388</a>             <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> n1C = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(n - 1, 0);
+<a class="jxr_linenumber" name="L389" href="#L389">389</a> 
+<a class="jxr_linenumber" name="L390" href="#L390">390</a>             <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> z = initial;
+<a class="jxr_linenumber" name="L391" href="#L391">391</a>             <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> oldz = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(Double.POSITIVE_INFINITY,
+<a class="jxr_linenumber" name="L392" href="#L392">392</a>                                        Double.POSITIVE_INFINITY);
+<a class="jxr_linenumber" name="L393" href="#L393">393</a>             <strong class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>) {
+<a class="jxr_linenumber" name="L394" href="#L394">394</a>                 <em class="jxr_comment">// Compute pv (polynomial value), dv (derivative value), and</em>
+<a class="jxr_linenumber" name="L395" href="#L395">395</a>                 <em class="jxr_comment">// d2v (second derivative value) simultaneously.</em>
+<a class="jxr_linenumber" name="L396" href="#L396">396</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> pv = coefficients[n];
+<a class="jxr_linenumber" name="L397" href="#L397">397</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dv = Complex.ZERO;
+<a class="jxr_linenumber" name="L398" href="#L398">398</a>                 <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> d2v = Complex.ZERO;
+<a class="jxr_linenumber" name="L399" href="#L399">399</a>                 <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = n-1; j &gt;= 0; j--) {
+<a class="jxr_linenumber" name="L400" href="#L400">400</a>                     d2v = dv.add(z.multiply(d2v));
+<a class="jxr_linenumber" name="L401" href="#L401">401</a>                     dv = pv.add(z.multiply(dv));
+<a class="jxr_linenumber" name="L402" href="#L402">402</a>                     pv = coefficients[j].add(z.multiply(pv));
+<a class="jxr_linenumber" name="L403" href="#L403">403</a>                 }
+<a class="jxr_linenumber" name="L404" href="#L404">404</a>                 d2v = d2v.multiply(<strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(2.0, 0.0));
+<a class="jxr_linenumber" name="L405" href="#L405">405</a> 
+<a class="jxr_linenumber" name="L406" href="#L406">406</a>                 <em class="jxr_comment">// Check for convergence.</em>
+<a class="jxr_linenumber" name="L407" href="#L407">407</a>                 <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tolerance = FastMath.max(relativeAccuracy * z.abs(),
+<a class="jxr_linenumber" name="L408" href="#L408">408</a>                                                       absoluteAccuracy);
+<a class="jxr_linenumber" name="L409" href="#L409">409</a>                 <strong class="jxr_keyword">if</strong> ((z.subtract(oldz)).abs() &lt;= tolerance) {
+<a class="jxr_linenumber" name="L410" href="#L410">410</a>                     <strong class="jxr_keyword">return</strong> z;
+<a class="jxr_linenumber" name="L411" href="#L411">411</a>                 }
+<a class="jxr_linenumber" name="L412" href="#L412">412</a>                 <strong class="jxr_keyword">if</strong> (pv.abs() &lt;= functionValueAccuracy) {
+<a class="jxr_linenumber" name="L413" href="#L413">413</a>                     <strong class="jxr_keyword">return</strong> z;
+<a class="jxr_linenumber" name="L414" href="#L414">414</a>                 }
+<a class="jxr_linenumber" name="L415" href="#L415">415</a> 
+<a class="jxr_linenumber" name="L416" href="#L416">416</a>                 <em class="jxr_comment">// Now pv != 0, calculate the new approximation.</em>
+<a class="jxr_linenumber" name="L417" href="#L417">417</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> G = dv.divide(pv);
+<a class="jxr_linenumber" name="L418" href="#L418">418</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> G2 = G.multiply(G);
+<a class="jxr_linenumber" name="L419" href="#L419">419</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> H = G2.subtract(d2v.divide(pv));
+<a class="jxr_linenumber" name="L420" href="#L420">420</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> delta = n1C.multiply((nC.multiply(H)).subtract(G2));
+<a class="jxr_linenumber" name="L421" href="#L421">421</a>                 <em class="jxr_comment">// Choose a denominator larger in magnitude.</em>
+<a class="jxr_linenumber" name="L422" href="#L422">422</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> deltaSqrt = delta.sqrt();
+<a class="jxr_linenumber" name="L423" href="#L423">423</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dplus = G.add(deltaSqrt);
+<a class="jxr_linenumber" name="L424" href="#L424">424</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> dminus = G.subtract(deltaSqrt);
+<a class="jxr_linenumber" name="L425" href="#L425">425</a>                 <strong class="jxr_keyword">final</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a> denominator = dplus.abs() &gt; dminus.abs() ? dplus : dminus;
+<a class="jxr_linenumber" name="L426" href="#L426">426</a>                 <em class="jxr_comment">// Perturb z if denominator is zero, for instance,</em>
+<a class="jxr_linenumber" name="L427" href="#L427">427</a>                 <em class="jxr_comment">// p(x) = x^3 + 1, z = 0.</em>
+<a class="jxr_linenumber" name="L428" href="#L428">428</a>                 <strong class="jxr_keyword">if</strong> (denominator.equals(<strong class="jxr_keyword">new</strong> Complex(0.0, 0.0))) {
+<a class="jxr_linenumber" name="L429" href="#L429">429</a>                     z = z.add(<strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(absoluteAccuracy, absoluteAccuracy));
+<a class="jxr_linenumber" name="L430" href="#L430">430</a>                     oldz = <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/math3/complex/Complex.html">Complex</a>(Double.POSITIVE_INFINITY,
+<a class="jxr_linenumber" name="L431" href="#L431">431</a>                                        Double.POSITIVE_INFINITY);
+<a class="jxr_linenumber" name="L432" href="#L432">432</a>                 } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L433" href="#L433">433</a>                     oldz = z;
+<a class="jxr_linenumber" name="L434" href="#L434">434</a>                     z = z.subtract(nC.divide(denominator));
+<a class="jxr_linenumber" name="L435" href="#L435">435</a>                 }
+<a class="jxr_linenumber" name="L436" href="#L436">436</a>                 incrementEvaluationCount();
+<a class="jxr_linenumber" name="L437" href="#L437">437</a>             }
+<a class="jxr_linenumber" name="L438" href="#L438">438</a>         }
+<a class="jxr_linenumber" name="L439" href="#L439">439</a>     }
+<a class="jxr_linenumber" name="L440" href="#L440">440</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>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/UnivariateSolverUtils.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/UnivariateSolverUtils.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/UnivariateSolverUtils.html Fri Apr 17 19:17:13 2015
@@ -251,7 +251,7 @@
 <a class="jxr_linenumber" name="L243" href="#L243">243</a> <em class="jxr_javadoccomment">     * \( \delta_{k+1} = r \delta_k + q, \delta_0 = 0\) and starting search with \( k=1 \).</em>
 <a class="jxr_linenumber" name="L244" href="#L244">244</a> <em class="jxr_javadoccomment">     * The algorithm stops when one of the following happens: &lt;ul&gt;</em>
 <a class="jxr_linenumber" name="L245" href="#L245">245</a> <em class="jxr_javadoccomment">     * &lt;li&gt; at least one positive and one negative value have been found --  success!&lt;/li&gt;</em>
-<a class="jxr_linenumber" name="L246" href="#L246">246</a> <em class="jxr_javadoccomment">     * &lt;li&gt; both endpoints have reached their respective limites -- NoBracketingException &lt;/li&gt;</em>
+<a class="jxr_linenumber" name="L246" href="#L246">246</a> <em class="jxr_javadoccomment">     * &lt;li&gt; both endpoints have reached their respective limits -- NoBracketingException &lt;/li&gt;</em>
 <a class="jxr_linenumber" name="L247" href="#L247">247</a> <em class="jxr_javadoccomment">     * &lt;li&gt; {@code maximumIterations} iterations elapse -- NoBracketingException &lt;/li&gt;&lt;/ul&gt;&lt;/p&gt;</em>
 <a class="jxr_linenumber" name="L248" href="#L248">248</a> <em class="jxr_javadoccomment">     * &lt;p&gt;</em>
 <a class="jxr_linenumber" name="L249" href="#L249">249</a> <em class="jxr_javadoccomment">     * If different signs are found at first iteration ({@code k=1}), then the returned</em>
@@ -265,7 +265,7 @@
 <a class="jxr_linenumber" name="L257" href="#L257">257</a> <em class="jxr_javadoccomment">     * Interval expansion rate is tuned by changing the recurrence parameters {@code r} and</em>
 <a class="jxr_linenumber" name="L258" href="#L258">258</a> <em class="jxr_javadoccomment">     * {@code q}. When the multiplicative factor {@code r} is set to 1, the sequence is a</em>
 <a class="jxr_linenumber" name="L259" href="#L259">259</a> <em class="jxr_javadoccomment">     * simple arithmetic sequence with linear increase. When the multiplicative factor {@code r}</em>
-<a class="jxr_linenumber" name="L260" href="#L260">260</a> <em class="jxr_javadoccomment">     * is larger than 1, the sequence has an asymtotically exponential rate. Note than the</em>
+<a class="jxr_linenumber" name="L260" href="#L260">260</a> <em class="jxr_javadoccomment">     * is larger than 1, the sequence has an asymptotically exponential rate. Note than the</em>
 <a class="jxr_linenumber" name="L261" href="#L261">261</a> <em class="jxr_javadoccomment">     * additive parameter {@code q} should never be set to zero, otherwise the interval would</em>
 <a class="jxr_linenumber" name="L262" href="#L262">262</a> <em class="jxr_javadoccomment">     * degenerate to the single initial point for all values of {@code k}.</em>
 <a class="jxr_linenumber" name="L263" href="#L263">263</a> <em class="jxr_javadoccomment">     * &lt;/p&gt;</em>
@@ -322,7 +322,7 @@
 <a class="jxr_linenumber" name="L314" href="#L314">314</a>         <strong class="jxr_keyword">double</strong> delta = 0;
 <a class="jxr_linenumber" name="L315" href="#L315">315</a> 
 <a class="jxr_linenumber" name="L316" href="#L316">316</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> numIterations = 0;
-<a class="jxr_linenumber" name="L317" href="#L317">317</a>              (numIterations &lt; maximumIterations) &amp;&amp; (a &gt; lowerBound || b &gt; upperBound);
+<a class="jxr_linenumber" name="L317" href="#L317">317</a>              (numIterations &lt; maximumIterations) &amp;&amp; (a &gt; lowerBound || b &lt; upperBound);
 <a class="jxr_linenumber" name="L318" href="#L318">318</a>              ++numIterations) {
 <a class="jxr_linenumber" name="L319" href="#L319">319</a> 
 <a class="jxr_linenumber" name="L320" href="#L320">320</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> previousA  = a;

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.analysis.solvers</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.analysis.solvers</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/analysis/solvers/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.analysis.solvers</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.analysis.solvers</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.complex</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.complex</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/complex/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.complex</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.complex</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.dfp</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.dfp</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/dfp/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.dfp</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.dfp</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/PoissonDistribution.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/PoissonDistribution.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/PoissonDistribution.html Fri Apr 17 19:17:13 2015
@@ -291,14 +291,14 @@
 <a class="jxr_linenumber" name="L283" href="#L283">283</a> <em class="jxr_javadoccomment">     * &lt;ul&gt;</em>
 <a class="jxr_linenumber" name="L284" href="#L284">284</a> <em class="jxr_javadoccomment">     *  &lt;li&gt;For small means, uses simulation of a Poisson process</em>
 <a class="jxr_linenumber" name="L285" href="#L285">285</a> <em class="jxr_javadoccomment">     *   using Uniform deviates, as described</em>
-<a class="jxr_linenumber" name="L286" href="#L286">286</a> <em class="jxr_javadoccomment">     *   &lt;a href="<a href="http://irmi.epfl.ch/cmos/Pmmi/interactive/rng7.htm" target="alexandria_uri">http://irmi.epfl.ch/cmos/Pmmi/interactive/rng7.htm</a>"&gt; here&lt;/a&gt;.</em>
+<a class="jxr_linenumber" name="L286" href="#L286">286</a> <em class="jxr_javadoccomment">     *   &lt;a href="<a href="http://mathaa.epfl.ch/cours/PMMI2001/interactive/rng7.htm" target="alexandria_uri">http://mathaa.epfl.ch/cours/PMMI2001/interactive/rng7.htm</a>"&gt; here&lt;/a&gt;.</em>
 <a class="jxr_linenumber" name="L287" href="#L287">287</a> <em class="jxr_javadoccomment">     *   The Poisson process (and hence value returned) is bounded by 1000 * mean.</em>
 <a class="jxr_linenumber" name="L288" href="#L288">288</a> <em class="jxr_javadoccomment">     *  &lt;/li&gt;</em>
 <a class="jxr_linenumber" name="L289" href="#L289">289</a> <em class="jxr_javadoccomment">     *  &lt;li&gt;For large means, uses the rejection algorithm described in</em>
-<a class="jxr_linenumber" name="L290" href="#L290">290</a> <em class="jxr_javadoccomment">     *   &lt;quote&gt;</em>
-<a class="jxr_linenumber" name="L291" href="#L291">291</a> <em class="jxr_javadoccomment">     *    Devroye, Luc. (1981).&lt;i&gt;The Computer Generation of Poisson Random Variables&lt;/i&gt;</em>
-<a class="jxr_linenumber" name="L292" href="#L292">292</a> <em class="jxr_javadoccomment">     *    &lt;strong&gt;Computing&lt;/strong&gt; vol. 26 pp. 197-207.</em>
-<a class="jxr_linenumber" name="L293" href="#L293">293</a> <em class="jxr_javadoccomment">     *   &lt;/quote&gt;</em>
+<a class="jxr_linenumber" name="L290" href="#L290">290</a> <em class="jxr_javadoccomment">     *   &lt;blockquote&gt;</em>
+<a class="jxr_linenumber" name="L291" href="#L291">291</a> <em class="jxr_javadoccomment">     *    Devroye, Luc. (1981).&lt;i&gt;The Computer Generation of Poisson Random Variables&lt;/i&gt;&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L292" href="#L292">292</a> <em class="jxr_javadoccomment">     *    &lt;strong&gt;Computing&lt;/strong&gt; vol. 26 pp. 197-207.&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L293" href="#L293">293</a> <em class="jxr_javadoccomment">     *   &lt;/blockquote&gt;</em>
 <a class="jxr_linenumber" name="L294" href="#L294">294</a> <em class="jxr_javadoccomment">     *  &lt;/li&gt;</em>
 <a class="jxr_linenumber" name="L295" href="#L295">295</a> <em class="jxr_javadoccomment">     * &lt;/ul&gt;</em>
 <a class="jxr_linenumber" name="L296" href="#L296">296</a> <em class="jxr_javadoccomment">     * &lt;/p&gt;</em>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.distribution.fitting</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.distribution.fitting</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/fitting/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.distribution.fitting</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.distribution.fitting</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.distribution</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.distribution</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/distribution/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.distribution</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.distribution</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-frame.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.exception</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.exception</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-summary.html (original)
+++ websites/production/commons/content/proper/commons-math/xref/org/apache/commons/math3/exception/package-summary.html Fri Apr 17 19:17:13 2015
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Math 3.4.1 Reference Package org.apache.commons.math3.exception</title>
+		<title>Apache Commons Math 3.5 Reference Package org.apache.commons.math3.exception</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>



Mime
View raw message