<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://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.12"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>DFNgen: DFNmain.cpp Source File</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>
<script type="text/javascript">
  $(document).ready(initResizable);
</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>
<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">DFNgen
   &#160;<span id="projectnumber">2.0</span>
   </div>
   <div id="projectbrief">DFN Model Generator</div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.12 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
$(function() {
  initMenu('',true,false,'search.php','Search');
  $(document).ready(function() { init_search(); });
});
</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">
$(document).ready(function(){initNavTree('_d_f_nmain_8cpp_source.html','');});
</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="headertitle">
<div class="title">DFNmain.cpp</div>  </div>
</div><!--header-->
<div class="contents">
<a href="_d_f_nmain_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="preprocessor">#include &lt;iostream&gt;</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="preprocessor">#include &lt;random&gt;</span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="preprocessor">#include &lt;algorithm&gt;</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="preprocessor">#include &lt;iterator&gt;</span> <span class="comment">// std::copy</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="preprocessor">#include &lt;iomanip&gt;</span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="structures_8h.html">structures.h</a>&quot;</span></div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="insert_shape_8h.html">insertShape.h</a>&quot;</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="preprocessor">#include &lt;vector&gt;</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="input_8h.html">input.h</a>&quot;</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="output_8h.html">output.h</a>&quot;</span></div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="read_input_functions_8h.html">readInputFunctions.h</a>&quot;</span></div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="cluster_groups_8h.html">clusterGroups.h</a>&quot;</span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="math_functions_8h.html">mathFunctions.h</a>&quot;</span></div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="domain_8h.html">domain.h</a>&quot;</span> </div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="computational_geometry_8h.html">computationalGeometry.h</a>&quot;</span></div><div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="generating_points_8h.html">generatingPoints.h</a>&quot;</span> </div><div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="distributions_8h.html">distributions.h</a>&quot;</span></div><div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="fracture_estimating_8h.html">fractureEstimating.h</a>&quot;</span></div><div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="debug_functions_8h.html">debugFunctions.h</a>&quot;</span></div><div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="remove_fractures_8h.html">removeFractures.h</a>&quot;</span></div><div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;</div><div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment">// Used for automated python testing</span></div><div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="testing_8h.html">testing.h</a>&quot;</span> </div><div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;</div><div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;<span class="comment">// Below includes are used for hotkey, </span></div><div class="line"><a name="l00026"></a><span class="lineno">   26</span>&#160;<span class="comment">// Stops inserting fractures and goes directly to output</span></div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;<span class="preprocessor">#include &lt;stdlib.h&gt;</span></div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;<span class="preprocessor">#include &lt;string.h&gt;</span></div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;<span class="preprocessor">#include &lt;unistd.h&gt;</span></div><div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;<span class="preprocessor">#include &lt;sys/select.h&gt;</span></div><div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;<span class="preprocessor">#include &lt;termios.h&gt;</span></div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="hotkey_8h.html">hotkey.h</a>&quot;</span></div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;</div><div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;</div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;<span class="comment">// DO NOT CHANGE VAR NAME</span></div><div class="line"><a name="l00036"></a><span class="lineno"><a class="line" href="hotkey_8h.html#a1705dee0e4e822859804bec377506fb5">   36</a></span>&#160;<span class="keyword">struct </span>termios <a class="code" href="_d_f_nmain_8cpp.html#a1705dee0e4e822859804bec377506fb5">orig_termios</a>; <span class="comment">//used for custom console and hotkey</span></div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;</div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;<span class="comment">// Global eps</span></div><div class="line"><a name="l00039"></a><span class="lineno"><a class="line" href="input_8h.html#a974b448c3d0a13d483d36ef15de369b5">   39</a></span>&#160;<span class="keywordtype">double</span> <a class="code" href="_d_f_nmain_8cpp.html#a974b448c3d0a13d483d36ef15de369b5">eps</a>;</div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;</div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;<span class="comment">// SEE STRUCTURES.H FOR ALL STRUCTURE DEFINITIONS</span></div><div class="line"><a name="l00042"></a><span class="lineno"><a class="line" href="_d_f_nmain_8cpp.html#a3c04138a5bfe5d72780bb7e82a18e627">   42</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_d_f_nmain_8cpp.html#a3c04138a5bfe5d72780bb7e82a18e627">main</a> (<span class="keywordtype">int</span> argc, <span class="keywordtype">char</span> **argv) {</div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;    </div><div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;    <span class="comment">// Error check on cmd line input:</span></div><div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;    <span class="comment">// 1st argument = input file path </span></div><div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;    <span class="comment">// 2nd argument = output folder path    </span></div><div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;    <span class="keywordflow">if</span> (argc != 3) {</div><div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;        <span class="keywordflow">if</span> (argc == 1 ) {</div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;ERROR: DFNWorks input and output file &quot;</span> </div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;paths were not included on command line.\n&quot;</span>;</div><div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;            <span class="keywordflow">return</span> 1;</div><div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;        }</div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;        <span class="keywordflow">else</span> <span class="keywordflow">if</span> (argc == 2) {</div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;ERROR: DFNWorks output file path &quot;</span> </div><div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;was not included on command line.\n&quot;</span>;</div><div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;            <span class="keywordflow">return</span> 1;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;        }</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;    }</div><div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;         </div><div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;<span class="comment">/************* Initialize Arrays/Vectors and Structures **************/</span></div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;</div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;    <span class="comment">// Vector to store accepted polygons/fractures</span></div><div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;    std::vector&lt;Poly&gt; acceptedPoly;</div><div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;    acceptedPoly.reserve(500); </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;    <span class="comment">// Vector for storing intersections</span></div><div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;    std::vector&lt;IntPoints&gt; intPts;</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;    intPts.reserve(250);</div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;    <span class="comment">// Vector for storing triple intersection points</span></div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;    std::vector&lt;Point&gt; triplePoints;</div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;    </div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;    <span class="comment">// Vector for shape families/ stochastic families</span></div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;    std::vector&lt;Shape&gt; shapeFamilies;</div><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="comment">// Statistics structure:</span></div><div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;    <span class="comment">// Keeps track of DFN statistics (see definition in structures.h)</span></div><div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;    <a class="code" href="struct_stats.html">Stats</a> pstats;</div><div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;</div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;<span class="comment">/*********************************************************************/</span></div><div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;</div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;    <span class="comment">// Read Input File</span></div><div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;    <span class="comment">// Initialize input variables. Most input variables are global</span></div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;    <a class="code" href="read_input_8cpp.html#a9af5f928fea3353fc4009118630b3f52">getInput</a>(argv[1], shapeFamilies);    </div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;</div><div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;    <span class="comment">// Set epsilon</span></div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;    <a class="code" href="_d_f_nmain_8cpp.html#a974b448c3d0a13d483d36ef15de369b5">eps</a> = <a class="code" href="input_8h.html#a8ee9be1b5aa75abae556de3088cba6d9">h</a> * 1e-8;</div><div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;</div><div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nh: &quot;</span> &lt;&lt; <a class="code" href="input_8h.html#a8ee9be1b5aa75abae556de3088cba6d9">h</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><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="keywordtype">int</span> totalFamilies = <a class="code" href="input_8h.html#aaf09f7cf249c6fd6247821bc01ead90c">nFamEll</a> + <a class="code" href="input_8h.html#a71b31a03f80cb26d46b69c9e96fe032b">nFamRect</a>;</div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;</div><div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;    <span class="comment">// Print shape families to screen</span></div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0) {</div><div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;        <a class="code" href="debug_functions_8cpp.html#a0d86e9931b24c2b392a6a307d7725c69">printShapeFams</a>(shapeFamilies);</div><div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;    }</div><div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160;</div><div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;    <span class="comment">// Initialize random generator with seed ( see c++ &lt;random&gt; )</span></div><div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;    <span class="comment">// Mersene Twister 19937 generator (64 bit)</span></div><div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ae21f357c223957d36046a0d71cc6aed7">seed</a> == 0) {</div><div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;        <a class="code" href="input_8h.html#ae21f357c223957d36046a0d71cc6aed7">seed</a> = <a class="code" href="read_input_functions_8cpp.html#af648b2c8b7622d441c1bc50f0f9c83b0">getTimeBasedSeed</a>();</div><div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160;    }</div><div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;    std::mt19937_64 generator(<a class="code" href="input_8h.html#ae21f357c223957d36046a0d71cc6aed7">seed</a>);</div><div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;</div><div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;    <span class="comment">// Init distributions class</span></div><div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;    <span class="comment">// Currenlty used only for exponential distribution</span></div><div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;    <a class="code" href="class_distributions.html">Distributions</a> distributions(generator, shapeFamilies);</div><div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;</div><div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;    <span class="keywordtype">float</span> domVol = <a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]*<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]*<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2];</div><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="keywordflow">if</span> (totalFamilies &gt; 0 ) { </div><div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;              </div><div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 0) { <span class="comment">// Npoly Option</span></div><div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;            <span class="comment">// Estimate fractures, generate radii lists for nPoly option</span></div><div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;            <a class="code" href="fracture_estimating_8cpp.html#acb37fb0878da2fa7f0648216fb4189e4">generateRadiiLists_nPolyOption</a>(shapeFamilies, <a class="code" href="input_8h.html#aeb1032cfdc6f748b280a4a0a59375d11">famProb</a>, generator, distributions);         </div><div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160;        }</div><div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// P32 Option</span></div><div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;            <span class="comment">// ESTIMATE # FRACTURES NEEDED</span></div><div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a0071fcf088e13e14d3cfa0c15e33e363">disableFram</a> == <span class="keyword">false</span>) {</div><div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;                std::cout &lt;&lt; <span class="stringliteral">&quot;\nEstimating number of fractures needed...\n&quot;</span>;</div><div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;                <a class="code" href="fracture_estimating_8cpp.html#a7cf21b594fd24b2879af467e23bdac0a">dryRun</a>(shapeFamilies, <a class="code" href="input_8h.html#aeb1032cfdc6f748b280a4a0a59375d11">famProb</a>, generator, distributions);</div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;            }</div><div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;        }</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;        <span class="comment">// Add a percentage more radii to each radii</span></div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;        <span class="comment">// list using families&#39; distribution.</span></div><div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;        <span class="comment">// First arg is percentage, eg: 0.1 will add 10% more fractures</span></div><div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;        <span class="comment">// to the radii list for each family</span></div><div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a0071fcf088e13e14d3cfa0c15e33e363">disableFram</a> == <span class="keyword">false</span>) {</div><div class="line"><a name="l00131"></a><span class="lineno">  131</span>&#160;                <a class="code" href="fracture_estimating_8cpp.html#a26080c2634b53083b47b763b66bd4911">addRadiiToLists</a>(<a class="code" href="input_8h.html#abd9fdedaba98f63d68fe3bb9d9098f63">radiiListIncrease</a>, shapeFamilies, generator, distributions); </div><div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;                        </div><div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> j = 0; j &lt; shapeFamilies.size(); j++) {</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;                <span class="keywordflow">if</span> (shapeFamilies[j].distributionType == 4) {</div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;                    <span class="comment">// Constant size</span></div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;                    std::cout &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[j]) </div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;                              &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(j,shapeFamilies[j].shapeFamily) </div><div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;                              &lt;&lt; <span class="stringliteral">&quot; using constant size\n&quot;</span>;</div><div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;                }</div><div class="line"><a name="l00141"></a><span class="lineno">  141</span>&#160;                <span class="keywordflow">else</span> {</div><div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;                    std::cout &lt;&lt; <span class="stringliteral">&quot;Estimated &quot;</span> &lt;&lt; shapeFamilies[j].radiiList.size()</div><div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;                              &lt;&lt; <span class="stringliteral">&quot; fractures for &quot;</span> &lt;&lt;  <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[j]) </div><div class="line"><a name="l00144"></a><span class="lineno">  144</span>&#160;                              &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(j,shapeFamilies[j].shapeFamily) </div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;                              &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;                }</div><div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;            }</div><div class="line"><a name="l00148"></a><span class="lineno">  148</span>&#160;            <a class="code" href="fracture_estimating_8cpp.html#ae3e14f5bf8ea7b1470e64a74e0354947">sortRadii</a>(shapeFamilies);</div><div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;        }</div><div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;        <span class="comment">// Keep count of accepted &amp; rejected fractures by family</span></div><div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;        pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[totalFamilies];</div><div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;        pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[totalFamilies];</div><div class="line"><a name="l00153"></a><span class="lineno">  153</span>&#160;        </div><div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;        <span class="comment">// Save sizes of pre-generated radii lists per family.</span></div><div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;        <span class="comment">// Print as part of statistics to user</span></div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;        pstats.<a class="code" href="struct_stats.html#a1041c9121838585a366175f3ddd23e91">expectedFromFam</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[totalFamilies];</div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;        </div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;        <span class="comment">// Zero arrays, init expectedFromFam array</span></div><div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) {</div><div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;            pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] = 0;</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;            pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] = 0;</div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;            pstats.<a class="code" href="struct_stats.html#a1041c9121838585a366175f3ddd23e91">expectedFromFam</a>[i] = shapeFamilies[i].radiiList.size();</div><div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;        }</div><div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;        <span class="comment">// Init first rejects per insertion attempt counter</span></div><div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;        pstats.<a class="code" href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">rejectsPerAttempt</a>.push_back(0);</div><div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;</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;    <span class="comment">/*********** SETUP HOT KEY *************/</span></div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;    <span class="keywordtype">char</span> key = 0;</div><div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;<span class="preprocessor">    #ifndef TESTING </span></div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;    <span class="comment">// Set custom terminal settings for hotkey functionality</span></div><div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;    <a class="code" href="hotkey_8cpp.html#a49d261c1ea17ddedcc836585c758e3b2">set_conio_terminal_mode</a>();</div><div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;    atexit(<a class="code" href="hotkey_8cpp.html#ac8373d63f264c2e956b071b5bd2588a1">reset_terminal_mode</a>);</div><div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;<span class="preprocessor">    #endif</span></div><div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;    <span class="comment">/*********  END SETUP HOT KEY **********/</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="comment">/********************* User Defined Shapes Insertion ************************/</span>    </div><div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;    </div><div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a92d138143c536b0bc8ca70ddc1479677">insertUserRectanglesFirst</a> == 1) {</div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;        <span class="comment">// Insert user rects first</span></div><div class="line"><a name="l00182"></a><span class="lineno">  182</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ab963e6ee0512d7651d002947e0bfd1df">userRectanglesOnOff</a> != 0) {</div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;            <a class="code" href="insert_shape_8h.html#af69365c2e0148366aabe6751ffc54a79">insertUserRects</a>(acceptedPoly, intPts, pstats, triplePoints);</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;        <span class="comment">// Insert all user rectangles by coordinates</span></div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a7ece519c4cc6ecd5512562594084e8aa">userRecByCoord</a> !=0 ) {</div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;            <a class="code" href="insert_shape_8h.html#a244f461634679dc91d08da3e183f9b8c">insertUserRectsByCoord</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;        }    </div><div class="line"><a name="l00189"></a><span class="lineno">  189</span>&#160;        <span class="comment">// Insert all user ellipses </span></div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#af93cb3f9ec0bb2490e0c2219b15aa1cf">userEllipsesOnOff</a> != 0) {</div><div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;            <a class="code" href="insert_shape_8h.html#a9841b356b0460e297748d6648c301558">insertUserEll</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;        }</div><div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;        <span class="comment">// Insert all user ellipses by coordinates</span></div><div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a3718da9decab4e36a1011af51ab11bd7">userEllByCoord</a> != 0) {</div><div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;            <a class="code" href="insert_shape_8h.html#a49e6b0be012ed44c6ef83ea66ed6dd86">insertUserEllByCoord</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;        }</div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;    }</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;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;        <span class="comment">// Insert all user ellipses first</span></div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#af93cb3f9ec0bb2490e0c2219b15aa1cf">userEllipsesOnOff</a> != 0) {</div><div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;            <a class="code" href="insert_shape_8h.html#a9841b356b0460e297748d6648c301558">insertUserEll</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00203"></a><span class="lineno">  203</span>&#160;        }</div><div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;        <span class="comment">// Insert all user ellipses by coordinates</span></div><div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a3718da9decab4e36a1011af51ab11bd7">userEllByCoord</a> != 0) {</div><div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;            <a class="code" href="insert_shape_8h.html#a49e6b0be012ed44c6ef83ea66ed6dd86">insertUserEllByCoord</a>(acceptedPoly, intPts, pstats, triplePoints);</div><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="comment">// Insert user rects</span></div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ab963e6ee0512d7651d002947e0bfd1df">userRectanglesOnOff</a> != 0) {</div><div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;            <a class="code" href="insert_shape_8h.html#af69365c2e0148366aabe6751ffc54a79">insertUserRects</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;        }</div><div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;        <span class="comment">// Insert all user rectangles by coordinates</span></div><div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a7ece519c4cc6ecd5512562594084e8aa">userRecByCoord</a> !=0 ) {</div><div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;            <a class="code" href="insert_shape_8h.html#a244f461634679dc91d08da3e183f9b8c">insertUserRectsByCoord</a>(acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00215"></a><span class="lineno">  215</span>&#160;        }    </div><div class="line"><a name="l00216"></a><span class="lineno">  216</span>&#160;    }        </div><div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;    </div><div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;    <span class="comment">/*********  Probabilities (famProb) setup, CDF init  *****************/</span></div><div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;    </div><div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;    <span class="comment">// &#39;CDF&#39; size will shrink along when used with fracture intensity (P32) option</span></div><div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;    <span class="keywordtype">float</span> *CDF;</div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;    <span class="keywordtype">int</span> cdfSize = totalFamilies; </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">if</span> (totalFamilies &gt; 0) {</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;        <span class="comment">// Convert famProb to CDF</span></div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;        CDF = <a class="code" href="math_functions_8cpp.html#a9422106ad945a110b0cf696050962556">createCDF</a>(<a class="code" href="input_8h.html#aeb1032cfdc6f748b280a4a0a59375d11">famProb</a>, cdfSize);</div><div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;    }</div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;</div><div class="line"><a name="l00229"></a><span class="lineno">  229</span>&#160;    std::ofstream radiiAll;</div><div class="line"><a name="l00230"></a><span class="lineno">  230</span>&#160;    std::string output = argv[2];</div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;    std::string radiiFolder = output + <span class="stringliteral">&quot;/radii&quot;</span>;</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;    <span class="comment">// Create output folder if not already created</span></div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;    <span class="comment">//makeDIR(argv[2]); </span></div><div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;</div><div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;    <span class="comment">// Create radii folder</span></div><div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;    <span class="comment">//makeDIR(radiiFolder.c_str()); </span></div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;</div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;    <span class="comment">// Create polys folder</span></div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;    std::string  polyFolder = output + <span class="stringliteral">&quot;/polys&quot;</span>;</div><div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;    <span class="comment">//makeDIR(polyFolder.c_str());</span></div><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="keywordflow">if</span> (<a class="code" href="input_8h.html#a6fc4e503b3c7fbcfbb51832ae2201e70">outputAllRadii</a> == 1) {</div><div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;        <span class="comment">// Option to include all radii in output (accepted and rejected)</span></div><div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;        std::string file = radiiFolder + <span class="stringliteral">&quot;/radii_All.dat&quot;</span>;       </div><div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;        radiiAll.open(file.c_str(), std::ofstream::out | std::ofstream::trunc);</div><div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;        radiiAll &lt;&lt; <span class="stringliteral">&quot;Format: xRadius yRadius Family# (-1 = userRectangle, 0 = userEllipse, &gt; 0 is family in order of famProb)\n&quot;</span>;</div><div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;    }</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;    <span class="comment">// Initialize uniform distribution on [0,1]</span></div><div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;    std::uniform_real_distribution&lt;double&gt; uniformDist(0,1);</div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;    </div><div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0) {</div><div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;        </div><div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;        <span class="comment">// Holds index to current &#39;shapeFamily&#39; being inserted </span></div><div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;        <span class="keywordtype">int</span> familyIndex;</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="comment">/******************************************************************************************/</span></div><div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;        <span class="comment">/**************************             MAIN LOOP            ******************************/</span></div><div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;        </div><div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;        <span class="comment">// NOTE: p32Complete() works on global array &#39;p32Status&#39;</span></div><div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;        <span class="comment">// p32Complete() only needs argument of the number of defined shape families </span></div><div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;        <span class="comment">// ********* Begin stochastic fracture insertion ***********</span></div><div class="line"><a name="l00264"></a><span class="lineno">  264</span>&#160;        <span class="keywordflow">while</span> (((<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 0 &amp;&amp; pstats.<a class="code" href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">acceptedPolyCount</a> &lt; <a class="code" href="input_8h.html#af336784d384b9e2ab2d990f7855ecdcf">nPoly</a>) || (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 1 &amp;&amp; <a class="code" href="insert_shape_8cpp.html#a262d50919e2ca1c4e8b32ec027a52b57">p32Complete</a>(totalFamilies) == 0)) &amp;&amp; key != <span class="charliteral">&#39;~&#39;</span> ) {</div><div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;</div><div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;            <span class="comment">// cdfIdx holds the index to the CDF array for the current shape family being inserted</span></div><div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;            <span class="keywordtype">int</span> cdfIdx; </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;            <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 0 ) { <span class="comment">// nPoly Option</span></div><div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;                <span class="comment">// Choose a family based purely on famProb probabilities</span></div><div class="line"><a name="l00271"></a><span class="lineno">  271</span>&#160;                familyIndex = <a class="code" href="math_functions_8cpp.html#a5eae1db9cbad8e807a0d677bcd805d36">indexFromProb</a>(CDF, uniformDist(generator), totalFamilies);</div><div class="line"><a name="l00272"></a><span class="lineno">  272</span>&#160;            }</div><div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;            </div><div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;            <span class="comment">// Choose a family based on probabiliyis AND their target p32 completion status.</span></div><div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;            <span class="comment">// If a family has already met is fracture intinisty req. (p32) don&#39;t choose that family anymore</span></div><div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;            <span class="keywordflow">else</span> { <span class="comment">// P32 Option</span></div><div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;                familyIndex = <a class="code" href="math_functions_8cpp.html#a2540d3ded7bd488a924ce2e4c2443261">indexFromProb_and_P32Status</a>(CDF, uniformDist(generator), totalFamilies, cdfSize, cdfIdx); </div><div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;            }</div><div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;        </div><div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;            <span class="keyword">struct </span><a class="code" href="struct_poly.html">Poly</a> newPoly = <a class="code" href="insert_shape_8cpp.html#a4f91304a14fefba4d4f93f2f3e54dee4">generatePoly</a>(shapeFamilies[familyIndex], generator, distributions, familyIndex, <span class="keyword">true</span>);</div><div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;</div><div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a6fc4e503b3c7fbcfbb51832ae2201e70">outputAllRadii</a> == 1) {</div><div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;                <span class="comment">// Output all radii </span></div><div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;                radiiAll &lt;&lt; std::setprecision(8) &lt;&lt;  newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a> </div><div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;                         &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; newPoly.<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> + 1 &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;    </div><div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;            }</div><div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;            </div><div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;            <span class="keywordtype">int</span> rejectCode = 1;                                 </div><div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;            <span class="keywordtype">int</span> rejectCounter = 0; </div><div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;            </div><div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;            <span class="keywordflow">while</span> (rejectCode != 0) { <span class="comment">// Loop used to reinsert same poly with different translation</span></div><div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;            </div><div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;                <span class="comment">// HOT KEY: check for keyboard input</span></div><div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="hotkey_8cpp.html#ad5451da499ab9d3907da8dd7060ab677">kbhit</a>()) {</div><div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;                    key = <a class="code" href="hotkey_8cpp.html#a160f6d2893bb6b89ab3ad5d863c20e3d">getch</a>();</div><div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;                }</div><div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;                </div><div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;                <span class="comment">// Truncate poly if needed</span></div><div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;                <span class="comment">// Reutnrs 1 if poly is outside of domain or has less than 3 vertices</span></div><div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;                <span class="keywordflow">if</span> ( <a class="code" href="domain_8cpp.html#aed3d77c17f33d2024704b479d0885332">domainTruncation</a>(newPoly, domainSize) == 1) {</div><div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;                    <span class="comment">// Poly was completely outside domain, or was truncated to less than </span></div><div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;                    <span class="comment">// 3 vertices due to vertices being too close together   </span></div><div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;                    pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a64f92b5fcce897497d067d072443dc60">outside</a>++;</div><div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160;                    <span class="comment">// Test if newPoly has reached its limit of insertion attempts</span></div><div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;                    <span class="keywordflow">if</span> (rejectCounter &gt;= <a class="code" href="input_8h.html#ad29409829ae78f5f928f417a2a7b265b">rejectsPerFracture</a>) {</div><div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;                        <span class="keyword">delete</span>[] newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>; <span class="comment">// Created with new, delete manually</span></div><div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;                        rejectCounter++; </div><div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;                        <span class="keywordflow">break</span>; <span class="comment">// Reject poly, generate new polygon</span></div><div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;                    }</div><div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;                    <span class="keywordflow">else</span> { <span class="comment">// Retranslate poly and try again, preserving normal, size, and shape</span></div><div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;                        <a class="code" href="insert_shape_8cpp.html#afa95ede56e6ecb8918b61c5a2815d40d">reTranslatePoly</a>(newPoly, shapeFamilies[familyIndex], generator);</div><div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;                        <span class="keywordflow">continue</span>; <span class="comment">// Go to next iteration of while loop, test new translation </span></div><div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;                    }                            </div><div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;                }</div><div class="line"><a name="l00315"></a><span class="lineno">  315</span>&#160;</div><div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;                <span class="comment">// Create/assign bounding box </span></div><div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160;                <a class="code" href="computational_geometry_8cpp.html#ace23b489a120bc6d6407f7e08e1e009d">createBoundingBox</a>(newPoly);</div><div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;           </div><div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;                <span class="comment">// Find line of intersection and FRAM check</span></div><div class="line"><a name="l00320"></a><span class="lineno">  320</span>&#160;                rejectCode = <a class="code" href="computational_geometry_8cpp.html#a38a8c322dba7d8518a4c3206cd55a5f6">intersectionChecking</a>(newPoly, acceptedPoly, intPts, pstats, triplePoints);</div><div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;<span class="preprocessor">                #ifdef TESTING</span></div><div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160;                    <span class="keywordflow">if</span> (rejectCode != 0) {</div><div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;                       <span class="keywordflow">return</span> 1;</div><div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;                    }</div><div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;<span class="preprocessor">                #endif</span></div><div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;</div><div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;                <span class="comment">// IF POLY ACCEPTED:</span></div><div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;                <span class="keywordflow">if</span>(rejectCode == 0) { <span class="comment">// Intersections are ok</span></div><div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;                    <span class="comment">// Incriment counter of accepted polys</span></div><div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;                    pstats.<a class="code" href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">acceptedPolyCount</a>++;</div><div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;                    pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[familyIndex]++;</div><div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;</div><div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;                    <span class="comment">// Make new rejection counter for next fracture attempt</span></div><div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;                    pstats.<a class="code" href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">rejectsPerAttempt</a>.push_back(0);</div><div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;                    </div><div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;                    <span class="keywordflow">if</span> (newPoly.<a class="code" href="struct_poly.html#a341a2ec25f71c2f58c1af119b3990387">truncated</a> == 1) { </div><div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;                        pstats.<a class="code" href="struct_stats.html#a0467673e7585331b1f609c9f7fc1de7a">truncated</a>++;</div><div class="line"><a name="l00338"></a><span class="lineno">  338</span>&#160;                    }</div><div class="line"><a name="l00339"></a><span class="lineno">  339</span>&#160;                    </div><div class="line"><a name="l00340"></a><span class="lineno">  340</span>&#160;                    <span class="comment">// Calculate poly&#39;s area</span></div><div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;                    newPoly.<a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a> = <a class="code" href="math_functions_8cpp.html#a78276693aa55fab858a41366e5b4aa73">getArea</a>(newPoly);</div><div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;    </div><div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;                    <span class="comment">// Update P32 </span></div><div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;                    <span class="keywordflow">if</span> (shapeFamilies[familyIndex].layer == 0) { <span class="comment">// Whole domain</span></div><div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;                        shapeFamilies[familyIndex].currentP32 += newPoly.<a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>*2/domVol;</div><div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;                    }</div><div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;                    <span class="keywordflow">else</span> { <span class="comment">// Layer</span></div><div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;                        shapeFamilies[familyIndex].currentP32 += newPoly.<a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>*2/<a class="code" href="input_8h.html#a8a68f69e78f6f9c011b7688ef3c48166">layerVol</a>[shapeFamilies[familyIndex].layer-1];</div><div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;                    }</div><div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;                        </div><div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;                    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 1) {    </div><div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;                        <span class="comment">// If the last inserted pologon met the p32 reqirement, set that familiy to no longer</span></div><div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;                        <span class="comment">// insert any more fractures. ajust the CDF and familiy probabilites to account for this    </span></div><div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;                        <span class="keywordflow">if</span> (shapeFamilies[familyIndex].currentP32 &gt;= shapeFamilies[familyIndex].p32Target ) {</div><div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;                        </div><div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;                            <a class="code" href="input_8h.html#a9d32f63bd011a629cf644c120e07eda8">p32Status</a>[familyIndex] = 1; <span class="comment">// Mark family as having its p32 requirement met</span></div><div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;                            std::cout &lt;&lt; <span class="stringliteral">&quot;\nP32 For Family &quot;</span> &lt;&lt; familyIndex+1 &lt;&lt; <span class="stringliteral">&quot; Completed\n\n&quot;</span>;</div><div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;</div><div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;                            <span class="comment">// Adjust CDF, PDF. Reduce their size by 1.</span></div><div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;                            <span class="comment">// Remove the completed family&#39;s element in &#39;CDF[]&#39; and &#39;famProb[]&#39;</span></div><div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;                            <span class="comment">// Distribute the removed family probability evenly among the others and rebuild the CDF</span></div><div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;                            <span class="comment">// familyIndex = index of family&#39;s probability</span></div><div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;                            <span class="comment">// cdfIdx = index of the completed family&#39;s correspongding CDF index</span></div><div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;                            <span class="keywordflow">if</span> (cdfSize &gt; 1 ) { <span class="comment">// If there are still more families to insert</span></div><div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;                                <span class="comment">// Remove completed family from CDF and famProb</span></div><div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;                                <a class="code" href="math_functions_8cpp.html#a55646cd87f15106134e9cb832492e1f6">adjustCDF_and_famProb</a>(CDF, <a class="code" href="input_8h.html#aeb1032cfdc6f748b280a4a0a59375d11">famProb</a>, cdfSize, cdfIdx);</div><div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;                            }    </div><div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;                        }</div><div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;                    }</div><div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;                    </div><div class="line"><a name="l00371"></a><span class="lineno">  371</span>&#160;                    </div><div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;                    <span class="comment">// Output to user: print running program status to user</span></div><div class="line"><a name="l00373"></a><span class="lineno">  373</span>&#160;                    <span class="keywordflow">if</span> (pstats.<a class="code" href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">acceptedPolyCount</a> % 200 == 0) {</div><div class="line"><a name="l00374"></a><span class="lineno">  374</span>&#160;                          std::cout &lt;&lt; <span class="stringliteral">&quot;\nAccepted &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">acceptedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot; fractures\n&quot;</span>;</div><div class="line"><a name="l00375"></a><span class="lineno">  375</span>&#160;                          std::cout &lt;&lt; <span class="stringliteral">&quot;Rejected &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ae78598e8e15f944d358722868f763e32">rejectedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot; fractures\n&quot;</span>;</div><div class="line"><a name="l00376"></a><span class="lineno">  376</span>&#160;                          std::cout &lt;&lt; <span class="stringliteral">&quot;Re-translated &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a0fe2c18e0646b84364712352595751ec">retranslatedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot; fractures\n\n&quot;</span>;</div><div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160;</div><div class="line"><a name="l00378"></a><span class="lineno">  378</span>&#160;                          std::cout &lt;&lt; <span class="stringliteral">&quot;Current p32 values per family:\n&quot;</span>;</div><div class="line"><a name="l00379"></a><span class="lineno">  379</span>&#160;                          <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) {</div><div class="line"><a name="l00380"></a><span class="lineno">  380</span>&#160;                              </div><div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;                              <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 0) {</div><div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;                                  std::cout &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> </div><div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;                                            &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i,shapeFamilies[i].shapeFamily) </div><div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;                                            &lt;&lt; <span class="stringliteral">&quot; Current P32 = &quot;</span> &lt;&lt; std::setprecision(8) </div><div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;                                            &lt;&lt; shapeFamilies[i].currentP32;</div><div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;                              }</div><div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;                              <span class="keywordflow">else</span> {</div><div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;                                  std::cout &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> </div><div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;                                            &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i,shapeFamilies[i].shapeFamily) </div><div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;                                            &lt;&lt; <span class="stringliteral">&quot; target P32 = &quot;</span> &lt;&lt; std::setprecision(8) </div><div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;                                            &lt;&lt; shapeFamilies[i].p32Target </div><div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160;                                            &lt;&lt; <span class="stringliteral">&quot;, &quot;</span> &lt;&lt; <span class="stringliteral">&quot;Current P32 = &quot;</span> &lt;&lt; shapeFamilies[i].currentP32;</div><div class="line"><a name="l00393"></a><span class="lineno">  393</span>&#160;                              }</div><div class="line"><a name="l00394"></a><span class="lineno">  394</span>&#160;                              </div><div class="line"><a name="l00395"></a><span class="lineno">  395</span>&#160;                              <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 1 &amp;&amp; shapeFamilies[i].p32Target &lt;= shapeFamilies[i].currentP32) {</div><div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;                                  std::cout &lt;&lt; <span class="stringliteral">&quot;...Done\n&quot;</span>;                          </div><div class="line"><a name="l00397"></a><span class="lineno">  397</span>&#160;                              } </div><div class="line"><a name="l00398"></a><span class="lineno">  398</span>&#160;                              <span class="keywordflow">else</span> {</div><div class="line"><a name="l00399"></a><span class="lineno">  399</span>&#160;                                  std::cout &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00400"></a><span class="lineno">  400</span>&#160;                              }    </div><div class="line"><a name="l00401"></a><span class="lineno">  401</span>&#160;                          }                      </div><div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160;                    }        </div><div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;                </div><div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;                    <span class="comment">// SAVING POLYGON (intersection and triple points saved witchin intersectionChecking())        </span></div><div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;                    acceptedPoly.push_back(newPoly); <span class="comment">// SAVE newPoly to accepted polys list</span></div><div class="line"><a name="l00406"></a><span class="lineno">  406</span>&#160;</div><div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;                }</div><div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;            </div><div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;                <span class="keywordflow">else</span> { <span class="comment">// Poly rejected</span></div><div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;</div><div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;                    <span class="comment">// Inc reject counter for current poly</span></div><div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;                    rejectCounter++;</div><div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;                    </div><div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;                    <span class="comment">// Inc reject counter for current attempt</span></div><div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;                    <span class="comment">// (number of rejects until next fracture accepted)</span></div><div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;                    pstats.<a class="code" href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">rejectsPerAttempt</a>[pstats.<a class="code" href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">acceptedPolyCount</a>]++;</div><div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;                                </div><div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;                                        </div><div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;                    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a9b979283d8d40a1417171c28a1b1d23f">printRejectReasons</a> != 0) {</div><div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;                        <a class="code" href="insert_shape_8cpp.html#ae9670e906ee5dc43ce3fbc0ce7290dc7">printRejectReason</a>(rejectCode, newPoly);</div><div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;                    }</div><div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;                    </div><div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;                    <span class="keywordflow">if</span> (rejectCounter &gt;= <a class="code" href="input_8h.html#ad29409829ae78f5f928f417a2a7b265b">rejectsPerFracture</a>) {</div><div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;                        <span class="keyword">delete</span>[] newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>; <span class="comment">// Delete manually, created with new[]</span></div><div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;                        pstats.<a class="code" href="struct_stats.html#ae78598e8e15f944d358722868f763e32">rejectedPolyCount</a>++;</div><div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;                        pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[familyIndex]++;</div><div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;                        <span class="comment">// Stop retranslating polygon if its reached its reject limit</span></div><div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;                        <span class="keywordflow">break</span>; <span class="comment">// Break will cause code to go to next poly</span></div><div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;                    }</div><div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;                    <span class="keywordflow">else</span> { </div><div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;                        <span class="comment">// Translate poly to new position</span></div><div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;                        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a9b979283d8d40a1417171c28a1b1d23f">printRejectReasons</a> != 0) {</div><div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;                            std::cout &lt;&lt; <span class="stringliteral">&quot;Translating rejected fracture to new position\n&quot;</span>;</div><div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;                        }</div><div class="line"><a name="l00435"></a><span class="lineno">  435</span>&#160;                        pstats.<a class="code" href="struct_stats.html#a0fe2c18e0646b84364712352595751ec">retranslatedPolyCount</a>++;</div><div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;                        <a class="code" href="insert_shape_8cpp.html#afa95ede56e6ecb8918b61c5a2815d40d">reTranslatePoly</a>(newPoly, shapeFamilies[familyIndex], generator);</div><div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160;                    }</div><div class="line"><a name="l00438"></a><span class="lineno">  438</span>&#160;                } <span class="comment">// End else poly rejected</span></div><div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;            } <span class="comment">// End loop while for re-translating polys option (reject == 1)</span></div><div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;        } <span class="comment">// !!!!  END MAIN LOOP !!!! end while loop for inserting polyons</span></div><div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;</div><div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;        <span class="comment">// Remove last element off the rejects per attempt counter. </span></div><div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;        <span class="comment">// It will have one extra item due to how eachelement is initialized. </span></div><div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160;        <span class="keywordflow">if</span> (!pstats.<a class="code" href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">rejectsPerAttempt</a>.empty()) {</div><div class="line"><a name="l00445"></a><span class="lineno">  445</span>&#160;            pstats.<a class="code" href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">rejectsPerAttempt</a>.pop_back();</div><div class="line"><a name="l00446"></a><span class="lineno">  446</span>&#160;        }</div><div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;</div><div class="line"><a name="l00448"></a><span class="lineno">  448</span>&#160;    } <span class="comment">// End if totalFamilies != 0</span></div><div class="line"><a name="l00449"></a><span class="lineno">  449</span>&#160; </div><div class="line"><a name="l00450"></a><span class="lineno">  450</span>&#160;    <span class="comment">/************************** DFN GENERATION IS COMPLETE ***************************/</span></div><div class="line"><a name="l00451"></a><span class="lineno">  451</span>&#160;<span class="comment">//    printIntersectionData(intPts);</span></div><div class="line"><a name="l00452"></a><span class="lineno">  452</span>&#160;<span class="comment">//    printGroupData(pstats,acceptedPoly);</span></div><div class="line"><a name="l00453"></a><span class="lineno">  453</span>&#160;    </div><div class="line"><a name="l00454"></a><span class="lineno">  454</span>&#160;    <span class="comment">// The close to node check is inside of the close to edge check function</span></div><div class="line"><a name="l00455"></a><span class="lineno">  455</span>&#160;    <span class="comment">// for optimization (only need do intersection close to node on one condition)</span></div><div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;    <span class="comment">// On close to node rejections, close to edge is counted as well.</span></div><div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;    <span class="comment">// To get the correct number we must subtract the close to node count </span></div><div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;    <span class="comment">// (they were counted in closeToEdge AND closeToNode)</span></div><div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;    pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#adbb08b81bcca27734f30fc8314a48719">closeToEdge</a> -= pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#ac581f8584a6bc05181a18c9145f2a26e">closeToNode</a>;</div><div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;   </div><div class="line"><a name="l00461"></a><span class="lineno">  461</span>&#160;    </div><div class="line"><a name="l00462"></a><span class="lineno">  462</span>&#160;</div><div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;<span class="preprocessor">    #ifndef TESTING    </span></div><div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;    <a class="code" href="hotkey_8cpp.html#ac8373d63f264c2e956b071b5bd2588a1">reset_terminal_mode</a>();</div><div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160;<span class="preprocessor">    #endif</span></div><div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;</div><div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a6fc4e503b3c7fbcfbb51832ae2201e70">outputAllRadii</a> == 1) {</div><div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160;        radiiAll.close();    </div><div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;    }</div><div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160;    </div><div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;    <span class="comment">// Assign apertures and permiability to accepted polygons</span></div><div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i &lt; acceptedPoly.size(); i++) {</div><div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;        <a class="code" href="insert_shape_8cpp.html#a19e9a86f45510bc87d7fe84a11caffcf">assignAperture</a>(acceptedPoly[i], generator);</div><div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;        <span class="comment">// NOTE: must assign aperture before permeability</span></div><div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;        <a class="code" href="insert_shape_8cpp.html#ab66dd8460f2f45feae57425d7c219bf5">assignPermeability</a>(acceptedPoly[i]);</div><div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;    }</div><div class="line"><a name="l00477"></a><span class="lineno">  477</span>&#160;</div><div class="line"><a name="l00478"></a><span class="lineno">  478</span>&#160;    <span class="comment">// Copy end of DFN generation stats to file, as well as print to screen</span></div><div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;    std::ofstream file;</div><div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;    std::string fileName = output + <span class="stringliteral">&quot;/DFN_output.txt&quot;</span>;</div><div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;    file.open(fileName.c_str(), std::ofstream::out | std::ofstream::trunc);</div><div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n========================================================\n&quot;</span>;</div><div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;            Network Generation Complete\n&quot;</span>;</div><div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;========================================================\n&quot;</span>; </div><div class="line"><a name="l00485"></a><span class="lineno">  485</span>&#160;   </div><div class="line"><a name="l00486"></a><span class="lineno">  486</span>&#160;</div><div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n========================================================\n&quot;</span>;</div><div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;            Network Generation Complete\n&quot;</span>;</div><div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;========================================================\n&quot;</span>; </div><div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160;</div><div class="line"><a name="l00491"></a><span class="lineno">  491</span>&#160;   </div><div class="line"><a name="l00492"></a><span class="lineno">  492</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a> == 1 ) {</div><div class="line"><a name="l00493"></a><span class="lineno">  493</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\nFinal p32 values per family:\n&quot;</span>;</div><div class="line"><a name="l00494"></a><span class="lineno">  494</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) {</div><div class="line"><a name="l00495"></a><span class="lineno">  495</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot; target P32 = &quot;</span> &lt;&lt; shapeFamilies[i].p32Target </div><div class="line"><a name="l00496"></a><span class="lineno">  496</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;, &quot;</span> &lt;&lt; <span class="stringliteral">&quot;Final P32 = &quot;</span> &lt;&lt; shapeFamilies[i].currentP32 &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>; </div><div class="line"><a name="l00497"></a><span class="lineno">  497</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot; target P32 = &quot;</span> &lt;&lt; shapeFamilies[i].p32Target </div><div class="line"><a name="l00498"></a><span class="lineno">  498</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;, &quot;</span> &lt;&lt; <span class="stringliteral">&quot;Final P32 = &quot;</span> &lt;&lt; shapeFamilies[i].currentP32 &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;              </div><div class="line"><a name="l00499"></a><span class="lineno">  499</span>&#160;        }</div><div class="line"><a name="l00500"></a><span class="lineno">  500</span>&#160;        </div><div class="line"><a name="l00501"></a><span class="lineno">  501</span>&#160;    }</div><div class="line"><a name="l00502"></a><span class="lineno">  502</span>&#160;    </div><div class="line"><a name="l00503"></a><span class="lineno">  503</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n&quot;</span>;</div><div class="line"><a name="l00504"></a><span class="lineno">  504</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n&quot;</span>;</div><div class="line"><a name="l00505"></a><span class="lineno">  505</span>&#160;</div><div class="line"><a name="l00506"></a><span class="lineno">  506</span>&#160;    <span class="comment">// Calculate total area, volume</span></div><div class="line"><a name="l00507"></a><span class="lineno">  507</span>&#160;    <span class="keywordtype">double</span> userDefinedShapesArea = 0;</div><div class="line"><a name="l00508"></a><span class="lineno">  508</span>&#160;    <span class="keywordtype">double</span> userDefinedVol = 0;</div><div class="line"><a name="l00509"></a><span class="lineno">  509</span>&#160;    <span class="keywordtype">double</span> *familyArea = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00510"></a><span class="lineno">  510</span>&#160;    <span class="keywordtype">double</span> *familyVol = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00511"></a><span class="lineno">  511</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0 ) {</div><div class="line"><a name="l00512"></a><span class="lineno">  512</span>&#160;        familyArea = <span class="keyword">new</span> <span class="keywordtype">double</span>[totalFamilies]; <span class="comment">// Holds fracture area per family</span></div><div class="line"><a name="l00513"></a><span class="lineno">  513</span>&#160;        familyVol = <span class="keyword">new</span> <span class="keywordtype">double</span>[totalFamilies];</div><div class="line"><a name="l00514"></a><span class="lineno">  514</span>&#160;        </div><div class="line"><a name="l00515"></a><span class="lineno">  515</span>&#160;        <span class="comment">// Zero array</span></div><div class="line"><a name="l00516"></a><span class="lineno">  516</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) { </div><div class="line"><a name="l00517"></a><span class="lineno">  517</span>&#160;            familyArea[i] = 0;</div><div class="line"><a name="l00518"></a><span class="lineno">  518</span>&#160;            familyVol[i] = 0; </div><div class="line"><a name="l00519"></a><span class="lineno">  519</span>&#160;        }</div><div class="line"><a name="l00520"></a><span class="lineno">  520</span>&#160;        </div><div class="line"><a name="l00521"></a><span class="lineno">  521</span>&#160;    }</div><div class="line"><a name="l00522"></a><span class="lineno">  522</span>&#160;    </div><div class="line"><a name="l00523"></a><span class="lineno">  523</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nStatistics Before Isolated Fractures Removed:\n\n&quot;</span>;</div><div class="line"><a name="l00524"></a><span class="lineno">  524</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Fractures: &quot;</span> &lt;&lt; acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00525"></a><span class="lineno">  525</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Truncated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a0467673e7585331b1f609c9f7fc1de7a">truncated</a> &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00526"></a><span class="lineno">  526</span>&#160;</div><div class="line"><a name="l00527"></a><span class="lineno">  527</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\nStatistics Before Isolated Fractures Removed:\n\n&quot;</span>;</div><div class="line"><a name="l00528"></a><span class="lineno">  528</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Fractures: &quot;</span> &lt;&lt; acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00529"></a><span class="lineno">  529</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Truncated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a0467673e7585331b1f609c9f7fc1de7a">truncated</a> &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00530"></a><span class="lineno">  530</span>&#160;      </div><div class="line"><a name="l00531"></a><span class="lineno">  531</span>&#160;    <span class="comment">// Calculate total fracture area, and area per family</span></div><div class="line"><a name="l00532"></a><span class="lineno">  532</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i &lt; acceptedPoly.size(); i++) {</div><div class="line"><a name="l00533"></a><span class="lineno">  533</span>&#160;        <span class="keywordtype">double</span> <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a> = acceptedPoly[i].area;</div><div class="line"><a name="l00534"></a><span class="lineno">  534</span>&#160;        <span class="keywordtype">double</span> vol = area * acceptedPoly[i].aperture;</div><div class="line"><a name="l00535"></a><span class="lineno">  535</span>&#160;</div><div class="line"><a name="l00536"></a><span class="lineno">  536</span>&#160;        pstats.<a class="code" href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">areaBeforeRemoval</a> += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00537"></a><span class="lineno">  537</span>&#160;        pstats.<a class="code" href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">volBeforeRemoval</a> += vol;</div><div class="line"><a name="l00538"></a><span class="lineno">  538</span>&#160;    </div><div class="line"><a name="l00539"></a><span class="lineno">  539</span>&#160;        <span class="keywordflow">if</span> (acceptedPoly[i].<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> &gt;= 0) {</div><div class="line"><a name="l00540"></a><span class="lineno">  540</span>&#160;                familyArea[acceptedPoly[i].familyNum] += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00541"></a><span class="lineno">  541</span>&#160;                familyVol[acceptedPoly[i].familyNum] += vol;</div><div class="line"><a name="l00542"></a><span class="lineno">  542</span>&#160;            }</div><div class="line"><a name="l00543"></a><span class="lineno">  543</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// User-defined polygon</span></div><div class="line"><a name="l00544"></a><span class="lineno">  544</span>&#160;            userDefinedShapesArea += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00545"></a><span class="lineno">  545</span>&#160;            userDefinedVol += vol;</div><div class="line"><a name="l00546"></a><span class="lineno">  546</span>&#160;        }    </div><div class="line"><a name="l00547"></a><span class="lineno">  547</span>&#160;    }</div><div class="line"><a name="l00548"></a><span class="lineno">  548</span>&#160;</div><div class="line"><a name="l00549"></a><span class="lineno">  549</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Surface Area:     &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">areaBeforeRemoval</a> * 2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;  </div><div class="line"><a name="l00550"></a><span class="lineno">  550</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fractures Volume: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">volBeforeRemoval</a> &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00551"></a><span class="lineno">  551</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Density   (P30): &quot;</span> &lt;&lt;acceptedPoly.size()/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00552"></a><span class="lineno">  552</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Intensity (P32): &quot;</span> &lt;&lt; (pstats.<a class="code" href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">areaBeforeRemoval</a> * 2)/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00553"></a><span class="lineno">  553</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Porosity  (P33): &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">volBeforeRemoval</a>/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00554"></a><span class="lineno">  554</span>&#160;    </div><div class="line"><a name="l00555"></a><span class="lineno">  555</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Surface Area:     &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">areaBeforeRemoval</a> * 2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;  </div><div class="line"><a name="l00556"></a><span class="lineno">  556</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fractures Volume: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">volBeforeRemoval</a> &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00557"></a><span class="lineno">  557</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Density   (P30): &quot;</span> &lt;&lt;acceptedPoly.size()/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00558"></a><span class="lineno">  558</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Intensity (P32): &quot;</span> &lt;&lt; (pstats.<a class="code" href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">areaBeforeRemoval</a> * 2)/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00559"></a><span class="lineno">  559</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Porosity  (P33): &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">volBeforeRemoval</a>/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00560"></a><span class="lineno">  560</span>&#160;        </div><div class="line"><a name="l00561"></a><span class="lineno">  561</span>&#160;    <span class="comment">// Print family stats to user</span></div><div class="line"><a name="l00562"></a><span class="lineno">  562</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) { </div><div class="line"><a name="l00563"></a><span class="lineno">  563</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot;:\n&quot;</span>;</div><div class="line"><a name="l00564"></a><span class="lineno">  564</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Accepted: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00565"></a><span class="lineno">  565</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00566"></a><span class="lineno">  566</span>&#160;</div><div class="line"><a name="l00567"></a><span class="lineno">  567</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot;:\n&quot;</span>;</div><div class="line"><a name="l00568"></a><span class="lineno">  568</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Accepted: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00569"></a><span class="lineno">  569</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00570"></a><span class="lineno">  570</span>&#160;</div><div class="line"><a name="l00571"></a><span class="lineno">  571</span>&#160;        <span class="keywordflow">if</span> ( shapeFamilies[i].layer &gt; 0) {</div><div class="line"><a name="l00572"></a><span class="lineno">  572</span>&#160;            <span class="keywordtype">int</span> idx = (shapeFamilies[i].layer-1)*2;</div><div class="line"><a name="l00573"></a><span class="lineno">  573</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer: &quot;</span> &lt;&lt; shapeFamilies[i].layer &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00574"></a><span class="lineno">  574</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer {-z, +z}: {&quot;</span> &lt;&lt;<a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[idx]&lt;&lt; <span class="stringliteral">&quot;, &quot;</span> &lt;&lt;<a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[idx+1]&lt;&lt; <span class="stringliteral">&quot;}\n&quot;</span>;</div><div class="line"><a name="l00575"></a><span class="lineno">  575</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer: &quot;</span> &lt;&lt; shapeFamilies[i].layer &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00576"></a><span class="lineno">  576</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer {-z, +z}: {&quot;</span> &lt;&lt;layers[idx]&lt;&lt; <span class="stringliteral">&quot;, &quot;</span> &lt;&lt;layers[idx+1]&lt;&lt; <span class="stringliteral">&quot;}\n&quot;</span>;</div><div class="line"><a name="l00577"></a><span class="lineno">  577</span>&#160;        }</div><div class="line"><a name="l00578"></a><span class="lineno">  578</span>&#160;        <span class="keywordflow">else</span> {</div><div class="line"><a name="l00579"></a><span class="lineno">  579</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer: Whole Domain \n&quot;</span>;</div><div class="line"><a name="l00580"></a><span class="lineno">  580</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer: Whole Domain \n&quot;</span>;</div><div class="line"><a name="l00581"></a><span class="lineno">  581</span>&#160;        }</div><div class="line"><a name="l00582"></a><span class="lineno">  582</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; familyArea[i]*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00583"></a><span class="lineno">  583</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; familyVol[i] &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00584"></a><span class="lineno">  584</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; shapeFamilies[i].currentP32 &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00585"></a><span class="lineno">  585</span>&#160;</div><div class="line"><a name="l00586"></a><span class="lineno">  586</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; familyArea[i]*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00587"></a><span class="lineno">  587</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; familyVol[i] &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00588"></a><span class="lineno">  588</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; shapeFamilies[i].currentP32 &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00589"></a><span class="lineno">  589</span>&#160;</div><div class="line"><a name="l00590"></a><span class="lineno">  590</span>&#160;      }    </div><div class="line"><a name="l00591"></a><span class="lineno">  591</span>&#160;     </div><div class="line"><a name="l00592"></a><span class="lineno">  592</span>&#160;      <span class="keywordflow">if</span> (userDefinedShapesArea &gt; 0) {</div><div class="line"><a name="l00593"></a><span class="lineno">  593</span>&#160;          std::cout &lt;&lt; <span class="stringliteral">&quot;User Defined: \n&quot;</span>;</div><div class="line"><a name="l00594"></a><span class="lineno">  594</span>&#160;          std::cout &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; userDefinedShapesArea*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00595"></a><span class="lineno">  595</span>&#160;          std::cout &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; userDefinedVol &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00596"></a><span class="lineno">  596</span>&#160;          std::cout &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; userDefinedShapesArea*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00597"></a><span class="lineno">  597</span>&#160;</div><div class="line"><a name="l00598"></a><span class="lineno">  598</span>&#160;          file &lt;&lt; <span class="stringliteral">&quot;User Defined: \n&quot;</span>;</div><div class="line"><a name="l00599"></a><span class="lineno">  599</span>&#160;          file &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; userDefinedShapesArea*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00600"></a><span class="lineno">  600</span>&#160;          file &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; userDefinedVol &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00601"></a><span class="lineno">  601</span>&#160;          file &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; userDefinedShapesArea*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00602"></a><span class="lineno">  602</span>&#160;</div><div class="line"><a name="l00603"></a><span class="lineno">  603</span>&#160;      }</div><div class="line"><a name="l00604"></a><span class="lineno">  604</span>&#160;  </div><div class="line"><a name="l00605"></a><span class="lineno">  605</span>&#160;</div><div class="line"><a name="l00606"></a><span class="lineno">  606</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a764fda4f24066351261eb6f777d2ff6c">removeFracturesLessThan</a> &gt; 0) {</div><div class="line"><a name="l00607"></a><span class="lineno">  607</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\nRemoving fractures with radius less than &quot;</span> &lt;&lt; <a class="code" href="input_8h.html#a764fda4f24066351261eb6f777d2ff6c">removeFracturesLessThan</a> &lt;&lt; <span class="stringliteral">&quot; and rebuilding DFN\n&quot;</span>;</div><div class="line"><a name="l00608"></a><span class="lineno">  608</span>&#160;        file      &lt;&lt; <span class="stringliteral">&quot;\nRemoving fractures with radius less than &quot;</span> &lt;&lt; removeFracturesLessThan &lt;&lt; <span class="stringliteral">&quot; and rebuilding DFN\n&quot;</span>;</div><div class="line"><a name="l00609"></a><span class="lineno">  609</span>&#160;        <span class="keywordtype">int</span> size = acceptedPoly.size();</div><div class="line"><a name="l00610"></a><span class="lineno">  610</span>&#160;        <a class="code" href="remove_fractures_8cpp.html#a016a4a3e0104523549279e6f6cf6c6d8">removeFractures</a>(removeFracturesLessThan, acceptedPoly, intPts, triplePoints, pstats);</div><div class="line"><a name="l00611"></a><span class="lineno">  611</span>&#160;        </div><div class="line"><a name="l00612"></a><span class="lineno">  612</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;Removed &quot;</span> &lt;&lt; size - acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot; fractures with radius less than &quot;</span> &lt;&lt; removeFracturesLessThan &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00613"></a><span class="lineno">  613</span>&#160;        file      &lt;&lt; <span class="stringliteral">&quot;Removed &quot;</span> &lt;&lt; size - acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot; fractures with radius less than &quot;</span> &lt;&lt; removeFracturesLessThan &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00614"></a><span class="lineno">  614</span>&#160;    }</div><div class="line"><a name="l00615"></a><span class="lineno">  615</span>&#160;</div><div class="line"><a name="l00616"></a><span class="lineno">  616</span>&#160;</div><div class="line"><a name="l00617"></a><span class="lineno">  617</span>&#160;<span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; acceptedPoly.size(); i++) {</div><div class="line"><a name="l00618"></a><span class="lineno">  618</span>&#160;    <a class="code" href="debug_functions_8cpp.html#a77317c191140e58c1ca12f2badadbc8b">printPolyData</a>(acceptedPoly[i]);</div><div class="line"><a name="l00619"></a><span class="lineno">  619</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00620"></a><span class="lineno">  620</span>&#160;</div><div class="line"><a name="l00621"></a><span class="lineno">  621</span>&#160;}</div><div class="line"><a name="l00622"></a><span class="lineno">  622</span>&#160;</div><div class="line"><a name="l00623"></a><span class="lineno">  623</span>&#160;</div><div class="line"><a name="l00624"></a><span class="lineno">  624</span>&#160;</div><div class="line"><a name="l00625"></a><span class="lineno">  625</span>&#160;<span class="comment">/*  Remove any isolated fractures and return </span></div><div class="line"><a name="l00626"></a><span class="lineno">  626</span>&#160;<span class="comment">    a list of polygon indices matching the users</span></div><div class="line"><a name="l00627"></a><span class="lineno">  627</span>&#160;<span class="comment">    boundaryFaces option. If input option </span></div><div class="line"><a name="l00628"></a><span class="lineno">  628</span>&#160;<span class="comment">    keepOnlyLargestCluster = 1, return largest </span></div><div class="line"><a name="l00629"></a><span class="lineno">  629</span>&#160;<span class="comment">    cluster matching users boundaryFaces option</span></div><div class="line"><a name="l00630"></a><span class="lineno">  630</span>&#160;<span class="comment">    If ignoreBoundaryFaces input option is on, </span></div><div class="line"><a name="l00631"></a><span class="lineno">  631</span>&#160;<span class="comment">    DFN will keep all fractures with intersections.</span></div><div class="line"><a name="l00632"></a><span class="lineno">  632</span>&#160;<span class="comment">*/</span>     </div><div class="line"><a name="l00633"></a><span class="lineno">  633</span>&#160;    std::vector&lt;unsigned int&gt; finalFractures =  <a class="code" href="cluster_groups_8cpp.html#aa3f6bef1533919464df1f55a65fcbad3">getCluster</a>(pstats);</div><div class="line"><a name="l00634"></a><span class="lineno">  634</span>&#160;    <span class="comment">// Sort fracture indecies to retain order by acceptance</span></div><div class="line"><a name="l00635"></a><span class="lineno">  635</span>&#160;    std::sort (finalFractures.begin(), finalFractures.end()); </div><div class="line"><a name="l00636"></a><span class="lineno">  636</span>&#160;</div><div class="line"><a name="l00637"></a><span class="lineno">  637</span>&#160;    <span class="comment">// Error check for no boundary connection </span></div><div class="line"><a name="l00638"></a><span class="lineno">  638</span>&#160;    <span class="keywordtype">bool</span> printConnectivityError = 0;</div><div class="line"><a name="l00639"></a><span class="lineno">  639</span>&#160;    <span class="keywordflow">if</span> (finalFractures.size() == 0 &amp;&amp; <a class="code" href="input_8h.html#ada6852d4da13aa222e5b959bf71e89d4">ignoreBoundaryFaces</a> == 0 ) {</div><div class="line"><a name="l00640"></a><span class="lineno">  640</span>&#160;        printConnectivityError = 1;</div><div class="line"><a name="l00641"></a><span class="lineno">  641</span>&#160;        <span class="comment">//if there is no fracture network connected useres liststed boundary faces</span></div><div class="line"><a name="l00642"></a><span class="lineno">  642</span>&#160;        <span class="comment">//switch to ignore boundary faces option with notice to user that there is no connectivity</span></div><div class="line"><a name="l00643"></a><span class="lineno">  643</span>&#160;        <a class="code" href="input_8h.html#ada6852d4da13aa222e5b959bf71e89d4">ignoreBoundaryFaces</a> = 1; <span class="comment">// NOW EXITING</span></div><div class="line"><a name="l00644"></a><span class="lineno">  644</span>&#160;        finalFractures =  <a class="code" href="cluster_groups_8cpp.html#aa3f6bef1533919464df1f55a65fcbad3">getCluster</a>(pstats);</div><div class="line"><a name="l00645"></a><span class="lineno">  645</span>&#160;        <span class="comment">//if still no fractures, there is no fracture network</span></div><div class="line"><a name="l00646"></a><span class="lineno">  646</span>&#160;    }</div><div class="line"><a name="l00647"></a><span class="lineno">  647</span>&#160;    </div><div class="line"><a name="l00648"></a><span class="lineno">  648</span>&#160;    <span class="keywordflow">if</span> (finalFractures.size() == 0) {</div><div class="line"><a name="l00649"></a><span class="lineno">  649</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\nERROR: DFN Generation has finished, however&quot;</span> </div><div class="line"><a name="l00650"></a><span class="lineno">  650</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot; there are no intersecting fractures.&quot;</span> </div><div class="line"><a name="l00651"></a><span class="lineno">  651</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot; Please adjust input parameters.\n&quot;</span>;</div><div class="line"><a name="l00652"></a><span class="lineno">  652</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;Try increasing the fracture density, or shrinking the domain.\n&quot;</span>;</div><div class="line"><a name="l00653"></a><span class="lineno">  653</span>&#160;</div><div class="line"><a name="l00654"></a><span class="lineno">  654</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;\nERROR: DFN Generation has finished, however&quot;</span> </div><div class="line"><a name="l00655"></a><span class="lineno">  655</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot; there are no intersecting fractures.&quot;</span> </div><div class="line"><a name="l00656"></a><span class="lineno">  656</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot; Please adjust input parameters.\n&quot;</span>;</div><div class="line"><a name="l00657"></a><span class="lineno">  657</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;Try increasing the fracture density, or shrinking the domain.\n&quot;</span>;</div><div class="line"><a name="l00658"></a><span class="lineno">  658</span>&#160;        file.close();</div><div class="line"><a name="l00659"></a><span class="lineno">  659</span>&#160;        exit(1);</div><div class="line"><a name="l00660"></a><span class="lineno">  660</span>&#160;    }</div><div class="line"><a name="l00661"></a><span class="lineno">  661</span>&#160;    </div><div class="line"><a name="l00662"></a><span class="lineno">  662</span>&#160;    <span class="comment">/************************* Print Statistics to User ***********************************/</span></div><div class="line"><a name="l00663"></a><span class="lineno">  663</span>&#160;     </div><div class="line"><a name="l00664"></a><span class="lineno">  664</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;</div><div class="line"><a name="l00665"></a><span class="lineno">  665</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Statistics After Isolated Fractures Removed:\n&quot;</span>;</div><div class="line"><a name="l00666"></a><span class="lineno">  666</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Final Number of Fractures: &quot;</span> &lt;&lt; finalFractures.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00667"></a><span class="lineno">  667</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Isolated Fractures Removed: &quot;</span> &lt;&lt; acceptedPoly.size() - finalFractures.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00668"></a><span class="lineno">  668</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Fractures before isolated fractures removed:: &quot;</span> &lt;&lt; acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00669"></a><span class="lineno">  669</span>&#160; </div><div class="line"><a name="l00670"></a><span class="lineno">  670</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;</div><div class="line"><a name="l00671"></a><span class="lineno">  671</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Statistics After Isolated Fractures Removed:\n&quot;</span>;</div><div class="line"><a name="l00672"></a><span class="lineno">  672</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Final Number of Fractures: &quot;</span> &lt;&lt; finalFractures.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00673"></a><span class="lineno">  673</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Isolated Fractures Removed: &quot;</span> &lt;&lt; acceptedPoly.size() - finalFractures.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00674"></a><span class="lineno">  674</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Fractures before isolated fractures removed:: &quot;</span> &lt;&lt; acceptedPoly.size() &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00675"></a><span class="lineno">  675</span>&#160;</div><div class="line"><a name="l00676"></a><span class="lineno">  676</span>&#160;   </div><div class="line"><a name="l00677"></a><span class="lineno">  677</span>&#160;    <span class="comment">// Reset totalVolume and totalArea to 0</span></div><div class="line"><a name="l00678"></a><span class="lineno">  678</span>&#160;    userDefinedShapesArea = 0;</div><div class="line"><a name="l00679"></a><span class="lineno">  679</span>&#160;    userDefinedVol = 0;</div><div class="line"><a name="l00680"></a><span class="lineno">  680</span>&#160;    </div><div class="line"><a name="l00681"></a><span class="lineno">  681</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0 ) {</div><div class="line"><a name="l00682"></a><span class="lineno">  682</span>&#160;        <span class="comment">//zero out array.</span></div><div class="line"><a name="l00683"></a><span class="lineno">  683</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) { </div><div class="line"><a name="l00684"></a><span class="lineno">  684</span>&#160;            familyArea[i] = 0; </div><div class="line"><a name="l00685"></a><span class="lineno">  685</span>&#160;            familyVol[i] = 0;</div><div class="line"><a name="l00686"></a><span class="lineno">  686</span>&#160;        }        </div><div class="line"><a name="l00687"></a><span class="lineno">  687</span>&#160;    }    </div><div class="line"><a name="l00688"></a><span class="lineno">  688</span>&#160;    </div><div class="line"><a name="l00689"></a><span class="lineno">  689</span>&#160;</div><div class="line"><a name="l00690"></a><span class="lineno">  690</span>&#160;    <span class="comment">// Calculate total fracture area, and area per family</span></div><div class="line"><a name="l00691"></a><span class="lineno">  691</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i &lt; finalFractures.size(); i++) {</div><div class="line"><a name="l00692"></a><span class="lineno">  692</span>&#160;        <span class="keywordtype">double</span> <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a> = acceptedPoly[finalFractures[i]].area;</div><div class="line"><a name="l00693"></a><span class="lineno">  693</span>&#160;        <span class="keywordtype">double</span> vol = area * acceptedPoly[finalFractures[i]].aperture;</div><div class="line"><a name="l00694"></a><span class="lineno">  694</span>&#160;               </div><div class="line"><a name="l00695"></a><span class="lineno">  695</span>&#160;        pstats.<a class="code" href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">areaAfterRemoval</a> += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00696"></a><span class="lineno">  696</span>&#160;        pstats.<a class="code" href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">volAfterRemoval</a> += vol;</div><div class="line"><a name="l00697"></a><span class="lineno">  697</span>&#160;    </div><div class="line"><a name="l00698"></a><span class="lineno">  698</span>&#160;        <span class="keywordflow">if</span> (acceptedPoly[finalFractures[i]].<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> &gt;= 0) {</div><div class="line"><a name="l00699"></a><span class="lineno">  699</span>&#160;                familyArea[acceptedPoly[finalFractures[i]].familyNum] += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00700"></a><span class="lineno">  700</span>&#160;                familyVol[acceptedPoly[finalFractures[i]].familyNum] += vol;</div><div class="line"><a name="l00701"></a><span class="lineno">  701</span>&#160;            }</div><div class="line"><a name="l00702"></a><span class="lineno">  702</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// User-defined polygon</span></div><div class="line"><a name="l00703"></a><span class="lineno">  703</span>&#160;            userDefinedShapesArea += <a class="code" href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">area</a>;</div><div class="line"><a name="l00704"></a><span class="lineno">  704</span>&#160;            userDefinedVol += vol;</div><div class="line"><a name="l00705"></a><span class="lineno">  705</span>&#160;        }    </div><div class="line"><a name="l00706"></a><span class="lineno">  706</span>&#160;    }</div><div class="line"><a name="l00707"></a><span class="lineno">  707</span>&#160;</div><div class="line"><a name="l00708"></a><span class="lineno">  708</span>&#160;    <span class="comment">// Re-count number of accepted fracture per family after isloated fractures were removed</span></div><div class="line"><a name="l00709"></a><span class="lineno">  709</span>&#160;    <span class="keywordtype">int</span> *acceptedFromFamCounters = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00710"></a><span class="lineno">  710</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0) {</div><div class="line"><a name="l00711"></a><span class="lineno">  711</span>&#160;        acceptedFromFamCounters = <span class="keyword">new</span> <span class="keywordtype">int</span>[totalFamilies];</div><div class="line"><a name="l00712"></a><span class="lineno">  712</span>&#160;</div><div class="line"><a name="l00713"></a><span class="lineno">  713</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) {</div><div class="line"><a name="l00714"></a><span class="lineno">  714</span>&#160;            <span class="comment">// zero counters</span></div><div class="line"><a name="l00715"></a><span class="lineno">  715</span>&#160;            acceptedFromFamCounters[i] = 0;</div><div class="line"><a name="l00716"></a><span class="lineno">  716</span>&#160;        }</div><div class="line"><a name="l00717"></a><span class="lineno">  717</span>&#160;</div><div class="line"><a name="l00718"></a><span class="lineno">  718</span>&#160;        <span class="keywordtype">int</span> size = finalFractures.size();</div><div class="line"><a name="l00719"></a><span class="lineno">  719</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; size; i++) {</div><div class="line"><a name="l00720"></a><span class="lineno">  720</span>&#160;            <span class="keywordflow">if</span> (acceptedPoly[finalFractures[i]].<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> &gt;= 0) {</div><div class="line"><a name="l00721"></a><span class="lineno">  721</span>&#160;                acceptedFromFamCounters[acceptedPoly[finalFractures[i]].familyNum]++;</div><div class="line"><a name="l00722"></a><span class="lineno">  722</span>&#160;            }</div><div class="line"><a name="l00723"></a><span class="lineno">  723</span>&#160;        }</div><div class="line"><a name="l00724"></a><span class="lineno">  724</span>&#160;    }</div><div class="line"><a name="l00725"></a><span class="lineno">  725</span>&#160;</div><div class="line"><a name="l00726"></a><span class="lineno">  726</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Surface Area:     &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">areaAfterRemoval</a> * 2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00727"></a><span class="lineno">  727</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fractures Volume: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">volAfterRemoval</a> &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00728"></a><span class="lineno">  728</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Density   (P30): &quot;</span> &lt;&lt; finalFractures.size()/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00729"></a><span class="lineno">  729</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Intensity (P32): &quot;</span> &lt;&lt; (pstats.<a class="code" href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">areaAfterRemoval</a> * 2)/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00730"></a><span class="lineno">  730</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fracture Porosity  (P33): &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">volAfterRemoval</a>/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00731"></a><span class="lineno">  731</span>&#160;</div><div class="line"><a name="l00732"></a><span class="lineno">  732</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Surface Area:     &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">areaAfterRemoval</a> * 2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00733"></a><span class="lineno">  733</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fractures Volume: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">volAfterRemoval</a> &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00734"></a><span class="lineno">  734</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Density   (P30): &quot;</span> &lt;&lt; finalFractures.size()/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00735"></a><span class="lineno">  735</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Intensity (P32): &quot;</span> &lt;&lt; (pstats.<a class="code" href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">areaAfterRemoval</a> * 2)/domVol &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00736"></a><span class="lineno">  736</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fracture Porosity  (P33): &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">volAfterRemoval</a>/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00737"></a><span class="lineno">  737</span>&#160;</div><div class="line"><a name="l00738"></a><span class="lineno">  738</span>&#160;    <span class="comment">// Print family stats to user</span></div><div class="line"><a name="l00739"></a><span class="lineno">  739</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) { </div><div class="line"><a name="l00740"></a><span class="lineno">  740</span>&#160;        </div><div class="line"><a name="l00741"></a><span class="lineno">  741</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot;:\n&quot;</span>;</div><div class="line"><a name="l00742"></a><span class="lineno">  742</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Fractures After Isloated Fracture Removal: &quot;</span> &lt;&lt; acceptedFromFamCounters[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00743"></a><span class="lineno">  743</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Isolated Fractures Removed: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] - acceptedFromFamCounters[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>; </div><div class="line"><a name="l00744"></a><span class="lineno">  744</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Accepted: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00745"></a><span class="lineno">  745</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00746"></a><span class="lineno">  746</span>&#160;</div><div class="line"><a name="l00747"></a><span class="lineno">  747</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;Family &quot;</span> &lt;&lt; i+1 &lt;&lt; <span class="stringliteral">&quot;:\n&quot;</span>;</div><div class="line"><a name="l00748"></a><span class="lineno">  748</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Fractures After Isloated Fracture Removal: &quot;</span> &lt;&lt; acceptedFromFamCounters[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00749"></a><span class="lineno">  749</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Isolated Fractures Removed: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] - acceptedFromFamCounters[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>; </div><div class="line"><a name="l00750"></a><span class="lineno">  750</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Accepted: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00751"></a><span class="lineno">  751</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00752"></a><span class="lineno">  752</span>&#160;</div><div class="line"><a name="l00753"></a><span class="lineno">  753</span>&#160;        <span class="keywordflow">if</span> ( shapeFamilies[i].layer &gt; 0) {</div><div class="line"><a name="l00754"></a><span class="lineno">  754</span>&#160;            <span class="keywordtype">int</span> idx = (shapeFamilies[i].layer-1)*2;</div><div class="line"><a name="l00755"></a><span class="lineno">  755</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer: &quot;</span> &lt;&lt; shapeFamilies[i].layer &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00756"></a><span class="lineno">  756</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer {-z, +z}: {&quot;</span> &lt;&lt;<a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[idx]&lt;&lt; <span class="stringliteral">&quot;,&quot;</span> &lt;&lt;<a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[idx+1]&lt;&lt; <span class="stringliteral">&quot;}\n&quot;</span>;</div><div class="line"><a name="l00757"></a><span class="lineno">  757</span>&#160;</div><div class="line"><a name="l00758"></a><span class="lineno">  758</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer: &quot;</span> &lt;&lt; shapeFamilies[i].layer &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00759"></a><span class="lineno">  759</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer {-z, +z}: {&quot;</span> &lt;&lt;layers[idx]&lt;&lt; <span class="stringliteral">&quot;,&quot;</span> &lt;&lt;layers[idx+1]&lt;&lt; <span class="stringliteral">&quot;}\n&quot;</span>;</div><div class="line"><a name="l00760"></a><span class="lineno">  760</span>&#160;        }</div><div class="line"><a name="l00761"></a><span class="lineno">  761</span>&#160;        <span class="keywordflow">else</span> {</div><div class="line"><a name="l00762"></a><span class="lineno">  762</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;    Layer: Whole Domain \n&quot;</span>;</div><div class="line"><a name="l00763"></a><span class="lineno">  763</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;    Layer: Whole Domain \n&quot;</span>;</div><div class="line"><a name="l00764"></a><span class="lineno">  764</span>&#160;        }</div><div class="line"><a name="l00765"></a><span class="lineno">  765</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; familyArea[i]*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00766"></a><span class="lineno">  766</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; familyVol[i] &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00767"></a><span class="lineno">  767</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Fracture Intesnsity (P32): &quot;</span> &lt;&lt; familyArea[i]*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00768"></a><span class="lineno">  768</span>&#160;</div><div class="line"><a name="l00769"></a><span class="lineno">  769</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; familyArea[i]*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00770"></a><span class="lineno">  770</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; familyVol[i] &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00771"></a><span class="lineno">  771</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Fracture Intesnsity (P32): &quot;</span> &lt;&lt; familyArea[i]*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00772"></a><span class="lineno">  772</span>&#160;    }    </div><div class="line"><a name="l00773"></a><span class="lineno">  773</span>&#160;      </div><div class="line"><a name="l00774"></a><span class="lineno">  774</span>&#160;    <span class="keywordflow">if</span> (userDefinedShapesArea &gt; 0) {</div><div class="line"><a name="l00775"></a><span class="lineno">  775</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;User Defined Shapes: \n&quot;</span>;</div><div class="line"><a name="l00776"></a><span class="lineno">  776</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; userDefinedShapesArea*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00777"></a><span class="lineno">  777</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; userDefinedVol &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00778"></a><span class="lineno">  778</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; userDefinedShapesArea*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00779"></a><span class="lineno">  779</span>&#160;</div><div class="line"><a name="l00780"></a><span class="lineno">  780</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;User Defined Shapes: \n&quot;</span>;</div><div class="line"><a name="l00781"></a><span class="lineno">  781</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Surface Area: &quot;</span> &lt;&lt; userDefinedShapesArea*2 &lt;&lt; <span class="stringliteral">&quot; m^2\n&quot;</span>;</div><div class="line"><a name="l00782"></a><span class="lineno">  782</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Volume: &quot;</span> &lt;&lt; userDefinedVol &lt;&lt; <span class="stringliteral">&quot; m^3\n&quot;</span>;</div><div class="line"><a name="l00783"></a><span class="lineno">  783</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;    Fracture Intensity (P32): &quot;</span> &lt;&lt; userDefinedShapesArea*2/domVol &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00784"></a><span class="lineno">  784</span>&#160;    }</div><div class="line"><a name="l00785"></a><span class="lineno">  785</span>&#160; </div><div class="line"><a name="l00786"></a><span class="lineno">  786</span>&#160;    <span class="keywordflow">if</span> (acceptedFromFamCounters != <span class="keyword">nullptr</span>) {</div><div class="line"><a name="l00787"></a><span class="lineno">  787</span>&#160;        <span class="keyword">delete</span>[] acceptedFromFamCounters;</div><div class="line"><a name="l00788"></a><span class="lineno">  788</span>&#160;    }</div><div class="line"><a name="l00789"></a><span class="lineno">  789</span>&#160;     </div><div class="line"><a name="l00790"></a><span class="lineno">  790</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00791"></a><span class="lineno">  791</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00792"></a><span class="lineno">  792</span>&#160;    </div><div class="line"><a name="l00793"></a><span class="lineno">  793</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span> &lt;&lt; acceptedPoly.size()&lt;&lt; <span class="stringliteral">&quot; Fractures Accepted (Before Isolated Fracture Removal)\n&quot;</span>;</div><div class="line"><a name="l00794"></a><span class="lineno">  794</span>&#160;    std::cout &lt;&lt; finalFractures.size()&lt;&lt; <span class="stringliteral">&quot; Final Fractures (After Isolated Fracture Removal)\n\n&quot;</span>;    </div><div class="line"><a name="l00795"></a><span class="lineno">  795</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fractures Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ae78598e8e15f944d358722868f763e32">rejectedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;    </div><div class="line"><a name="l00796"></a><span class="lineno">  796</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Total Fractures Re-translated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a0fe2c18e0646b84364712352595751ec">retranslatedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00797"></a><span class="lineno">  797</span>&#160; </div><div class="line"><a name="l00798"></a><span class="lineno">  798</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span> &lt;&lt; acceptedPoly.size()&lt;&lt; <span class="stringliteral">&quot; Fractures Accepted (Before Isolated Fracture Removal)\n&quot;</span>;</div><div class="line"><a name="l00799"></a><span class="lineno">  799</span>&#160;    file &lt;&lt; finalFractures.size()&lt;&lt; <span class="stringliteral">&quot; Final Fractures (After Isolated Fracture Removal)\n\n&quot;</span>;    </div><div class="line"><a name="l00800"></a><span class="lineno">  800</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fractures Rejected: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ae78598e8e15f944d358722868f763e32">rejectedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;    </div><div class="line"><a name="l00801"></a><span class="lineno">  801</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Total Fractures Re-translated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a0fe2c18e0646b84364712352595751ec">retranslatedPolyCount</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00802"></a><span class="lineno">  802</span>&#160;  </div><div class="line"><a name="l00803"></a><span class="lineno">  803</span>&#160;    <span class="keywordflow">if</span> (printConnectivityError == 1) {</div><div class="line"><a name="l00804"></a><span class="lineno">  804</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\nERROR: DFN generation has finished but the formed\n&quot;</span> </div><div class="line"><a name="l00805"></a><span class="lineno">  805</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;fracture network does not make a connection between\n&quot;</span> </div><div class="line"><a name="l00806"></a><span class="lineno">  806</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;the user&#39;s specifed boundary faces.\n&quot;</span>;</div><div class="line"><a name="l00807"></a><span class="lineno">  807</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;Try increasing the fracture density, shrinking the domain\n&quot;</span></div><div class="line"><a name="l00808"></a><span class="lineno">  808</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;or consider using the &#39;ignoreBoundaryFaces&#39; option.\n&quot;</span>;</div><div class="line"><a name="l00809"></a><span class="lineno">  809</span>&#160;</div><div class="line"><a name="l00810"></a><span class="lineno">  810</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;\nERROR: DFN generation has finished but the formed\n&quot;</span> </div><div class="line"><a name="l00811"></a><span class="lineno">  811</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;fracture network does not make a connection between\n&quot;</span> </div><div class="line"><a name="l00812"></a><span class="lineno">  812</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;the user&#39;s specifed boundary faces.\n&quot;</span>;</div><div class="line"><a name="l00813"></a><span class="lineno">  813</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;Try increasing the fracture density, shrinking the domain\n&quot;</span></div><div class="line"><a name="l00814"></a><span class="lineno">  814</span>&#160;                  &lt;&lt; <span class="stringliteral">&quot;or consider using the &#39;ignoreBoundaryFaces&#39; option.\n&quot;</span>;</div><div class="line"><a name="l00815"></a><span class="lineno">  815</span>&#160;</div><div class="line"><a name="l00816"></a><span class="lineno">  816</span>&#160;        file.close();</div><div class="line"><a name="l00817"></a><span class="lineno">  817</span>&#160;        exit(1);</div><div class="line"><a name="l00818"></a><span class="lineno">  818</span>&#160;    }</div><div class="line"><a name="l00819"></a><span class="lineno">  819</span>&#160;    </div><div class="line"><a name="l00820"></a><span class="lineno">  820</span>&#160;    <span class="comment">//************ Intersection Stats ***************</span></div><div class="line"><a name="l00821"></a><span class="lineno">  821</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nNumber of Triple Intersection Points (Before Isolated Fracture Removal): &quot;</span> &lt;&lt; triplePoints.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;  </div><div class="line"><a name="l00822"></a><span class="lineno">  822</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\nNumber of Triple Intersection Points (Before Isolated Fracture Removal): &quot;</span> &lt;&lt; triplePoints.size() &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;   </div><div class="line"><a name="l00823"></a><span class="lineno">  823</span>&#160;    <span class="comment">// Shrink intersection stats</span></div><div class="line"><a name="l00824"></a><span class="lineno">  824</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nIntersection Statistics:\n&quot;</span>;</div><div class="line"><a name="l00825"></a><span class="lineno">  825</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ae27fc3d500d31275f8d97fde166b055c">intersectionsShortened</a> &lt;&lt; <span class="stringliteral">&quot; of &quot;</span> &lt;&lt; intPts.size() &lt;&lt; <span class="stringliteral">&quot; Intersections Shortened\n&quot;</span>;</div><div class="line"><a name="l00826"></a><span class="lineno">  826</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    Original Intersection (Before Intersection Shrinking) Length: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a8cf48c5b0b96005ef7dbb03cceaf9003">originalLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00827"></a><span class="lineno">  827</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    Intersection Length Discarded: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a5cd7cf21a73c6d359b63732ef1d7eebb">discardedLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00828"></a><span class="lineno">  828</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    Final Intersection Length: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a8cf48c5b0b96005ef7dbb03cceaf9003">originalLength</a> - pstats.<a class="code" href="struct_stats.html#a5cd7cf21a73c6d359b63732ef1d7eebb">discardedLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00829"></a><span class="lineno">  829</span>&#160;</div><div class="line"><a name="l00830"></a><span class="lineno">  830</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\nIntersection Statistics:\n&quot;</span>;</div><div class="line"><a name="l00831"></a><span class="lineno">  831</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#ae27fc3d500d31275f8d97fde166b055c">intersectionsShortened</a> &lt;&lt; <span class="stringliteral">&quot; of &quot;</span> &lt;&lt; intPts.size() &lt;&lt; <span class="stringliteral">&quot; Intersections Shortened\n&quot;</span>;</div><div class="line"><a name="l00832"></a><span class="lineno">  832</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    Original Intersection (Before Intersection Shrinking) Length: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a8cf48c5b0b96005ef7dbb03cceaf9003">originalLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00833"></a><span class="lineno">  833</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    Intersection Length Discarded: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a5cd7cf21a73c6d359b63732ef1d7eebb">discardedLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00834"></a><span class="lineno">  834</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    Final Intersection Length: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a8cf48c5b0b96005ef7dbb03cceaf9003">originalLength</a> - pstats.<a class="code" href="struct_stats.html#a5cd7cf21a73c6d359b63732ef1d7eebb">discardedLength</a> &lt;&lt; <span class="stringliteral">&quot;m\n&quot;</span>;</div><div class="line"><a name="l00835"></a><span class="lineno">  835</span>&#160;   </div><div class="line"><a name="l00836"></a><span class="lineno">  836</span>&#160;    <span class="comment">//*********** Rejection Stats *******************</span></div><div class="line"><a name="l00837"></a><span class="lineno">  837</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nRejection Statistics: \n&quot;</span>;  </div><div class="line"><a name="l00838"></a><span class="lineno">  838</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#ab935b00ed9bf36d47d435a3f4fbd328d">shortIntersection</a> &lt;&lt; <span class="stringliteral">&quot; Short Intersections \n&quot;</span>;</div><div class="line"><a name="l00839"></a><span class="lineno">  839</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#ac581f8584a6bc05181a18c9145f2a26e">closeToNode</a> &lt;&lt; <span class="stringliteral">&quot; Close to Node\n&quot;</span>;</div><div class="line"><a name="l00840"></a><span class="lineno">  840</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#adbb08b81bcca27734f30fc8314a48719">closeToEdge</a> &lt;&lt; <span class="stringliteral">&quot; Close to Edge\n&quot;</span>;</div><div class="line"><a name="l00841"></a><span class="lineno">  841</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a97d5e4f50262ee7f1142fa7a61c81749">closePointToEdge</a> &lt;&lt; <span class="stringliteral">&quot; Vertice Close to Edge\n&quot;</span>;</div><div class="line"><a name="l00842"></a><span class="lineno">  842</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a64f92b5fcce897497d067d072443dc60">outside</a> &lt;&lt; <span class="stringliteral">&quot; Outside of Domain\n&quot;</span>;</div><div class="line"><a name="l00843"></a><span class="lineno">  843</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a71e3bcee9eb7d5eb99f60ec8832a8754">triple</a> &lt;&lt; <span class="stringliteral">&quot; Triple intersection Rejections\n&quot;</span>;  </div><div class="line"><a name="l00844"></a><span class="lineno">  844</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a8937dddd46dbdb1731e980500c02e1f3">interCloseToInter</a> &lt;&lt; <span class="stringliteral">&quot; Intersections Close to Other Intersections\n\n&quot;</span>;</div><div class="line"><a name="l00845"></a><span class="lineno">  845</span>&#160;</div><div class="line"><a name="l00846"></a><span class="lineno">  846</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\nRejection Statistics: \n&quot;</span>;  </div><div class="line"><a name="l00847"></a><span class="lineno">  847</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#ab935b00ed9bf36d47d435a3f4fbd328d">shortIntersection</a> &lt;&lt; <span class="stringliteral">&quot; Short Intersections \n&quot;</span>;</div><div class="line"><a name="l00848"></a><span class="lineno">  848</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#ac581f8584a6bc05181a18c9145f2a26e">closeToNode</a> &lt;&lt; <span class="stringliteral">&quot; Close to Node\n&quot;</span>;</div><div class="line"><a name="l00849"></a><span class="lineno">  849</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#adbb08b81bcca27734f30fc8314a48719">closeToEdge</a> &lt;&lt; <span class="stringliteral">&quot; Close to Edge\n&quot;</span>;</div><div class="line"><a name="l00850"></a><span class="lineno">  850</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a97d5e4f50262ee7f1142fa7a61c81749">closePointToEdge</a> &lt;&lt; <span class="stringliteral">&quot; Vertice Close to Edge\n&quot;</span>;</div><div class="line"><a name="l00851"></a><span class="lineno">  851</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a64f92b5fcce897497d067d072443dc60">outside</a> &lt;&lt; <span class="stringliteral">&quot; Outside of Domain\n&quot;</span>;</div><div class="line"><a name="l00852"></a><span class="lineno">  852</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a71e3bcee9eb7d5eb99f60ec8832a8754">triple</a> &lt;&lt; <span class="stringliteral">&quot; Triple intersection Rejections\n&quot;</span>;  </div><div class="line"><a name="l00853"></a><span class="lineno">  853</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">rejectionReasons</a>.<a class="code" href="struct_rejection_reasons.html#a8937dddd46dbdb1731e980500c02e1f3">interCloseToInter</a> &lt;&lt; <span class="stringliteral">&quot; Intersections Close to Other Intersections\n\n&quot;</span>;</div><div class="line"><a name="l00854"></a><span class="lineno">  854</span>&#160;</div><div class="line"><a name="l00855"></a><span class="lineno">  855</span>&#160;</div><div class="line"><a name="l00856"></a><span class="lineno">  856</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00857"></a><span class="lineno">  857</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00858"></a><span class="lineno">  858</span>&#160;</div><div class="line"><a name="l00859"></a><span class="lineno">  859</span>&#160;    <span class="keywordflow">if</span> (totalFamilies &gt; 0) {</div><div class="line"><a name="l00860"></a><span class="lineno">  860</span>&#160;    </div><div class="line"><a name="l00861"></a><span class="lineno">  861</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;Fracture Estimation statistics:\n&quot;</span>;</div><div class="line"><a name="l00862"></a><span class="lineno">  862</span>&#160;            file &lt;&lt; <span class="stringliteral">&quot;Fracture Estimation statistics:\n&quot;</span>;</div><div class="line"><a name="l00863"></a><span class="lineno">  863</span>&#160; </div><div class="line"><a name="l00864"></a><span class="lineno">  864</span>&#160;           std::cout &lt;&lt; <span class="stringliteral">&quot;NOTE: If estimation and actual are very different, \nexpected family distributions might &quot;</span> </div><div class="line"><a name="l00865"></a><span class="lineno">  865</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;not be accurate. \nIf this is the case, try increasing or decreasing \nthe &#39;radiiListIncrease&#39; option &quot;</span></div><div class="line"><a name="l00866"></a><span class="lineno">  866</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;in the input file.\n\n&quot;</span>;</div><div class="line"><a name="l00867"></a><span class="lineno">  867</span>&#160;           file &lt;&lt; <span class="stringliteral">&quot;NOTE: If estimation and actual are very different, \nexpected family distributions might &quot;</span> </div><div class="line"><a name="l00868"></a><span class="lineno">  868</span>&#160;                &lt;&lt; <span class="stringliteral">&quot;not be accurate. \nIf this is the case, try increasing or decreasing \nthe &#39;radiiListIncrease&#39; option &quot;</span></div><div class="line"><a name="l00869"></a><span class="lineno">  869</span>&#160;                &lt;&lt; <span class="stringliteral">&quot;in the input file.\n\n&quot;</span>;</div><div class="line"><a name="l00870"></a><span class="lineno">  870</span>&#160;</div><div class="line"><a name="l00871"></a><span class="lineno">  871</span>&#160;        <span class="comment">// Compare expected radii/poly size and actual </span></div><div class="line"><a name="l00872"></a><span class="lineno">  872</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; totalFamilies; i++) {</div><div class="line"><a name="l00873"></a><span class="lineno">  873</span>&#160;            <span class="keywordflow">if</span> (shapeFamilies[i].distributionType == 4) { <span class="comment">// Constant</span></div><div class="line"><a name="l00874"></a><span class="lineno">  874</span>&#160;               std::cout &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; Family &quot;</span></div><div class="line"><a name="l00875"></a><span class="lineno">  875</span>&#160;                          &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i, shapeFamilies[i].shapeFamily) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span> </div><div class="line"><a name="l00876"></a><span class="lineno">  876</span>&#160;                          &lt;&lt; <span class="stringliteral">&quot;Using constant size\n\n&quot;</span>;</div><div class="line"><a name="l00877"></a><span class="lineno">  877</span>&#160;</div><div class="line"><a name="l00878"></a><span class="lineno">  878</span>&#160;                file &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; Family &quot;</span></div><div class="line"><a name="l00879"></a><span class="lineno">  879</span>&#160;                     &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i, shapeFamilies[i].shapeFamily) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span></div><div class="line"><a name="l00880"></a><span class="lineno">  880</span>&#160;                     &lt;&lt; <span class="stringliteral">&quot;Using constant size\n\n&quot;</span>;</div><div class="line"><a name="l00881"></a><span class="lineno">  881</span>&#160;</div><div class="line"><a name="l00882"></a><span class="lineno">  882</span>&#160;            }</div><div class="line"><a name="l00883"></a><span class="lineno">  883</span>&#160;            <span class="keywordflow">else</span> {</div><div class="line"><a name="l00884"></a><span class="lineno">  884</span>&#160;                std::cout &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; Family &quot;</span></div><div class="line"><a name="l00885"></a><span class="lineno">  885</span>&#160;                          &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i, shapeFamilies[i].shapeFamily) &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span> </div><div class="line"><a name="l00886"></a><span class="lineno">  886</span>&#160;                          &lt;&lt; <span class="stringliteral">&quot;Estimated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a1041c9121838585a366175f3ddd23e91">expectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00887"></a><span class="lineno">  887</span>&#160;                std::cout &lt;&lt; <span class="stringliteral">&quot;Actual:    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] + pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00888"></a><span class="lineno">  888</span>&#160;</div><div class="line"><a name="l00889"></a><span class="lineno">  889</span>&#160;                file &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFamilies[i]) &lt;&lt; <span class="stringliteral">&quot; Family &quot;</span></div><div class="line"><a name="l00890"></a><span class="lineno">  890</span>&#160;                     &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(i, shapeFamilies[i].shapeFamily) </div><div class="line"><a name="l00891"></a><span class="lineno">  891</span>&#160;                     &lt;&lt; <span class="stringliteral">&quot;Estimated: &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a1041c9121838585a366175f3ddd23e91">expectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00892"></a><span class="lineno">  892</span>&#160;                file &lt;&lt; <span class="stringliteral">&quot;Actual:    &quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">acceptedFromFam</a>[i] + pstats.<a class="code" href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">rejectedFromFam</a>[i] &lt;&lt; <span class="stringliteral">&quot;\n\n&quot;</span>;</div><div class="line"><a name="l00893"></a><span class="lineno">  893</span>&#160;            }</div><div class="line"><a name="l00894"></a><span class="lineno">  894</span>&#160;</div><div class="line"><a name="l00895"></a><span class="lineno">  895</span>&#160;        }</div><div class="line"><a name="l00896"></a><span class="lineno">  896</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00897"></a><span class="lineno">  897</span>&#160;        file &lt;&lt; <span class="stringliteral">&quot;\n________________________________________________________\n\n&quot;</span>;  </div><div class="line"><a name="l00898"></a><span class="lineno">  898</span>&#160;    }</div><div class="line"><a name="l00899"></a><span class="lineno">  899</span>&#160;</div><div class="line"><a name="l00900"></a><span class="lineno">  900</span>&#160;</div><div class="line"><a name="l00901"></a><span class="lineno">  901</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;Seed: &quot;</span> &lt;&lt; <a class="code" href="input_8h.html#ae21f357c223957d36046a0d71cc6aed7">seed</a> &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00902"></a><span class="lineno">  902</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;Seed: &quot;</span> &lt;&lt; seed &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00903"></a><span class="lineno">  903</span>&#160;</div><div class="line"><a name="l00904"></a><span class="lineno">  904</span>&#160;    <span class="comment">// Write all output files</span></div><div class="line"><a name="l00905"></a><span class="lineno">  905</span>&#160;    <a class="code" href="output_8cpp.html#ad84815b625587d2194321694a13f2477">writeOutput</a>(argv[2], acceptedPoly, intPts, triplePoints, pstats, finalFractures, shapeFamilies);</div><div class="line"><a name="l00906"></a><span class="lineno">  906</span>&#160;</div><div class="line"><a name="l00907"></a><span class="lineno">  907</span>&#160;    <span class="comment">// Duplicate node counters are set in writeOutput(). Write output must happen before</span></div><div class="line"><a name="l00908"></a><span class="lineno">  908</span>&#160;    <span class="comment">// duplicate node prints</span></div><div class="line"><a name="l00909"></a><span class="lineno">  909</span>&#160;    <span class="comment">// Print number of duplicate nodes (pstats.intersectionsNodeCount is set in writeOutpu() )</span></div><div class="line"><a name="l00910"></a><span class="lineno">  910</span>&#160;    std::cout &lt;&lt; <span class="stringliteral">&quot;\nLagrit Should Remove &quot;</span> </div><div class="line"><a name="l00911"></a><span class="lineno">  911</span>&#160;              &lt;&lt; pstats.<a class="code" href="struct_stats.html#a84c9fecbd7e80b7cf93b185952618062">intersectionNodeCount</a>/2 - pstats.<a class="code" href="struct_stats.html#a6d521913f1a280406a6536cb632c4356">tripleNodeCount</a> </div><div class="line"><a name="l00912"></a><span class="lineno">  912</span>&#160;              &lt;&lt; <span class="stringliteral">&quot; Nodes (&quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a84c9fecbd7e80b7cf93b185952618062">intersectionNodeCount</a> &lt;&lt; <span class="stringliteral">&quot;/2 - &quot;</span> </div><div class="line"><a name="l00913"></a><span class="lineno">  913</span>&#160;              &lt;&lt; pstats.<a class="code" href="struct_stats.html#a6d521913f1a280406a6536cb632c4356">tripleNodeCount</a> &lt;&lt; <span class="stringliteral">&quot;)\n&quot;</span>;</div><div class="line"><a name="l00914"></a><span class="lineno">  914</span>&#160;</div><div class="line"><a name="l00915"></a><span class="lineno">  915</span>&#160;    file &lt;&lt; <span class="stringliteral">&quot;\nLagrit Should Remove &quot;</span> </div><div class="line"><a name="l00916"></a><span class="lineno">  916</span>&#160;              &lt;&lt; pstats.<a class="code" href="struct_stats.html#a84c9fecbd7e80b7cf93b185952618062">intersectionNodeCount</a>/2 - pstats.<a class="code" href="struct_stats.html#a6d521913f1a280406a6536cb632c4356">tripleNodeCount</a></div><div class="line"><a name="l00917"></a><span class="lineno">  917</span>&#160;              &lt;&lt; <span class="stringliteral">&quot; Nodes (&quot;</span> &lt;&lt; pstats.<a class="code" href="struct_stats.html#a84c9fecbd7e80b7cf93b185952618062">intersectionNodeCount</a> &lt;&lt; <span class="stringliteral">&quot;/2 - &quot;</span> </div><div class="line"><a name="l00918"></a><span class="lineno">  918</span>&#160;              &lt;&lt; pstats.<a class="code" href="struct_stats.html#a6d521913f1a280406a6536cb632c4356">tripleNodeCount</a> &lt;&lt; <span class="stringliteral">&quot;)\n&quot;</span>;</div><div class="line"><a name="l00919"></a><span class="lineno">  919</span>&#160;</div><div class="line"><a name="l00920"></a><span class="lineno">  920</span>&#160;    file.close();</div><div class="line"><a name="l00921"></a><span class="lineno">  921</span>&#160;    <span class="keywordflow">return</span> 0;</div><div class="line"><a name="l00922"></a><span class="lineno">  922</span>&#160;}</div><div class="line"><a name="l00923"></a><span class="lineno">  923</span>&#160;</div><div class="line"><a name="l00924"></a><span class="lineno">  924</span>&#160;<span class="comment">/******************************** END MAIN ***********************************/</span></div><div class="line"><a name="l00925"></a><span class="lineno">  925</span>&#160;<span class="comment">/*****************************************************************************/</span></div><div class="line"><a name="l00926"></a><span class="lineno">  926</span>&#160;</div><div class="line"><a name="l00927"></a><span class="lineno">  927</span>&#160;</div><div class="ttc" id="hotkey_8cpp_html_ac8373d63f264c2e956b071b5bd2588a1"><div class="ttname"><a href="hotkey_8cpp.html#ac8373d63f264c2e956b071b5bd2588a1">reset_terminal_mode</a></div><div class="ttdeci">void reset_terminal_mode()</div><div class="ttdef"><b>Definition:</b> <a href="hotkey_8cpp_source.html#l00012">hotkey.cpp:12</a></div></div>
<div class="ttc" id="computational_geometry_8h_html"><div class="ttname"><a href="computational_geometry_8h.html">computationalGeometry.h</a></div></div>
<div class="ttc" id="distributions_8h_html"><div class="ttname"><a href="distributions_8h.html">distributions.h</a></div></div>
<div class="ttc" id="output_8h_html"><div class="ttname"><a href="output_8h.html">output.h</a></div></div>
<div class="ttc" id="math_functions_8cpp_html_a55646cd87f15106134e9cb832492e1f6"><div class="ttname"><a href="math_functions_8cpp.html#a55646cd87f15106134e9cb832492e1f6">adjustCDF_and_famProb</a></div><div class="ttdeci">void adjustCDF_and_famProb(float *&amp;CDF, float *&amp;famProbability, int &amp;cdfSize, int idx2Remove)</div><div class="ttdef"><b>Definition:</b> <a href="math_functions_8cpp_source.html#l00269">mathFunctions.cpp:269</a></div></div>
<div class="ttc" id="input_8h_html"><div class="ttname"><a href="input_8h.html">input.h</a></div></div>
<div class="ttc" id="fracture_estimating_8cpp_html_a7cf21b594fd24b2879af467e23bdac0a"><div class="ttname"><a href="fracture_estimating_8cpp.html#a7cf21b594fd24b2879af467e23bdac0a">dryRun</a></div><div class="ttdeci">void dryRun(std::vector&lt; Shape &gt; &amp;shapeFamilies, float *shapeProb, std::mt19937_64 &amp;generator, Distributions &amp;distributions)</div><div class="ttdef"><b>Definition:</b> <a href="fracture_estimating_8cpp_source.html#l00174">fractureEstimating.cpp:174</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_a262d50919e2ca1c4e8b32ec027a52b57"><div class="ttname"><a href="insert_shape_8cpp.html#a262d50919e2ca1c4e8b32ec027a52b57">p32Complete</a></div><div class="ttdeci">bool p32Complete(int size)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00598">insertShape.cpp:598</a></div></div>
<div class="ttc" id="insert_shape_8h_html_a244f461634679dc91d08da3e183f9b8c"><div class="ttname"><a href="insert_shape_8h.html#a244f461634679dc91d08da3e183f9b8c">insertUserRectsByCoord</a></div><div class="ttdeci">void insertUserRectsByCoord(std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intpts, struct Stats &amp;pstats, std::vector&lt; Point &gt; &amp;triplePoints)</div><div class="ttdef"><b>Definition:</b> <a href="insert_user_rects_by_coord_8cpp_source.html#l00021">insertUserRectsByCoord.cpp:21</a></div></div>
<div class="ttc" id="struct_poly_html_a341a2ec25f71c2f58c1af119b3990387"><div class="ttname"><a href="struct_poly.html#a341a2ec25f71c2f58c1af119b3990387">Poly::truncated</a></div><div class="ttdeci">bool truncated</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00095">structures.h:95</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_ab935b00ed9bf36d47d435a3f4fbd328d"><div class="ttname"><a href="struct_rejection_reasons.html#ab935b00ed9bf36d47d435a3f4fbd328d">RejectionReasons::shortIntersection</a></div><div class="ttdeci">unsigned long long int shortIntersection</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00252">structures.h:252</a></div></div>
<div class="ttc" id="struct_stats_html_ab9e2ed63a3c7d89d4c2f128c918f3c24"><div class="ttname"><a href="struct_stats.html#ab9e2ed63a3c7d89d4c2f128c918f3c24">Stats::areaAfterRemoval</a></div><div class="ttdeci">double areaAfterRemoval</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00324">structures.h:324</a></div></div>
<div class="ttc" id="struct_poly_html"><div class="ttname"><a href="struct_poly.html">Poly</a></div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00012">structures.h:12</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_ac581f8584a6bc05181a18c9145f2a26e"><div class="ttname"><a href="struct_rejection_reasons.html#ac581f8584a6bc05181a18c9145f2a26e">RejectionReasons::closeToNode</a></div><div class="ttdeci">unsigned long long int closeToNode</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00255">structures.h:255</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_adbb08b81bcca27734f30fc8314a48719"><div class="ttname"><a href="struct_rejection_reasons.html#adbb08b81bcca27734f30fc8314a48719">RejectionReasons::closeToEdge</a></div><div class="ttdeci">unsigned long long int closeToEdge</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00258">structures.h:258</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_a71e3bcee9eb7d5eb99f60ec8832a8754"><div class="ttname"><a href="struct_rejection_reasons.html#a71e3bcee9eb7d5eb99f60ec8832a8754">RejectionReasons::triple</a></div><div class="ttdeci">unsigned long long int triple</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00264">structures.h:264</a></div></div>
<div class="ttc" id="hotkey_8cpp_html_a49d261c1ea17ddedcc836585c758e3b2"><div class="ttname"><a href="hotkey_8cpp.html#a49d261c1ea17ddedcc836585c758e3b2">set_conio_terminal_mode</a></div><div class="ttdeci">void set_conio_terminal_mode()</div><div class="ttdef"><b>Definition:</b> <a href="hotkey_8cpp_source.html#l00023">hotkey.cpp:23</a></div></div>
<div class="ttc" id="struct_stats_html_ae27fc3d500d31275f8d97fde166b055c"><div class="ttname"><a href="struct_stats.html#ae27fc3d500d31275f8d97fde166b055c">Stats::intersectionsShortened</a></div><div class="ttdeci">unsigned int intersectionsShortened</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00344">structures.h:344</a></div></div>
<div class="ttc" id="input_8h_html_af336784d384b9e2ab2d990f7855ecdcf"><div class="ttname"><a href="input_8h.html#af336784d384b9e2ab2d990f7855ecdcf">nPoly</a></div><div class="ttdeci">unsigned int nPoly</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00015">readInput.cpp:15</a></div></div>
<div class="ttc" id="struct_stats_html_a8cf48c5b0b96005ef7dbb03cceaf9003"><div class="ttname"><a href="struct_stats.html#a8cf48c5b0b96005ef7dbb03cceaf9003">Stats::originalLength</a></div><div class="ttdeci">double originalLength</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00347">structures.h:347</a></div></div>
<div class="ttc" id="_d_f_nmain_8cpp_html_a974b448c3d0a13d483d36ef15de369b5"><div class="ttname"><a href="_d_f_nmain_8cpp.html#a974b448c3d0a13d483d36ef15de369b5">eps</a></div><div class="ttdeci">double eps</div><div class="ttdef"><b>Definition:</b> <a href="_d_f_nmain_8cpp_source.html#l00039">DFNmain.cpp:39</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_a8937dddd46dbdb1731e980500c02e1f3"><div class="ttname"><a href="struct_rejection_reasons.html#a8937dddd46dbdb1731e980500c02e1f3">RejectionReasons::interCloseToInter</a></div><div class="ttdeci">unsigned long long int interCloseToInter</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00267">structures.h:267</a></div></div>
<div class="ttc" id="input_8h_html_a7ece519c4cc6ecd5512562594084e8aa"><div class="ttname"><a href="input_8h.html#a7ece519c4cc6ecd5512562594084e8aa">userRecByCoord</a></div><div class="ttdeci">bool userRecByCoord</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00340">readInput.cpp:340</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_ae4fa66ac2a60ae519ee05168b44c9917"><div class="ttname"><a href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a></div><div class="ttdeci">std::string shapeType(struct Shape &amp;shapeFam)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00700">insertShape.cpp:700</a></div></div>
<div class="ttc" id="structures_8h_html"><div class="ttname"><a href="structures_8h.html">structures.h</a></div></div>
<div class="ttc" id="read_input_8cpp_html_a9af5f928fea3353fc4009118630b3f52"><div class="ttname"><a href="read_input_8cpp.html#a9af5f928fea3353fc4009118630b3f52">getInput</a></div><div class="ttdeci">void getInput(char *input, std::vector&lt; Shape &gt; &amp;shapeFamily)</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00485">readInput.cpp:485</a></div></div>
<div class="ttc" id="input_8h_html_ab963e6ee0512d7651d002947e0bfd1df"><div class="ttname"><a href="input_8h.html#ab963e6ee0512d7651d002947e0bfd1df">userRectanglesOnOff</a></div><div class="ttdeci">bool userRectanglesOnOff</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00336">readInput.cpp:336</a></div></div>
<div class="ttc" id="input_8h_html_a6fc4e503b3c7fbcfbb51832ae2201e70"><div class="ttname"><a href="input_8h.html#a6fc4e503b3c7fbcfbb51832ae2201e70">outputAllRadii</a></div><div class="ttdeci">bool outputAllRadii</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00357">readInput.cpp:357</a></div></div>
<div class="ttc" id="math_functions_8cpp_html_a78276693aa55fab858a41366e5b4aa73"><div class="ttname"><a href="math_functions_8cpp.html#a78276693aa55fab858a41366e5b4aa73">getArea</a></div><div class="ttdeci">double getArea(struct Poly &amp;poly)</div><div class="ttdef"><b>Definition:</b> <a href="math_functions_8cpp_source.html#l00110">mathFunctions.cpp:110</a></div></div>
<div class="ttc" id="struct_stats_html_a0467673e7585331b1f609c9f7fc1de7a"><div class="ttname"><a href="struct_stats.html#a0467673e7585331b1f609c9f7fc1de7a">Stats::truncated</a></div><div class="ttdeci">unsigned int truncated</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00311">structures.h:311</a></div></div>
<div class="ttc" id="math_functions_8cpp_html_a9422106ad945a110b0cf696050962556"><div class="ttname"><a href="math_functions_8cpp.html#a9422106ad945a110b0cf696050962556">createCDF</a></div><div class="ttdeci">float * createCDF(float *famProb, int size)</div><div class="ttdef"><b>Definition:</b> <a href="math_functions_8cpp_source.html#l00244">mathFunctions.cpp:244</a></div></div>
<div class="ttc" id="input_8h_html_ae21f357c223957d36046a0d71cc6aed7"><div class="ttname"><a href="input_8h.html#ae21f357c223957d36046a0d71cc6aed7">seed</a></div><div class="ttdeci">unsigned int seed</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00105">readInput.cpp:105</a></div></div>
<div class="ttc" id="struct_stats_html_a3c0ec96761b92e70bbb96de0df58cce5"><div class="ttname"><a href="struct_stats.html#a3c0ec96761b92e70bbb96de0df58cce5">Stats::rejectionReasons</a></div><div class="ttdeci">struct RejectionReasons rejectionReasons</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00339">structures.h:339</a></div></div>
<div class="ttc" id="input_8h_html_a03ac7830eb87def8d67cc0c95b24ec82"><div class="ttname"><a href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a></div><div class="ttdeci">float * layers</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00443">readInput.cpp:443</a></div></div>
<div class="ttc" id="hotkey_8cpp_html_a160f6d2893bb6b89ab3ad5d863c20e3d"><div class="ttname"><a href="hotkey_8cpp.html#a160f6d2893bb6b89ab3ad5d863c20e3d">getch</a></div><div class="ttdeci">int getch()</div><div class="ttdef"><b>Definition:</b> <a href="hotkey_8cpp_source.html#l00057">hotkey.cpp:57</a></div></div>
<div class="ttc" id="hotkey_8h_html"><div class="ttname"><a href="hotkey_8h.html">hotkey.h</a></div></div>
<div class="ttc" id="math_functions_8h_html"><div class="ttname"><a href="math_functions_8h.html">mathFunctions.h</a></div></div>
<div class="ttc" id="input_8h_html_af93cb3f9ec0bb2490e0c2219b15aa1cf"><div class="ttname"><a href="input_8h.html#af93cb3f9ec0bb2490e0c2219b15aa1cf">userEllipsesOnOff</a></div><div class="ttdeci">bool userEllipsesOnOff</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00306">readInput.cpp:306</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_a97d5e4f50262ee7f1142fa7a61c81749"><div class="ttname"><a href="struct_rejection_reasons.html#a97d5e4f50262ee7f1142fa7a61c81749">RejectionReasons::closePointToEdge</a></div><div class="ttdeci">unsigned long long int closePointToEdge</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00260">structures.h:260</a></div></div>
<div class="ttc" id="struct_stats_html_a84c9fecbd7e80b7cf93b185952618062"><div class="ttname"><a href="struct_stats.html#a84c9fecbd7e80b7cf93b185952618062">Stats::intersectionNodeCount</a></div><div class="ttdeci">unsigned int intersectionNodeCount</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00358">structures.h:358</a></div></div>
<div class="ttc" id="struct_stats_html_a1041c9121838585a366175f3ddd23e91"><div class="ttname"><a href="struct_stats.html#a1041c9121838585a366175f3ddd23e91">Stats::expectedFromFam</a></div><div class="ttdeci">int * expectedFromFam</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00296">structures.h:296</a></div></div>
<div class="ttc" id="input_8h_html_a92d138143c536b0bc8ca70ddc1479677"><div class="ttname"><a href="input_8h.html#a92d138143c536b0bc8ca70ddc1479677">insertUserRectanglesFirst</a></div><div class="ttdeci">bool insertUserRectanglesFirst</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00096">readInput.cpp:96</a></div></div>
<div class="ttc" id="struct_stats_html_a5cd7cf21a73c6d359b63732ef1d7eebb"><div class="ttname"><a href="struct_stats.html#a5cd7cf21a73c6d359b63732ef1d7eebb">Stats::discardedLength</a></div><div class="ttdeci">double discardedLength</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00352">structures.h:352</a></div></div>
<div class="ttc" id="insert_shape_8h_html_a49e6b0be012ed44c6ef83ea66ed6dd86"><div class="ttname"><a href="insert_shape_8h.html#a49e6b0be012ed44c6ef83ea66ed6dd86">insertUserEllByCoord</a></div><div class="ttdeci">void insertUserEllByCoord(std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intpts, struct Stats &amp;pstats, std::vector&lt; Point &gt; &amp;triplePoints)</div><div class="ttdef"><b>Definition:</b> <a href="insert_user_ell_by_coord_8cpp_source.html#l00021">insertUserEllByCoord.cpp:21</a></div></div>
<div class="ttc" id="input_8h_html_ad29409829ae78f5f928f417a2a7b265b"><div class="ttname"><a href="input_8h.html#ad29409829ae78f5f928f417a2a7b265b">rejectsPerFracture</a></div><div class="ttdeci">int rejectsPerFracture</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00439">readInput.cpp:439</a></div></div>
<div class="ttc" id="struct_rejection_reasons_html_a64f92b5fcce897497d067d072443dc60"><div class="ttname"><a href="struct_rejection_reasons.html#a64f92b5fcce897497d067d072443dc60">RejectionReasons::outside</a></div><div class="ttdeci">unsigned long long int outside</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00262">structures.h:262</a></div></div>
<div class="ttc" id="input_8h_html_abd9fdedaba98f63d68fe3bb9d9098f63"><div class="ttname"><a href="input_8h.html#abd9fdedaba98f63d68fe3bb9d9098f63">radiiListIncrease</a></div><div class="ttdeci">float radiiListIncrease</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00026">readInput.cpp:26</a></div></div>
<div class="ttc" id="cluster_groups_8cpp_html_aa3f6bef1533919464df1f55a65fcbad3"><div class="ttname"><a href="cluster_groups_8cpp.html#aa3f6bef1533919464df1f55a65fcbad3">getCluster</a></div><div class="ttdeci">std::vector&lt; unsigned int &gt; getCluster(Stats &amp;pstats)</div><div class="ttdef"><b>Definition:</b> <a href="cluster_groups_8cpp_source.html#l00045">clusterGroups.cpp:45</a></div></div>
<div class="ttc" id="struct_stats_html_adf7f38474bfb4c7ab6e35877c5c174bb"><div class="ttname"><a href="struct_stats.html#adf7f38474bfb4c7ab6e35877c5c174bb">Stats::areaBeforeRemoval</a></div><div class="ttdeci">double areaBeforeRemoval</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00320">structures.h:320</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_a9b3f632508ac90991bd6ebee7bcfbc2f"><div class="ttname"><a href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a></div><div class="ttdeci">int getFamilyNumber(int familyIndex, int familyShape)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00688">insertShape.cpp:688</a></div></div>
<div class="ttc" id="struct_stats_html_ae2c4ebee5152b67185fc4932fa6cca22"><div class="ttname"><a href="struct_stats.html#ae2c4ebee5152b67185fc4932fa6cca22">Stats::rejectsPerAttempt</a></div><div class="ttdeci">std::vector&lt; unsigned int &gt; rejectsPerAttempt</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00371">structures.h:371</a></div></div>
<div class="ttc" id="fracture_estimating_8cpp_html_a26080c2634b53083b47b763b66bd4911"><div class="ttname"><a href="fracture_estimating_8cpp.html#a26080c2634b53083b47b763b66bd4911">addRadiiToLists</a></div><div class="ttdeci">void addRadiiToLists(float percent, std::vector&lt; Shape &gt; &amp;shapeFamilies, std::mt19937_64 &amp;generator, Distributions &amp;distributions)</div><div class="ttdef"><b>Definition:</b> <a href="fracture_estimating_8cpp_source.html#l00082">fractureEstimating.cpp:82</a></div></div>
<div class="ttc" id="generating_points_8h_html"><div class="ttname"><a href="generating_points_8h.html">generatingPoints.h</a></div></div>
<div class="ttc" id="insert_shape_8h_html_af69365c2e0148366aabe6751ffc54a79"><div class="ttname"><a href="insert_shape_8h.html#af69365c2e0148366aabe6751ffc54a79">insertUserRects</a></div><div class="ttdeci">void insertUserRects(std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intpts, struct Stats &amp;pstats, std::vector&lt; Point &gt; &amp;triplePoints)</div><div class="ttdef"><b>Definition:</b> <a href="insert_user_rects_8cpp_source.html#l00022">insertUserRects.cpp:22</a></div></div>
<div class="ttc" id="output_8cpp_html_ad84815b625587d2194321694a13f2477"><div class="ttname"><a href="output_8cpp.html#ad84815b625587d2194321694a13f2477">writeOutput</a></div><div class="ttdeci">void writeOutput(char *outputFolder, std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intPts, std::vector&lt; Point &gt; &amp;triplePoints, struct Stats &amp;pstats, std::vector&lt; unsigned int &gt; &amp;finalFractures, std::vector&lt; Shape &gt; &amp;shapeFamilies)</div><div class="ttdef"><b>Definition:</b> <a href="output_8cpp_source.html#l00031">output.cpp:31</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_a19e9a86f45510bc87d7fe84a11caffcf"><div class="ttname"><a href="insert_shape_8cpp.html#a19e9a86f45510bc87d7fe84a11caffcf">assignAperture</a></div><div class="ttdeci">void assignAperture(struct Poly &amp;newPoly, std::mt19937_64 &amp;generator)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00330">insertShape.cpp:330</a></div></div>
<div class="ttc" id="_d_f_nmain_8cpp_html_a3c04138a5bfe5d72780bb7e82a18e627"><div class="ttname"><a href="_d_f_nmain_8cpp.html#a3c04138a5bfe5d72780bb7e82a18e627">main</a></div><div class="ttdeci">int main(int argc, char **argv)</div><div class="ttdef"><b>Definition:</b> <a href="_d_f_nmain_8cpp_source.html#l00042">DFNmain.cpp:42</a></div></div>
<div class="ttc" id="struct_stats_html_a04e3ecdf0010372b9dae313b7dc6ecc3"><div class="ttname"><a href="struct_stats.html#a04e3ecdf0010372b9dae313b7dc6ecc3">Stats::volBeforeRemoval</a></div><div class="ttdeci">double volBeforeRemoval</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00328">structures.h:328</a></div></div>
<div class="ttc" id="struct_poly_html_a2cb519bcf0553215de907c8fc878a255"><div class="ttname"><a href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">Poly::vertices</a></div><div class="ttdeci">double * vertices</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00070">structures.h:70</a></div></div>
<div class="ttc" id="struct_stats_html_a334fa0009d90e69ec0013a17ff81e7f3"><div class="ttname"><a href="struct_stats.html#a334fa0009d90e69ec0013a17ff81e7f3">Stats::acceptedFromFam</a></div><div class="ttdeci">int * acceptedFromFam</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00283">structures.h:283</a></div></div>
<div class="ttc" id="insert_shape_8h_html"><div class="ttname"><a href="insert_shape_8h.html">insertShape.h</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_ab66dd8460f2f45feae57425d7c219bf5"><div class="ttname"><a href="insert_shape_8cpp.html#ab66dd8460f2f45feae57425d7c219bf5">assignPermeability</a></div><div class="ttdeci">void assignPermeability(struct Poly &amp;newPoly)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00380">insertShape.cpp:380</a></div></div>
<div class="ttc" id="struct_stats_html_a7e6c4e52e898eb48d5df350955dda7a8"><div class="ttname"><a href="struct_stats.html#a7e6c4e52e898eb48d5df350955dda7a8">Stats::volAfterRemoval</a></div><div class="ttdeci">double volAfterRemoval</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00332">structures.h:332</a></div></div>
<div class="ttc" id="insert_shape_8h_html_a9841b356b0460e297748d6648c301558"><div class="ttname"><a href="insert_shape_8h.html#a9841b356b0460e297748d6648c301558">insertUserEll</a></div><div class="ttdeci">void insertUserEll(std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intpts, struct Stats &amp;pstats, std::vector&lt; Point &gt; &amp;triplePoints)</div><div class="ttdef"><b>Definition:</b> <a href="insert_user_ell_8cpp_source.html#l00021">insertUserEll.cpp:21</a></div></div>
<div class="ttc" id="debug_functions_8cpp_html_a77317c191140e58c1ca12f2badadbc8b"><div class="ttname"><a href="debug_functions_8cpp.html#a77317c191140e58c1ca12f2badadbc8b">printPolyData</a></div><div class="ttdeci">void printPolyData(struct Poly &amp;poly)</div><div class="ttdef"><b>Definition:</b> <a href="debug_functions_8cpp_source.html#l00079">debugFunctions.cpp:79</a></div></div>
<div class="ttc" id="struct_poly_html_ad1acec5af38b346429eada06f05078c6"><div class="ttname"><a href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">Poly::yradius</a></div><div class="ttdeci">double yradius</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00046">structures.h:46</a></div></div>
<div class="ttc" id="_d_f_nmain_8cpp_html_a1705dee0e4e822859804bec377506fb5"><div class="ttname"><a href="_d_f_nmain_8cpp.html#a1705dee0e4e822859804bec377506fb5">orig_termios</a></div><div class="ttdeci">struct termios orig_termios</div><div class="ttdef"><b>Definition:</b> <a href="_d_f_nmain_8cpp_source.html#l00036">DFNmain.cpp:36</a></div></div>
<div class="ttc" id="computational_geometry_8cpp_html_ace23b489a120bc6d6407f7e08e1e009d"><div class="ttname"><a href="computational_geometry_8cpp.html#ace23b489a120bc6d6407f7e08e1e009d">createBoundingBox</a></div><div class="ttdeci">void createBoundingBox(struct Poly &amp;newPoly)</div><div class="ttdef"><b>Definition:</b> <a href="computational_geometry_8cpp_source.html#l00322">computationalGeometry.cpp:322</a></div></div>
<div class="ttc" id="struct_stats_html_a0fe2c18e0646b84364712352595751ec"><div class="ttname"><a href="struct_stats.html#a0fe2c18e0646b84364712352595751ec">Stats::retranslatedPolyCount</a></div><div class="ttdeci">unsigned int retranslatedPolyCount</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00308">structures.h:308</a></div></div>
<div class="ttc" id="fracture_estimating_8cpp_html_acb37fb0878da2fa7f0648216fb4189e4"><div class="ttname"><a href="fracture_estimating_8cpp.html#acb37fb0878da2fa7f0648216fb4189e4">generateRadiiLists_nPolyOption</a></div><div class="ttdeci">void generateRadiiLists_nPolyOption(std::vector&lt; Shape &gt; &amp;shapeFamilies, float *famProb, std::mt19937_64 &amp;generator, Distributions &amp;distributions)</div><div class="ttdef"><b>Definition:</b> <a href="fracture_estimating_8cpp_source.html#l00031">fractureEstimating.cpp:31</a></div></div>
<div class="ttc" id="input_8h_html_aaf09f7cf249c6fd6247821bc01ead90c"><div class="ttname"><a href="input_8h.html#aaf09f7cf249c6fd6247821bc01ead90c">nFamEll</a></div><div class="ttdeci">int nFamEll</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00120">readInput.cpp:120</a></div></div>
<div class="ttc" id="hotkey_8cpp_html_ad5451da499ab9d3907da8dd7060ab677"><div class="ttname"><a href="hotkey_8cpp.html#ad5451da499ab9d3907da8dd7060ab677">kbhit</a></div><div class="ttdeci">int kbhit()</div><div class="ttdef"><b>Definition:</b> <a href="hotkey_8cpp_source.html#l00044">hotkey.cpp:44</a></div></div>
<div class="ttc" id="input_8h_html_aeb1032cfdc6f748b280a4a0a59375d11"><div class="ttname"><a href="input_8h.html#aeb1032cfdc6f748b280a4a0a59375d11">famProb</a></div><div class="ttdeci">float * famProb</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00133">readInput.cpp:133</a></div></div>
<div class="ttc" id="struct_stats_html_a6d521913f1a280406a6536cb632c4356"><div class="ttname"><a href="struct_stats.html#a6d521913f1a280406a6536cb632c4356">Stats::tripleNodeCount</a></div><div class="ttdeci">unsigned int tripleNodeCount</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00364">structures.h:364</a></div></div>
<div class="ttc" id="testing_8h_html"><div class="ttname"><a href="testing_8h.html">testing.h</a></div></div>
<div class="ttc" id="domain_8cpp_html_aed3d77c17f33d2024704b479d0885332"><div class="ttname"><a href="domain_8cpp.html#aed3d77c17f33d2024704b479d0885332">domainTruncation</a></div><div class="ttdeci">bool domainTruncation(Poly &amp;newPoly, double *domainSize)</div><div class="ttdef"><b>Definition:</b> <a href="domain_8cpp_source.html#l00019">domain.cpp:19</a></div></div>
<div class="ttc" id="struct_poly_html_a2e10824b90d60b54522b7ce475ed2a17"><div class="ttname"><a href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">Poly::xradius</a></div><div class="ttdeci">double xradius</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00040">structures.h:40</a></div></div>
<div class="ttc" id="math_functions_8cpp_html_a5eae1db9cbad8e807a0d677bcd805d36"><div class="ttname"><a href="math_functions_8cpp.html#a5eae1db9cbad8e807a0d677bcd805d36">indexFromProb</a></div><div class="ttdeci">int indexFromProb(float *CDF, double roll, int size)</div><div class="ttdef"><b>Definition:</b> <a href="math_functions_8cpp_source.html#l00166">mathFunctions.cpp:166</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_afa95ede56e6ecb8918b61c5a2815d40d"><div class="ttname"><a href="insert_shape_8cpp.html#afa95ede56e6ecb8918b61c5a2815d40d">reTranslatePoly</a></div><div class="ttdeci">void reTranslatePoly(struct Poly &amp;newPoly, struct Shape &amp;shapeFam, std::mt19937_64 &amp;generator)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00455">insertShape.cpp:455</a></div></div>
<div class="ttc" id="struct_stats_html"><div class="ttname"><a href="struct_stats.html">Stats</a></div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00277">structures.h:277</a></div></div>
<div class="ttc" id="input_8h_html_a71b31a03f80cb26d46b69c9e96fe032b"><div class="ttname"><a href="input_8h.html#a71b31a03f80cb26d46b69c9e96fe032b">nFamRect</a></div><div class="ttdeci">int nFamRect</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00116">readInput.cpp:116</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_a4f91304a14fefba4d4f93f2f3e54dee4"><div class="ttname"><a href="insert_shape_8cpp.html#a4f91304a14fefba4d4f93f2f3e54dee4">generatePoly</a></div><div class="ttdeci">struct Poly generatePoly(struct Shape &amp;shapeFam, std::mt19937_64 &amp;generator, Distributions &amp;distributions, int familyIndex, bool useList)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00026">insertShape.cpp:26</a></div></div>
<div class="ttc" id="input_8h_html_a9b979283d8d40a1417171c28a1b1d23f"><div class="ttname"><a href="input_8h.html#a9b979283d8d40a1417171c28a1b1d23f">printRejectReasons</a></div><div class="ttdeci">bool printRejectReasons</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00067">readInput.cpp:67</a></div></div>
<div class="ttc" id="struct_stats_html_aa8289ad2f4e3480362ff9d9f13adb1a2"><div class="ttname"><a href="struct_stats.html#aa8289ad2f4e3480362ff9d9f13adb1a2">Stats::rejectedFromFam</a></div><div class="ttdeci">int * rejectedFromFam</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00289">structures.h:289</a></div></div>
<div class="ttc" id="domain_8h_html"><div class="ttname"><a href="domain_8h.html">domain.h</a></div></div>
<div class="ttc" id="input_8h_html_a9d32f63bd011a629cf644c120e07eda8"><div class="ttname"><a href="input_8h.html#a9d32f63bd011a629cf644c120e07eda8">p32Status</a></div><div class="ttdeci">bool * p32Status</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00468">readInput.cpp:468</a></div></div>
<div class="ttc" id="class_distributions_html"><div class="ttname"><a href="class_distributions.html">Distributions</a></div><div class="ttdef"><b>Definition:</b> <a href="distributions_8h_source.html#l00015">distributions.h:15</a></div></div>
<div class="ttc" id="remove_fractures_8cpp_html_a016a4a3e0104523549279e6f6cf6c6d8"><div class="ttname"><a href="remove_fractures_8cpp.html#a016a4a3e0104523549279e6f6cf6c6d8">removeFractures</a></div><div class="ttdeci">void removeFractures(double minSize, std::vector&lt; Poly &gt; &amp;acceptedPolys, std::vector&lt; IntPoints &gt; &amp;intPts, std::vector&lt; Point &gt; triplePoints, Stats &amp;pstats)</div><div class="ttdef"><b>Definition:</b> <a href="remove_fractures_8cpp_source.html#l00031">removeFractures.cpp:31</a></div></div>
<div class="ttc" id="input_8h_html_a3718da9decab4e36a1011af51ab11bd7"><div class="ttname"><a href="input_8h.html#a3718da9decab4e36a1011af51ab11bd7">userEllByCoord</a></div><div class="ttdeci">bool userEllByCoord</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00344">readInput.cpp:344</a></div></div>
<div class="ttc" id="input_8h_html_ada6852d4da13aa222e5b959bf71e89d4"><div class="ttname"><a href="input_8h.html#ada6852d4da13aa222e5b959bf71e89d4">ignoreBoundaryFaces</a></div><div class="ttdeci">bool ignoreBoundaryFaces</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00473">readInput.cpp:473</a></div></div>
<div class="ttc" id="struct_poly_html_a26a509e53f5fe73722b1d82443ce1841"><div class="ttname"><a href="struct_poly.html#a26a509e53f5fe73722b1d82443ce1841">Poly::area</a></div><div class="ttdeci">float area</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00033">structures.h:33</a></div></div>
<div class="ttc" id="computational_geometry_8cpp_html_a38a8c322dba7d8518a4c3206cd55a5f6"><div class="ttname"><a href="computational_geometry_8cpp.html#a38a8c322dba7d8518a4c3206cd55a5f6">intersectionChecking</a></div><div class="ttdeci">int intersectionChecking(struct Poly &amp;newPoly, std::vector&lt; Poly &gt; &amp;acceptedPoly, std::vector&lt; IntPoints &gt; &amp;intPtsList, struct Stats &amp;pstats, std::vector&lt; Point &gt; &amp;triplePoints)</div><div class="ttdef"><b>Definition:</b> <a href="computational_geometry_8cpp_source.html#l00870">computationalGeometry.cpp:870</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_ae9670e906ee5dc43ce3fbc0ce7290dc7"><div class="ttname"><a href="insert_shape_8cpp.html#ae9670e906ee5dc43ce3fbc0ce7290dc7">printRejectReason</a></div><div class="ttdeci">void printRejectReason(int rejectCode, struct Poly newPoly)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00617">insertShape.cpp:617</a></div></div>
<div class="ttc" id="fracture_estimating_8h_html"><div class="ttname"><a href="fracture_estimating_8h.html">fractureEstimating.h</a></div></div>
<div class="ttc" id="debug_functions_8h_html"><div class="ttname"><a href="debug_functions_8h.html">debugFunctions.h</a></div></div>
<div class="ttc" id="input_8h_html_a764fda4f24066351261eb6f777d2ff6c"><div class="ttname"><a href="input_8h.html#a764fda4f24066351261eb6f777d2ff6c">removeFracturesLessThan</a></div><div class="ttdeci">float removeFracturesLessThan</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00260">readInput.cpp:260</a></div></div>
<div class="ttc" id="input_8h_html_a5cd62663cb067d79183d8de4dd0224db"><div class="ttname"><a href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a></div><div class="ttdeci">double domainSize[3]</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00018">readInput.cpp:18</a></div></div>
<div class="ttc" id="debug_functions_8cpp_html_a0d86e9931b24c2b392a6a307d7725c69"><div class="ttname"><a href="debug_functions_8cpp.html#a0d86e9931b24c2b392a6a307d7725c69">printShapeFams</a></div><div class="ttdeci">void printShapeFams(std::vector&lt; Shape &gt; &amp;shapeFamilies)</div><div class="ttdef"><b>Definition:</b> <a href="debug_functions_8cpp_source.html#l00119">debugFunctions.cpp:119</a></div></div>
<div class="ttc" id="struct_poly_html_ae211dd5368d3b023abc943dc940837af"><div class="ttname"><a href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">Poly::familyNum</a></div><div class="ttdeci">int familyNum</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00022">structures.h:22</a></div></div>
<div class="ttc" id="math_functions_8cpp_html_a2540d3ded7bd488a924ce2e4c2443261"><div class="ttname"><a href="math_functions_8cpp.html#a2540d3ded7bd488a924ce2e4c2443261">indexFromProb_and_P32Status</a></div><div class="ttdeci">int indexFromProb_and_P32Status(float *CDF, double roll, int famSize, int cdfSize, int &amp;cdfIdx)</div><div class="ttdef"><b>Definition:</b> <a href="math_functions_8cpp_source.html#l00187">mathFunctions.cpp:187</a></div></div>
<div class="ttc" id="input_8h_html_a8ee9be1b5aa75abae556de3088cba6d9"><div class="ttname"><a href="input_8h.html#a8ee9be1b5aa75abae556de3088cba6d9">h</a></div><div class="ttdeci">double h</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00021">readInput.cpp:21</a></div></div>
<div class="ttc" id="input_8h_html_a0071fcf088e13e14d3cfa0c15e33e363"><div class="ttname"><a href="input_8h.html#a0071fcf088e13e14d3cfa0c15e33e363">disableFram</a></div><div class="ttdeci">bool disableFram</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00030">readInput.cpp:30</a></div></div>
<div class="ttc" id="fracture_estimating_8cpp_html_ae3e14f5bf8ea7b1470e64a74e0354947"><div class="ttname"><a href="fracture_estimating_8cpp.html#ae3e14f5bf8ea7b1470e64a74e0354947">sortRadii</a></div><div class="ttdeci">void sortRadii(std::vector&lt; Shape &gt; &amp;shapeFam)</div><div class="ttdef"><b>Definition:</b> <a href="fracture_estimating_8cpp_source.html#l00016">fractureEstimating.cpp:16</a></div></div>
<div class="ttc" id="read_input_functions_8cpp_html_af648b2c8b7622d441c1bc50f0f9c83b0"><div class="ttname"><a href="read_input_functions_8cpp.html#af648b2c8b7622d441c1bc50f0f9c83b0">getTimeBasedSeed</a></div><div class="ttdeci">unsigned int getTimeBasedSeed()</div><div class="ttdef"><b>Definition:</b> <a href="read_input_functions_8cpp_source.html#l00112">readInputFunctions.cpp:112</a></div></div>
<div class="ttc" id="cluster_groups_8h_html"><div class="ttname"><a href="cluster_groups_8h.html">clusterGroups.h</a></div></div>
<div class="ttc" id="struct_stats_html_aff2b4f1ca734a438e89ad561b82dd7ea"><div class="ttname"><a href="struct_stats.html#aff2b4f1ca734a438e89ad561b82dd7ea">Stats::acceptedPolyCount</a></div><div class="ttdeci">unsigned int acceptedPolyCount</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00300">structures.h:300</a></div></div>
<div class="ttc" id="input_8h_html_a8a68f69e78f6f9c011b7688ef3c48166"><div class="ttname"><a href="input_8h.html#a8a68f69e78f6f9c011b7688ef3c48166">layerVol</a></div><div class="ttdeci">float * layerVol</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00447">readInput.cpp:447</a></div></div>
<div class="ttc" id="remove_fractures_8h_html"><div class="ttname"><a href="remove_fractures_8h.html">removeFractures.h</a></div></div>
<div class="ttc" id="input_8h_html_ad47617e2251da4eb87c56365e217c8e3"><div class="ttname"><a href="input_8h.html#ad47617e2251da4eb87c56365e217c8e3">stopCondition</a></div><div class="ttdeci">short stopCondition</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00012">readInput.cpp:12</a></div></div>
<div class="ttc" id="struct_stats_html_ae78598e8e15f944d358722868f763e32"><div class="ttname"><a href="struct_stats.html#ae78598e8e15f944d358722868f763e32">Stats::rejectedPolyCount</a></div><div class="ttdeci">unsigned long long int rejectedPolyCount</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00304">structures.h:304</a></div></div>
<div class="ttc" id="read_input_functions_8h_html"><div class="ttname"><a href="read_input_functions_8h.html">readInputFunctions.h</a></div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
  <ul>
    <li class="navelem"><a class="el" href="_d_f_nmain_8cpp.html">DFNmain.cpp</a></li>
    <li class="footer">Generated by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.12 </li>
  </ul>
</div>
</body>
</html>
