/usr/share/doc/siscone-doc-html/html/devel/quadtree_8cpp_source.html is in siscone-doc-html 2.0.6-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 | <!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.3.1"/>
<title>SISCone: siscone/quadtree.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="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 style="padding-left: 0.5em;">
<div id="projectname">SISCone
 <span id="projectnumber">2.0.6</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.3.1 -->
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_1331b0450b50de117660f97bbea14878.html">siscone</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">quadtree.cpp</div> </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> </div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment">// File: quadtree.cpp //</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">// Description: source file for quadtree management (Cquadtree class) //</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment">// This file is part of the SISCone project. //</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment">// For more details, see http://projects.hepforge.org/siscone //</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment">// //</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment">// Copyright (c) 2006 Gavin Salam and Gregory Soyez //</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment">// //</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment">// This program is free software; you can redistribute it and/or modify //</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment">// it under the terms of the GNU General Public License as published by //</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment">// the Free Software Foundation; either version 2 of the License, or //</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">// (at your option) any later version. //</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment">// //</span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment">// This program is distributed in the hope that it will be useful, //</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment">// but WITHOUT ANY WARRANTY; without even the implied warranty of //</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment">// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the //</span></div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment">// GNU General Public License for more details. //</span></div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment">// //</span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment">// You should have received a copy of the GNU General Public License //</span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment">// along with this program; if not, write to the Free Software //</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment">// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA //</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment">// //</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment">// $Revision:: 320 $//</span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="comment">// $Date:: 2011-11-15 09:54:50 +0100 (Tue, 15 Nov 2011) $//</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="comment"></span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="preprocessor">#include "quadtree.h"</span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor">#include <math.h></span></div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="preprocessor">#include <stdio.h></span></div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor">#include <iostream></span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> </div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="keyword">namespace </span>siscone{</div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span> </div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="keyword">using namespace </span>std;</div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> </div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment">/*******************************************************************</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment"> * Cquadtree implementation *</span></div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment"> * Implementation of a 2D quadtree. *</span></div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment"> * This class implements the traditional two-dimensional quadtree. *</span></div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment"> * The elements at each node are of 'Cmomentum' type. *</span></div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"> *******************************************************************/</span></div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> </div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="comment">// default ctor</span></div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment">//--------------</span></div>
<div class="line"><a name="l00045"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a74f087ccb4668bf0da9930a783dcdb31"> 45</a></span> <a class="code" href="classsiscone_1_1Cquadtree.html#a74f087ccb4668bf0da9930a783dcdb31" title="default ctor">Cquadtree::Cquadtree</a>(){</div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  v = NULL;</div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> </div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  children[0][0] = children[0][1] = children[1][0] = children[1][1] = NULL;</div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  has_child = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span> }</div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span> </div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> </div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="comment">// ctor with initialisation (see init for details)</span></div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="comment">//--------------------------</span></div>
<div class="line"><a name="l00055"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a8787e345646364ee4acfa33290ab32b7"> 55</a></span> <a class="code" href="classsiscone_1_1Cquadtree.html#a74f087ccb4668bf0da9930a783dcdb31" title="default ctor">Cquadtree::Cquadtree</a>(<span class="keywordtype">double</span> _x, <span class="keywordtype">double</span> _y, <span class="keywordtype">double</span> _half_size_x, <span class="keywordtype">double</span> _half_size_y){</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  v = NULL;</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span> </div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  children[0][0] = children[0][1] = children[1][0] = children[1][1] = NULL;</div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  has_child = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span> </div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  init(_x, _y, _half_size_x, _half_size_y);</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span> }</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span> </div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span> </div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span> <span class="comment">// default destructor</span></div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span> <span class="comment">// at destruction, everything is destroyed except </span></div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span> <span class="comment">// physical values at the leaves</span></div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span> <span class="comment">//------------------------------------------------</span></div>
<div class="line"><a name="l00069"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a54a4da1cf2fb1dd9eb838ef940acf380"> 69</a></span> <a class="code" href="classsiscone_1_1Cquadtree.html#a54a4da1cf2fb1dd9eb838ef940acf380" title="default destructor at destruction, everything is destroyed except physical values at the leaves...">Cquadtree::~Cquadtree</a>(){</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  <span class="keywordflow">if</span> (has_child){</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  <span class="keywordflow">if</span> (v!=NULL) <span class="keyword">delete</span> v;</div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="keyword">delete</span> children[0][0];</div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="keyword">delete</span> children[0][1];</div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  <span class="keyword">delete</span> children[1][0];</div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="keyword">delete</span> children[1][1];</div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  }</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span> }</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span> </div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span> </div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span> <span class="comment"> * init the tree.</span></div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span> <span class="comment"> * By initializing the tree, we mean setting the cell parameters</span></div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span> <span class="comment"> * and preparing the object to act as a seed for a new tree.</span></div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span> <span class="comment"> * - _x x-position of the center</span></div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span> <span class="comment"> * - _y y-position of the center</span></div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span> <span class="comment"> * - half_size_x half x-size of the cell</span></div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span> <span class="comment"> * - half_size_y half y-size of the cell</span></div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span> <span class="comment"> * return 0 on success, 1 on error. Note that if the cell</span></div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> <span class="comment"> * is already filled, we return an error.</span></div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span> <span class="comment"> ******************************************************************/</span></div>
<div class="line"><a name="l00091"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#ad24a504ba563607b7f2c8b33cfd6f404"> 91</a></span> <span class="keywordtype">int</span> <a class="code" href="classsiscone_1_1Cquadtree.html#ad24a504ba563607b7f2c8b33cfd6f404" title="init the tree.">Cquadtree::init</a>(<span class="keywordtype">double</span> _x, <span class="keywordtype">double</span> _y, <span class="keywordtype">double</span> _half_size_x, <span class="keywordtype">double</span> _half_size_y){</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="keywordflow">if</span> (v!=NULL)</div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="keywordflow">return</span> 1;</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span> </div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  centre_x = _x;</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  centre_y = _y;</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  half_size_x = _half_size_x;</div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  half_size_y = _half_size_y;</div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span> </div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span> }</div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span> </div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span> </div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span> <span class="comment"> * adding a particle to the tree.</span></div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span> <span class="comment"> * This method adds one vector to the quadtree structure which </span></div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span> <span class="comment"> * is updated consequently.</span></div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span> <span class="comment"> * - v vector to add</span></div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> <span class="comment"> * return 0 on success 1 on error</span></div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span> <span class="comment"> ******************************************************************/</span></div>
<div class="line"><a name="l00111"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a5eaf290c99a3cad3685f17766c3e459e"> 111</a></span> <span class="keywordtype">int</span> <a class="code" href="classsiscone_1_1Cquadtree.html#a5eaf290c99a3cad3685f17766c3e459e" title="adding a particle to the tree.">Cquadtree::add</a>(<a class="code" href="classsiscone_1_1Cmomentum.html" title="base class for dynamic coordinates management">Cmomentum</a> *v_add){</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <span class="comment">// Description of the method:</span></div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="comment">// --------------------------</span></div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  <span class="comment">// the addition process goes as follows:</span></div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="comment">// 1. check if the cell is empty, in which case, add the particle </span></div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="comment">// here and leave.</span></div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <span class="comment">// 2. If there is a unique particle already inside,</span></div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="comment">// (a) create children</span></div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  <span class="comment">// (b) forward the existing particle to the appropriate child</span></div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <span class="comment">// 3. Add current particle to this cell and forward to the </span></div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  <span class="comment">// adequate child</span></div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  <span class="comment">// NOTE: we assume in the whole procedure that the particle is </span></div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="comment">// indeed inside the cell !</span></div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  <span class="comment">// step 1: the case of empty cells</span></div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="keywordflow">if</span> (v==NULL){</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  v = v_add;</div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  }</div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span> </div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="comment">// step 2: additional work if 1! particle already present</span></div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <span class="comment">// we use the fact that only 1-particle systems have no child</span></div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  <span class="keywordflow">if</span> (!has_child){</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="keywordtype">double</span> new_half_size_x = 0.5*half_size_x;</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="keywordtype">double</span> new_half_size_y = 0.5*half_size_y;</div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  <span class="comment">// create children</span></div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  children[0][0] = <span class="keyword">new</span> <a class="code" href="classsiscone_1_1Cquadtree.html" title="Implementation of a 2D quadtree.">Cquadtree</a>(centre_x-new_half_size_x, centre_y-new_half_size_y,</div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  new_half_size_x, new_half_size_y);</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  children[0][1] = <span class="keyword">new</span> <a class="code" href="classsiscone_1_1Cquadtree.html" title="Implementation of a 2D quadtree.">Cquadtree</a>(centre_x-new_half_size_x, centre_y+new_half_size_y,</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  new_half_size_x, new_half_size_y);</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  children[1][0] = <span class="keyword">new</span> <a class="code" href="classsiscone_1_1Cquadtree.html" title="Implementation of a 2D quadtree.">Cquadtree</a>(centre_x+new_half_size_x, centre_y-new_half_size_y,</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  new_half_size_x, new_half_size_y);</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  children[1][1] = <span class="keyword">new</span> <a class="code" href="classsiscone_1_1Cquadtree.html" title="Implementation of a 2D quadtree.">Cquadtree</a>(centre_x+new_half_size_x, centre_y+new_half_size_y,</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  new_half_size_x, new_half_size_y);</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span> </div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  has_child = <span class="keyword">true</span>;</div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span> </div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  <span class="comment">// forward to child</span></div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <span class="comment">//? The following line assumes 'true'==1 and 'false'==0</span></div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="comment">// Note: v being a single particle, eta and phi are correct</span></div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  children[v->eta>centre_x][v->phi>centre_y]->add(v);</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span> </div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  <span class="comment">// copy physical params</span></div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  v = <span class="keyword">new</span> <a class="code" href="classsiscone_1_1Cmomentum.html" title="base class for dynamic coordinates management">Cmomentum</a>(*v);</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  }</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> </div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="comment">// step 3: add new particle</span></div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  <span class="comment">// Note: v_add being a single particle, eta and phi are correct</span></div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  children[v_add-><a class="code" href="classsiscone_1_1Cmomentum.html#a5ca61b8cdd1e8225ac5da7dd1814b7b3" title="particle pseudo-rapidity">eta</a>>centre_x][v_add-><a class="code" href="classsiscone_1_1Cmomentum.html#ae16f3eee4a19315fc0ced5cde2c265a6" title="particle azimuthal angle">phi</a>>centre_y]->add(v_add);</div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  *v+=*v_add;</div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span> </div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span> }</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span> </div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span> </div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span> <span class="comment"> * circle intersection.</span></div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span> <span class="comment"> * computes the intersection with a circle of given centre and radius.</span></div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span> <span class="comment"> * The output takes the form of a quadtree with all squares included </span></div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span> <span class="comment"> * in the circle.</span></div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span> <span class="comment"> * - cx circle centre x coordinate</span></div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span> <span class="comment"> * - cy circle centre y coordinate</span></div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span> <span class="comment"> * - cR2 circle radius SQUARED</span></div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span> <span class="comment"> * return the checksum for the intersection</span></div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span> <span class="comment"> ******************************************************************/</span></div>
<div class="line"><a name="l00176"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a4533e29e559e72ed5fd0fd44d5015dc9"> 176</a></span> <a class="code" href="classsiscone_1_1Creference.html" title="references used for checksums.">Creference</a> <a class="code" href="classsiscone_1_1Cquadtree.html#a4533e29e559e72ed5fd0fd44d5015dc9" title="circle intersection.">Cquadtree::circle_intersect</a>(<span class="keywordtype">double</span> cx, <span class="keywordtype">double</span> cy, <span class="keywordtype">double</span> cR2){</div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="comment">// Description of the method:</span></div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="comment">// --------------------------</span></div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="comment">// 1. check if cell is empty => no intersection</span></div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="comment">// 2. if cell has 1! particle, check if it is inside the circle.</span></div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  <span class="comment">// If yes, add it and return, if not simply return.</span></div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="comment">// 3. check if the circle intersects the square. If not, return.</span></div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  <span class="comment">// 4. check if the square is inside the circle. </span></div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  <span class="comment">// If yes, add it to qt and return.</span></div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <span class="comment">// 5. check intersections with children.</span></div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span> </div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="comment">// step 1: if there is no particle inside te square, no reason to go further</span></div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  <span class="keywordflow">if</span> (v==NULL)</div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="keywordflow">return</span> <a class="code" href="classsiscone_1_1Creference.html" title="references used for checksums.">Creference</a>();</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span> </div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  <span class="keywordtype">double</span> dx, dy;</div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span> </div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  <span class="comment">// step 2: if there is only one particle inside the square, test if it is in</span></div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <span class="comment">// the circle, in which case return associated reference</span></div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="keywordflow">if</span> (!has_child){</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="comment">// compute the distance</span></div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="comment">// Note: v has only one particle => eta and phi are defined</span></div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  dx = cx - v->eta;</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  dy = fabs(cy - v->phi);</div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  <span class="keywordflow">if</span> (dy>M_PI) </div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  dy -= 2.0*M_PI;</div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span> </div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  <span class="comment">// test distance</span></div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="keywordflow">if</span> (dx*dx+dy*dy<cR2){</div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  <span class="keywordflow">return</span> v->ref;</div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  }</div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span> </div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <span class="keywordflow">return</span> <a class="code" href="classsiscone_1_1Creference.html" title="references used for checksums.">Creference</a>();</div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  }</div>
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span> </div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  <span class="comment">// step 3: check if there is an intersection</span></div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="comment">//double ryp, rym;</span></div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <span class="keywordtype">double</span> dx_c, dy_c;</div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span> </div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  <span class="comment">// store distance with the centre of the square</span></div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  dx_c = fabs(cx-centre_x);</div>
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  dy_c = fabs(cy-centre_y);</div>
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  <span class="keywordflow">if</span> (dy_c>M_PI) dy_c = 2.0*M_PI-dy_c;</div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span> </div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  <span class="comment">// compute (minimal) the distance (pay attention to the periodicity in phi).</span></div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  dx = dx_c-half_size_x;</div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  <span class="keywordflow">if</span> (dx<0) dx=0;</div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  dy = dy_c-half_size_y;</div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  <span class="keywordflow">if</span> (dy<0) dy=0;</div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span> </div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  <span class="comment">// check the distance </span></div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  <span class="keywordflow">if</span> (dx*dx+dy*dy>=cR2){</div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  <span class="comment">// no intersection</span></div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span>  <span class="keywordflow">return</span> <a class="code" href="classsiscone_1_1Creference.html" title="references used for checksums.">Creference</a>();</div>
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span>  }</div>
<div class="line"><a name="l00231"></a><span class="lineno"> 231</span> </div>
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span>  <span class="comment">// step 4: check if included</span></div>
<div class="line"><a name="l00233"></a><span class="lineno"> 233</span> </div>
<div class="line"><a name="l00234"></a><span class="lineno"> 234</span>  <span class="comment">// compute the (maximal) distance</span></div>
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span>  dx = dx_c+half_size_x;</div>
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  dy = dy_c+half_size_y;</div>
<div class="line"><a name="l00237"></a><span class="lineno"> 237</span>  <span class="keywordflow">if</span> (dy>M_PI) dy = M_PI;</div>
<div class="line"><a name="l00238"></a><span class="lineno"> 238</span> </div>
<div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <span class="comment">// compute the distance</span></div>
<div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  <span class="keywordflow">if</span> (dx*dx+dy*dy<cR2){</div>
<div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  <span class="keywordflow">return</span> v->ref;</div>
<div class="line"><a name="l00242"></a><span class="lineno"> 242</span>  }</div>
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span> </div>
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  <span class="comment">// step 5: the square is not fully in. Recurse to children</span></div>
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  <span class="keywordflow">return</span> children[0][0]->circle_intersect(cx, cy, cR2)</div>
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span>  + children[0][1]->circle_intersect(cx, cy, cR2)</div>
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  + children[1][0]->circle_intersect(cx, cy, cR2)</div>
<div class="line"><a name="l00248"></a><span class="lineno"> 248</span>  + children[1][1]->circle_intersect(cx, cy, cR2);</div>
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span> }</div>
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span> </div>
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span> </div>
<div class="line"><a name="l00252"></a><span class="lineno"> 252</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00253"></a><span class="lineno"> 253</span> <span class="comment"> * output a data file for drawing the grid.</span></div>
<div class="line"><a name="l00254"></a><span class="lineno"> 254</span> <span class="comment"> * This can be used to output a data file containing all the</span></div>
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span> <span class="comment"> * grid subdivisions. The file contents is as follows:</span></div>
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span> <span class="comment"> * first and second columns give center of the cell, the third </span></div>
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span> <span class="comment"> * gives the size.</span></div>
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span> <span class="comment"> * - flux opened stream to write to</span></div>
<div class="line"><a name="l00259"></a><span class="lineno"> 259</span> <span class="comment"> * return 0 on success, 1 on error</span></div>
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span> <span class="comment"> ******************************************************************/</span></div>
<div class="line"><a name="l00261"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a772c03817c41f09e2df712c94d8b74bc"> 261</a></span> <span class="keywordtype">int</span> <a class="code" href="classsiscone_1_1Cquadtree.html#a772c03817c41f09e2df712c94d8b74bc" title="output a data file for drawing the grid.">Cquadtree::save</a>(FILE *flux){</div>
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span> </div>
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  <span class="keywordflow">if</span> (flux==NULL)</div>
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  <span class="keywordflow">return</span> 1;</div>
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span> </div>
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span>  <span class="keywordflow">if</span> (has_child){</div>
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span>  fprintf(flux, <span class="stringliteral">"%e\t%e\t%e\t%e\n"</span>, centre_x, centre_y, half_size_x, half_size_y);</div>
<div class="line"><a name="l00268"></a><span class="lineno"> 268</span>  children[0][0]->save(flux);</div>
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  children[0][1]->save(flux);</div>
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  children[1][0]->save(flux);</div>
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  children[1][1]->save(flux);</div>
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  }</div>
<div class="line"><a name="l00273"></a><span class="lineno"> 273</span> </div>
<div class="line"><a name="l00274"></a><span class="lineno"> 274</span>  <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00275"></a><span class="lineno"> 275</span> }</div>
<div class="line"><a name="l00276"></a><span class="lineno"> 276</span> </div>
<div class="line"><a name="l00277"></a><span class="lineno"> 277</span> </div>
<div class="line"><a name="l00278"></a><span class="lineno"> 278</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00279"></a><span class="lineno"> 279</span> <span class="comment"> * output a data file for drawing the tree leaves.</span></div>
<div class="line"><a name="l00280"></a><span class="lineno"> 280</span> <span class="comment"> * This can be used to output a data file containing all the</span></div>
<div class="line"><a name="l00281"></a><span class="lineno"> 281</span> <span class="comment"> * tree leaves. The file contents is as follows:</span></div>
<div class="line"><a name="l00282"></a><span class="lineno"> 282</span> <span class="comment"> * first and second columns give center of the cell, the third </span></div>
<div class="line"><a name="l00283"></a><span class="lineno"> 283</span> <span class="comment"> * gives the size.</span></div>
<div class="line"><a name="l00284"></a><span class="lineno"> 284</span> <span class="comment"> * - flux opened stream to write to</span></div>
<div class="line"><a name="l00285"></a><span class="lineno"> 285</span> <span class="comment"> * return 0 on success, 1 on error</span></div>
<div class="line"><a name="l00286"></a><span class="lineno"> 286</span> <span class="comment"> ******************************************************************/</span></div>
<div class="line"><a name="l00287"></a><span class="lineno"><a class="code" href="classsiscone_1_1Cquadtree.html#a8ba89ad2c9f2d1f8a5adc150b5ea56b8"> 287</a></span> <span class="keywordtype">int</span> <a class="code" href="classsiscone_1_1Cquadtree.html#a8ba89ad2c9f2d1f8a5adc150b5ea56b8" title="output a data file for drawing the tree leaves.">Cquadtree::save_leaves</a>(FILE *flux){</div>
<div class="line"><a name="l00288"></a><span class="lineno"> 288</span> </div>
<div class="line"><a name="l00289"></a><span class="lineno"> 289</span>  <span class="keywordflow">if</span> (flux==NULL)</div>
<div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  <span class="keywordflow">return</span> 1;</div>
<div class="line"><a name="l00291"></a><span class="lineno"> 291</span> </div>
<div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  <span class="keywordflow">if</span> (has_child){</div>
<div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  <span class="keywordflow">if</span> (children[0][0]!=NULL) children[0][0]->save_leaves(flux);</div>
<div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  <span class="keywordflow">if</span> (children[0][1]!=NULL) children[0][1]->save_leaves(flux);</div>
<div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  <span class="keywordflow">if</span> (children[1][0]!=NULL) children[1][0]->save_leaves(flux);</div>
<div class="line"><a name="l00296"></a><span class="lineno"> 296</span>  <span class="keywordflow">if</span> (children[1][1]!=NULL) children[1][1]->save_leaves(flux);</div>
<div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  } <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00298"></a><span class="lineno"> 298</span>  fprintf(flux, <span class="stringliteral">"%e\t%e\t%e\t%e\n"</span>, centre_x, centre_y, half_size_x, half_size_y);</div>
<div class="line"><a name="l00299"></a><span class="lineno"> 299</span>  }</div>
<div class="line"><a name="l00300"></a><span class="lineno"> 300</span> </div>
<div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00302"></a><span class="lineno"> 302</span> }</div>
<div class="line"><a name="l00303"></a><span class="lineno"> 303</span> </div>
<div class="line"><a name="l00304"></a><span class="lineno"> 304</span> }</div>
</div><!-- fragment --></div><!-- contents -->
<address style="align: right;"><small>
The <a href="http://projects.hepforge.org/siscone" target="_top">
SISCone</a> project has been developed by
<a href="http://www.lpthe.jussieu.fr/~salam" target="_top">
Gavin Salam</a> and
<a href="http://www.theo.phys.ulg.ac.be/~soyez" target="_top">
Gregory Soyez</a><br>
Documentation generated on Mon May 6 2013 11:30:35 for SISCone by
<a href="http://www.doxygen.org/index.html">Doxygen</a> 1.8.3.1
</small>
</address>
</body>
</html>
|