<!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: insertShape.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('insert_shape_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">insertShape.cpp</div>  </div>
</div><!--header-->
<div class="contents">
<a href="insert_shape_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 &quot;<a class="code" href="structures_8h.html">structures.h</a>&quot;</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</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="l00005"></a><span class="lineno">    5</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="l00006"></a><span class="lineno">    6</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="l00007"></a><span class="lineno">    7</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="l00008"></a><span class="lineno">    8</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="vector_functions_8h.html">vectorFunctions.h</a>&quot;</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="preprocessor">#include &lt;string&gt;</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;</div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;</div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="comment">/**************************************************************************/</span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="comment">/*****************  Generate Polygon/Fracture  ****************************/</span></div><div class="line"><a name="l00026"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#a4f91304a14fefba4d4f93f2f3e54dee4">   26</a></span>&#160;<span class="keyword">struct </span><a class="code" href="struct_poly.html">Poly</a> <a class="code" href="insert_shape_8cpp.html#a4f91304a14fefba4d4f93f2f3e54dee4">generatePoly</a>(struct <a class="code" href="struct_shape.html">Shape</a> &amp;shapeFam, std::mt19937_64 &amp;generator, <a class="code" href="class_distributions.html">Distributions</a> &amp;distributions, <span class="keywordtype">int</span> familyIndex, <span class="keywordtype">bool</span> useList) {</div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;</div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;    <span class="comment">// New polygon to build</span></div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;    <span class="keyword">struct </span><a class="code" href="struct_poly.html">Poly</a> newPoly;</div><div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;    <span class="comment">// Initialize normal to {0,0,1}. ( All polys start on x-y plane )</span></div><div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = 0; <span class="comment">// x</span></div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = 0; <span class="comment">// y</span></div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] =    1; <span class="comment">// z</span></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">// Assign number of nodes</span></div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a> = shapeFam.numPoints; </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;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a> = <span class="keyword">new</span> <span class="keywordtype">double</span>[3*newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a>]; <span class="comment">//numPoints*{x,y,z}</span></div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;    </div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;    <span class="comment">// Assign family number (index of array)</span></div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;    newPoly.<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> = familyIndex;</div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;     </div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;    <span class="comment">// Switch based on distribution type</span></div><div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;    <span class="keywordflow">switch</span> (shapeFam.distributionType) {</div><div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;</div><div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;        <span class="keywordflow">case</span> 1: { <span class="comment">// Lognormal</span></div><div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;            <span class="keywordtype">double</span> radius;</div><div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;            <span class="keywordtype">int</span> count = 1;</div><div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;            <span class="keywordflow">if</span> (shapeFam.radiiIdx &gt;= shapeFam.radiiList.size() || useList == <span class="keyword">false</span>) { </div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;            <span class="comment">// If out of radii from list, insert random radius</span></div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;                std::lognormal_distribution&lt;double&gt; logDistribution(shapeFam.mean, shapeFam.sd);</div><div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;                <span class="keywordflow">do</span> {                    </div><div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;                    radius = logDistribution(generator);</div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;                    </div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;                    <span class="keywordflow">if</span> (count % 1000 == 0) {</div><div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;                        std::cout &lt;&lt; <span class="stringliteral">&quot;\nWARNING: Lognormal distribution for &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFam) </div><div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;                                  &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(familyIndex, shapeFam.shapeFamily) </div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;                                  &lt;&lt; <span class="stringliteral">&quot; has been  unable to generate a fracture with radius within set parameters after &quot;</span> </div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;                                  &lt;&lt; count &lt;&lt; <span class="stringliteral">&quot; consecutive tries.\n&quot;</span>; </div><div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;                        std::cout &lt;&lt; <span class="stringliteral">&quot;Consider adjusting the lognormal paramerters for this family in the input file.\n&quot;</span>;</div><div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;                        <span class="keywordflow">break</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;                    count++;    </div><div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;                } <span class="keywordflow">while</span> (radius &lt; <a class="code" href="input_8h.html#a8ee9be1b5aa75abae556de3088cba6d9">h</a> || radius &lt; shapeFam.logMin || radius &gt; shapeFam.logMax);</div><div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;            }</div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;           <span class="keywordflow">else</span> { <span class="comment">// Insert radius from list</span></div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;                radius = shapeFam.radiiList[shapeFam.radiiIdx];</div><div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;                 shapeFam.radiiIdx++;</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;            }        </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="keywordflow">if</span> (shapeFam.shapeFamily == 1) { <span class="comment">// Rectangle</span></div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;                <span class="comment">// Initialize rectangles vertices using lognormal dist.</span></div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;                <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>);</div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;            } </div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;            <span class="keywordflow">else</span> {<span class="comment">// Ellipse</span></div><div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;                <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>, shapeFam.thetaList, shapeFam.numPoints);</div><div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;            }        </div><div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;            </div><div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;            <span class="keywordflow">break</span>;</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;</div><div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;        <span class="keywordflow">case</span> 2: { <span class="comment">// Truncated power-law</span></div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;            </div><div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;            <span class="keywordtype">double</span> radius;</div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;            </div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;            <span class="keywordflow">if</span> (shapeFam.radiiIdx &gt;= shapeFam.radiiList.size() || useList == <span class="keyword">false</span>) { </div><div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;            <span class="comment">// If out of radii from list, generate random radius</span></div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;                std::uniform_real_distribution&lt;double&gt; uniformDist(0.0,1.0);</div><div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;                 radius = <a class="code" href="generating_points_8cpp.html#a91e58d3a0d7a5e9d65cf15471f02fab0">truncatedPowerLaw</a>( uniformDist(generator), shapeFam.min, shapeFam.max, shapeFam.alpha);</div><div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;             }</div><div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;             <span class="keywordflow">else</span> { <span class="comment">// Pull radius from list</span></div><div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;                 radius = shapeFam.radiiList[shapeFam.radiiIdx];</div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;                 shapeFam.radiiIdx++;</div><div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;             }</div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;             </div><div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;            <span class="keywordflow">if</span> (shapeFam.shapeFamily == 1) {</div><div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;                <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>);</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="keywordflow">else</span> {</div><div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;                <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>, shapeFam.thetaList, shapeFam.numPoints);</div><div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;            }    </div><div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;            <span class="keywordflow">break</span>;</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;            </div><div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;        <span class="keywordflow">case</span> 3: { <span class="comment">// Exponential</span></div><div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;            <span class="keywordtype">double</span> radius;</div><div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;            <span class="keywordtype">int</span> count = 1;</div><div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;            <span class="keywordflow">if</span> (shapeFam.radiiIdx &gt;= shapeFam.radiiList.size() || useList == <span class="keyword">false</span>) { </div><div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;            <span class="comment">// If out of radii from list, generate random radius</span></div><div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;                <span class="keywordflow">do</span> {</div><div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;                    radius = distributions.expDist-&gt;getValue(shapeFam.expLambda, shapeFam.minDistInput, shapeFam.maxDistInput);</div><div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;                    <span class="keywordflow">if</span> (count % 1000 == 0) {</div><div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;                        std::cout &lt;&lt; <span class="stringliteral">&quot;\nWARNING: Exponential distribution for &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(shapeFam) </div><div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;                                  &lt;&lt; <span class="stringliteral">&quot; family &quot;</span> &lt;&lt; <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(familyIndex, shapeFam.shapeFamily) </div><div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;                                  &lt;&lt; <span class="stringliteral">&quot; has been  unable to generate a fracture with radius within set parameters after &quot;</span> </div><div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;                                  &lt;&lt; count &lt;&lt; <span class="stringliteral">&quot; consecutive tries.\n&quot;</span>; </div><div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160;                        std::cout &lt;&lt; <span class="stringliteral">&quot;Consider adjusting the exponential parameters for this family in the input file.\n&quot;</span>; </div><div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;                        <span class="keywordflow">break</span>;</div><div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;                    }</div><div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;                            </div><div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;                    count++;</div><div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;</div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;                } <span class="keywordflow">while</span> (radius &lt; <a class="code" href="input_8h.html#a8ee9be1b5aa75abae556de3088cba6d9">h</a> || radius &lt; shapeFam.expMin || radius &gt; shapeFam.expMax);</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;            <span class="keywordflow">else</span> { <span class="comment">// Insert radius from list</span></div><div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;                    radius = shapeFam.radiiList[shapeFam.radiiIdx];</div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;                     shapeFam.radiiIdx++;</div><div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;            }</div><div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;              </div><div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;            <span class="keywordflow">if</span> (shapeFam.shapeFamily == 1) { <span class="comment">// Rectangle</span></div><div class="line"><a name="l00131"></a><span class="lineno">  131</span>&#160;                <span class="comment">// Initialize rectangles vertices using exp. dist.</span></div><div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;                <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>);        </div><div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;            } </div><div class="line"><a name="l00134"></a><span class="lineno">  134</span>&#160;            <span class="keywordflow">else</span> { <span class="comment">// Ellipse</span></div><div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160;                <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, radius , shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>, shapeFam.thetaList, shapeFam.numPoints);                    </div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;            }    </div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;        }</div><div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;            </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">case</span> 4: { <span class="comment">// Constant</span></div><div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;            <span class="keywordflow">if</span> (shapeFam.shapeFamily == 1) { <span class="comment">// Rectangle</span></div><div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;                <span class="comment">// Initialize rectangles vertices</span></div><div class="line"><a name="l00144"></a><span class="lineno">  144</span>&#160;                <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(newPoly, shapeFam.constRadi, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>);</div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;            }</div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;            <span class="keywordflow">else</span> { <span class="comment">// Ellipse</span></div><div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;                <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, shapeFam.constRadi, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>,shapeFam.thetaList, shapeFam.numPoints);</div><div class="line"><a name="l00148"></a><span class="lineno">  148</span>&#160;            }</div><div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;        }</div><div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;    } </div><div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;    </div><div class="line"><a name="l00153"></a><span class="lineno">  153</span>&#160;    </div><div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;    <span class="keywordtype">double</span> beta;</div><div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;    <span class="comment">// Initialize beta based on distrubution type: 0 = unifrom on [0,2PI], 1 = constant</span></div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;    <span class="keywordflow">if</span> (shapeFam.betaDistribution == 0) { <span class="comment">//uniform distribution</span></div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;        std::uniform_real_distribution&lt;double&gt; uniformDist (0, 2*M_PI);</div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;        beta = uniformDist(generator);</div><div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;    }</div><div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;        beta = shapeFam.beta;</div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;    }</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">// Apply 2d rotation matrix, twist around origin</span></div><div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;    <span class="comment">// Assumes polygon on x-y plane</span></div><div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;    <span class="comment">// Angle must be in rad</span></div><div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#a878eeb1b4990cde9d7d12a64f647af4d">applyRotation2D</a>(newPoly, beta);</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">// Fisher distribution / get normal vector</span></div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;    <span class="keywordtype">double</span> *norm = <a class="code" href="generating_points_8cpp.html#a2f7b981bd314f565407a2db0358bf81b">fisherDistribution</a>(shapeFam.theta, shapeFam.phi, shapeFam.kappa, generator);</div><div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;    <span class="keywordtype">double</span> mag = <a class="code" href="vector_functions_8h.html#a92da2058f3ed2a6db16d2fe0e445794a">magnitude</a>(norm[0], norm[1], norm[2]); </div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;    <span class="keywordflow">if</span> (mag &lt; 1-eps || mag &gt; 1+<a class="code" href="_d_f_nmain_8cpp.html#a974b448c3d0a13d483d36ef15de369b5">eps</a>) { </div><div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;        <a class="code" href="vector_functions_8h.html#a48bb00fc4195377b229efe044f3bb67e">normalize</a>(norm); <span class="comment">// Ensure norm is normalized</span></div><div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;    }</div><div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;</div><div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#af7f1ba8a57908229078780a2695e4e46">applyRotation3D</a>(newPoly, norm); <span class="comment">// Rotate vertices to norm (new normal)    </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">// Save newPoly&#39;s new normal vector</span></div><div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = norm[0];</div><div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = norm[1];</div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] = norm[2];</div><div class="line"><a name="l00182"></a><span class="lineno">  182</span>&#160;    <span class="keyword">delete</span>[] norm;</div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;</div><div class="line"><a name="l00184"></a><span class="lineno">  184</span>&#160;    <span class="keywordtype">double</span> *t;</div><div class="line"><a name="l00185"></a><span class="lineno">  185</span>&#160;    </div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;    <span class="keywordflow">if</span> (shapeFam.layer == 0 ) { <span class="comment">// The family layer is the whole domain</span></div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2);</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;    }</div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;    <span class="keywordflow">else</span> { <span class="comment">// Family belongs to a certain layer, shapeFam.layer is &gt; zero</span></div><div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;    </div><div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;        <span class="comment">// Layers start at 1, but the array of layers start at 0 hence sub of 1</span></div><div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;        <span class="comment">// Layer 0 is reservered to denote the entire domain</span></div><div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;        <span class="keywordtype">int</span> layerIdx = (shapeFam.layer - 1) * 2; </div><div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;        </div><div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;        <span class="comment">// Layers only apply to z coordinates</span></div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;        t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx] , <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx+1]);</div><div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;    } <span class="comment">// End else</span></div><div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;</div><div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;    <span class="comment">// Translate - will also set translation vector in poly structure</span></div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#a3098cccf29c7a47d2eae661b0f545224">translate</a>(newPoly, t);</div><div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;    <span class="keyword">delete</span>[] t;</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="keywordflow">return</span> newPoly;</div><div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;}</div><div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;</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">/**************************************************************************/</span></div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;<span class="comment">/*************  Generate Polygon/Fracture With Given Radius  **************/</span></div><div class="line"><a name="l00221"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#a4a4d307b0a091a6159358f312371ae95">  221</a></span>&#160;<span class="keyword">struct </span><a class="code" href="struct_poly.html">Poly</a> <a class="code" href="insert_shape_8cpp.html#a4a4d307b0a091a6159358f312371ae95">generatePoly_withRadius</a>(double radius,struct <a class="code" href="struct_shape.html">Shape</a> &amp;shapeFam, std::mt19937_64 &amp;generator, <a class="code" href="class_distributions.html">Distributions</a> &amp;distributions, <span class="keywordtype">int</span> familyIndex) {</div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;</div><div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;    <span class="comment">// New polygon to build</span></div><div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160;    <span class="keyword">struct </span><a class="code" href="struct_poly.html">Poly</a> newPoly;</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;    <span class="comment">// Initialize normal to {0,0,1}. ( All polys start on x-y plane )</span></div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = 0; <span class="comment">//x</span></div><div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = 0; <span class="comment">//y</span></div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] = 1; <span class="comment">//z</span></div><div class="line"><a name="l00229"></a><span class="lineno">  229</span>&#160;</div><div class="line"><a name="l00230"></a><span class="lineno">  230</span>&#160;    <span class="comment">// Assign number of nodes</span></div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a> = shapeFam.numPoints; </div><div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a> = <span class="keyword">new</span> <span class="keywordtype">double</span>[3*newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a>]; <span class="comment">//numPoints*{x,y,z}</span></div><div class="line"><a name="l00233"></a><span class="lineno">  233</span>&#160;</div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;    <span class="comment">// Assign family number (index of shapeFam array)</span></div><div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;    newPoly.<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> = familyIndex;</div><div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;    </div><div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;    <span class="comment">// TODO: Convert any degrees to rad </span></div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;    <span class="comment">// in readInput() to avoid continuous checking</span></div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;</div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;    <span class="keywordflow">if</span> (shapeFam.shapeFamily == 1) { <span class="comment">// If rectangle shape</span></div><div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;        <span class="comment">// Initialize rectangles vertices</span></div><div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;        <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>);</div><div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;    }</div><div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;    <span class="keywordflow">else</span> { <span class="comment">// Ellipse</span></div><div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;        <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, radius, shapeFam.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>,shapeFam.thetaList, shapeFam.numPoints);</div><div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;    }</div><div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;    </div><div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;    <span class="keywordtype">double</span> beta;</div><div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;    <span class="comment">// Initialize beta based on distrubution type: 0 = unifrom on [0,2PI], 1 = constant</span></div><div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;    <span class="keywordflow">if</span> (shapeFam.betaDistribution == 0) { <span class="comment">// Uniform distribution</span></div><div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;        std::uniform_real_distribution&lt;double&gt; uniformDist (0, 2*M_PI);</div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;        beta = uniformDist(generator);</div><div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;    }</div><div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;        beta = shapeFam.beta;</div><div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;    }</div><div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;</div><div class="line"><a name="l00258"></a><span class="lineno">  258</span>&#160;    <span class="comment">// Apply 2d rotation matrix, twist around origin</span></div><div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;    <span class="comment">// assumes polygon on x-y plane</span></div><div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;    <span class="comment">// Angle must be in rad</span></div><div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#a878eeb1b4990cde9d7d12a64f647af4d">applyRotation2D</a>(newPoly, beta);</div><div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;</div><div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;    <span class="comment">// Fisher distribution / get normal vector</span></div><div class="line"><a name="l00264"></a><span class="lineno">  264</span>&#160;    <span class="keywordtype">double</span> *norm = <a class="code" href="generating_points_8cpp.html#a2f7b981bd314f565407a2db0358bf81b">fisherDistribution</a>(shapeFam.theta, shapeFam.phi, shapeFam.kappa, generator);</div><div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;    <span class="keywordtype">double</span> mag = <a class="code" href="vector_functions_8h.html#a92da2058f3ed2a6db16d2fe0e445794a">magnitude</a>(norm[0], norm[1], norm[2]); </div><div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;    <span class="keywordflow">if</span> (mag &lt; 1-eps || mag &gt; 1+<a class="code" href="_d_f_nmain_8cpp.html#a974b448c3d0a13d483d36ef15de369b5">eps</a>) { </div><div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;        <a class="code" href="vector_functions_8h.html#a48bb00fc4195377b229efe044f3bb67e">normalize</a>(norm); <span class="comment">//ensure norm is normalized</span></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;</div><div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#af7f1ba8a57908229078780a2695e4e46">applyRotation3D</a>(newPoly, norm); <span class="comment">// Rotate vertices to norm (new normal)    </span></div><div class="line"><a name="l00271"></a><span class="lineno">  271</span>&#160;    </div><div class="line"><a name="l00272"></a><span class="lineno">  272</span>&#160;    <span class="comment">// Save newPoly&#39;s new normal vector</span></div><div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = norm[0];</div><div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = norm[1];</div><div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] = norm[2];</div><div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;    <span class="keyword">delete</span>[] norm;</div><div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;    </div><div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;    <span class="keywordtype">double</span> *t;</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="keywordflow">if</span> (shapeFam.layer == 0 ) { <span class="comment">// The family layer is the whole domain</span></div><div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;        t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, </div><div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;            (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, </div><div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;            (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2, </div><div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;            (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2);</div><div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;        </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;    <span class="keywordflow">else</span> { <span class="comment">// Family belongs to a certain layer, shapeFam.layer is &gt; zero</span></div><div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;    </div><div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;        <span class="comment">// Layers start at 1, but the array of layers start at 0, hence</span></div><div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;        <span class="comment">// the subtraction by 1</span></div><div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;        <span class="comment">// Layer 0 is reservered to be the entire domain</span></div><div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;        <span class="keywordtype">int</span> layerIdx = (shapeFam.layer - 1) * 2; </div><div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;        </div><div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160;        <span class="comment">// Layers only apply to z coordinates</span></div><div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;        t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, </div><div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;            (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, </div><div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;            (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx] , <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx+1]);</div><div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;    </div><div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;    } <span class="comment">// End else</span></div><div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;</div><div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;    <span class="comment">// Translate - will also set translation vector in poly structure</span></div><div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;    <a class="code" href="computational_geometry_8cpp.html#a3098cccf29c7a47d2eae661b0f545224">translate</a>(newPoly, t);</div><div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;    <span class="keyword">delete</span>[] t;</div><div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160;           </div><div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;    <span class="keywordflow">return</span> newPoly;    </div><div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;}</div><div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;</div><div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;</div><div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;<span class="comment">/*******************************************************************************/</span></div><div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;<span class="comment">/******************** Aperture  assignment function ****************************/</span></div><div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;<span class="comment">// water density = 997.7</span></div><div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;<span class="comment">// gravity = 9.8</span></div><div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;<span class="comment">// water Visc = 8.9e-4</span></div><div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;<span class="comment">// constant scalar = waterDesnsity*gravity/waterVisc = 48.3868</span></div><div class="line"><a name="l00329"></a><span class="lineno"><a class="line" href="insert_shape_8cpp.html#ab707f54136b034587f94b633476787a0">  329</a></span>&#160;<span class="preprocessor">#define _CONSTSCALAR 48.3868</span></div><div class="line"><a name="l00330"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#a19e9a86f45510bc87d7fe84a11caffcf">  330</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#a19e9a86f45510bc87d7fe84a11caffcf">assignAperture</a>(<span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> &amp;newPoly, std::mt19937_64 &amp;generator) {</div><div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;    </div><div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;    <span class="comment">// Most aperture variables are currently declared globaly</span></div><div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;    <span class="keywordflow">switch</span> (<a class="code" href="input_8h.html#abc39a0cf64740d1080e5c9867191fa0e">aperture</a>) {</div><div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;        <span class="keywordflow">case</span> 1: { <span class="comment">// Lognormal</span></div><div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;            std::lognormal_distribution&lt;double&gt; logDistribution(<a class="code" href="input_8h.html#a8ca71d341876667108c011c790ed33de">meanAperture</a>, <a class="code" href="input_8h.html#a0b60c5d860054f964ec05348b186e86e">stdAperture</a>);</div><div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a> = logDistribution(generator);</div><div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;            <span class="keywordflow">break</span>;</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="keywordflow">case</span> 2:{</div><div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;        <span class="comment">/*  Transmissivity is calculated as transmissivity = F*R^k,</span></div><div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;<span class="comment">            where F is a first element in aperturefromTransmissivity,</span></div><div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;<span class="comment">            k is a second element and R is a mean radius of a polygon. </span></div><div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;<span class="comment">            Aperture is calculated according to cubic law as </span></div><div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;<span class="comment">            b=(transmissivity*12)^1/3  */</span></div><div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;            <span class="keywordtype">float</span> radiAvg = (newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a> + newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>)*.5;</div><div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;            <span class="keywordtype">double</span> F = <a class="code" href="input_8h.html#aafcf8c3d5090c1cda444c5a79854b434">apertureFromTransmissivity</a>[0];</div><div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;            <span class="keywordtype">double</span> k = <a class="code" href="input_8h.html#aafcf8c3d5090c1cda444c5a79854b434">apertureFromTransmissivity</a>[1];</div><div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;            <span class="keywordtype">double</span> transmissivity = F * std::pow(radiAvg, k);                            </div><div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a> = std::cbrt((transmissivity*12/<a class="code" href="insert_shape_8cpp.html#ab707f54136b034587f94b633476787a0">_CONSTSCALAR</a>)); <span class="comment">//cube root</span></div><div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;        }</div><div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;    </div><div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;    </div><div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;        <span class="keywordflow">case</span> 3:{</div><div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a> = <a class="code" href="input_8h.html#a48936b0f8f54a1cd2c7e83ca2392ec0a">constantAperture</a>;            </div><div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;            <span class="keywordflow">break</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;    </div><div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;        <span class="keywordflow">case</span> 4:{</div><div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;            <span class="keywordtype">double</span> radiAvg = (newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a> + newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>)*.5;</div><div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;            <span class="keywordtype">double</span> apertureMeanF = <a class="code" href="input_8h.html#a96e535ecb2fde247beaba684f23b1c3c">lengthCorrelatedAperture</a>[0];</div><div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;            <span class="keywordtype">double</span> b = <a class="code" href="input_8h.html#a96e535ecb2fde247beaba684f23b1c3c">lengthCorrelatedAperture</a>[1];</div><div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a> = apertureMeanF * std::pow(radiAvg, b);</div><div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;        </div><div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;            <span class="keywordflow">break</span>;</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;<span class="comment">/**************************************************************************/</span></div><div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;<span class="comment">/**************** permiability  assignment function ***********************/</span></div><div class="line"><a name="l00380"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#ab66dd8460f2f45feae57425d7c219bf5">  380</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#ab66dd8460f2f45feae57425d7c219bf5">assignPermeability</a>(<span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> &amp;newPoly) {</div><div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;        </div><div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;        <span class="comment">// Aperture variables are declared globaly</span></div><div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;        <span class="keywordflow">switch</span> (<a class="code" href="input_8h.html#a6fda5f9c90f07130c1b0001301c6c945">permOption</a>) {</div><div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;            <span class="keywordflow">case</span> 0: { <span class="comment">// Perm as function of aperture</span></div><div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;                newPoly.<a class="code" href="struct_poly.html#a48c792ab99f5f7ca554a3f77062419db">permeability</a> = (newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a> * newPoly.<a class="code" href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">aperture</a>)/12;</div><div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;                <span class="keywordflow">break</span>;</div><div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;            }</div><div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;            </div><div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;            <span class="keywordflow">case</span> 1:{</div><div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;                newPoly.<a class="code" href="struct_poly.html#a48c792ab99f5f7ca554a3f77062419db">permeability</a> = <a class="code" href="input_8h.html#a903896dcf39ee23cbea885abd3bc926d">constantPermeability</a>;</div><div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;                <span class="keywordflow">break</span>;</div><div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160;            }</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;</div><div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;<span class="comment">/**************  Initialize Rectangular Vertices  *****************/</span></div><div class="line"><a name="l00403"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#aed01cfb18a8a6b4ec394783593d5765f">  403</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a>(<span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> &amp;newPoly, <span class="keywordtype">float</span> radius, <span class="keywordtype">float</span> <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>) {</div><div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;    <span class="keywordtype">float</span> x = radius;</div><div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;    <span class="keywordtype">float</span> y = radius * <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>;</div><div class="line"><a name="l00406"></a><span class="lineno">  406</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a> = x;</div><div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;    newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a> = y;</div><div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a> = <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>;</div><div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;</div><div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;    <span class="comment">// Initialize vertices</span></div><div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[0] = x;</div><div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[1] = y;</div><div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[2] = 0;</div><div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[3] = -x;</div><div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[4] = y;</div><div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[5] = 0;</div><div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[6] = -x;</div><div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[7] = -y;</div><div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[8] = 0;</div><div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[9] = x;</div><div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[10] = -y;</div><div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[11] = 0; </div><div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;}</div><div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;</div><div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;</div><div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;<span class="comment">/****************************************************************/</span></div><div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;<span class="comment">/************* Initialize Ellipse Vertices **********************/</span></div><div class="line"><a name="l00432"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#a8ed0d23da7ff9cd1f7c4109a4c478071">  432</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(<span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> &amp;newPoly, <span class="keywordtype">float</span> radius, <span class="keywordtype">float</span> <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>, <span class="keywordtype">float</span> *thetaList, <span class="keywordtype">int</span> numPoints) {</div><div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;    </div><div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a> = radius;</div><div class="line"><a name="l00435"></a><span class="lineno">  435</span>&#160;    newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a> = radius * <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>;</div><div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;    newPoly.<a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a> = <a class="code" href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">aspectRatio</a>;</div><div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; numPoints; i++ ) {</div><div class="line"><a name="l00438"></a><span class="lineno">  438</span>&#160;        <span class="keywordtype">int</span> idx = i*3;</div><div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx]   = radius * std::cos(thetaList[i]);</div><div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx+1] = radius * aspectRatio * std::sin(thetaList[i]);</div><div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx+2] = 0;</div><div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;    }</div><div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;}</div><div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160;</div><div class="line"><a name="l00445"></a><span class="lineno">  445</span>&#160;</div><div class="line"><a name="l00446"></a><span class="lineno">  446</span>&#160;<span class="comment">/**********************************************************************/</span></div><div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;<span class="comment">/********************  Retranslate Polygon  ***************************/</span></div><div class="line"><a name="l00455"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#afa95ede56e6ecb8918b61c5a2815d40d">  455</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#afa95ede56e6ecb8918b61c5a2815d40d">reTranslatePoly</a>(<span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> &amp;newPoly, <span class="keyword">struct</span> <a class="code" href="struct_shape.html">Shape</a> &amp;shapeFam, std::mt19937_64 &amp;generator) {</div><div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;</div><div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;</div><div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;    <span class="keywordflow">if</span> (newPoly.<a class="code" href="struct_poly.html#a341a2ec25f71c2f58c1af119b3990387">truncated</a> == 0) {  </div><div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;        <span class="comment">// If poly isn&#39;t truncated we can skip a lot of steps such </span></div><div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;        <span class="comment">// as reallocating vertice memory, rotations, etc.. </span></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;        newPoly.<a class="code" href="struct_poly.html#add8e39b7f48274757f46c80014aaeb6b">groupNum</a> = 0; <span class="comment">// Clear cluster group information</span></div><div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;        newPoly.<a class="code" href="struct_poly.html#aa6c373c49a4e0ab7c83a353a9f8cf518">intersectionIndex</a>.clear(); <span class="comment">// Clear any saved intersections</span></div><div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;    </div><div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160;        <span class="comment">// Move poly back to origin    </span></div><div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a>; i++) {</div><div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;            <span class="keywordtype">int</span> idx = 3*i;</div><div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx]   -= newPoly.<a class="code" href="struct_poly.html#ac93df10849f4916e195a899d2ba55485">translation</a>[0]; <span class="comment">// x</span></div><div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx+1] -= newPoly.<a class="code" href="struct_poly.html#ac93df10849f4916e195a899d2ba55485">translation</a>[1]; <span class="comment">// y</span></div><div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[idx+2] -= newPoly.<a class="code" href="struct_poly.html#ac93df10849f4916e195a899d2ba55485">translation</a>[2]; <span class="comment">// z</span></div><div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;        }</div><div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;                </div><div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;        <span class="comment">// Translate to new position</span></div><div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;        <span class="keywordtype">double</span> *t;</div><div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;        <span class="keywordflow">if</span> (shapeFam.<a class="code" href="struct_shape.html#a0ce256fc0e226ea6559786a2ae9d007f">layer</a> == 0 ) { <span class="comment">// The family layer is the whole domain</span></div><div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;            t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2);</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="keywordflow">else</span> { <span class="comment">// Family belongs to a certain layer</span></div><div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;</div><div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;            <span class="comment">// Layers start at 1, but the array of layers start at 0. Layer 0 is reservered to be the entire domain</span></div><div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;            <span class="keywordtype">int</span> layerIdx = (shapeFam.<a class="code" href="struct_shape.html#a0ce256fc0e226ea6559786a2ae9d007f">layer</a> - 1) * 2; </div><div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;</div><div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;            <span class="comment">// Layers only apply to z coordinates</span></div><div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;            t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx] , <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx+1]);</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;        } <span class="comment">// End else</span></div><div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;        </div><div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;        <span class="comment">// Translate - will also set translation vector in poly structure</span></div><div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;        <a class="code" href="computational_geometry_8cpp.html#a3098cccf29c7a47d2eae661b0f545224">translate</a>(newPoly, t);</div><div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160;        <span class="keyword">delete</span>[] t;    </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;        </div><div class="line"><a name="l00493"></a><span class="lineno">  493</span>&#160;    <span class="keywordflow">else</span> { <span class="comment">// Poly was truncated, need to rebuild the polygon</span></div><div class="line"><a name="l00494"></a><span class="lineno">  494</span>&#160;</div><div class="line"><a name="l00495"></a><span class="lineno">  495</span>&#160;        <span class="keyword">delete</span>[] newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>; <span class="comment">// Delete truncated vertices</span></div><div class="line"><a name="l00496"></a><span class="lineno">  496</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a> = <span class="keyword">new</span> <span class="keywordtype">double</span>[shapeFam.<a class="code" href="struct_shape.html#a463c3c0ba9465daaac854eaeccdb3341">numPoints</a>*3];</div><div class="line"><a name="l00497"></a><span class="lineno">  497</span>&#160;</div><div class="line"><a name="l00498"></a><span class="lineno">  498</span>&#160;        <span class="comment">// Reset boundary faces (0 means poly is no longer touching a boundary)</span></div><div class="line"><a name="l00499"></a><span class="lineno">  499</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[0] = 0;</div><div class="line"><a name="l00500"></a><span class="lineno">  500</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[1] = 0;</div><div class="line"><a name="l00501"></a><span class="lineno">  501</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[2] = 0;</div><div class="line"><a name="l00502"></a><span class="lineno">  502</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[3] = 0;</div><div class="line"><a name="l00503"></a><span class="lineno">  503</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[4] = 0;</div><div class="line"><a name="l00504"></a><span class="lineno">  504</span>&#160;        newPoly.<a class="code" href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">faces</a>[5] = 0;        </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;        newPoly.<a class="code" href="struct_poly.html#a341a2ec25f71c2f58c1af119b3990387">truncated</a> = 0; <span class="comment">// Set to 0 to mean not truncated</span></div><div class="line"><a name="l00507"></a><span class="lineno">  507</span>&#160;        newPoly.<a class="code" href="struct_poly.html#add8e39b7f48274757f46c80014aaeb6b">groupNum</a> = 0;  <span class="comment">// Clear cluster group information</span></div><div class="line"><a name="l00508"></a><span class="lineno">  508</span>&#160;        </div><div class="line"><a name="l00509"></a><span class="lineno">  509</span>&#160;        newPoly.<a class="code" href="struct_poly.html#aa6c373c49a4e0ab7c83a353a9f8cf518">intersectionIndex</a>.clear(); <span class="comment">// Clear any saved intersections</span></div><div class="line"><a name="l00510"></a><span class="lineno">  510</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">numberOfNodes</a> = shapeFam.<a class="code" href="struct_shape.html#a463c3c0ba9465daaac854eaeccdb3341">numPoints</a>;</div><div class="line"><a name="l00511"></a><span class="lineno">  511</span>&#160;                </div><div class="line"><a name="l00512"></a><span class="lineno">  512</span>&#160;        <span class="keywordflow">if</span> (shapeFam.<a class="code" href="struct_shape.html#a766d04a69991e931ba2176a927c747cb">shapeFamily</a> == 1) { <span class="comment">// 1 is rectanglular families. rebuild rectangle</span></div><div class="line"><a name="l00513"></a><span class="lineno">  513</span>&#160;            </div><div class="line"><a name="l00514"></a><span class="lineno">  514</span>&#160;            <span class="comment">// Rebuild poly at origin using previous size</span></div><div class="line"><a name="l00515"></a><span class="lineno">  515</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[0] = newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a>;</div><div class="line"><a name="l00516"></a><span class="lineno">  516</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[1] = newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>;</div><div class="line"><a name="l00517"></a><span class="lineno">  517</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[2] = 0;</div><div class="line"><a name="l00518"></a><span class="lineno">  518</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[3] = -newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a>;</div><div class="line"><a name="l00519"></a><span class="lineno">  519</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[4] = newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>;</div><div class="line"><a name="l00520"></a><span class="lineno">  520</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[5] = 0;</div><div class="line"><a name="l00521"></a><span class="lineno">  521</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[6] = -newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a>;</div><div class="line"><a name="l00522"></a><span class="lineno">  522</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[7] = -newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>;</div><div class="line"><a name="l00523"></a><span class="lineno">  523</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[8] = 0;</div><div class="line"><a name="l00524"></a><span class="lineno">  524</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[9] = newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a>;</div><div class="line"><a name="l00525"></a><span class="lineno">  525</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[10] = -newPoly.<a class="code" href="struct_poly.html#ad1acec5af38b346429eada06f05078c6">yradius</a>;</div><div class="line"><a name="l00526"></a><span class="lineno">  526</span>&#160;            newPoly.<a class="code" href="struct_poly.html#a2cb519bcf0553215de907c8fc878a255">vertices</a>[11] = 0; </div><div class="line"><a name="l00527"></a><span class="lineno">  527</span>&#160;        }</div><div class="line"><a name="l00528"></a><span class="lineno">  528</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// Rebuild ellipse</span></div><div class="line"><a name="l00529"></a><span class="lineno">  529</span>&#160;            <a class="code" href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a>(newPoly, newPoly.<a class="code" href="struct_poly.html#a2e10824b90d60b54522b7ce475ed2a17">xradius</a>, shapeFam.<a class="code" href="struct_shape.html#a6126c3f7cd57df03350d4b58b844bba1">aspectRatio</a>, shapeFam.<a class="code" href="struct_shape.html#a0d1fd1363bacd61560809dc379d17363">thetaList</a>, shapeFam.<a class="code" href="struct_shape.html#a463c3c0ba9465daaac854eaeccdb3341">numPoints</a>);        </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;                                        </div><div class="line"><a name="l00532"></a><span class="lineno">  532</span>&#160;        <span class="comment">// Save newPoly&#39;s previous normal vector and then reset poly normal to {0,0,1} for applyRotation3D function</span></div><div class="line"><a name="l00533"></a><span class="lineno">  533</span>&#160;        <span class="keywordtype">double</span> normalB[3] = {newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0], newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1], newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2]}; </div><div class="line"><a name="l00534"></a><span class="lineno">  534</span>&#160;         </div><div class="line"><a name="l00535"></a><span class="lineno">  535</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = 0; <span class="comment">//x</span></div><div class="line"><a name="l00536"></a><span class="lineno">  536</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = 0; <span class="comment">//y</span></div><div class="line"><a name="l00537"></a><span class="lineno">  537</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] = 1; <span class="comment">//z</span></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="keywordtype">double</span> beta;</div><div class="line"><a name="l00540"></a><span class="lineno">  540</span>&#160;        <span class="comment">// Initialize beta based on distrubution type: 0 = unifrom on [0,2PI], 1 = constant</span></div><div class="line"><a name="l00541"></a><span class="lineno">  541</span>&#160;        <span class="keywordflow">if</span> (shapeFam.<a class="code" href="struct_shape.html#a874e4756ad81f2e6f135e066bb0b4347">betaDistribution</a> == 0) { </div><div class="line"><a name="l00542"></a><span class="lineno">  542</span>&#160;            <span class="comment">// Uniform distribution</span></div><div class="line"><a name="l00543"></a><span class="lineno">  543</span>&#160;            std::uniform_real_distribution&lt;double&gt; uniformDist (0, 2*M_PI);</div><div class="line"><a name="l00544"></a><span class="lineno">  544</span>&#160;            beta = uniformDist(generator);</div><div class="line"><a name="l00545"></a><span class="lineno">  545</span>&#160;        }</div><div class="line"><a name="l00546"></a><span class="lineno">  546</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// Constant</span></div><div class="line"><a name="l00547"></a><span class="lineno">  547</span>&#160;            beta = shapeFam.<a class="code" href="struct_shape.html#ab48639160dae4dfe2e3a0a55ffa5a691">beta</a>;</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;</div><div class="line"><a name="l00550"></a><span class="lineno">  550</span>&#160;        <span class="comment">// Apply 2d rotation matrix, twist around origin</span></div><div class="line"><a name="l00551"></a><span class="lineno">  551</span>&#160;        <span class="comment">// Assumes polygon on x-y plane</span></div><div class="line"><a name="l00552"></a><span class="lineno">  552</span>&#160;        <span class="comment">// Angle must be in rad</span></div><div class="line"><a name="l00553"></a><span class="lineno">  553</span>&#160;        <a class="code" href="computational_geometry_8cpp.html#a878eeb1b4990cde9d7d12a64f647af4d">applyRotation2D</a>(newPoly, beta);</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;        <span class="comment">// Rotates poly from {0,0,1} to normalB, NEED to save normalB to newPoly.normal afterwards</span></div><div class="line"><a name="l00556"></a><span class="lineno">  556</span>&#160;        <a class="code" href="computational_geometry_8cpp.html#af7f1ba8a57908229078780a2695e4e46">applyRotation3D</a>(newPoly, normalB);</div><div class="line"><a name="l00557"></a><span class="lineno">  557</span>&#160;        </div><div class="line"><a name="l00558"></a><span class="lineno">  558</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[0] = normalB[0];</div><div class="line"><a name="l00559"></a><span class="lineno">  559</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[1] = normalB[1];</div><div class="line"><a name="l00560"></a><span class="lineno">  560</span>&#160;        newPoly.<a class="code" href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">normal</a>[2] = normalB[2];</div><div class="line"><a name="l00561"></a><span class="lineno">  561</span>&#160;</div><div class="line"><a name="l00562"></a><span class="lineno">  562</span>&#160;        <span class="comment">// Translate to new position</span></div><div class="line"><a name="l00563"></a><span class="lineno">  563</span>&#160;        <span class="comment">// Translate() will also set translation vector in poly structure</span></div><div class="line"><a name="l00564"></a><span class="lineno">  564</span>&#160;        <span class="keywordtype">double</span> *t;</div><div class="line"><a name="l00565"></a><span class="lineno">  565</span>&#160;        <span class="keywordflow">if</span> (shapeFam.<a class="code" href="struct_shape.html#a0ce256fc0e226ea6559786a2ae9d007f">layer</a> == 0 ) { <span class="comment">// The family layer is the whole domain</span></div><div class="line"><a name="l00566"></a><span class="lineno">  566</span>&#160;            t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[2]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[2])/2);</div><div class="line"><a name="l00567"></a><span class="lineno">  567</span>&#160;        }</div><div class="line"><a name="l00568"></a><span class="lineno">  568</span>&#160;        <span class="keywordflow">else</span> { <span class="comment">// Family belongs to a certain layer</span></div><div class="line"><a name="l00569"></a><span class="lineno">  569</span>&#160;</div><div class="line"><a name="l00570"></a><span class="lineno">  570</span>&#160;            <span class="comment">// Layers start at #1, but the array index of layers start at 0. ShapeFam.layer=0 is reservered to be the entire domain</span></div><div class="line"><a name="l00571"></a><span class="lineno">  571</span>&#160;            <span class="keywordtype">int</span> layerIdx = (shapeFam.<a class="code" href="struct_shape.html#a0ce256fc0e226ea6559786a2ae9d007f">layer</a> - 1) * 2; </div><div class="line"><a name="l00572"></a><span class="lineno">  572</span>&#160;    </div><div class="line"><a name="l00573"></a><span class="lineno">  573</span>&#160;            <span class="comment">// Layers only apply to z coordinates</span></div><div class="line"><a name="l00574"></a><span class="lineno">  574</span>&#160;            t = <a class="code" href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a>(generator, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[0]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[0])/2, (-<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]-<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, (<a class="code" href="input_8h.html#a5cd62663cb067d79183d8de4dd0224db">domainSize</a>[1]+<a class="code" href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a>[1])/2, <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx] , <a class="code" href="input_8h.html#a03ac7830eb87def8d67cc0c95b24ec82">layers</a>[layerIdx+1]);</div><div class="line"><a name="l00575"></a><span class="lineno">  575</span>&#160;</div><div class="line"><a name="l00576"></a><span class="lineno">  576</span>&#160;        } <span class="comment">// End else</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;        <a class="code" href="computational_geometry_8cpp.html#a3098cccf29c7a47d2eae661b0f545224">translate</a>(newPoly, t);</div><div class="line"><a name="l00579"></a><span class="lineno">  579</span>&#160;        <span class="keyword">delete</span>[] t;        </div><div class="line"><a name="l00580"></a><span class="lineno">  580</span>&#160;    }</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;</div><div class="line"><a name="l00583"></a><span class="lineno">  583</span>&#160;</div><div class="line"><a name="l00584"></a><span class="lineno">  584</span>&#160;<span class="comment">/**********************************************************************/</span></div><div class="line"><a name="l00585"></a><span class="lineno">  585</span>&#160;<span class="comment">/*************  Check if Target P32 Has Been Met  *********************/</span></div><div class="line"><a name="l00598"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#a262d50919e2ca1c4e8b32ec027a52b57">  598</a></span>&#160;<span class="keywordtype">bool</span> <a class="code" href="insert_shape_8cpp.html#a262d50919e2ca1c4e8b32ec027a52b57">p32Complete</a>(<span class="keywordtype">int</span> size) {</div><div class="line"><a name="l00599"></a><span class="lineno">  599</span>&#160;</div><div class="line"><a name="l00600"></a><span class="lineno">  600</span>&#160;    <span class="comment">// Check if p32Status array is all 1&#39;s, if not return 0</span></div><div class="line"><a name="l00601"></a><span class="lineno">  601</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; size; i++) {</div><div class="line"><a name="l00602"></a><span class="lineno">  602</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="input_8h.html#a9d32f63bd011a629cf644c120e07eda8">p32Status</a>[i] == 0) {</div><div class="line"><a name="l00603"></a><span class="lineno">  603</span>&#160;            <span class="keywordflow">return</span> 0;</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="comment">// If function has not returned yet, array is all 1&#39;s</span></div><div class="line"><a name="l00607"></a><span class="lineno">  607</span>&#160;    <span class="keywordflow">return</span> 1;</div><div class="line"><a name="l00608"></a><span class="lineno">  608</span>&#160;}</div><div class="line"><a name="l00609"></a><span class="lineno">  609</span>&#160;</div><div class="line"><a name="l00610"></a><span class="lineno">  610</span>&#160;</div><div class="line"><a name="l00611"></a><span class="lineno">  611</span>&#160;<span class="comment">/***************************************************************************/</span></div><div class="line"><a name="l00612"></a><span class="lineno">  612</span>&#160;<span class="comment">/**********************  Print Rejection Reson  ****************************/</span></div><div class="line"><a name="l00617"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#ae9670e906ee5dc43ce3fbc0ce7290dc7">  617</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="insert_shape_8cpp.html#ae9670e906ee5dc43ce3fbc0ce7290dc7">printRejectReason</a>(<span class="keywordtype">int</span> rejectCode, <span class="keyword">struct</span> <a class="code" href="struct_poly.html">Poly</a> newPoly) {</div><div class="line"><a name="l00618"></a><span class="lineno">  618</span>&#160;    </div><div class="line"><a name="l00619"></a><span class="lineno">  619</span>&#160;    <span class="keywordflow">if</span> (newPoly.<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> &gt;=0 ) {</div><div class="line"><a name="l00620"></a><span class="lineno">  620</span>&#160;        std::cout &lt;&lt; <span class="stringliteral">&quot;\nAttempted fracture from family &quot;</span> </div><div class="line"><a name="l00621"></a><span class="lineno">  621</span>&#160;                  &lt;&lt;  newPoly.<a class="code" href="struct_poly.html#ae211dd5368d3b023abc943dc940837af">familyNum</a> &lt;&lt; <span class="stringliteral">&quot; was rejected:\n&quot;</span>;</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;    <span class="keywordflow">switch</span> (rejectCode) {</div><div class="line"><a name="l00625"></a><span class="lineno">  625</span>&#160;        <span class="keywordflow">case</span> -2: </div><div class="line"><a name="l00626"></a><span class="lineno">  626</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -2: Intersection of length &lt; h.\n&quot;</span>;</div><div class="line"><a name="l00627"></a><span class="lineno">  627</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00628"></a><span class="lineno">  628</span>&#160;            </div><div class="line"><a name="l00629"></a><span class="lineno">  629</span>&#160;        <span class="keywordflow">case</span> -1:</div><div class="line"><a name="l00630"></a><span class="lineno">  630</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -1: Fracture too close to a node.\n&quot;</span>;            </div><div class="line"><a name="l00631"></a><span class="lineno">  631</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00632"></a><span class="lineno">  632</span>&#160;            </div><div class="line"><a name="l00633"></a><span class="lineno">  633</span>&#160;        <span class="keywordflow">case</span> -6:</div><div class="line"><a name="l00634"></a><span class="lineno">  634</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -6: Fracture too close &quot;</span></div><div class="line"><a name="l00635"></a><span class="lineno">  635</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;to another fracture&#39;s edge.\n&quot;</span>;</div><div class="line"><a name="l00636"></a><span class="lineno">  636</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00637"></a><span class="lineno">  637</span>&#160;        <span class="keywordflow">case</span> -7:</div><div class="line"><a name="l00638"></a><span class="lineno">  638</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -7: Fractures intersecting on same plane\n&quot;</span>;</div><div class="line"><a name="l00639"></a><span class="lineno">  639</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00640"></a><span class="lineno">  640</span>&#160;</div><div class="line"><a name="l00641"></a><span class="lineno">  641</span>&#160;        <span class="keywordflow">case</span> -10:</div><div class="line"><a name="l00642"></a><span class="lineno">  642</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -10: Rejected triple intersection &quot;</span> </div><div class="line"><a name="l00643"></a><span class="lineno">  643</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;due to triple intersections being turned off in input file.\n&quot;</span>;</div><div class="line"><a name="l00644"></a><span class="lineno">  644</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00645"></a><span class="lineno">  645</span>&#160;            </div><div class="line"><a name="l00646"></a><span class="lineno">  646</span>&#160;        <span class="keywordflow">case</span> -11:</div><div class="line"><a name="l00647"></a><span class="lineno">  647</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -11: Fracture&#39;s intersection &quot;</span> </div><div class="line"><a name="l00648"></a><span class="lineno">  648</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;landed too close to a previous intersection.\n&quot;</span>;</div><div class="line"><a name="l00649"></a><span class="lineno">  649</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00650"></a><span class="lineno">  650</span>&#160;        </div><div class="line"><a name="l00651"></a><span class="lineno">  651</span>&#160;        <span class="keywordflow">case</span> -12:</div><div class="line"><a name="l00652"></a><span class="lineno">  652</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -12: Fracture created a triple &quot;</span> </div><div class="line"><a name="l00653"></a><span class="lineno">  653</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;intersection with an angle too small.\n&quot;</span>;</div><div class="line"><a name="l00654"></a><span class="lineno">  654</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00655"></a><span class="lineno">  655</span>&#160;            </div><div class="line"><a name="l00656"></a><span class="lineno">  656</span>&#160;        <span class="keywordflow">case</span> -13:</div><div class="line"><a name="l00657"></a><span class="lineno">  657</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -13: Fracture created a triple intersection &quot;</span> </div><div class="line"><a name="l00658"></a><span class="lineno">  658</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;with the triple intersection point too close to an intersection&#39;s endpoint.\n&quot;</span>;</div><div class="line"><a name="l00659"></a><span class="lineno">  659</span>&#160;            <span class="keywordflow">break</span>;</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;        <span class="keywordflow">case</span> -14:</div><div class="line"><a name="l00662"></a><span class="lineno">  662</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = -14: Fracture created a triple intersection with &quot;</span> </div><div class="line"><a name="l00663"></a><span class="lineno">  663</span>&#160;                      &lt;&lt; <span class="stringliteral">&quot;the triple intersection point too close to another triple intersection point.\n&quot;</span>;</div><div class="line"><a name="l00664"></a><span class="lineno">  664</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00665"></a><span class="lineno">  665</span>&#160;            </div><div class="line"><a name="l00666"></a><span class="lineno">  666</span>&#160;        <span class="keywordflow">default</span>:</div><div class="line"><a name="l00667"></a><span class="lineno">  667</span>&#160;            std::cout &lt;&lt; <span class="stringliteral">&quot;\trejectCode = &quot;</span> &lt;&lt; rejectCode &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;</div><div class="line"><a name="l00668"></a><span class="lineno">  668</span>&#160;            <span class="keywordflow">break</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;}</div><div class="line"><a name="l00671"></a><span class="lineno">  671</span>&#160;</div><div class="line"><a name="l00672"></a><span class="lineno">  672</span>&#160;</div><div class="line"><a name="l00673"></a><span class="lineno">  673</span>&#160;<span class="comment">/******************************************************************/</span></div><div class="line"><a name="l00674"></a><span class="lineno">  674</span>&#160;<span class="comment">/*********************  Get Family Number *************************/</span></div><div class="line"><a name="l00688"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#ab6f3546e27a00cb14ce3c9ccf09f1f65">  688</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="insert_shape_8cpp.html#a9b3f632508ac90991bd6ebee7bcfbc2f">getFamilyNumber</a>(<span class="keywordtype">int</span> familyIndex, <span class="keywordtype">int</span> familyShape) {</div><div class="line"><a name="l00689"></a><span class="lineno">  689</span>&#160;    <span class="keywordflow">if</span> (familyShape != 0) { <span class="comment">// if not ellipse family</span></div><div class="line"><a name="l00690"></a><span class="lineno">  690</span>&#160;        <span class="keywordflow">return</span> familyIndex - <a class="code" href="input_8h.html#aaf09f7cf249c6fd6247821bc01ead90c">nFamEll</a> + 1;</div><div class="line"><a name="l00691"></a><span class="lineno">  691</span>&#160;    }</div><div class="line"><a name="l00692"></a><span class="lineno">  692</span>&#160;    <span class="keywordflow">else</span> <span class="keywordflow">return</span> familyIndex+1;</div><div class="line"><a name="l00693"></a><span class="lineno">  693</span>&#160;}</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;</div><div class="line"><a name="l00696"></a><span class="lineno">  696</span>&#160;<span class="comment">/******************************************************************/</span></div><div class="line"><a name="l00697"></a><span class="lineno">  697</span>&#160;<span class="comment">/********************  Print Shape Type  **************************/</span></div><div class="line"><a name="l00700"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#ae4fa66ac2a60ae519ee05168b44c9917">  700</a></span>&#160;std::string <a class="code" href="insert_shape_8cpp.html#ae4fa66ac2a60ae519ee05168b44c9917">shapeType</a>(<span class="keyword">struct</span> <a class="code" href="struct_shape.html">Shape</a> &amp;shapeFam) {</div><div class="line"><a name="l00701"></a><span class="lineno">  701</span>&#160;    <span class="keywordflow">if</span> (shapeFam.<a class="code" href="struct_shape.html#a766d04a69991e931ba2176a927c747cb">shapeFamily</a> == 0) {</div><div class="line"><a name="l00702"></a><span class="lineno">  702</span>&#160;        <span class="keywordflow">return</span> <span class="stringliteral">&quot;Ellipse&quot;</span>;</div><div class="line"><a name="l00703"></a><span class="lineno">  703</span>&#160;    }</div><div class="line"><a name="l00704"></a><span class="lineno">  704</span>&#160;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00705"></a><span class="lineno">  705</span>&#160;        <span class="keywordflow">return</span> <span class="stringliteral">&quot;Rectangular&quot;</span>;</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;</div><div class="line"><a name="l00709"></a><span class="lineno">  709</span>&#160;<span class="comment">/******************************************************************/</span></div><div class="line"><a name="l00710"></a><span class="lineno">  710</span>&#160;<span class="comment">/************  Get Max Fracture Radii From Family  ****************/</span></div><div class="line"><a name="l00717"></a><span class="lineno"><a class="line" href="insert_shape_8h.html#aa49097ed959c129cc9c8d11c6053ec71">  717</a></span>&#160;<span class="keywordtype">double</span> <a class="code" href="insert_shape_8cpp.html#aa49097ed959c129cc9c8d11c6053ec71">getLargestFractureRadius</a>(<a class="code" href="struct_shape.html">Shape</a> &amp;shapeFam) {</div><div class="line"><a name="l00718"></a><span class="lineno">  718</span>&#160;    <span class="keywordflow">switch</span> (shapeFam.<a class="code" href="struct_shape.html#a01b965e10318731a32fa68f0c87c9f20">distributionType</a>) {</div><div class="line"><a name="l00719"></a><span class="lineno">  719</span>&#160;        <span class="keywordflow">case</span> 1: <span class="comment">// Log-normal</span></div><div class="line"><a name="l00720"></a><span class="lineno">  720</span>&#160;                <span class="keywordflow">return</span> shapeFam.<a class="code" href="struct_shape.html#a58884d0543597d129cd6c118d2c2ad61">logMax</a>;</div><div class="line"><a name="l00721"></a><span class="lineno">  721</span>&#160;            <span class="keywordflow">break</span>;</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;        <span class="keywordflow">case</span> 2: <span class="comment">// Power-law</span></div><div class="line"><a name="l00724"></a><span class="lineno">  724</span>&#160;                <span class="keywordflow">return</span> shapeFam.<a class="code" href="struct_shape.html#ac674a8084da0711a93fb48d976467697">max</a>;</div><div class="line"><a name="l00725"></a><span class="lineno">  725</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00726"></a><span class="lineno">  726</span>&#160;</div><div class="line"><a name="l00727"></a><span class="lineno">  727</span>&#160;        <span class="keywordflow">case</span> 3: <span class="comment">// Exponential</span></div><div class="line"><a name="l00728"></a><span class="lineno">  728</span>&#160;                <span class="keywordflow">return</span> shapeFam.<a class="code" href="struct_shape.html#adcf523f838d1c624bd9d13931a342e29">expMax</a>;</div><div class="line"><a name="l00729"></a><span class="lineno">  729</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00730"></a><span class="lineno">  730</span>&#160;        <span class="keywordflow">default</span>: <span class="comment">// Constant</span></div><div class="line"><a name="l00731"></a><span class="lineno">  731</span>&#160;                <span class="keywordflow">return</span> shapeFam.<a class="code" href="struct_shape.html#ae1d2c2e5ad0273aed7474f3c4f7b8b7c">constRadi</a>;</div><div class="line"><a name="l00732"></a><span class="lineno">  732</span>&#160;</div><div class="line"><a name="l00733"></a><span class="lineno">  733</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00734"></a><span class="lineno">  734</span>&#160;    }</div><div class="line"><a name="l00735"></a><span class="lineno">  735</span>&#160;}</div><div class="line"><a name="l00736"></a><span class="lineno">  736</span>&#160;</div><div class="line"><a name="l00737"></a><span class="lineno">  737</span>&#160;</div><div class="ttc" id="input_8h_html_aafcf8c3d5090c1cda444c5a79854b434"><div class="ttname"><a href="input_8h.html#aafcf8c3d5090c1cda444c5a79854b434">apertureFromTransmissivity</a></div><div class="ttdeci">float apertureFromTransmissivity[2]</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00420">readInput.cpp:420</a></div></div>
<div class="ttc" id="computational_geometry_8cpp_html_a3098cccf29c7a47d2eae661b0f545224"><div class="ttname"><a href="computational_geometry_8cpp.html#a3098cccf29c7a47d2eae661b0f545224">translate</a></div><div class="ttdeci">void translate(Poly &amp;newPoly, double *translation)</div><div class="ttdef"><b>Definition:</b> <a href="computational_geometry_8cpp_source.html#l00047">computationalGeometry.cpp:47</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="struct_poly_html_add8e39b7f48274757f46c80014aaeb6b"><div class="ttname"><a href="struct_poly.html#add8e39b7f48274757f46c80014aaeb6b">Poly::groupNum</a></div><div class="ttdeci">unsigned int groupNum</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00030">structures.h:30</a></div></div>
<div class="ttc" id="vector_functions_8h_html_a92da2058f3ed2a6db16d2fe0e445794a"><div class="ttname"><a href="vector_functions_8h.html#a92da2058f3ed2a6db16d2fe0e445794a">magnitude</a></div><div class="ttdeci">T magnitude(T x, T y, T z)</div><div class="ttdef"><b>Definition:</b> <a href="vector_functions_8h_source.html#l00086">vectorFunctions.h:86</a></div></div>
<div class="ttc" id="vector_functions_8h_html_a48bb00fc4195377b229efe044f3bb67e"><div class="ttname"><a href="vector_functions_8h.html#a48bb00fc4195377b229efe044f3bb67e">normalize</a></div><div class="ttdeci">void normalize(T *vec)</div><div class="ttdef"><b>Definition:</b> <a href="vector_functions_8h_source.html#l00036">vectorFunctions.h:36</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="input_8h_html_aaa9154067587782d6e6eefc9652559a6"><div class="ttname"><a href="input_8h.html#aaa9154067587782d6e6eefc9652559a6">domainSizeIncrease</a></div><div class="ttdeci">float domainSizeIncrease[3]</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00111">readInput.cpp:111</a></div></div>
<div class="ttc" id="struct_shape_html_ac674a8084da0711a93fb48d976467697"><div class="ttname"><a href="struct_shape.html#ac674a8084da0711a93fb48d976467697">Shape::max</a></div><div class="ttdeci">float max</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00502">structures.h:502</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="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_poly_html_a8cc8d6067e940d33a1b29844b3e8c6ec"><div class="ttname"><a href="struct_poly.html#a8cc8d6067e940d33a1b29844b3e8c6ec">Poly::numberOfNodes</a></div><div class="ttdeci">int numberOfNodes</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00014">structures.h:14</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_shape_html_ae1d2c2e5ad0273aed7474f3c4f7b8b7c"><div class="ttname"><a href="struct_shape.html#ae1d2c2e5ad0273aed7474f3c4f7b8b7c">Shape::constRadi</a></div><div class="ttdeci">float constRadi</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00496">structures.h:496</a></div></div>
<div class="ttc" id="struct_shape_html_a874e4756ad81f2e6f135e066bb0b4347"><div class="ttname"><a href="struct_shape.html#a874e4756ad81f2e6f135e066bb0b4347">Shape::betaDistribution</a></div><div class="ttdeci">bool betaDistribution</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00431">structures.h:431</a></div></div>
<div class="ttc" id="struct_shape_html_a0d1fd1363bacd61560809dc379d17363"><div class="ttname"><a href="struct_shape.html#a0d1fd1363bacd61560809dc379d17363">Shape::thetaList</a></div><div class="ttdeci">float * thetaList</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00402">structures.h:402</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="vector_functions_8h_html"><div class="ttname"><a href="vector_functions_8h.html">vectorFunctions.h</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="struct_shape_html_a0ce256fc0e226ea6559786a2ae9d007f"><div class="ttname"><a href="struct_shape.html#a0ce256fc0e226ea6559786a2ae9d007f">Shape::layer</a></div><div class="ttdeci">short layer</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00412">structures.h:412</a></div></div>
<div class="ttc" id="generating_points_8cpp_html_a91e58d3a0d7a5e9d65cf15471f02fab0"><div class="ttname"><a href="generating_points_8cpp.html#a91e58d3a0d7a5e9d65cf15471f02fab0">truncatedPowerLaw</a></div><div class="ttdeci">float truncatedPowerLaw(float randomNum, float min, float max, float alpha)</div><div class="ttdef"><b>Definition:</b> <a href="generating_points_8cpp_source.html#l00208">generatingPoints.cpp:208</a></div></div>
<div class="ttc" id="insert_shape_8cpp_html_aed01cfb18a8a6b4ec394783593d5765f"><div class="ttname"><a href="insert_shape_8cpp.html#aed01cfb18a8a6b4ec394783593d5765f">initializeRectVertices</a></div><div class="ttdeci">void initializeRectVertices(struct Poly &amp;newPoly, float radius, float aspectRatio)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00403">insertShape.cpp:403</a></div></div>
<div class="ttc" id="input_8h_html_a6fda5f9c90f07130c1b0001301c6c945"><div class="ttname"><a href="input_8h.html#a6fda5f9c90f07130c1b0001301c6c945">permOption</a></div><div class="ttdeci">bool permOption</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00349">readInput.cpp:349</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="input_8h_html_abc39a0cf64740d1080e5c9867191fa0e"><div class="ttname"><a href="input_8h.html#abc39a0cf64740d1080e5c9867191fa0e">aperture</a></div><div class="ttdeci">int aperture</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00413">readInput.cpp:413</a></div></div>
<div class="ttc" id="generating_points_8cpp_html_a152eb5ef32d8e1dd69785a1e05d9c8c0"><div class="ttname"><a href="generating_points_8cpp.html#a152eb5ef32d8e1dd69785a1e05d9c8c0">randomTranslation</a></div><div class="ttdeci">double * randomTranslation(std::mt19937_64 &amp;generator, float xMin, float xMax, float yMin, float yMax, float zMin, float zMax)</div><div class="ttdef"><b>Definition:</b> <a href="generating_points_8cpp_source.html#l00177">generatingPoints.cpp:177</a></div></div>
<div class="ttc" id="struct_poly_html_a38a501ec1855a34b0dbe769bf0319481"><div class="ttname"><a href="struct_poly.html#a38a501ec1855a34b0dbe769bf0319481">Poly::aspectRatio</a></div><div class="ttdeci">float aspectRatio</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00049">structures.h:49</a></div></div>
<div class="ttc" id="input_8h_html_a8ca71d341876667108c011c790ed33de"><div class="ttname"><a href="input_8h.html#a8ca71d341876667108c011c790ed33de">meanAperture</a></div><div class="ttdeci">float meanAperture</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00401">readInput.cpp:401</a></div></div>
<div class="ttc" id="input_8h_html_a0b60c5d860054f964ec05348b186e86e"><div class="ttname"><a href="input_8h.html#a0b60c5d860054f964ec05348b186e86e">stdAperture</a></div><div class="ttdeci">float stdAperture</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00405">readInput.cpp:405</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_shape_html_adcf523f838d1c624bd9d13931a342e29"><div class="ttname"><a href="struct_shape.html#adcf523f838d1c624bd9d13931a342e29">Shape::expMax</a></div><div class="ttdeci">float expMax</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00478">structures.h:478</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="computational_geometry_8cpp_html_af7f1ba8a57908229078780a2695e4e46"><div class="ttname"><a href="computational_geometry_8cpp.html#af7f1ba8a57908229078780a2695e4e46">applyRotation3D</a></div><div class="ttdeci">void applyRotation3D(Poly &amp;newPoly, double *normalB)</div><div class="ttdef"><b>Definition:</b> <a href="computational_geometry_8cpp_source.html#l00147">computationalGeometry.cpp:147</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="insert_shape_8cpp_html_a4a4d307b0a091a6159358f312371ae95"><div class="ttname"><a href="insert_shape_8cpp.html#a4a4d307b0a091a6159358f312371ae95">generatePoly_withRadius</a></div><div class="ttdeci">struct Poly generatePoly_withRadius(double radius, struct Shape &amp;shapeFam, std::mt19937_64 &amp;generator, Distributions &amp;distributions, int familyIndex)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00221">insertShape.cpp:221</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="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_ab707f54136b034587f94b633476787a0"><div class="ttname"><a href="insert_shape_8cpp.html#ab707f54136b034587f94b633476787a0">_CONSTSCALAR</a></div><div class="ttdeci">#define _CONSTSCALAR</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00329">insertShape.cpp:329</a></div></div>
<div class="ttc" id="struct_shape_html_ab48639160dae4dfe2e3a0a55ffa5a691"><div class="ttname"><a href="struct_shape.html#ab48639160dae4dfe2e3a0a55ffa5a691">Shape::beta</a></div><div class="ttdeci">float beta</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00435">structures.h:435</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_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="struct_poly_html_a48c792ab99f5f7ca554a3f77062419db"><div class="ttname"><a href="struct_poly.html#a48c792ab99f5f7ca554a3f77062419db">Poly::permeability</a></div><div class="ttdeci">double permeability</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00074">structures.h:74</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="struct_shape_html_a463c3c0ba9465daaac854eaeccdb3341"><div class="ttname"><a href="struct_shape.html#a463c3c0ba9465daaac854eaeccdb3341">Shape::numPoints</a></div><div class="ttdeci">short numPoints</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00399">structures.h:399</a></div></div>
<div class="ttc" id="struct_shape_html_a58884d0543597d129cd6c118d2c2ad61"><div class="ttname"><a href="struct_shape.html#a58884d0543597d129cd6c118d2c2ad61">Shape::logMax</a></div><div class="ttdeci">float logMax</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00493">structures.h:493</a></div></div>
<div class="ttc" id="struct_shape_html_a01b965e10318731a32fa68f0c87c9f20"><div class="ttname"><a href="struct_shape.html#a01b965e10318731a32fa68f0c87c9f20">Shape::distributionType</a></div><div class="ttdeci">short distributionType</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00396">structures.h:396</a></div></div>
<div class="ttc" id="struct_shape_html_a6126c3f7cd57df03350d4b58b844bba1"><div class="ttname"><a href="struct_shape.html#a6126c3f7cd57df03350d4b58b844bba1">Shape::aspectRatio</a></div><div class="ttdeci">float aspectRatio</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00415">structures.h:415</a></div></div>
<div class="ttc" id="input_8h_html_a48936b0f8f54a1cd2c7e83ca2392ec0a"><div class="ttname"><a href="input_8h.html#a48936b0f8f54a1cd2c7e83ca2392ec0a">constantAperture</a></div><div class="ttdeci">double constantAperture</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00423">readInput.cpp:423</a></div></div>
<div class="ttc" id="generating_points_8cpp_html_a2f7b981bd314f565407a2db0358bf81b"><div class="ttname"><a href="generating_points_8cpp.html#a2f7b981bd314f565407a2db0358bf81b">fisherDistribution</a></div><div class="ttdeci">double * fisherDistribution(double theta, double phi, double kappa, std::mt19937_64 &amp;generator)</div><div class="ttdef"><b>Definition:</b> <a href="generating_points_8cpp_source.html#l00093">generatingPoints.cpp:93</a></div></div>
<div class="ttc" id="input_8h_html_a96e535ecb2fde247beaba684f23b1c3c"><div class="ttname"><a href="input_8h.html#a96e535ecb2fde247beaba684f23b1c3c">lengthCorrelatedAperture</a></div><div class="ttdeci">double lengthCorrelatedAperture[2]</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00429">readInput.cpp:429</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="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_poly_html_aa6c373c49a4e0ab7c83a353a9f8cf518"><div class="ttname"><a href="struct_poly.html#aa6c373c49a4e0ab7c83a353a9f8cf518">Poly::intersectionIndex</a></div><div class="ttdeci">std::vector&lt; unsigned int &gt; intersectionIndex</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00098">structures.h:98</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="insert_shape_8cpp_html_aa49097ed959c129cc9c8d11c6053ec71"><div class="ttname"><a href="insert_shape_8cpp.html#aa49097ed959c129cc9c8d11c6053ec71">getLargestFractureRadius</a></div><div class="ttdeci">double getLargestFractureRadius(Shape &amp;shapeFam)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00717">insertShape.cpp:717</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="struct_poly_html_a4b92c8fd93db07ccb928298defab1651"><div class="ttname"><a href="struct_poly.html#a4b92c8fd93db07ccb928298defab1651">Poly::normal</a></div><div class="ttdeci">double normal[3]</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00055">structures.h:55</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="input_8h_html_a903896dcf39ee23cbea885abd3bc926d"><div class="ttname"><a href="input_8h.html#a903896dcf39ee23cbea885abd3bc926d">constantPermeability</a></div><div class="ttdeci">double constantPermeability</div><div class="ttdef"><b>Definition:</b> <a href="read_input_8cpp_source.html#l00432">readInput.cpp:432</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="struct_shape_html_a766d04a69991e931ba2176a927c747cb"><div class="ttname"><a href="struct_shape.html#a766d04a69991e931ba2176a927c747cb">Shape::shapeFamily</a></div><div class="ttdeci">short shapeFamily</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00393">structures.h:393</a></div></div>
<div class="ttc" id="struct_poly_html_a5b78228013be48bf6b3eb6debce59f50"><div class="ttname"><a href="struct_poly.html#a5b78228013be48bf6b3eb6debce59f50">Poly::aperture</a></div><div class="ttdeci">double aperture</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00066">structures.h:66</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="insert_shape_8cpp_html_a8ed0d23da7ff9cd1f7c4109a4c478071"><div class="ttname"><a href="insert_shape_8cpp.html#a8ed0d23da7ff9cd1f7c4109a4c478071">initializeEllVertices</a></div><div class="ttdeci">void initializeEllVertices(struct Poly &amp;newPoly, float radius, float aspectRatio, float *thetaList, int numPoints)</div><div class="ttdef"><b>Definition:</b> <a href="insert_shape_8cpp_source.html#l00432">insertShape.cpp:432</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="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="struct_shape_html"><div class="ttname"><a href="struct_shape.html">Shape</a></div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00390">structures.h:390</a></div></div>
<div class="ttc" id="struct_poly_html_ae110548a9a60aacab55ecce02d759d8d"><div class="ttname"><a href="struct_poly.html#ae110548a9a60aacab55ecce02d759d8d">Poly::faces</a></div><div class="ttdeci">bool faces[6]</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00084">structures.h:84</a></div></div>
<div class="ttc" id="computational_geometry_8cpp_html_a878eeb1b4990cde9d7d12a64f647af4d"><div class="ttname"><a href="computational_geometry_8cpp.html#a878eeb1b4990cde9d7d12a64f647af4d">applyRotation2D</a></div><div class="ttdeci">void applyRotation2D(Poly &amp;newPoly, float angle)</div><div class="ttdef"><b>Definition:</b> <a href="computational_geometry_8cpp_source.html#l00024">computationalGeometry.cpp:24</a></div></div>
<div class="ttc" id="struct_poly_html_ac93df10849f4916e195a899d2ba55485"><div class="ttname"><a href="struct_poly.html#ac93df10849f4916e195a899d2ba55485">Poly::translation</a></div><div class="ttdeci">double translation[3]</div><div class="ttdef"><b>Definition:</b> <a href="structures_8h_source.html#l00052">structures.h:52</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="insert_shape_8cpp.html">insertShape.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>
