backup/C-Plus-Plus
1
0
Fork 0
mirror of https://github.com/TheAlgorithms/C-Plus-Plus.git synced 2026-04-30 13:33:58 +08:00
Code Issues Projects Releases Wiki Activity

Documentation for 2dadbf73f4

Browse Source
This commit is contained in:
github-actions
2024-02-12 20:07:08 +00:00
parent b5b4d78981
commit a1eb9bd385
1261 changed files with 30404 additions and 29974 deletions
Download Patch File Download Diff File Expand all files Collapse all files

112
da/df2/durand__kerner__roots_8cpp.html
View File

@@ -3,16 +3,18 @@
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=11"/>
<meta name="generator" content="Doxygen 1.9.8"/>
<meta name="generator" content="Doxygen 1.10.0"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Algorithms_in_C++: numerical_methods/durand_kerner_roots.cpp File Reference</title>
<link href="../../tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="../../jquery.js"></script>
<script type="text/javascript" src="../../dynsections.js"></script>
<script type="text/javascript" src="../../clipboard.js"></script>
<link href="../../navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="../../resize.js"></script>
<script type="text/javascript" src="../../navtreedata.js"></script>
<script type="text/javascript" src="../../navtree.js"></script>
<script type="text/javascript" src="../../cookie.js"></script>
<link href="../../search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="../../search/searchdata.js"></script>
<script type="text/javascript" src="../../search/search.js"></script>
@@ -41,7 +43,7 @@ MathJax.Hub.Config({
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.8 -->
<!-- Generated by Doxygen 1.10.0 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
var searchBox = new SearchBox("searchBox", "../../search/",'.html');
@@ -53,7 +55,7 @@ var searchBox = new SearchBox("searchBox", "../../search/",'.html');
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(function() {
initMenu('../../',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
$(function() { init_search(); });
});
/* @license-end */
</script>
@@ -71,7 +73,7 @@ $(function() {
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(document).ready(function(){initNavTree('da/df2/durand__kerner__roots_8cpp.html','../../'); initResizable(); });
$(function(){initNavTree('da/df2/durand__kerner__roots_8cpp.html','../../'); initResizable(); });
/* @license-end */
</script>
<div id="doc-content">
@@ -123,7 +125,7 @@ Include dependency graph for durand_kerner_roots.cpp:</div>
</div><table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a id="define-members" name="define-members"></a>
Macros</h2></td></tr>
<tr class="memitem:af270a96662132d0385cb6b4637c5a689" id="r_af270a96662132d0385cb6b4637c5a689"><td class="memItemLeft" align="right" valign="top">#define&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a>&#160;&#160;&#160;1e-10</td></tr>
<tr class="memitem:af270a96662132d0385cb6b4637c5a689" id="r_af270a96662132d0385cb6b4637c5a689"><td class="memItemLeft" align="right" valign="top">#define&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#af270a96662132d0385cb6b4637c5a689">ACCURACY</a>&#160;&#160;&#160;1e-10</td></tr>
<tr class="separator:af270a96662132d0385cb6b4637c5a689"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a212979f5fc9dcf0a420e90b9b21d780a" id="r_a212979f5fc9dcf0a420e90b9b21d780a"><td class="memItemLeft" align="right" valign="top"><a id="a212979f5fc9dcf0a420e90b9b21d780a" name="a212979f5fc9dcf0a420e90b9b21d780a"></a>
#define&#160;</td><td class="memItemRight" valign="bottom"><b>MAX_BUFF_SIZE</b>&#160;&#160;&#160;50</td></tr>
@@ -131,19 +133,19 @@ Macros</h2></td></tr>
</table><table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a id="func-members" name="func-members"></a>
Functions</h2></td></tr>
<tr class="memitem:ae2ab49ba1969d7a5988f9bc2c8924569" id="r_ae2ab49ba1969d7a5988f9bc2c8924569"><td class="memItemLeft" align="right" valign="top"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;coeffs, <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; x)</td></tr>
<tr class="memitem:ae2ab49ba1969d7a5988f9bc2c8924569" id="r_ae2ab49ba1969d7a5988f9bc2c8924569"><td class="memItemLeft" align="right" valign="top"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;coeffs, <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; x)</td></tr>
<tr class="separator:ae2ab49ba1969d7a5988f9bc2c8924569"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a90219e35062007d1f1b68e9af071ab5c" id="r_a90219e35062007d1f1b68e9af071ab5c"><td class="memItemLeft" align="right" valign="top">const char *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &amp;x)</td></tr>
<tr class="memitem:a90219e35062007d1f1b68e9af071ab5c" id="r_a90219e35062007d1f1b68e9af071ab5c"><td class="memItemLeft" align="right" valign="top">const char *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &amp;x)</td></tr>
<tr class="separator:a90219e35062007d1f1b68e9af071ab5c"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a024b8bc4755863315456d573a6732377" id="r_a024b8bc4755863315456d573a6732377"><td class="memItemLeft" align="right" valign="top">bool&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377">check_termination</a> (long double delta)</td></tr>
<tr class="memitem:a024b8bc4755863315456d573a6732377" id="r_a024b8bc4755863315456d573a6732377"><td class="memItemLeft" align="right" valign="top">bool&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a024b8bc4755863315456d573a6732377">check_termination</a> (long double delta)</td></tr>
<tr class="separator:a024b8bc4755863315456d573a6732377"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a2c35b320ace8677f9b331faf94f8b2fd" id="r_a2c35b320ace8677f9b331faf94f8b2fd"><td class="memItemLeft" align="right" valign="top"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/utility/pair.html">std::pair</a>&lt; uint32_t, double &gt;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;coeffs, <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &gt; *roots, bool write_log=false)</td></tr>
<tr class="memitem:a2c35b320ace8677f9b331faf94f8b2fd" id="r_a2c35b320ace8677f9b331faf94f8b2fd"><td class="memItemLeft" align="right" valign="top"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/utility/pair.html">std::pair</a>&lt; uint32_t, double &gt;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a> (const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;coeffs, <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &gt; *roots, bool write_log=false)</td></tr>
<tr class="separator:a2c35b320ace8677f9b331faf94f8b2fd"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a1440a7779ac56f47a3f355ce4a8c7da0" id="r_a1440a7779ac56f47a3f355ce4a8c7da0"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a> ()</td></tr>
<tr class="memitem:a1440a7779ac56f47a3f355ce4a8c7da0" id="r_a1440a7779ac56f47a3f355ce4a8c7da0"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a> ()</td></tr>
<tr class="separator:a1440a7779ac56f47a3f355ce4a8c7da0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a0283886819c7c140a023582b7269e2d0" id="r_a0283886819c7c140a023582b7269e2d0"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a0283886819c7c140a023582b7269e2d0">test2</a> ()</td></tr>
<tr class="memitem:a0283886819c7c140a023582b7269e2d0" id="r_a0283886819c7c140a023582b7269e2d0"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a0283886819c7c140a023582b7269e2d0">test2</a> ()</td></tr>
<tr class="separator:a0283886819c7c140a023582b7269e2d0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a3c04138a5bfe5d72780bb7e82a18e627" id="r_a3c04138a5bfe5d72780bb7e82a18e627"><td class="memItemLeft" align="right" valign="top">int&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html#a3c04138a5bfe5d72780bb7e82a18e627">main</a> (int argc, char **argv)</td></tr>
<tr class="memitem:a3c04138a5bfe5d72780bb7e82a18e627" id="r_a3c04138a5bfe5d72780bb7e82a18e627"><td class="memItemLeft" align="right" valign="top">int&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="#a3c04138a5bfe5d72780bb7e82a18e627">main</a> (int argc, char **argv)</td></tr>
<tr class="separator:a3c04138a5bfe5d72780bb7e82a18e627"><td class="memSeparator" colspan="2">&#160;</td></tr>
</table>
<a name="details" id="details"></a><h2 class="groupheader">Detailed Description</h2>
@@ -188,8 +190,7 @@ approximations" style="pointer-events: none;" width="400" class="inline"/> </p>
<tr>
<td class="memname">bool check_termination </td>
<td>(</td>
<td class="paramtype">long double&#160;</td>
<td class="paramname"><em>delta</em></td><td>)</td>
<td class="paramtype">long double</td> <td class="paramname"><span class="paramname"><em>delta</em></span></td><td>)</td>
<td></td>
</tr>
</table>
@@ -205,12 +206,12 @@ approximations" style="pointer-events: none;" width="400" class="inline"/> </p>
<code>true</code> if termination reached </dd></dl>
<div class="fragment"><div class="line"><span class="lineno"> 91</span> {</div>
<div class="line"><span class="lineno"> 92</span> <span class="keyword">static</span> <span class="keywordtype">long</span> <span class="keywordtype">double</span> past_delta = INFINITY;</div>
<div class="line"><span class="lineno"> 93</span> <span class="keywordflow">if</span> (std::abs(past_delta - delta) &lt;= <a class="code hl_define" href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a> || delta &lt; <a class="code hl_define" href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a>)</div>
<div class="line"><span class="lineno"> 93</span> <span class="keywordflow">if</span> (std::abs(past_delta - delta) &lt;= <a class="code hl_define" href="#af270a96662132d0385cb6b4637c5a689">ACCURACY</a> || delta &lt; <a class="code hl_define" href="#af270a96662132d0385cb6b4637c5a689">ACCURACY</a>)</div>
<div class="line"><span class="lineno"> 94</span> <span class="keywordflow">return</span> <span class="keyword">true</span>;</div>
<div class="line"><span class="lineno"> 95</span> past_delta = delta;</div>
<div class="line"><span class="lineno"> 96</span> <span class="keywordflow">return</span> <span class="keyword">false</span>;</div>
<div class="line"><span class="lineno"> 97</span>}</div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_af270a96662132d0385cb6b4637c5a689"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a></div><div class="ttdeci">#define ACCURACY</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:45</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_af270a96662132d0385cb6b4637c5a689"><div class="ttname"><a href="#af270a96662132d0385cb6b4637c5a689">ACCURACY</a></div><div class="ttdeci">#define ACCURACY</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:45</div></div>
</div><!-- fragment -->
</div>
</div>
@@ -223,8 +224,7 @@ approximations" style="pointer-events: none;" width="400" class="inline"/> </p>
<tr>
<td class="memname">const char * complex_str </td>
<td>(</td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &amp;&#160;</td>
<td class="paramname"><em>x</em></td><td>)</td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &amp;</td> <td class="paramname"><span class="paramname"><em>x</em></span></td><td>)</td>
<td></td>
</tr>
</table>
@@ -264,25 +264,17 @@ Here is the call graph for this function:</div>
<tr>
<td class="memname"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/utility/pair.html">std::pair</a>&lt; uint32_t, double &gt; durand_kerner_algo </td>
<td>(</td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;&#160;</td>
<td class="paramname"><em>coeffs</em>, </td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;</td> <td class="paramname"><span class="paramname"><em>coeffs</em>, </span></td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &gt; *&#160;</td>
<td class="paramname"><em>roots</em>, </td>
<td class="paramtype"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; &gt; *</td> <td class="paramname"><span class="paramname"><em>roots</em>, </span></td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">bool&#160;</td>
<td class="paramname"><em>write_log</em> = <code>false</code>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
<td class="paramtype">bool</td> <td class="paramname"><span class="paramname"><em>write_log</em><span class="paramdefsep"> = </span><span class="paramdefval">false</span></span>&#160;)</td>
</tr>
</table>
</div><div class="memdoc">
@@ -318,11 +310,11 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 130</span> log_file &lt;&lt; <span class="stringliteral">&quot;avg. correction&quot;</span>;</div>
<div class="line"><span class="lineno"> 131</span> log_file &lt;&lt; <span class="stringliteral">&quot;\n0,&quot;</span>;</div>
<div class="line"><span class="lineno"> 132</span> <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><span class="lineno"> 133</span> log_file &lt;&lt; <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><span class="lineno"> 133</span> log_file &lt;&lt; <a class="code hl_function" href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><span class="lineno"> 134</span> }</div>
<div class="line"><span class="lineno"> 135</span> </div>
<div class="line"><span class="lineno"> 136</span> <span class="keywordtype">bool</span> break_loop = <span class="keyword">false</span>;</div>
<div class="line"><span class="lineno"> 137</span> <span class="keywordflow">while</span> (!<a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377">check_termination</a>(tol_condition) &amp;&amp; iter &lt; INT16_MAX &amp;&amp;</div>
<div class="line"><span class="lineno"> 137</span> <span class="keywordflow">while</span> (!<a class="code hl_function" href="#a024b8bc4755863315456d573a6732377">check_termination</a>(tol_condition) &amp;&amp; iter &lt; INT16_MAX &amp;&amp;</div>
<div class="line"><span class="lineno"> 138</span> !break_loop) {</div>
<div class="line"><span class="lineno"> 139</span> tol_condition = 0;</div>
<div class="line"><span class="lineno"> 140</span> iter++;</div>
@@ -339,7 +331,7 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 151</span> <span class="keywordflow">continue</span>;</div>
<div class="line"><span class="lineno"> 152</span> </div>
<div class="line"><span class="lineno"> 153</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a> numerator, denominator;</div>
<div class="line"><span class="lineno"> 154</span> numerator = <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a>(coeffs, (*roots)[n]);</div>
<div class="line"><span class="lineno"> 154</span> numerator = <a class="code hl_function" href="#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a>(coeffs, (*roots)[n]);</div>
<div class="line"><span class="lineno"> 155</span> denominator = 1.0;</div>
<div class="line"><span class="lineno"> 156</span> <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots-&gt;size(); i++)</div>
<div class="line"><span class="lineno"> 157</span> <span class="keywordflow">if</span> (i != n)</div>
@@ -368,14 +360,14 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 180</span> </div>
<div class="line"><span class="lineno"> 181</span> <span class="keywordflow">if</span> (log_file.<a class="code hl_functionRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">is_open</a>()) {</div>
<div class="line"><span class="lineno"> 182</span> <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><span class="lineno"> 183</span> log_file &lt;&lt; <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><span class="lineno"> 183</span> log_file &lt;&lt; <a class="code hl_function" href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><span class="lineno"> 184</span> }</div>
<div class="line"><span class="lineno"> 185</span> </div>
<div class="line"><span class="lineno"> 186</span><span class="preprocessor">#if defined(DEBUG) || !defined(NDEBUG)</span></div>
<div class="line"><span class="lineno"> 187</span> <span class="keywordflow">if</span> (iter % 500 == 0) {</div>
<div class="line"><span class="lineno"> 188</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;Iter: &quot;</span> &lt;&lt; iter &lt;&lt; <span class="stringliteral">&quot;\t&quot;</span>;</div>
<div class="line"><span class="lineno"> 189</span> <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><span class="lineno"> 190</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\t&quot;</span> &lt;&lt; <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]);</div>
<div class="line"><span class="lineno"> 190</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\t&quot;</span> &lt;&lt; <a class="code hl_function" href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>((*roots)[n]);</div>
<div class="line"><span class="lineno"> 191</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\t\tabsolute average change: &quot;</span> &lt;&lt; tol_condition</div>
<div class="line"><span class="lineno"> 192</span> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><span class="lineno"> 193</span> }</div>
@@ -390,9 +382,9 @@ Here is the call graph for this function:</div>
<div class="ttc" id="abasic_ofstream_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/basic_ofstream.html">std::ofstream</a></div></div>
<div class="ttc" id="abasic_ostream_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cerr</a></div></div>
<div class="ttc" id="acomplex_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a></div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a024b8bc4755863315456d573a6732377"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377">check_termination</a></div><div class="ttdeci">bool check_termination(long double delta)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:91</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a90219e35062007d1f1b68e9af071ab5c"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a></div><div class="ttdeci">const char * complex_str(const std::complex&lt; double &gt; &amp;x)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:76</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_ae2ab49ba1969d7a5988f9bc2c8924569"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a></div><div class="ttdeci">std::complex&lt; double &gt; poly_function(const std::valarray&lt; double &gt; &amp;coeffs, std::complex&lt; double &gt; x)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:53</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a024b8bc4755863315456d573a6732377"><div class="ttname"><a href="#a024b8bc4755863315456d573a6732377">check_termination</a></div><div class="ttdeci">bool check_termination(long double delta)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:91</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a90219e35062007d1f1b68e9af071ab5c"><div class="ttname"><a href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a></div><div class="ttdeci">const char * complex_str(const std::complex&lt; double &gt; &amp;x)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:76</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_ae2ab49ba1969d7a5988f9bc2c8924569"><div class="ttname"><a href="#ae2ab49ba1969d7a5988f9bc2c8924569">poly_function</a></div><div class="ttdeci">std::complex&lt; double &gt; poly_function(const std::valarray&lt; double &gt; &amp;coeffs, std::complex&lt; double &gt; x)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:53</div></div>
<div class="ttc" id="aexit_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/utility/program/exit.html">std::exit</a></div><div class="ttdeci">T exit(T... args)</div></div>
<div class="ttc" id="ais_open_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">std::ofstream::is_open</a></div><div class="ttdeci">T is_open(T... args)</div></div>
<div class="ttc" id="aisinf_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/math/isinf.html">std::isinf</a></div><div class="ttdeci">T isinf(T... args)</div></div>
@@ -418,19 +410,12 @@ Here is the call graph for this function:</div>
<tr>
<td class="memname">int main </td>
<td>(</td>
<td class="paramtype">int&#160;</td>
<td class="paramname"><em>argc</em>, </td>
<td class="paramtype">int</td> <td class="paramname"><span class="paramname"><em>argc</em>, </span></td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">char **&#160;</td>
<td class="paramname"><em>argv</em>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
<td class="paramtype">char **</td> <td class="paramname"><span class="paramname"><em>argv</em></span>&#160;)</td>
</tr>
</table>
</div><div class="memdoc">
@@ -439,8 +424,8 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 285</span> <a class="code hl_functionRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/random/srand.html">std::srand</a>(<a class="code hl_functionRef" target="_blank" href="http://en.cppreference.com/w/cpp/chrono/c/time.html">std::time</a>(<span class="keyword">nullptr</span>));</div>
<div class="line"><span class="lineno"> 286</span> </div>
<div class="line"><span class="lineno"> 287</span> <span class="keywordflow">if</span> (argc &lt; 2) {</div>
<div class="line"><span class="lineno"> 288</span> <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a>(); <span class="comment">// run tests when no input is provided</span></div>
<div class="line"><span class="lineno"> 289</span> <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a0283886819c7c140a023582b7269e2d0">test2</a>(); <span class="comment">// and skip tests when input polynomial is provided</span></div>
<div class="line"><span class="lineno"> 288</span> <a class="code hl_function" href="#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a>(); <span class="comment">// run tests when no input is provided</span></div>
<div class="line"><span class="lineno"> 289</span> <a class="code hl_function" href="#a0283886819c7c140a023582b7269e2d0">test2</a>(); <span class="comment">// and skip tests when input polynomial is provided</span></div>
<div class="line"><span class="lineno"> 290</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;Please pass the coefficients of the polynomial as &quot;</span></div>
<div class="line"><span class="lineno"> 291</span> <span class="stringliteral">&quot;commandline &quot;</span></div>
<div class="line"><span class="lineno"> 292</span> <span class="stringliteral">&quot;arguments.\n&quot;</span>;</div>
@@ -478,12 +463,12 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 324</span> }</div>
<div class="line"><span class="lineno"> 325</span> </div>
<div class="line"><span class="lineno"> 326</span> clock_t end_time, start_time = <a class="code hl_functionRef" target="_blank" href="http://en.cppreference.com/w/cpp/chrono/c/clock.html">clock</a>();</div>
<div class="line"><span class="lineno"> 327</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;s0, <span class="keyword">true</span>);</div>
<div class="line"><span class="lineno"> 327</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;s0, <span class="keyword">true</span>);</div>
<div class="line"><span class="lineno"> 328</span> end_time = <a class="code hl_functionRef" target="_blank" href="http://en.cppreference.com/w/cpp/chrono/c/clock.html">clock</a>();</div>
<div class="line"><span class="lineno"> 329</span> </div>
<div class="line"><span class="lineno"> 330</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\nIterations: &quot;</span> &lt;&lt; <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a>.first &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><span class="lineno"> 331</span> <span class="keywordflow">for</span> (n = 0; n &lt; degree - 1; n++)</div>
<div class="line"><span class="lineno"> 332</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\t&quot;</span> &lt;&lt; <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>(s0[n]) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><span class="lineno"> 332</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;\t&quot;</span> &lt;&lt; <a class="code hl_function" href="#a90219e35062007d1f1b68e9af071ab5c">complex_str</a>(s0[n]) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><span class="lineno"> 333</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;absolute average change: &quot;</span> &lt;&lt; <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a>.second &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><span class="lineno"> 334</span> <a class="code hl_classRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;Time taken: &quot;</span></div>
<div class="line"><span class="lineno"> 335</span> &lt;&lt; <span class="keyword">static_cast&lt;</span><span class="keywordtype">double</span><span class="keyword">&gt;</span>(end_time - start_time) / CLOCKS_PER_SEC</div>
@@ -492,9 +477,9 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 338</span> <span class="keywordflow">return</span> 0;</div>
<div class="line"><span class="lineno"> 339</span>}</div>
<div class="ttc" id="aclock_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/chrono/c/clock.html">std::clock</a></div><div class="ttdeci">T clock(T... args)</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a0283886819c7c140a023582b7269e2d0"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a0283886819c7c140a023582b7269e2d0">test2</a></div><div class="ttdeci">void test2()</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:242</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a1440a7779ac56f47a3f355ce4a8c7da0"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a></div><div class="ttdeci">void test1()</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:207</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a2c35b320ace8677f9b331faf94f8b2fd"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a></div><div class="ttdeci">std::pair&lt; uint32_t, double &gt; durand_kerner_algo(const std::valarray&lt; double &gt; &amp;coeffs, std::valarray&lt; std::complex&lt; double &gt; &gt; *roots, bool write_log=false)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:109</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a0283886819c7c140a023582b7269e2d0"><div class="ttname"><a href="#a0283886819c7c140a023582b7269e2d0">test2</a></div><div class="ttdeci">void test2()</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:242</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a1440a7779ac56f47a3f355ce4a8c7da0"><div class="ttname"><a href="#a1440a7779ac56f47a3f355ce4a8c7da0">test1</a></div><div class="ttdeci">void test1()</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:207</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a2c35b320ace8677f9b331faf94f8b2fd"><div class="ttname"><a href="#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a></div><div class="ttdeci">std::pair&lt; uint32_t, double &gt; durand_kerner_algo(const std::valarray&lt; double &gt; &amp;coeffs, std::valarray&lt; std::complex&lt; double &gt; &gt; *roots, bool write_log=false)</div><div class="ttdef"><b>Definition</b> durand_kerner_roots.cpp:109</div></div>
<div class="ttc" id="afibonacci__sum_8cpp_html_aadb40ac4c74a7efc0680b83eeee138aa"><div class="ttname"><a href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">math::fibonacci_sum::result</a></div><div class="ttdeci">uint64_t result(uint64_t n)</div><div class="ttdef"><b>Definition</b> fibonacci_sum.cpp:76</div></div>
<div class="ttc" id="arand_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/random/rand.html">std::rand</a></div><div class="ttdeci">T rand(T... args)</div></div>
<div class="ttc" id="asrand_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/random/srand.html">std::srand</a></div><div class="ttdeci">T srand(T... args)</div></div>
@@ -513,19 +498,12 @@ Here is the call graph for this function:</div>
<tr>
<td class="memname"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt; poly_function </td>
<td>(</td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;&#160;</td>
<td class="paramname"><em>coeffs</em>, </td>
<td class="paramtype">const <a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a>&lt; double &gt; &amp;</td> <td class="paramname"><span class="paramname"><em>coeffs</em>, </span></td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt;&#160;</td>
<td class="paramname"><em>x</em>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
<td class="paramtype"><a class="elRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a>&lt; double &gt;</td> <td class="paramname"><span class="paramname"><em>x</em></span>&#160;)</td>
</tr>
</table>
</div><div class="memdoc">
@@ -571,7 +549,7 @@ Here is the call graph for this function:</div>
<tr>
<td class="memname">void test1 </td>
<td>(</td>
<td class="paramname"></td><td>)</td>
<td class="paramname"><span class="paramname"></span></td><td>)</td>
<td></td>
</tr>
</table>
@@ -592,7 +570,7 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 219</span> roots[n] /= 25.f;</div>
<div class="line"><span class="lineno"> 220</span> }</div>
<div class="line"><span class="lineno"> 221</span> </div>
<div class="line"><span class="lineno"> 222</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><span class="lineno"> 222</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><span class="lineno"> 223</span> </div>
<div class="line"><span class="lineno"> 224</span> <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots.size(); i++) {</div>
<div class="line"><span class="lineno"> 225</span> <span class="comment">// check if approximations are have &lt; 0.1% error with one of the</span></div>
@@ -624,7 +602,7 @@ Here is the call graph for this function:</div>
<tr>
<td class="memname">void test2 </td>
<td>(</td>
<td class="paramname"></td><td>)</td>
<td class="paramname"><span class="paramname"></span></td><td>)</td>
<td></td>
</tr>
</table>
@@ -646,7 +624,7 @@ Here is the call graph for this function:</div>
<div class="line"><span class="lineno"> 255</span> roots[n] /= 25.f;</div>
<div class="line"><span class="lineno"> 256</span> }</div>
<div class="line"><span class="lineno"> 257</span> </div>
<div class="line"><span class="lineno"> 258</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="../../da/df2/durand__kerner__roots_8cpp.html#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><span class="lineno"> 258</span> <span class="keyword">auto</span> <a class="code hl_function" href="../../de/dc3/fibonacci__sum_8cpp.html#aadb40ac4c74a7efc0680b83eeee138aa">result</a> = <a class="code hl_function" href="#a2c35b320ace8677f9b331faf94f8b2fd">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><span class="lineno"> 259</span> </div>
<div class="line"><span class="lineno"> 260</span> <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots.size(); i++) {</div>
<div class="line"><span class="lineno"> 261</span> <span class="comment">// check if approximations are have &lt; 0.1% error with one of the</span></div>
@@ -675,7 +653,7 @@ Here is the call graph for this function:</div>
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
<ul>
<li class="navelem"><a class="el" href="../../dir_9c6faab82c22511b50177aa2e38e2780.html">numerical_methods</a></li><li class="navelem"><a class="el" href="../../da/df2/durand__kerner__roots_8cpp.html">durand_kerner_roots.cpp</a></li>
<li class="footer">Generated by <a href="https://www.doxygen.org/index.html"><img class="footer" src="../../doxygen.svg" width="104" height="31" alt="doxygen"/></a> 1.9.8 </li>
<li class="footer">Generated by <a href="https://www.doxygen.org/index.html"><img class="footer" src="../../doxygen.svg" width="104" height="31" alt="doxygen"/></a> 1.10.0 </li>
</ul>
</div>
</body>