Documentation for aaa08b0150

This commit is contained in:
github-actions
2020-06-19 16:07:44 +00:00
parent 4c261744f6
commit 93021161c2
1714 changed files with 114526 additions and 0 deletions

View File

@@ -0,0 +1,586 @@
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.18"/>
<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>
<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>
<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>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
extensions: ["tex2jax.js", "TeX/AMSmath.js", "TeX/AMSsymbols.js"],
jax: ["input/TeX","output/HTML-CSS"],
});
</script>
<script type="text/javascript" async="async" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.7/MathJax.js?config=TeX-MML-AM_CHTML/MathJax.js"></script>
<link href="../../doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectalign" style="padding-left: 0.5em;">
<div id="projectname">Algorithms_in_C++
&#160;<span id="projectnumber">1.0.0</span>
</div>
<div id="projectbrief">Set of algorithms implemented in C++.</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.18 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "../../search",false,'Search');
/* @license-end */
</script>
<script type="text/javascript" src="../../menudata.js"></script>
<script type="text/javascript" src="../../menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(function() {
initMenu('../../',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
/* @license-end */</script>
<div id="main-nav"></div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
<div id="nav-tree">
<div id="nav-tree-contents">
<div id="nav-sync" class="sync"></div>
</div>
</div>
<div id="splitbar" style="-moz-user-select:none;"
class="ui-resizable-handle">
</div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function(){initNavTree('da/df2/durand__kerner__roots_8cpp.html','../../'); initResizable(); });
/* @license-end */
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div class="header">
<div class="summary">
<a href="#define-members">Macros</a> &#124;
<a href="#func-members">Functions</a> </div>
<div class="headertitle">
<div class="title">durand_kerner_roots.cpp File Reference</div> </div>
</div><!--header-->
<div class="contents">
<p>Compute all possible approximate roots of any given polynomial using <a href="https://en.wikipedia.org/wiki/Durand%E2%80%93Kerner_method">Durand Kerner algorithm</a>
<a href="#details">More...</a></p>
<div class="textblock"><code>#include &lt;algorithm&gt;</code><br />
<code>#include &lt;cassert&gt;</code><br />
<code>#include &lt;cmath&gt;</code><br />
<code>#include &lt;complex&gt;</code><br />
<code>#include &lt;cstdlib&gt;</code><br />
<code>#include &lt;ctime&gt;</code><br />
<code>#include &lt;fstream&gt;</code><br />
<code>#include &lt;iostream&gt;</code><br />
<code>#include &lt;valarray&gt;</code><br />
</div><div class="textblock"><div class="dynheader">
Include dependency graph for durand_kerner_roots.cpp:</div>
<div class="dyncontent">
<div class="center"><div class="zoom"><iframe scrolling="no" frameborder="0" src="../../d4/d26/durand__kerner__roots_8cpp__incl.svg" width="100%" height="395"><p><b>This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.</b></p></iframe></div>
</div>
</div>
</div><table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a name="define-members"></a>
Macros</h2></td></tr>
<tr class="memitem: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="separator:af270a96662132d0385cb6b4637c5a689"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a212979f5fc9dcf0a420e90b9b21d780a"><td class="memItemLeft" align="right" valign="top"><a id="a212979f5fc9dcf0a420e90b9b21d780a"></a>
#define&#160;</td><td class="memItemRight" valign="bottom"><b>MAX_BUFF_SIZE</b>&#160;&#160;&#160;50</td></tr>
<tr class="separator:a212979f5fc9dcf0a420e90b9b21d780a"><td class="memSeparator" colspan="2">&#160;</td></tr>
</table><table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a name="func-members"></a>
Functions</h2></td></tr>
<tr class="memitem:a51637ff89e55ca749993d3e84b95d153"><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#a51637ff89e55ca749993d3e84b95d153">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:a51637ff89e55ca749993d3e84b95d153"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a3bada3e58ad8d1b00918db2ffe83a740"><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#a3bada3e58ad8d1b00918db2ffe83a740">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:a3bada3e58ad8d1b00918db2ffe83a740"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem: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="separator:a024b8bc4755863315456d573a6732377"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ac38b9cfd2a39d856e0feba3b5d241ded"><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#ac38b9cfd2a39d856e0feba3b5d241ded">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:ac38b9cfd2a39d856e0feba3b5d241ded"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem: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="separator:a1440a7779ac56f47a3f355ce4a8c7da0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem: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="separator:a0283886819c7c140a023582b7269e2d0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a3c04138a5bfe5d72780bb7e82a18e627"><td class="memItemLeft" align="right" valign="top"><a id="a3c04138a5bfe5d72780bb7e82a18e627"></a>
int&#160;</td><td class="memItemRight" valign="bottom"><b>main</b> (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>
<div class="textblock"><p>Compute all possible approximate roots of any given polynomial using <a href="https://en.wikipedia.org/wiki/Durand%E2%80%93Kerner_method">Durand Kerner algorithm</a> </p>
<dl class="section author"><dt>Author</dt><dd><a href="https://github.com/kvedala">Krishna Vedala</a></dd></dl>
<p>Test the algorithm online: <a href="https://gist.github.com/kvedala/27f1b0b6502af935f6917673ec43bcd7">https://gist.github.com/kvedala/27f1b0b6502af935f6917673ec43bcd7</a></p>
<p>Try the highly unstable Wilkinson's polynomial: </p><div class="fragment"><div class="line">./numerical_methods/durand_kerner_roots 1 -210 20615 -1256850 53327946</div>
<div class="line">-1672280820 40171771630 -756111184500 11310276995381 -135585182899530</div>
<div class="line">1307535010540395 -10142299865511450 63030812099294896 -311333643161390640</div>
<div class="line">1206647803780373360 -3599979517947607200 8037811822645051776</div>
<div class="line">-12870931245150988800 13803759753640704000 -8752948036761600000</div>
<div class="line">2432902008176640000</div>
</div><!-- fragment --><p>Sample implementation results to compute approximate roots of the equation \(x^4-1=0\):<br />
<img src="https://raw.githubusercontent.com/TheAlgorithms/C-Plus-Plus/docs/images/numerical_methods/durand_kerner_error.svg" alt="Error evolution during root approximations computed every
iteration." style="pointer-events: none;" width="400" class="inline"/> <img src="https://raw.githubusercontent.com/TheAlgorithms/C-Plus-Plus/docs/images/numerical_methods/durand_kerner_roots.svg" alt="Roots evolution - shows the initial approximation of the
roots and their convergence to a final approximation along with the iterative
approximations" style="pointer-events: none;" width="400" class="inline"/> </p>
</div><h2 class="groupheader">Macro Definition Documentation</h2>
<a id="af270a96662132d0385cb6b4637c5a689"></a>
<h2 class="memtitle"><span class="permalink"><a href="#af270a96662132d0385cb6b4637c5a689">&#9670;&nbsp;</a></span>ACCURACY</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname">#define ACCURACY&#160;&#160;&#160;1e-10</td>
</tr>
</table>
</div><div class="memdoc">
<p>maximum accuracy limit </p>
</div>
</div>
<h2 class="groupheader">Function Documentation</h2>
<a id="a024b8bc4755863315456d573a6732377"></a>
<h2 class="memtitle"><span class="permalink"><a href="#a024b8bc4755863315456d573a6732377">&#9670;&nbsp;</a></span>check_termination()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<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></td>
</tr>
</table>
</div><div class="memdoc">
<p>check for termination condition </p><dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">delta</td><td>point at which to evaluate the polynomial </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd><code>false</code> if termination not reached </dd>
<dd>
<code>true</code> if termination reached </dd></dl>
<div class="fragment"><div class="line"><a name="l00091"></a><span class="lineno"> 91</span>&#160; {</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>&#160; <span class="keyword">static</span> <span class="keywordtype">long</span> <span class="keywordtype">double</span> past_delta = INFINITY;</div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>&#160; <span class="keywordflow">if</span> (std::abs(past_delta - delta) &lt;= <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a> || delta &lt; <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689">ACCURACY</a>)</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>&#160; <span class="keywordflow">return</span> <span class="keyword">true</span>;</div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>&#160; past_delta = delta;</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>&#160; <span class="keywordflow">return</span> <span class="keyword">false</span>;</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>&#160;}</div>
</div><!-- fragment -->
</div>
</div>
<a id="a3bada3e58ad8d1b00918db2ffe83a740"></a>
<h2 class="memtitle"><span class="permalink"><a href="#a3bada3e58ad8d1b00918db2ffe83a740">&#9670;&nbsp;</a></span>complex_str()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<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></td>
</tr>
</table>
</div><div class="memdoc">
<p>create a textual form of complex number </p><dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">x</td><td>point at which to evaluate the polynomial </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>pointer to converted string </dd></dl>
<div class="fragment"><div class="line"><a name="l00076"></a><span class="lineno"> 76</span>&#160; {</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>&#160;<span class="preprocessor">#define MAX_BUFF_SIZE 50</span></div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>&#160; <span class="keyword">static</span> <span class="keywordtype">char</span> msg[MAX_BUFF_SIZE];</div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>&#160; </div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/c/fprintf.html">std::snprintf</a>(msg, MAX_BUFF_SIZE, <span class="stringliteral">&quot;% 7.04g%+7.04gj&quot;</span>, x.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex/real.html">real</a>(), x.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html">imag</a>());</div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>&#160; </div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>&#160; <span class="keywordflow">return</span> msg;</div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>&#160;}</div>
</div><!-- fragment -->
</div>
</div>
<a id="ac38b9cfd2a39d856e0feba3b5d241ded"></a>
<h2 class="memtitle"><span class="permalink"><a href="#ac38b9cfd2a39d856e0feba3b5d241ded">&#9670;&nbsp;</a></span>durand_kerner_algo()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<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>
</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>
</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>
</tr>
</table>
</div><div class="memdoc">
<p>Implements Durand Kerner iterative algorithm to compute all roots of a polynomial.</p>
<dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">coeffs</td><td>coefficients of the polynomial </td></tr>
<tr><td class="paramdir">[out]</td><td class="paramname">roots</td><td>the computed roots of the polynomial </td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">write_log</td><td>flag whether to save the log file (default = <code>false</code>) </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>pair of values - number of iterations taken and final accuracy achieved </dd></dl>
<div class="fragment"><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>&#160; {</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>&#160; <span class="keywordtype">long</span> <span class="keywordtype">double</span> tol_condition = 1;</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>&#160; uint32_t iter = 0;</div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>&#160; <span class="keywordtype">int</span> n;</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream.html">std::ofstream</a> log_file;</div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>&#160; </div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>&#160; <span class="keywordflow">if</span> (write_log) {</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>&#160; <span class="comment">/*</span></div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>&#160;<span class="comment"> * store intermediate values to a CSV file</span></div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>&#160; log_file.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/open.html">open</a>(<span class="stringliteral">&quot;durand_kerner.log.csv&quot;</span>);</div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>&#160; <span class="keywordflow">if</span> (!log_file.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">is_open</a>()) {</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/c/perror.html">perror</a>(<span class="stringliteral">&quot;Unable to create a storage log file!&quot;</span>);</div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/utility/program/exit.html">std::exit</a>(EXIT_FAILURE);</div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>&#160; }</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>&#160; log_file &lt;&lt; <span class="stringliteral">&quot;iter#,&quot;</span>;</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>&#160; </div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++) log_file &lt;&lt; <span class="stringliteral">&quot;root_&quot;</span> &lt;&lt; n &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>&#160; </div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>&#160; log_file &lt;&lt; <span class="stringliteral">&quot;avg. correction&quot;</span>;</div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>&#160; log_file &lt;&lt; <span class="stringliteral">&quot;\n0,&quot;</span>;</div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>&#160; log_file &lt;&lt; <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>&#160; }</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>&#160; </div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>&#160; <span class="keywordtype">bool</span> break_loop = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>&#160; <span class="keywordflow">while</span> (!<a class="code" 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"><a name="l00138"></a><span class="lineno"> 138</span>&#160; !break_loop) {</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>&#160; tol_condition = 0;</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>&#160; iter++;</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>&#160; break_loop = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>&#160; </div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>&#160; <span class="keywordflow">if</span> (log_file.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">is_open</a>())</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>&#160; log_file &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span> &lt;&lt; iter &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>&#160; </div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>&#160;<span class="preprocessor">#ifdef _OPENMP</span></div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>&#160;<span class="preprocessor">#pragma omp parallel for shared(break_loop, tol_condition)</span></div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++) {</div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>&#160; <span class="keywordflow">if</span> (break_loop)</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>&#160; <span class="keywordflow">continue</span>;</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>&#160; </div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a> numerator, denominator;</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>&#160; numerator = <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153">poly_function</a>(coeffs, (*roots)[n]);</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>&#160; denominator = 1.0;</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots-&gt;size(); i++)</div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>&#160; <span class="keywordflow">if</span> (i != n)</div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>&#160; denominator *= (*roots)[n] - (*roots)[i];</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>&#160; </div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;long double&gt;</a> delta = numerator / denominator;</div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>&#160; </div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>&#160; <span class="keywordflow">if</span> (<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/math/isnan.html">std::isnan</a>(std::abs(delta)) || <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/math/isinf.html">std::isinf</a>(std::abs(delta))) {</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cerr</a> &lt;&lt; <span class="stringliteral">&quot;\n\nOverflow/underrun error - got value = &quot;</span></div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>&#160; &lt;&lt; std::abs(delta) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>&#160; <span class="comment">// return std::pair&lt;uint32_t, double&gt;(iter, tol_condition);</span></div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>&#160; break_loop = <span class="keyword">true</span>;</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>&#160; }</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>&#160; </div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>&#160; (*roots)[n] -= delta;</div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>&#160; </div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>&#160;<span class="preprocessor">#ifdef _OPENMP</span></div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>&#160;<span class="preprocessor">#pragma omp critical</span></div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>&#160; tol_condition = <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/algorithm/max.html">std::max</a>(tol_condition, std::abs(std::abs(delta)));</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>&#160; }</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>&#160; <span class="comment">// tol_condition /= (degree - 1);</span></div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>&#160; </div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>&#160; <span class="keywordflow">if</span> (break_loop)</div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>&#160; <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>&#160; </div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>&#160; <span class="keywordflow">if</span> (log_file.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">is_open</a>()) {</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>&#160; log_file &lt;&lt; <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740">complex_str</a>((*roots)[n]) &lt;&lt; <span class="stringliteral">&quot;,&quot;</span>;</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>&#160; }</div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>&#160; </div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>&#160;<span class="preprocessor">#if defined(DEBUG) || !defined(NDEBUG)</span></div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>&#160; <span class="keywordflow">if</span> (iter % 500 == 0) {</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>&#160; <a class="codeRef" 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"><a name="l00189"></a><span class="lineno"> 189</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; roots-&gt;size(); n++)</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>&#160; <a class="codeRef" 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" href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740">complex_str</a>((*roots)[n]);</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>&#160; <a class="codeRef" 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"><a name="l00192"></a><span class="lineno"> 192</span>&#160; &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>&#160; }</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>&#160; </div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>&#160; <span class="keywordflow">if</span> (log_file.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html">is_open</a>())</div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>&#160; log_file &lt;&lt; tol_condition;</div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>&#160; }</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>&#160; </div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>&#160; <span class="keywordflow">return</span> <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/utility/pair.html">std::pair&lt;uint32_t, long double&gt;</a>(iter, tol_condition);</div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>&#160;}</div>
</div><!-- fragment --><div class="dynheader">
Here is the call graph for this function:</div>
<div class="dyncontent">
<div class="center"><iframe scrolling="no" frameborder="0" src="../../da/df2/durand__kerner__roots_8cpp_ac38b9cfd2a39d856e0feba3b5d241ded_cgraph.svg" width="527" height="494"><p><b>This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.</b></p></iframe>
</div>
</div>
</div>
</div>
<a id="a51637ff89e55ca749993d3e84b95d153"></a>
<h2 class="memtitle"><span class="permalink"><a href="#a51637ff89e55ca749993d3e84b95d153">&#9670;&nbsp;</a></span>poly_function()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<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>
</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>
</tr>
</table>
</div><div class="memdoc">
<p>Evaluate the value of a polynomial with given coefficients </p><dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">coeffs</td><td>coefficients of the polynomial </td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">x</td><td>point at which to evaluate the polynomial </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>\(f(x)\) </dd></dl>
<div class="fragment"><div class="line"><a name="l00054"></a><span class="lineno"> 54</span>&#160; {</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>&#160; <span class="keywordtype">double</span> real = 0.f, imag = 0.f;</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>&#160; <span class="keywordtype">int</span> n;</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>&#160; </div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>&#160; <span class="comment">// #ifdef _OPENMP</span></div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>&#160; <span class="comment">// #pragma omp target teams distribute reduction(+ : real, imag)</span></div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>&#160; <span class="comment">// #endif</span></div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>&#160; <span class="keywordflow">for</span> (n = 0; n &lt; coeffs.size(); n++) {</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a> tmp =</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>&#160; coeffs[n] * <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/math/pow.html">std::pow</a>(x, coeffs.size() - n - 1);</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>&#160; real += tmp.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex/real.html">real</a>();</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>&#160; imag += tmp.<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html">imag</a>();</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>&#160; }</div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>&#160; </div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>&#160; <span class="keywordflow">return</span> <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(real, imag);</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>&#160;}</div>
</div><!-- fragment --><div class="dynheader">
Here is the call graph for this function:</div>
<div class="dyncontent">
<div class="center"><iframe scrolling="no" frameborder="0" src="../../da/df2/durand__kerner__roots_8cpp_a51637ff89e55ca749993d3e84b95d153_cgraph.svg" width="291" height="139"><p><b>This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.</b></p></iframe>
</div>
</div>
</div>
</div>
<a id="a1440a7779ac56f47a3f355ce4a8c7da0"></a>
<h2 class="memtitle"><span class="permalink"><a href="#a1440a7779ac56f47a3f355ce4a8c7da0">&#9670;&nbsp;</a></span>test1()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname">void test1 </td>
<td>(</td>
<td class="paramname"></td><td>)</td>
<td></td>
</tr>
</table>
</div><div class="memdoc">
<p>Self test the algorithm by checking the roots for \(x^2+4=0\) to which the roots are \(0 \pm 2i\) </p>
<div class="fragment"><div class="line"><a name="l00207"></a><span class="lineno"> 207</span>&#160; {</div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>&#160; <span class="keyword">const</span> <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;double&gt;</a> coeffs = {1, 0, 4}; <span class="comment">// x^2 - 2 = 0</span></div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;std::complex&lt;double&gt;</a>&gt; roots(2);</div>
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;std::complex&lt;double&gt;</a>&gt; expected = {</div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(0., 2.),</div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(0., -2.) <span class="comment">// known expected roots</span></div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>&#160; };</div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>&#160; </div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>&#160; <span class="comment">/* initialize root approximations with random values */</span></div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> n = 0; n &lt; roots.size(); n++) {</div>
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>&#160; roots[n] = <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/random/rand.html">std::rand</a>() % 100, <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/random/rand.html">std::rand</a>() % 100);</div>
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>&#160; roots[n] -= 50.f;</div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span>&#160; roots[n] /= 25.f;</div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>&#160; }</div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>&#160; </div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span>&#160; <span class="keyword">auto</span> result = <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>&#160; </div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots.size(); i++) {</div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span>&#160; <span class="comment">// check if approximations are have &lt; 0.1% error with one of the</span></div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>&#160; <span class="comment">// expected roots</span></div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span>&#160; <span class="keywordtype">bool</span> err1 = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j = 0; j &lt; roots.size(); j++)</div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span>&#160; err1 |= std::abs(std::abs(roots[i] - expected[j])) &lt; 1e-3;</div>
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span>&#160; assert(err1);</div>
<div class="line"><a name="l00231"></a><span class="lineno"> 231</span>&#160; }</div>
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span>&#160; </div>
<div class="line"><a name="l00233"></a><span class="lineno"> 233</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;Test 1 passed! - &quot;</span> &lt;&lt; result.first &lt;&lt; <span class="stringliteral">&quot; iterations, &quot;</span></div>
<div class="line"><a name="l00234"></a><span class="lineno"> 234</span>&#160; &lt;&lt; result.second &lt;&lt; <span class="stringliteral">&quot; accuracy&quot;</span></div>
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span>&#160; &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span>&#160;}</div>
</div><!-- fragment --><div class="dynheader">
Here is the call graph for this function:</div>
<div class="dyncontent">
<div class="center"><iframe scrolling="no" frameborder="0" src="../../da/df2/durand__kerner__roots_8cpp_a1440a7779ac56f47a3f355ce4a8c7da0_cgraph.svg" width="626" height="494"><p><b>This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.</b></p></iframe>
</div>
</div>
</div>
</div>
<a id="a0283886819c7c140a023582b7269e2d0"></a>
<h2 class="memtitle"><span class="permalink"><a href="#a0283886819c7c140a023582b7269e2d0">&#9670;&nbsp;</a></span>test2()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname">void test2 </td>
<td>(</td>
<td class="paramname"></td><td>)</td>
<td></td>
</tr>
</table>
</div><div class="memdoc">
<p>Self test the algorithm by checking the roots for \(0.015625x^3-1=0\) to which the roots are \((4+0i),\,(-2\pm3.464i)\) </p>
<div class="fragment"><div class="line"><a name="l00242"></a><span class="lineno"> 242</span>&#160; {</div>
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span>&#160; <span class="keyword">const</span> <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;double&gt;</a> coeffs = {<span class="comment">// 0.015625 x^3 - 1 = 0</span></div>
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>&#160; 1. / 64., 0., 0., -1.};</div>
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;std::complex&lt;double&gt;</a>&gt; roots(3);</div>
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span>&#160; <span class="keyword">const</span> <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray&lt;std::complex&lt;double&gt;</a>&gt; expected = {</div>
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(4., 0.), <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(-2., 3.46410162),</div>
<div class="line"><a name="l00248"></a><span class="lineno"> 248</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(-2., -3.46410162) <span class="comment">// known expected roots</span></div>
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span>&#160; };</div>
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span>&#160; </div>
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span>&#160; <span class="comment">/* initialize root approximations with random values */</span></div>
<div class="line"><a name="l00252"></a><span class="lineno"> 252</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> n = 0; n &lt; roots.size(); n++) {</div>
<div class="line"><a name="l00253"></a><span class="lineno"> 253</span>&#160; roots[n] = <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex&lt;double&gt;</a>(<a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/random/rand.html">std::rand</a>() % 100, <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/numeric/random/rand.html">std::rand</a>() % 100);</div>
<div class="line"><a name="l00254"></a><span class="lineno"> 254</span>&#160; roots[n] -= 50.f;</div>
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span>&#160; roots[n] /= 25.f;</div>
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span>&#160; }</div>
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span>&#160; </div>
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span>&#160; <span class="keyword">auto</span> result = <a class="code" href="../../da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded">durand_kerner_algo</a>(coeffs, &amp;roots, <span class="keyword">false</span>);</div>
<div class="line"><a name="l00259"></a><span class="lineno"> 259</span>&#160; </div>
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; roots.size(); i++) {</div>
<div class="line"><a name="l00261"></a><span class="lineno"> 261</span>&#160; <span class="comment">// check if approximations are have &lt; 0.1% error with one of the</span></div>
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span>&#160; <span class="comment">// expected roots</span></div>
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span>&#160; <span class="keywordtype">bool</span> err1 = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span>&#160; <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j = 0; j &lt; roots.size(); j++)</div>
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span>&#160; err1 |= std::abs(std::abs(roots[i] - expected[j])) &lt; 1e-3;</div>
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span>&#160; assert(err1);</div>
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span>&#160; }</div>
<div class="line"><a name="l00268"></a><span class="lineno"> 268</span>&#160; </div>
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>&#160; <a class="codeRef" target="_blank" href="http://en.cppreference.com/w/cpp/io/basic_ostream.html">std::cout</a> &lt;&lt; <span class="stringliteral">&quot;Test 2 passed! - &quot;</span> &lt;&lt; result.first &lt;&lt; <span class="stringliteral">&quot; iterations, &quot;</span></div>
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span>&#160; &lt;&lt; result.second &lt;&lt; <span class="stringliteral">&quot; accuracy&quot;</span></div>
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span>&#160; &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div>
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span>&#160;}</div>
</div><!-- fragment --><div class="dynheader">
Here is the call graph for this function:</div>
<div class="dyncontent">
<div class="center"><iframe scrolling="no" frameborder="0" src="../../da/df2/durand__kerner__roots_8cpp_a0283886819c7c140a023582b7269e2d0_cgraph.svg" width="626" height="494"><p><b>This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.</b></p></iframe>
</div>
</div>
</div>
</div>
</div><!-- contents -->
</div><!-- doc-content -->
<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="apair_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/utility/pair.html">std::pair</a></div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a3bada3e58ad8d1b00918db2ffe83a740"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740">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_ac38b9cfd2a39d856e0feba3b5d241ded"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded">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="aisnan_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/math/isnan.html">std::isnan</a></div><div class="ttdeci">T isnan(T... args)</div></div>
<div class="ttc" id="afprintf_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/c/fprintf.html">std::snprintf</a></div><div class="ttdeci">T snprintf(T... args)</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="abasic_ofstream_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/basic_ofstream.html">std::ofstream</a></div><div class="ttdoc">STL class.</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>
<div class="ttc" id="aperror_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/c/perror.html">std::perror</a></div><div class="ttdeci">T perror(T... args)</div></div>
<div class="ttc" id="adurand__kerner__roots_8cpp_html_a51637ff89e55ca749993d3e84b95d153"><div class="ttname"><a href="../../da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153">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="avalarray_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/valarray.html">std::valarray</a></div><div class="ttdoc">STL class.</div></div>
<div class="ttc" id="areal_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/complex/real.html">std::complex::real</a></div><div class="ttdeci">T real(T... args)</div></div>
<div class="ttc" id="aopen_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/io/basic_ofstream/open.html">std::ofstream::open</a></div><div class="ttdeci">T open(T... args)</div></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="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="acomplex_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/complex.html">std::complex</a></div><div class="ttdoc">STL class.</div></div>
<div class="ttc" id="aimag_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html">std::complex::imag</a></div><div class="ttdeci">T imag(T... args)</div></div>
<div class="ttc" id="amax_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/algorithm/max.html">std::max</a></div><div class="ttdeci">T max(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="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="apow_html"><div class="ttname"><a href="http://en.cppreference.com/w/cpp/numeric/math/pow.html">std::pow</a></div><div class="ttdeci">T pow(T... args)</div></div>
<!-- start footer part -->
<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 on Fri Jun 19 2020 16:06:44 for Algorithms_in_C++ by
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="../../doxygen.png" alt="doxygen"/></a> 1.8.18 </li>
</ul>
</div>
</body>
</html>

View File

@@ -0,0 +1,12 @@
var durand__kerner__roots_8cpp =
[
[ "ACCURACY", "da/df2/durand__kerner__roots_8cpp.html#af270a96662132d0385cb6b4637c5a689", null ],
[ "MAX_BUFF_SIZE", "da/df2/durand__kerner__roots_8cpp.html#a212979f5fc9dcf0a420e90b9b21d780a", null ],
[ "check_termination", "da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377", null ],
[ "complex_str", "da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740", null ],
[ "durand_kerner_algo", "da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded", null ],
[ "main", "da/df2/durand__kerner__roots_8cpp.html#a3c04138a5bfe5d72780bb7e82a18e627", null ],
[ "poly_function", "da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153", null ],
[ "test1", "da/df2/durand__kerner__roots_8cpp.html#a1440a7779ac56f47a3f355ce4a8c7da0", null ],
[ "test2", "da/df2/durand__kerner__roots_8cpp.html#a0283886819c7c140a023582b7269e2d0", null ]
];

View File

@@ -0,0 +1,17 @@
<map id="test2" name="test2">
<area shape="rect" id="node1" title=" " alt="" coords="5,233,56,260"/>
<area shape="rect" id="node2" href="$da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded" title=" " alt="" coords="104,208,239,235"/>
<area shape="rect" id="node15" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/random/rand.html#" title=" " alt="" coords="135,259,208,285"/>
<area shape="rect" id="node3" href="$da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" title=" " alt="" coords="299,5,427,32"/>
<area shape="rect" id="node4" href="$da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" title=" " alt="" coords="316,56,409,83"/>
<area shape="rect" id="node5" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/utility/program/exit.html#" title=" " alt="" coords="328,107,397,133"/>
<area shape="rect" id="node6" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/is_open.html#" title=" " alt="" coords="287,157,439,184"/>
<area shape="rect" id="node7" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isinf.html#" title=" " alt="" coords="327,208,399,235"/>
<area shape="rect" id="node8" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isnan.html#" title=" " alt="" coords="323,259,402,285"/>
<area shape="rect" id="node9" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/algorithm/max.html#" title=" " alt="" coords="326,309,399,336"/>
<area shape="rect" id="node10" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/open.html#" title=" " alt="" coords="295,360,430,387"/>
<area shape="rect" id="node11" href="$da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" title=" " alt="" coords="313,411,412,437"/>
<area shape="rect" id="node12" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/imag.html#" title=" " alt="" coords="487,360,620,387"/>
<area shape="rect" id="node13" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/pow.html#" title=" " alt="" coords="517,411,589,437"/>
<area shape="rect" id="node14" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/real.html#" title=" " alt="" coords="490,461,617,488"/>
</map>

View File

@@ -0,0 +1 @@
d2f637fcf014076ba2cbd3f9b1b93fe4

View File

@@ -0,0 +1,232 @@
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
<!-- Generated by graphviz version 2.44.0 (20200408.0750)
-->
<!-- Title: test2 Pages: 1 -->
<svg width="469pt" height="370pt"
viewBox="0.00 0.00 469.00 370.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 366)">
<title>test2</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-366 465,-366 465,4 -4,4"/>
<!-- Node1 -->
<g id="node1" class="node">
<title>Node1</title>
<g id="a_node1"><a xlink:title=" ">
<polygon fill="#bfbfbf" stroke="black" points="0,-171.5 0,-190.5 38,-190.5 38,-171.5 0,-171.5"/>
<text text-anchor="middle" x="19" y="-178.5" font-family="Helvetica,sans-Serif" font-size="10.00">test2</text>
</a>
</g>
</g>
<!-- Node2 -->
<g id="node2" class="node">
<title>Node2</title>
<g id="a_node2"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="74,-190.5 74,-209.5 175,-209.5 175,-190.5 74,-190.5"/>
<text text-anchor="middle" x="124.5" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">durand_kerner_algo</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node2 -->
<g id="edge1" class="edge">
<title>Node1&#45;&gt;Node2</title>
<path fill="none" stroke="midnightblue" d="M38.03,-184.31C45.42,-185.67 54.44,-187.32 63.82,-189.04"/>
<polygon fill="midnightblue" stroke="midnightblue" points="63.2,-192.49 73.66,-190.85 64.46,-185.6 63.2,-192.49"/>
</g>
<!-- Node15 -->
<g id="node15" class="node">
<title>Node15</title>
<g id="a_node15"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/random/rand.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="97,-152.5 97,-171.5 152,-171.5 152,-152.5 97,-152.5"/>
<text text-anchor="middle" x="124.5" y="-159.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::rand</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node15 -->
<g id="edge14" class="edge">
<title>Node1&#45;&gt;Node15</title>
<path fill="none" stroke="midnightblue" d="M38.03,-177.69C51.46,-175.23 70.25,-171.77 86.8,-168.74"/>
<polygon fill="midnightblue" stroke="midnightblue" points="87.54,-172.16 96.75,-166.91 86.28,-165.27 87.54,-172.16"/>
</g>
<!-- Node3 -->
<g id="node3" class="node">
<title>Node3</title>
<g id="a_node3"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="220,-342.5 220,-361.5 316,-361.5 316,-342.5 220,-342.5"/>
<text text-anchor="middle" x="268" y="-349.5" font-family="Helvetica,sans-Serif" font-size="10.00">check_termination</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node3 -->
<g id="edge2" class="edge">
<title>Node2&#45;&gt;Node3</title>
<path fill="none" stroke="midnightblue" d="M129.05,-209.61C138.51,-233.91 166.66,-298.35 211,-333 213.31,-334.81 215.81,-336.45 218.42,-337.94"/>
<polygon fill="midnightblue" stroke="midnightblue" points="217.09,-341.19 227.61,-342.43 220.17,-334.9 217.09,-341.19"/>
</g>
<!-- Node4 -->
<g id="node4" class="node">
<title>Node4</title>
<g id="a_node4"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="233,-304.5 233,-323.5 303,-323.5 303,-304.5 233,-304.5"/>
<text text-anchor="middle" x="268" y="-311.5" font-family="Helvetica,sans-Serif" font-size="10.00">complex_str</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node4 -->
<g id="edge3" class="edge">
<title>Node2&#45;&gt;Node4</title>
<path fill="none" stroke="midnightblue" d="M131.6,-209.56C144.17,-228.54 174.87,-271.2 211,-295 214.82,-297.51 219,-299.73 223.3,-301.67"/>
<polygon fill="midnightblue" stroke="midnightblue" points="222.25,-305.02 232.83,-305.5 224.86,-298.52 222.25,-305.02"/>
</g>
<!-- Node5 -->
<g id="node5" class="node">
<title>Node5</title>
<g id="a_node5"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/utility/program/exit.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="242,-266.5 242,-285.5 294,-285.5 294,-266.5 242,-266.5"/>
<text text-anchor="middle" x="268" y="-273.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::exit</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node5 -->
<g id="edge4" class="edge">
<title>Node2&#45;&gt;Node5</title>
<path fill="none" stroke="midnightblue" d="M137.63,-209.58C153.66,-221.88 183.15,-243.29 211,-257 217.71,-260.31 225.14,-263.27 232.33,-265.82"/>
<polygon fill="midnightblue" stroke="midnightblue" points="231.31,-269.17 241.91,-269.03 233.54,-262.53 231.31,-269.17"/>
</g>
<!-- Node6 -->
<g id="node6" class="node">
<title>Node6</title>
<g id="a_node6"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="211,-228.5 211,-247.5 325,-247.5 325,-228.5 211,-228.5"/>
<text text-anchor="middle" x="268" y="-235.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::is_open</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node6 -->
<g id="edge5" class="edge">
<title>Node2&#45;&gt;Node6</title>
<path fill="none" stroke="midnightblue" d="M161.21,-209.59C179.47,-214.49 201.91,-220.52 221.41,-225.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="220.74,-229.2 231.3,-228.41 222.55,-222.44 220.74,-229.2"/>
</g>
<!-- Node7 -->
<g id="node7" class="node">
<title>Node7</title>
<g id="a_node7"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isinf.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="241,-190.5 241,-209.5 295,-209.5 295,-190.5 241,-190.5"/>
<text text-anchor="middle" x="268" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isinf</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node7 -->
<g id="edge6" class="edge">
<title>Node2&#45;&gt;Node7</title>
<path fill="none" stroke="midnightblue" d="M175.39,-200C193.71,-200 214.07,-200 230.94,-200"/>
<polygon fill="midnightblue" stroke="midnightblue" points="230.95,-203.5 240.94,-200 230.94,-196.5 230.95,-203.5"/>
</g>
<!-- Node8 -->
<g id="node8" class="node">
<title>Node8</title>
<g id="a_node8"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isnan.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="238.5,-152.5 238.5,-171.5 297.5,-171.5 297.5,-152.5 238.5,-152.5"/>
<text text-anchor="middle" x="268" y="-159.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isnan</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node8 -->
<g id="edge7" class="edge">
<title>Node2&#45;&gt;Node8</title>
<path fill="none" stroke="midnightblue" d="M161.21,-190.41C181.75,-184.89 207.57,-177.96 228.56,-172.32"/>
<polygon fill="midnightblue" stroke="midnightblue" points="229.68,-175.65 238.43,-169.67 227.86,-168.89 229.68,-175.65"/>
</g>
<!-- Node9 -->
<g id="node9" class="node">
<title>Node9</title>
<g id="a_node9"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/algorithm/max.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="240.5,-114.5 240.5,-133.5 295.5,-133.5 295.5,-114.5 240.5,-114.5"/>
<text text-anchor="middle" x="268" y="-121.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::max</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node9 -->
<g id="edge8" class="edge">
<title>Node2&#45;&gt;Node9</title>
<path fill="none" stroke="midnightblue" d="M156.66,-190.48C162.99,-187.89 169.41,-184.74 175,-181 194.34,-168.07 191.41,-155.55 211,-143 216.98,-139.17 223.82,-136.02 230.62,-133.46"/>
<polygon fill="midnightblue" stroke="midnightblue" points="232.09,-136.66 240.44,-130.14 229.85,-130.03 232.09,-136.66"/>
</g>
<!-- Node10 -->
<g id="node10" class="node">
<title>Node10</title>
<g id="a_node10"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="217.5,-76.5 217.5,-95.5 318.5,-95.5 318.5,-76.5 217.5,-76.5"/>
<text text-anchor="middle" x="268" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::open</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node10 -->
<g id="edge9" class="edge">
<title>Node2&#45;&gt;Node10</title>
<path fill="none" stroke="midnightblue" d="M160.48,-190.49C165.75,-187.99 170.81,-184.88 175,-181 202.41,-155.59 183.09,-129.86 211,-105 212.91,-103.3 214.98,-101.75 217.16,-100.33"/>
<polygon fill="midnightblue" stroke="midnightblue" points="218.89,-103.38 226.01,-95.54 215.55,-97.22 218.89,-103.38"/>
</g>
<!-- Node11 -->
<g id="node11" class="node">
<title>Node11</title>
<g id="a_node11"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="231,-38.5 231,-57.5 305,-57.5 305,-38.5 231,-38.5"/>
<text text-anchor="middle" x="268" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">poly_function</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node11 -->
<g id="edge10" class="edge">
<title>Node2&#45;&gt;Node11</title>
<path fill="none" stroke="midnightblue" d="M161.96,-190.46C166.8,-187.99 171.34,-184.89 175,-181 211.41,-142.31 173.83,-104.97 211,-67 214.16,-63.78 217.88,-61.1 221.88,-58.88"/>
<polygon fill="midnightblue" stroke="midnightblue" points="223.35,-62.06 230.97,-54.7 220.42,-55.69 223.35,-62.06"/>
</g>
<!-- Node12 -->
<g id="node12" class="node">
<title>Node12</title>
<g id="a_node12"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="361,-76.5 361,-95.5 461,-95.5 461,-76.5 361,-76.5"/>
<text text-anchor="middle" x="411" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::imag</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node12 -->
<g id="edge11" class="edge">
<title>Node11&#45;&gt;Node12</title>
<path fill="none" stroke="midnightblue" d="M304.58,-57.59C322.78,-62.49 345.14,-68.52 364.57,-73.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="363.86,-77.19 374.43,-76.41 365.69,-70.43 363.86,-77.19"/>
</g>
<!-- Node13 -->
<g id="node13" class="node">
<title>Node13</title>
<g id="a_node13"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/pow.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="384,-38.5 384,-57.5 438,-57.5 438,-38.5 384,-38.5"/>
<text text-anchor="middle" x="411" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::pow</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node13 -->
<g id="edge12" class="edge">
<title>Node11&#45;&gt;Node13</title>
<path fill="none" stroke="midnightblue" d="M305.31,-48C326.25,-48 352.5,-48 373.48,-48"/>
<polygon fill="midnightblue" stroke="midnightblue" points="373.61,-51.5 383.61,-48 373.61,-44.5 373.61,-51.5"/>
</g>
<!-- Node14 -->
<g id="node14" class="node">
<title>Node14</title>
<g id="a_node14"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/real.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="363.5,-0.5 363.5,-19.5 458.5,-19.5 458.5,-0.5 363.5,-0.5"/>
<text text-anchor="middle" x="411" y="-7.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::real</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node14 -->
<g id="edge13" class="edge">
<title>Node11&#45;&gt;Node14</title>
<path fill="none" stroke="midnightblue" d="M304.58,-38.41C322.78,-33.51 345.14,-27.48 364.57,-22.24"/>
<polygon fill="midnightblue" stroke="midnightblue" points="365.69,-25.57 374.43,-19.59 363.86,-18.81 365.69,-25.57"/>
</g>
</g>
</svg>

After

Width:  |  Height:  |  Size: 12 KiB

View File

@@ -0,0 +1,17 @@
<map id="test1" name="test1">
<area shape="rect" id="node1" title=" " alt="" coords="5,233,56,260"/>
<area shape="rect" id="node2" href="$da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded" title=" " alt="" coords="104,208,239,235"/>
<area shape="rect" id="node15" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/random/rand.html#" title=" " alt="" coords="135,259,208,285"/>
<area shape="rect" id="node3" href="$da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" title=" " alt="" coords="299,5,427,32"/>
<area shape="rect" id="node4" href="$da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" title=" " alt="" coords="316,56,409,83"/>
<area shape="rect" id="node5" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/utility/program/exit.html#" title=" " alt="" coords="328,107,397,133"/>
<area shape="rect" id="node6" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/is_open.html#" title=" " alt="" coords="287,157,439,184"/>
<area shape="rect" id="node7" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isinf.html#" title=" " alt="" coords="327,208,399,235"/>
<area shape="rect" id="node8" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isnan.html#" title=" " alt="" coords="323,259,402,285"/>
<area shape="rect" id="node9" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/algorithm/max.html#" title=" " alt="" coords="326,309,399,336"/>
<area shape="rect" id="node10" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/open.html#" title=" " alt="" coords="295,360,430,387"/>
<area shape="rect" id="node11" href="$da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" title=" " alt="" coords="313,411,412,437"/>
<area shape="rect" id="node12" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/imag.html#" title=" " alt="" coords="487,360,620,387"/>
<area shape="rect" id="node13" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/pow.html#" title=" " alt="" coords="517,411,589,437"/>
<area shape="rect" id="node14" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/real.html#" title=" " alt="" coords="490,461,617,488"/>
</map>

View File

@@ -0,0 +1 @@
a9cd90795b692401b9335f1a73a8038e

View File

@@ -0,0 +1,232 @@
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
<!-- Generated by graphviz version 2.44.0 (20200408.0750)
-->
<!-- Title: test1 Pages: 1 -->
<svg width="469pt" height="370pt"
viewBox="0.00 0.00 469.00 370.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 366)">
<title>test1</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-366 465,-366 465,4 -4,4"/>
<!-- Node1 -->
<g id="node1" class="node">
<title>Node1</title>
<g id="a_node1"><a xlink:title=" ">
<polygon fill="#bfbfbf" stroke="black" points="0,-171.5 0,-190.5 38,-190.5 38,-171.5 0,-171.5"/>
<text text-anchor="middle" x="19" y="-178.5" font-family="Helvetica,sans-Serif" font-size="10.00">test1</text>
</a>
</g>
</g>
<!-- Node2 -->
<g id="node2" class="node">
<title>Node2</title>
<g id="a_node2"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#ac38b9cfd2a39d856e0feba3b5d241ded" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="74,-190.5 74,-209.5 175,-209.5 175,-190.5 74,-190.5"/>
<text text-anchor="middle" x="124.5" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">durand_kerner_algo</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node2 -->
<g id="edge1" class="edge">
<title>Node1&#45;&gt;Node2</title>
<path fill="none" stroke="midnightblue" d="M38.03,-184.31C45.42,-185.67 54.44,-187.32 63.82,-189.04"/>
<polygon fill="midnightblue" stroke="midnightblue" points="63.2,-192.49 73.66,-190.85 64.46,-185.6 63.2,-192.49"/>
</g>
<!-- Node15 -->
<g id="node15" class="node">
<title>Node15</title>
<g id="a_node15"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/random/rand.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="97,-152.5 97,-171.5 152,-171.5 152,-152.5 97,-152.5"/>
<text text-anchor="middle" x="124.5" y="-159.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::rand</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node15 -->
<g id="edge14" class="edge">
<title>Node1&#45;&gt;Node15</title>
<path fill="none" stroke="midnightblue" d="M38.03,-177.69C51.46,-175.23 70.25,-171.77 86.8,-168.74"/>
<polygon fill="midnightblue" stroke="midnightblue" points="87.54,-172.16 96.75,-166.91 86.28,-165.27 87.54,-172.16"/>
</g>
<!-- Node3 -->
<g id="node3" class="node">
<title>Node3</title>
<g id="a_node3"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="220,-342.5 220,-361.5 316,-361.5 316,-342.5 220,-342.5"/>
<text text-anchor="middle" x="268" y="-349.5" font-family="Helvetica,sans-Serif" font-size="10.00">check_termination</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node3 -->
<g id="edge2" class="edge">
<title>Node2&#45;&gt;Node3</title>
<path fill="none" stroke="midnightblue" d="M129.05,-209.61C138.51,-233.91 166.66,-298.35 211,-333 213.31,-334.81 215.81,-336.45 218.42,-337.94"/>
<polygon fill="midnightblue" stroke="midnightblue" points="217.09,-341.19 227.61,-342.43 220.17,-334.9 217.09,-341.19"/>
</g>
<!-- Node4 -->
<g id="node4" class="node">
<title>Node4</title>
<g id="a_node4"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="233,-304.5 233,-323.5 303,-323.5 303,-304.5 233,-304.5"/>
<text text-anchor="middle" x="268" y="-311.5" font-family="Helvetica,sans-Serif" font-size="10.00">complex_str</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node4 -->
<g id="edge3" class="edge">
<title>Node2&#45;&gt;Node4</title>
<path fill="none" stroke="midnightblue" d="M131.6,-209.56C144.17,-228.54 174.87,-271.2 211,-295 214.82,-297.51 219,-299.73 223.3,-301.67"/>
<polygon fill="midnightblue" stroke="midnightblue" points="222.25,-305.02 232.83,-305.5 224.86,-298.52 222.25,-305.02"/>
</g>
<!-- Node5 -->
<g id="node5" class="node">
<title>Node5</title>
<g id="a_node5"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/utility/program/exit.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="242,-266.5 242,-285.5 294,-285.5 294,-266.5 242,-266.5"/>
<text text-anchor="middle" x="268" y="-273.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::exit</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node5 -->
<g id="edge4" class="edge">
<title>Node2&#45;&gt;Node5</title>
<path fill="none" stroke="midnightblue" d="M137.63,-209.58C153.66,-221.88 183.15,-243.29 211,-257 217.71,-260.31 225.14,-263.27 232.33,-265.82"/>
<polygon fill="midnightblue" stroke="midnightblue" points="231.31,-269.17 241.91,-269.03 233.54,-262.53 231.31,-269.17"/>
</g>
<!-- Node6 -->
<g id="node6" class="node">
<title>Node6</title>
<g id="a_node6"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="211,-228.5 211,-247.5 325,-247.5 325,-228.5 211,-228.5"/>
<text text-anchor="middle" x="268" y="-235.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::is_open</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node6 -->
<g id="edge5" class="edge">
<title>Node2&#45;&gt;Node6</title>
<path fill="none" stroke="midnightblue" d="M161.21,-209.59C179.47,-214.49 201.91,-220.52 221.41,-225.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="220.74,-229.2 231.3,-228.41 222.55,-222.44 220.74,-229.2"/>
</g>
<!-- Node7 -->
<g id="node7" class="node">
<title>Node7</title>
<g id="a_node7"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isinf.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="241,-190.5 241,-209.5 295,-209.5 295,-190.5 241,-190.5"/>
<text text-anchor="middle" x="268" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isinf</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node7 -->
<g id="edge6" class="edge">
<title>Node2&#45;&gt;Node7</title>
<path fill="none" stroke="midnightblue" d="M175.39,-200C193.71,-200 214.07,-200 230.94,-200"/>
<polygon fill="midnightblue" stroke="midnightblue" points="230.95,-203.5 240.94,-200 230.94,-196.5 230.95,-203.5"/>
</g>
<!-- Node8 -->
<g id="node8" class="node">
<title>Node8</title>
<g id="a_node8"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isnan.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="238.5,-152.5 238.5,-171.5 297.5,-171.5 297.5,-152.5 238.5,-152.5"/>
<text text-anchor="middle" x="268" y="-159.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isnan</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node8 -->
<g id="edge7" class="edge">
<title>Node2&#45;&gt;Node8</title>
<path fill="none" stroke="midnightblue" d="M161.21,-190.41C181.75,-184.89 207.57,-177.96 228.56,-172.32"/>
<polygon fill="midnightblue" stroke="midnightblue" points="229.68,-175.65 238.43,-169.67 227.86,-168.89 229.68,-175.65"/>
</g>
<!-- Node9 -->
<g id="node9" class="node">
<title>Node9</title>
<g id="a_node9"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/algorithm/max.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="240.5,-114.5 240.5,-133.5 295.5,-133.5 295.5,-114.5 240.5,-114.5"/>
<text text-anchor="middle" x="268" y="-121.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::max</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node9 -->
<g id="edge8" class="edge">
<title>Node2&#45;&gt;Node9</title>
<path fill="none" stroke="midnightblue" d="M156.66,-190.48C162.99,-187.89 169.41,-184.74 175,-181 194.34,-168.07 191.41,-155.55 211,-143 216.98,-139.17 223.82,-136.02 230.62,-133.46"/>
<polygon fill="midnightblue" stroke="midnightblue" points="232.09,-136.66 240.44,-130.14 229.85,-130.03 232.09,-136.66"/>
</g>
<!-- Node10 -->
<g id="node10" class="node">
<title>Node10</title>
<g id="a_node10"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="217.5,-76.5 217.5,-95.5 318.5,-95.5 318.5,-76.5 217.5,-76.5"/>
<text text-anchor="middle" x="268" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::open</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node10 -->
<g id="edge9" class="edge">
<title>Node2&#45;&gt;Node10</title>
<path fill="none" stroke="midnightblue" d="M160.48,-190.49C165.75,-187.99 170.81,-184.88 175,-181 202.41,-155.59 183.09,-129.86 211,-105 212.91,-103.3 214.98,-101.75 217.16,-100.33"/>
<polygon fill="midnightblue" stroke="midnightblue" points="218.89,-103.38 226.01,-95.54 215.55,-97.22 218.89,-103.38"/>
</g>
<!-- Node11 -->
<g id="node11" class="node">
<title>Node11</title>
<g id="a_node11"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="231,-38.5 231,-57.5 305,-57.5 305,-38.5 231,-38.5"/>
<text text-anchor="middle" x="268" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">poly_function</text>
</a>
</g>
</g>
<!-- Node2&#45;&gt;Node11 -->
<g id="edge10" class="edge">
<title>Node2&#45;&gt;Node11</title>
<path fill="none" stroke="midnightblue" d="M161.96,-190.46C166.8,-187.99 171.34,-184.89 175,-181 211.41,-142.31 173.83,-104.97 211,-67 214.16,-63.78 217.88,-61.1 221.88,-58.88"/>
<polygon fill="midnightblue" stroke="midnightblue" points="223.35,-62.06 230.97,-54.7 220.42,-55.69 223.35,-62.06"/>
</g>
<!-- Node12 -->
<g id="node12" class="node">
<title>Node12</title>
<g id="a_node12"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="361,-76.5 361,-95.5 461,-95.5 461,-76.5 361,-76.5"/>
<text text-anchor="middle" x="411" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::imag</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node12 -->
<g id="edge11" class="edge">
<title>Node11&#45;&gt;Node12</title>
<path fill="none" stroke="midnightblue" d="M304.58,-57.59C322.78,-62.49 345.14,-68.52 364.57,-73.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="363.86,-77.19 374.43,-76.41 365.69,-70.43 363.86,-77.19"/>
</g>
<!-- Node13 -->
<g id="node13" class="node">
<title>Node13</title>
<g id="a_node13"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/pow.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="384,-38.5 384,-57.5 438,-57.5 438,-38.5 384,-38.5"/>
<text text-anchor="middle" x="411" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::pow</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node13 -->
<g id="edge12" class="edge">
<title>Node11&#45;&gt;Node13</title>
<path fill="none" stroke="midnightblue" d="M305.31,-48C326.25,-48 352.5,-48 373.48,-48"/>
<polygon fill="midnightblue" stroke="midnightblue" points="373.61,-51.5 383.61,-48 373.61,-44.5 373.61,-51.5"/>
</g>
<!-- Node14 -->
<g id="node14" class="node">
<title>Node14</title>
<g id="a_node14"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/real.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="363.5,-0.5 363.5,-19.5 458.5,-19.5 458.5,-0.5 363.5,-0.5"/>
<text text-anchor="middle" x="411" y="-7.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::real</text>
</a>
</g>
</g>
<!-- Node11&#45;&gt;Node14 -->
<g id="edge13" class="edge">
<title>Node11&#45;&gt;Node14</title>
<path fill="none" stroke="midnightblue" d="M304.58,-38.41C322.78,-33.51 345.14,-27.48 364.57,-22.24"/>
<polygon fill="midnightblue" stroke="midnightblue" points="365.69,-25.57 374.43,-19.59 363.86,-18.81 365.69,-25.57"/>
</g>
</g>
</svg>

After

Width:  |  Height:  |  Size: 12 KiB

View File

@@ -0,0 +1,6 @@
<map id="poly_function" name="poly_function">
<area shape="rect" id="node1" title=" " alt="" coords="5,56,104,83"/>
<area shape="rect" id="node2" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/imag.html#" title=" " alt="" coords="152,5,285,32"/>
<area shape="rect" id="node3" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/pow.html#" title=" " alt="" coords="183,56,255,83"/>
<area shape="rect" id="node4" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/real.html#" title=" " alt="" coords="155,107,282,133"/>
</map>

View File

@@ -0,0 +1 @@
2c145b5ecf01642704ae315c2884e089

View File

@@ -0,0 +1,67 @@
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
<!-- Generated by graphviz version 2.44.0 (20200408.0750)
-->
<!-- Title: poly_function Pages: 1 -->
<svg width="218pt" height="104pt"
viewBox="0.00 0.00 218.00 104.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 100)">
<title>poly_function</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-100 214,-100 214,4 -4,4"/>
<!-- Node1 -->
<g id="node1" class="node">
<title>Node1</title>
<g id="a_node1"><a xlink:title=" ">
<polygon fill="#bfbfbf" stroke="black" points="0,-38.5 0,-57.5 74,-57.5 74,-38.5 0,-38.5"/>
<text text-anchor="middle" x="37" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">poly_function</text>
</a>
</g>
</g>
<!-- Node2 -->
<g id="node2" class="node">
<title>Node2</title>
<g id="a_node2"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="110,-76.5 110,-95.5 210,-95.5 210,-76.5 110,-76.5"/>
<text text-anchor="middle" x="160" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::imag</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node2 -->
<g id="edge1" class="edge">
<title>Node1&#45;&gt;Node2</title>
<path fill="none" stroke="midnightblue" d="M68.53,-57.59C83.72,-62.36 102.29,-68.19 118.63,-73.32"/>
<polygon fill="midnightblue" stroke="midnightblue" points="117.89,-76.76 128.48,-76.41 119.98,-70.08 117.89,-76.76"/>
</g>
<!-- Node3 -->
<g id="node3" class="node">
<title>Node3</title>
<g id="a_node3"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/pow.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="133,-38.5 133,-57.5 187,-57.5 187,-38.5 133,-38.5"/>
<text text-anchor="middle" x="160" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::pow</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node3 -->
<g id="edge2" class="edge">
<title>Node1&#45;&gt;Node3</title>
<path fill="none" stroke="midnightblue" d="M74.3,-48C89.58,-48 107.31,-48 122.59,-48"/>
<polygon fill="midnightblue" stroke="midnightblue" points="122.72,-51.5 132.72,-48 122.72,-44.5 122.72,-51.5"/>
</g>
<!-- Node4 -->
<g id="node4" class="node">
<title>Node4</title>
<g id="a_node4"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/real.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="112.5,-0.5 112.5,-19.5 207.5,-19.5 207.5,-0.5 112.5,-0.5"/>
<text text-anchor="middle" x="160" y="-7.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::real</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node4 -->
<g id="edge3" class="edge">
<title>Node1&#45;&gt;Node4</title>
<path fill="none" stroke="midnightblue" d="M68.53,-38.41C83.72,-33.64 102.29,-27.81 118.63,-22.68"/>
<polygon fill="midnightblue" stroke="midnightblue" points="119.98,-25.92 128.48,-19.59 117.89,-19.24 119.98,-25.92"/>
</g>
</g>
</svg>

After

Width:  |  Height:  |  Size: 3.1 KiB

View File

@@ -0,0 +1,15 @@
<map id="durand_kerner_algo" name="durand_kerner_algo">
<area shape="rect" id="node1" title=" " alt="" coords="5,208,140,235"/>
<area shape="rect" id="node2" href="$da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" title=" " alt="" coords="200,5,328,32"/>
<area shape="rect" id="node3" href="$da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" title=" " alt="" coords="217,56,311,83"/>
<area shape="rect" id="node4" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/utility/program/exit.html#" title=" " alt="" coords="229,107,299,133"/>
<area shape="rect" id="node5" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/is_open.html#" title=" " alt="" coords="188,157,340,184"/>
<area shape="rect" id="node6" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isinf.html#" title=" " alt="" coords="228,208,300,235"/>
<area shape="rect" id="node7" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/isnan.html#" title=" " alt="" coords="225,259,303,285"/>
<area shape="rect" id="node8" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/algorithm/max.html#" title=" " alt="" coords="227,309,301,336"/>
<area shape="rect" id="node9" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/io/basic_ofstream/open.html#" title=" " alt="" coords="197,360,331,387"/>
<area shape="rect" id="node10" href="$da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" title=" " alt="" coords="215,411,313,437"/>
<area shape="rect" id="node11" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/imag.html#" title=" " alt="" coords="388,360,521,387"/>
<area shape="rect" id="node12" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/math/pow.html#" title=" " alt="" coords="419,411,491,437"/>
<area shape="rect" id="node13" href="/Users/runner/runners/2.263.0/work/C-Plus-Plus/C-Plus-Plus/doc/cppreference-doxygen-web.tag.xml$cpp/numeric/complex/real.html#" title=" " alt="" coords="391,461,518,488"/>
</map>

View File

@@ -0,0 +1 @@
d975e88b18f81b99129d130f661e7bd8

View File

@@ -0,0 +1,202 @@
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
<!-- Generated by graphviz version 2.44.0 (20200408.0750)
-->
<!-- Title: durand_kerner_algo Pages: 1 -->
<svg width="395pt" height="370pt"
viewBox="0.00 0.00 395.00 370.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 366)">
<title>durand_kerner_algo</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-366 391,-366 391,4 -4,4"/>
<!-- Node1 -->
<g id="node1" class="node">
<title>Node1</title>
<g id="a_node1"><a xlink:title=" ">
<polygon fill="#bfbfbf" stroke="black" points="0,-190.5 0,-209.5 101,-209.5 101,-190.5 0,-190.5"/>
<text text-anchor="middle" x="50.5" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">durand_kerner_algo</text>
</a>
</g>
</g>
<!-- Node2 -->
<g id="node2" class="node">
<title>Node2</title>
<g id="a_node2"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a024b8bc4755863315456d573a6732377" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="146,-342.5 146,-361.5 242,-361.5 242,-342.5 146,-342.5"/>
<text text-anchor="middle" x="194" y="-349.5" font-family="Helvetica,sans-Serif" font-size="10.00">check_termination</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node2 -->
<g id="edge1" class="edge">
<title>Node1&#45;&gt;Node2</title>
<path fill="none" stroke="midnightblue" d="M55.05,-209.61C64.51,-233.91 92.66,-298.35 137,-333 139.31,-334.81 141.81,-336.45 144.42,-337.94"/>
<polygon fill="midnightblue" stroke="midnightblue" points="143.09,-341.19 153.61,-342.43 146.17,-334.9 143.09,-341.19"/>
</g>
<!-- Node3 -->
<g id="node3" class="node">
<title>Node3</title>
<g id="a_node3"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a3bada3e58ad8d1b00918db2ffe83a740" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="159,-304.5 159,-323.5 229,-323.5 229,-304.5 159,-304.5"/>
<text text-anchor="middle" x="194" y="-311.5" font-family="Helvetica,sans-Serif" font-size="10.00">complex_str</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node3 -->
<g id="edge2" class="edge">
<title>Node1&#45;&gt;Node3</title>
<path fill="none" stroke="midnightblue" d="M57.6,-209.56C70.17,-228.54 100.87,-271.2 137,-295 140.82,-297.51 145,-299.73 149.3,-301.67"/>
<polygon fill="midnightblue" stroke="midnightblue" points="148.25,-305.02 158.83,-305.5 150.86,-298.52 148.25,-305.02"/>
</g>
<!-- Node4 -->
<g id="node4" class="node">
<title>Node4</title>
<g id="a_node4"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/utility/program/exit.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="168,-266.5 168,-285.5 220,-285.5 220,-266.5 168,-266.5"/>
<text text-anchor="middle" x="194" y="-273.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::exit</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node4 -->
<g id="edge3" class="edge">
<title>Node1&#45;&gt;Node4</title>
<path fill="none" stroke="midnightblue" d="M63.63,-209.58C79.66,-221.88 109.15,-243.29 137,-257 143.71,-260.31 151.14,-263.27 158.33,-265.82"/>
<polygon fill="midnightblue" stroke="midnightblue" points="157.31,-269.17 167.91,-269.03 159.54,-262.53 157.31,-269.17"/>
</g>
<!-- Node5 -->
<g id="node5" class="node">
<title>Node5</title>
<g id="a_node5"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/is_open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="137,-228.5 137,-247.5 251,-247.5 251,-228.5 137,-228.5"/>
<text text-anchor="middle" x="194" y="-235.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::is_open</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node5 -->
<g id="edge4" class="edge">
<title>Node1&#45;&gt;Node5</title>
<path fill="none" stroke="midnightblue" d="M87.21,-209.59C105.47,-214.49 127.91,-220.52 147.41,-225.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="146.74,-229.2 157.3,-228.41 148.55,-222.44 146.74,-229.2"/>
</g>
<!-- Node6 -->
<g id="node6" class="node">
<title>Node6</title>
<g id="a_node6"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isinf.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="167,-190.5 167,-209.5 221,-209.5 221,-190.5 167,-190.5"/>
<text text-anchor="middle" x="194" y="-197.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isinf</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node6 -->
<g id="edge5" class="edge">
<title>Node1&#45;&gt;Node6</title>
<path fill="none" stroke="midnightblue" d="M101.39,-200C119.71,-200 140.07,-200 156.94,-200"/>
<polygon fill="midnightblue" stroke="midnightblue" points="156.95,-203.5 166.94,-200 156.94,-196.5 156.95,-203.5"/>
</g>
<!-- Node7 -->
<g id="node7" class="node">
<title>Node7</title>
<g id="a_node7"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/isnan.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="164.5,-152.5 164.5,-171.5 223.5,-171.5 223.5,-152.5 164.5,-152.5"/>
<text text-anchor="middle" x="194" y="-159.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::isnan</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node7 -->
<g id="edge6" class="edge">
<title>Node1&#45;&gt;Node7</title>
<path fill="none" stroke="midnightblue" d="M87.21,-190.41C107.75,-184.89 133.57,-177.96 154.56,-172.32"/>
<polygon fill="midnightblue" stroke="midnightblue" points="155.68,-175.65 164.43,-169.67 153.86,-168.89 155.68,-175.65"/>
</g>
<!-- Node8 -->
<g id="node8" class="node">
<title>Node8</title>
<g id="a_node8"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/algorithm/max.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="166.5,-114.5 166.5,-133.5 221.5,-133.5 221.5,-114.5 166.5,-114.5"/>
<text text-anchor="middle" x="194" y="-121.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::max</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node8 -->
<g id="edge7" class="edge">
<title>Node1&#45;&gt;Node8</title>
<path fill="none" stroke="midnightblue" d="M63.63,-190.42C79.66,-178.12 109.15,-156.71 137,-143 143.19,-139.95 149.99,-137.2 156.65,-134.79"/>
<polygon fill="midnightblue" stroke="midnightblue" points="157.89,-138.06 166.22,-131.52 155.63,-131.44 157.89,-138.06"/>
</g>
<!-- Node9 -->
<g id="node9" class="node">
<title>Node9</title>
<g id="a_node9"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/io/basic_ofstream/open.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="143.5,-76.5 143.5,-95.5 244.5,-95.5 244.5,-76.5 143.5,-76.5"/>
<text text-anchor="middle" x="194" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::ofstream::open</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node9 -->
<g id="edge8" class="edge">
<title>Node1&#45;&gt;Node9</title>
<path fill="none" stroke="midnightblue" d="M57.6,-190.44C70.17,-171.46 100.87,-128.8 137,-105 139.96,-103.05 143.14,-101.28 146.43,-99.67"/>
<polygon fill="midnightblue" stroke="midnightblue" points="148.1,-102.77 155.89,-95.59 145.32,-96.34 148.1,-102.77"/>
</g>
<!-- Node10 -->
<g id="node10" class="node">
<title>Node10</title>
<g id="a_node10"><a xlink:href="../../da/df2/durand__kerner__roots_8cpp.html#a51637ff89e55ca749993d3e84b95d153" target="_top" xlink:title=" ">
<polygon fill="white" stroke="black" points="157,-38.5 157,-57.5 231,-57.5 231,-38.5 157,-38.5"/>
<text text-anchor="middle" x="194" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">poly_function</text>
</a>
</g>
</g>
<!-- Node1&#45;&gt;Node10 -->
<g id="edge9" class="edge">
<title>Node1&#45;&gt;Node10</title>
<path fill="none" stroke="midnightblue" d="M55.05,-190.39C64.51,-166.09 92.66,-101.65 137,-67 140.19,-64.51 143.72,-62.34 147.41,-60.44"/>
<polygon fill="midnightblue" stroke="midnightblue" points="149.09,-63.53 156.82,-56.28 146.26,-57.12 149.09,-63.53"/>
</g>
<!-- Node11 -->
<g id="node11" class="node">
<title>Node11</title>
<g id="a_node11"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/imag.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="287,-76.5 287,-95.5 387,-95.5 387,-76.5 287,-76.5"/>
<text text-anchor="middle" x="337" y="-83.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::imag</text>
</a>
</g>
</g>
<!-- Node10&#45;&gt;Node11 -->
<g id="edge10" class="edge">
<title>Node10&#45;&gt;Node11</title>
<path fill="none" stroke="midnightblue" d="M230.58,-57.59C248.78,-62.49 271.14,-68.52 290.57,-73.76"/>
<polygon fill="midnightblue" stroke="midnightblue" points="289.86,-77.19 300.43,-76.41 291.69,-70.43 289.86,-77.19"/>
</g>
<!-- Node12 -->
<g id="node12" class="node">
<title>Node12</title>
<g id="a_node12"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/math/pow.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="310,-38.5 310,-57.5 364,-57.5 364,-38.5 310,-38.5"/>
<text text-anchor="middle" x="337" y="-45.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::pow</text>
</a>
</g>
</g>
<!-- Node10&#45;&gt;Node12 -->
<g id="edge11" class="edge">
<title>Node10&#45;&gt;Node12</title>
<path fill="none" stroke="midnightblue" d="M231.31,-48C252.25,-48 278.5,-48 299.48,-48"/>
<polygon fill="midnightblue" stroke="midnightblue" points="299.61,-51.5 309.61,-48 299.61,-44.5 299.61,-51.5"/>
</g>
<!-- Node13 -->
<g id="node13" class="node">
<title>Node13</title>
<g id="a_node13"><a target="_blank" xlink:href="http://en.cppreference.com/w/cpp/numeric/complex/real.html#" xlink:title=" ">
<polygon fill="white" stroke="black" points="289.5,-0.5 289.5,-19.5 384.5,-19.5 384.5,-0.5 289.5,-0.5"/>
<text text-anchor="middle" x="337" y="-7.5" font-family="Helvetica,sans-Serif" font-size="10.00">std::complex::real</text>
</a>
</g>
</g>
<!-- Node10&#45;&gt;Node13 -->
<g id="edge12" class="edge">
<title>Node10&#45;&gt;Node13</title>
<path fill="none" stroke="midnightblue" d="M230.58,-38.41C248.78,-33.51 271.14,-27.48 290.57,-22.24"/>
<polygon fill="midnightblue" stroke="midnightblue" points="291.69,-25.57 300.43,-19.59 289.86,-18.81 291.69,-25.57"/>
</g>
</g>
</svg>

After

Width:  |  Height:  |  Size: 10 KiB