report-1cpyQQ.html

Deil Christoph, 01/05/2014 11:57 PM

Download (118 KB)

 
1 
<!doctype html>
2 
<html>
3 
<head>
4 
<title>GHealpix.cpp</title>
5 
<style type="text/css">
6 
 body { color:#000000; background-color:#ffffff }
7 
 body { font-family:Helvetica, sans-serif; font-size:10pt }
8 
 h1 { font-size:14pt }
9 
 .code { border-collapse:collapse; width:100%; }
10 
 .code { font-family: "Monospace", monospace; font-size:10pt }
11 
 .code { line-height: 1.2em }
12 
 .comment { color: green; font-style: oblique }
13 
 .keyword { color: blue }
14 
 .string_literal { color: red }
15 
 .directive { color: darkmagenta }
16 
 .expansion { display: none; }
17 
 .macro:hover .expansion { display: block; border: 2px solid #FF0000; padding: 2px; background-color:#FFF0F0; font-weight: normal;   -webkit-border-radius:5px;  -webkit-box-shadow:1px 1px 7px #000; position: absolute; top: -1em; left:10em; z-index: 1 }
18 
 .macro { color: darkmagenta; background-color:LemonChiffon; position: relative }
19 
 .num { width:2.5em; padding-right:2ex; background-color:#eeeeee }
20 
 .num { text-align:right; font-size:8pt }
21 
 .num { color:#444444 }
22 
 .line { padding-left: 1ex; border-left: 3px solid #ccc }
23 
 .line { white-space: pre }
24 
 .msg { -webkit-box-shadow:1px 1px 7px #000 }
25 
 .msg { -webkit-border-radius:5px }
26 
 .msg { font-family:Helvetica, sans-serif; font-size:8pt }
27 
 .msg { float:left }
28 
 .msg { padding:0.25em 1ex 0.25em 1ex }
29 
 .msg { margin-top:10px; margin-bottom:10px }
30 
 .msg { font-weight:bold }
31 
 .msg { max-width:60em; word-wrap: break-word; white-space: pre-wrap }
32 
 .msgT { padding:0x; spacing:0x }
33 
 .msgEvent { background-color:#fff8b4; color:#000000 }
34 
 .msgControl { background-color:#bbbbbb; color:#000000 }
35 
 .mrange { background-color:#dfddf3 }
36 
 .mrange { border-bottom:1px solid #6F9DBE }
37 
 .PathIndex { font-weight: bold; padding:0px 5px; margin-right:5px; }
38 
 .PathIndex { -webkit-border-radius:8px }
39 
 .PathIndexEvent { background-color:#bfba87 }
40 
 .PathIndexControl { background-color:#8c8c8c }
41 
 .PathNav a { text-decoration:none; font-size: larger }
42 
 .CodeInsertionHint { font-weight: bold; background-color: #10dd10 }
43 
 .CodeRemovalHint { background-color:#de1010 }
44 
 .CodeRemovalHint { border-bottom:1px solid #6F9DBE }
45 
 table.simpletable {
46 
   padding: 5px;
47 
   font-size:12pt;
48 
   margin:20px;
49 
   border-collapse: collapse; border-spacing: 0px;
50 
 }
51 
 td.rowname {
52 
   text-align:right; font-weight:bold; color:#444444;
53 
   padding-right:2ex; }
54 
</style>
55 
</head>
56 
<body>
57 
<!-- BUGDESC Function call argument is an uninitialized value -->
58 

    		
  

  

    
      59
    
    
      
<!-- BUGTYPE Uninitialized argument value -->


    		
  

  

    
      60
    
    
      

    		
  

  

    
      61
    
    
      
<!-- BUGCATEGORY Logic error -->


    		
  

  

    
      62
    
    
      

    		
  

  

    
      63
    
    
      
<!-- BUGFILE /Users/deil/code/gammalib/src/sky/GHealpix.cpp -->


    		
  

  

    
      64
    
    
      

    		
  

  

    
      65
    
    
      
<!-- BUGLINE 421 -->


    		
  

  

    
      66
    
    
      

    		
  

  

    
      67
    
    
      
<!-- BUGPATHLENGTH 5 -->


    		
  

  

    
      68
    
    
      

    		
  

  

    
      69
    
    
      
<!-- BUGMETAEND -->


    		
  

  

    
      70
    
    
      
<!-- REPORTHEADER -->


    		
  

  

    
      71
    
    
      
<h3>Bug Summary</h3>


    		
  

  

    
      72
    
    
      
<table class="simpletable">


    		
  

  

    
      73
    
    
      
<tr><td class="rowname">File:</td><td>src/sky/GHealpix.cpp</td></tr>


    		
  

  

    
      74
    
    
      
<tr><td class="rowname">Location:</td><td><a href="#EndPath">line 421, column 17</a></td></tr>


    		
  

  

    
      75
    
    
      
<tr><td class="rowname">Description:</td><td>Function call argument is an uninitialized value</td></tr>


    		
  

  

    
      76
    
    
      
</table>


    		
  

  

    
      77
    
    
      
<!-- REPORTSUMMARYEXTRA -->


    		
  

  

    
      78
    
    
      
<h3>Annotated Source Code</h3>


    		
  

  

    
      79
    
    
      
<table class="code">


    		
  

  

    
      80
    
    
      
<tr><td class="num" id="LN1">1</td><td class="line"><span class='comment'>/***************************************************************************</span></td></tr>


    		
  

  

    
      81
    
    
      
<tr><td class="num" id="LN2">2</td><td class="line"> <span class='comment'>*                 GHealpix.cpp - Healpix projection class                 *</span></td></tr>


    		
  

  

    
      82
    
    
      
<tr><td class="num" id="LN3">3</td><td class="line"> <span class='comment'>* ----------------------------------------------------------------------- *</span></td></tr>


    		
  

  

    
      83
    
    
      
<tr><td class="num" id="LN4">4</td><td class="line"> <span class='comment'>*  copyright (C) 2010-2013 by Juergen Knoedlseder                         *</span></td></tr>


    		
  

  

    
      84
    
    
      
<tr><td class="num" id="LN5">5</td><td class="line"> <span class='comment'>* ----------------------------------------------------------------------- *</span></td></tr>


    		
  

  

    
      85
    
    
      
<tr><td class="num" id="LN6">6</td><td class="line"> <span class='comment'>*                                                                         *</span></td></tr>


    		
  

  

    
      86
    
    
      
<tr><td class="num" id="LN7">7</td><td class="line"> <span class='comment'>*  This program is free software: you can redistribute it and/or modify   *</span></td></tr>


    		
  

  

    
      87
    
    
      
<tr><td class="num" id="LN8">8</td><td class="line"> <span class='comment'>*  it under the terms of the GNU General Public License as published by   *</span></td></tr>


    		
  

  

    
      88
    
    
      
<tr><td class="num" id="LN9">9</td><td class="line"> <span class='comment'>*  the Free Software Foundation, either version 3 of the License, or      *</span></td></tr>


    		
  

  

    
      89
    
    
      
<tr><td class="num" id="LN10">10</td><td class="line"> <span class='comment'>*  (at your option) any later version.                                    *</span></td></tr>


    		
  

  

    
      90
    
    
      
<tr><td class="num" id="LN11">11</td><td class="line"> <span class='comment'>*                                                                         *</span></td></tr>


    		
  

  

    
      91
    
    
      
<tr><td class="num" id="LN12">12</td><td class="line"> <span class='comment'>*  This program is distributed in the hope that it will be useful,        *</span></td></tr>


    		
  

  

    
      92
    
    
      
<tr><td class="num" id="LN13">13</td><td class="line"> <span class='comment'>*  but WITHOUT ANY WARRANTY; without even the implied warranty of         *</span></td></tr>


    		
  

  

    
      93
    
    
      
<tr><td class="num" id="LN14">14</td><td class="line"> <span class='comment'>*  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          *</span></td></tr>


    		
  

  

    
      94
    
    
      
<tr><td class="num" id="LN15">15</td><td class="line"> <span class='comment'>*  GNU General Public License for more details.                           *</span></td></tr>


    		
  

  

    
      95
    
    
      
<tr><td class="num" id="LN16">16</td><td class="line"> <span class='comment'>*                                                                         *</span></td></tr>


    		
  

  

    
      96
    
    
      
<tr><td class="num" id="LN17">17</td><td class="line"> <span class='comment'>*  You should have received a copy of the GNU General Public License      *</span></td></tr>


    		
  

  

    
      97
    
    
      
<tr><td class="num" id="LN18">18</td><td class="line"> <span class='comment'>*  along with this program.  If not, see &lt;http://www.gnu.org/licenses/&gt;.  *</span></td></tr>


    		
  

  

    
      98
    
    
      
<tr><td class="num" id="LN19">19</td><td class="line"> <span class='comment'>*                                                                         *</span></td></tr>


    		
  

  

    
      99
    
    
      
<tr><td class="num" id="LN20">20</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      100
    
    
      
<tr><td class="num" id="LN21">21</td><td class="line"><span class='comment'>/**</span></td></tr>


    		
  

  

    
      101
    
    
      
<tr><td class="num" id="LN22">22</td><td class="line"> <span class='comment'>* @file GHealpix.cpp</span></td></tr>


    		
  

  

    
      102
    
    
      
<tr><td class="num" id="LN23">23</td><td class="line"> <span class='comment'>* @brief HealPix projection class implementation</span></td></tr>


    		
  

  

    
      103
    
    
      
<tr><td class="num" id="LN24">24</td><td class="line"> <span class='comment'>* @author Juergen Knoedlseder</span></td></tr>


    		
  

  

    
      104
    
    
      
<tr><td class="num" id="LN25">25</td><td class="line"> <span class='comment'>*/</span></td></tr>


    		
  

  

    
      105
    
    
      
<tr><td class="num" id="LN26">26</td><td class="line"> </td></tr>


    		
  

  

    
      106
    
    
      
<tr><td class="num" id="LN27">27</td><td class="line"><span class='comment'>/* __ Includes ___________________________________________________________ */</span></td></tr>


    		
  

  

    
      107
    
    
      
<tr><td class="num" id="LN28">28</td><td class="line"><span class='directive'>#ifdef <span class='macro'>HAVE_CONFIG_H<span class='expansion'>1</span></span></span></td></tr>


    		
  

  

    
      108
    
    
      
<tr><td class="num" id="LN29">29</td><td class="line"><span class='directive'>#include &lt;config.h&gt;</span></td></tr>


    		
  

  

    
      109
    
    
      
<tr><td class="num" id="LN30">30</td><td class="line"><span class='directive'>#endif</span></td></tr>


    		
  

  

    
      110
    
    
      
<tr><td class="num" id="LN31">31</td><td class="line"><span class='directive'>#include &lt;cmath&gt;</span></td></tr>


    		
  

  

    
      111
    
    
      
<tr><td class="num" id="LN32">32</td><td class="line"><span class='directive'>#include "GException.hpp"</span></td></tr>


    		
  

  

    
      112
    
    
      
<tr><td class="num" id="LN33">33</td><td class="line"><span class='directive'>#include "GTools.hpp"</span></td></tr>


    		
  

  

    
      113
    
    
      
<tr><td class="num" id="LN34">34</td><td class="line"><span class='directive'>#include "GMath.hpp"</span></td></tr>


    		
  

  

    
      114
    
    
      
<tr><td class="num" id="LN35">35</td><td class="line"><span class='directive'>#include "GHealpix.hpp"</span></td></tr>


    		
  

  

    
      115
    
    
      
<tr><td class="num" id="LN36">36</td><td class="line"> </td></tr>


    		
  

  

    
      116
    
    
      
<tr><td class="num" id="LN37">37</td><td class="line"><span class='comment'>/* __ Method name definitions ____________________________________________ */</span></td></tr>


    		
  

  

    
      117
    
    
      
<tr><td class="num" id="LN38">38</td><td class="line"><span class='directive'>#define <span class='macro'>G_CONSTRUCT<span class='expansion'>"GHealpix::GHealpix(int&amp; ,std::string&amp; ,std::string&amp;)"</span></span>    "GHealpix::GHealpix(int&amp; ,std::string&amp; ,std::string&amp;)"</span></td></tr>


    		
  

  

    
      118
    
    
      
<tr><td class="num" id="LN39">39</td><td class="line"><span class='directive'>#define <span class='macro'>G_READ<span class='expansion'>"GHealpix::read(GFitsHDU&amp;)"</span></span>                                    "GHealpix::read(GFitsHDU&amp;)"</span></td></tr>


    		
  

  

    
      119
    
    
      
<tr><td class="num" id="LN40">40</td><td class="line"><span class='directive'>#define <span class='macro'>G_XY2DIR<span class='expansion'>"GHealpix::xy2dir(GSkyPixel&amp;)"</span></span>                               "GHealpix::xy2dir(GSkyPixel&amp;)"</span></td></tr>


    		
  

  

    
      120
    
    
      
<tr><td class="num" id="LN41">41</td><td class="line"><span class='directive'>#define <span class='macro'>G_DIR2XY2<span class='expansion'>"GHealpix::dir2xy(GSkyDir&amp;)"</span></span>                                "GHealpix::dir2xy(GSkyDir&amp;)"</span></td></tr>


    		
  

  

    
      121
    
    
      
<tr><td class="num" id="LN42">42</td><td class="line"><span class='directive'>#define <span class='macro'>G_PIX2ANG_RING<span class='expansion'>"GHealpix::pix2ang_ring(int, double*, double*)"</span></span>        "GHealpix::pix2ang_ring(int, double*, double*)"</span></td></tr>


    		
  

  

    
      122
    
    
      
<tr><td class="num" id="LN43">43</td><td class="line"><span class='directive'>#define <span class='macro'>G_PIX2ANG_NEST<span class='expansion'>"GHealpix::pix2ang_nest(int, double*, double*)"</span></span>        "GHealpix::pix2ang_nest(int, double*, double*)"</span></td></tr>


    		
  

  

    
      123
    
    
      
<tr><td class="num" id="LN44">44</td><td class="line"><span class='directive'>#define <span class='macro'>G_ORDERING_SET<span class='expansion'>"GHealpix::ordering(std::string&amp;)"</span></span>                     "GHealpix::ordering(std::string&amp;)"</span></td></tr>


    		
  

  

    
      124
    
    
      
<tr><td class="num" id="LN45">45</td><td class="line"> </td></tr>


    		
  

  

    
      125
    
    
      
<tr><td class="num" id="LN46">46</td><td class="line"><span class='comment'>/* __ Macros _____________________________________________________________ */</span></td></tr>


    		
  

  

    
      126
    
    
      
<tr><td class="num" id="LN47">47</td><td class="line"> </td></tr>


    		
  

  

    
      127
    
    
      
<tr><td class="num" id="LN48">48</td><td class="line"><span class='comment'>/* __ Coding definitions _________________________________________________ */</span></td></tr>


    		
  

  

    
      128
    
    
      
<tr><td class="num" id="LN49">49</td><td class="line"> </td></tr>


    		
  

  

    
      129
    
    
      
<tr><td class="num" id="LN50">50</td><td class="line"><span class='comment'>/* __ Debug definitions __________________________________________________ */</span></td></tr>


    		
  

  

    
      130
    
    
      
<tr><td class="num" id="LN51">51</td><td class="line"> </td></tr>


    		
  

  

    
      131
    
    
      
<tr><td class="num" id="LN52">52</td><td class="line"><span class='comment'>/* __ Local prototypes ___________________________________________________ */</span></td></tr>


    		
  

  

    
      132
    
    
      
<tr><td class="num" id="LN53">53</td><td class="line"> </td></tr>


    		
  

  

    
      133
    
    
      
<tr><td class="num" id="LN54">54</td><td class="line"><span class='comment'>/* __ Constants __________________________________________________________ */</span></td></tr>


    		
  

  

    
      134
    
    
      
<tr><td class="num" id="LN55">55</td><td class="line"><span class='keyword'>const</span> <span class='keyword'>int</span> jrll[12]  = {2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4};</td></tr>


    		
  

  

    
      135
    
    
      
<tr><td class="num" id="LN56">56</td><td class="line"><span class='keyword'>const</span> <span class='keyword'>int</span> jpll[12]  = {1, 3, 5, 7, 0, 2, 4, 6, 1, 3, 5, 7};</td></tr>


    		
  

  

    
      136
    
    
      
<tr><td class="num" id="LN57">57</td><td class="line"><span class='keyword'>const</span> <span class='keyword'>int</span> order_max = 13;</td></tr>


    		
  

  

    
      137
    
    
      
<tr><td class="num" id="LN58">58</td><td class="line"><span class='keyword'>const</span> <span class='keyword'>int</span> ns_max    = 1 &lt;&lt; order_max;</td></tr>


    		
  

  

    
      138
    
    
      
<tr><td class="num" id="LN59">59</td><td class="line"> </td></tr>


    		
  

  

    
      139
    
    
      
<tr><td class="num" id="LN60">60</td><td class="line"><span class='comment'>/* __ Static conversion arrays ___________________________________________ */</span></td></tr>


    		
  

  

    
      140
    
    
      
<tr><td class="num" id="LN61">61</td><td class="line"><span class='keyword'>static</span> <span class='keyword'>short</span> ctab[0x100];</td></tr>


    		
  

  

    
      141
    
    
      
<tr><td class="num" id="LN62">62</td><td class="line"><span class='keyword'>static</span> <span class='keyword'>short</span> utab[0x100];</td></tr>


    		
  

  

    
      142
    
    
      
<tr><td class="num" id="LN63">63</td><td class="line"> </td></tr>


    		
  

  

    
      143
    
    
      
<tr><td class="num" id="LN64">64</td><td class="line"><span class='comment'>/* __ Globals ____________________________________________________________ */</span></td></tr>


    		
  

  

    
      144
    
    
      
<tr><td class="num" id="LN65">65</td><td class="line"> </td></tr>


    		
  

  

    
      145
    
    
      
<tr><td class="num" id="LN66">66</td><td class="line"> </td></tr>


    		
  

  

    
      146
    
    
      
<tr><td class="num" id="LN67">67</td><td class="line"><span class='comment'>/*==========================================================================</span></td></tr>


    		
  

  

    
      147
    
    
      
<tr><td class="num" id="LN68">68</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      148
    
    
      
<tr><td class="num" id="LN69">69</td><td class="line"> <span class='comment'>=                         Constructors/destructors                        =</span></td></tr>


    		
  

  

    
      149
    
    
      
<tr><td class="num" id="LN70">70</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      150
    
    
      
<tr><td class="num" id="LN71">71</td><td class="line"> <span class='comment'>==========================================================================*/</span></td></tr>


    		
  

  

    
      151
    
    
      
<tr><td class="num" id="LN72">72</td><td class="line"> </td></tr>


    		
  

  

    
      152
    
    
      
<tr><td class="num" id="LN73">73</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      153
    
    
      
<tr><td class="num" id="LN74">74</td><td class="line"> <span class='comment'>* @brief Void constructor</span></td></tr>


    		
  

  

    
      154
    
    
      
<tr><td class="num" id="LN75">75</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      155
    
    
      
<tr><td class="num" id="LN76">76</td><td class="line">GHealpix::GHealpix(<span class='keyword'>void</span>) : GSkyProjection()</td></tr>


    		
  

  

    
      156
    
    
      
<tr><td class="num" id="LN77">77</td><td class="line">{</td></tr>


    		
  

  

    
      157
    
    
      
<tr><td class="num" id="LN78">78</td><td class="line">    <span class='comment'>// Initialise class members</span></td></tr>


    		
  

  

    
      158
    
    
      
<tr><td class="num" id="LN79">79</td><td class="line">    init_members();</td></tr>


    		
  

  

    
      159
    
    
      
<tr><td class="num" id="LN80">80</td><td class="line"> </td></tr>


    		
  

  

    
      160
    
    
      
<tr><td class="num" id="LN81">81</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      161
    
    
      
<tr><td class="num" id="LN82">82</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      162
    
    
      
<tr><td class="num" id="LN83">83</td><td class="line">}</td></tr>


    		
  

  

    
      163
    
    
      
<tr><td class="num" id="LN84">84</td><td class="line"> </td></tr>


    		
  

  

    
      164
    
    
      
<tr><td class="num" id="LN85">85</td><td class="line"> </td></tr>


    		
  

  

    
      165
    
    
      
<tr><td class="num" id="LN86">86</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      166
    
    
      
<tr><td class="num" id="LN87">87</td><td class="line"> <span class='comment'>* @brief Constructor</span></td></tr>


    		
  

  

    
      167
    
    
      
<tr><td class="num" id="LN88">88</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      168
    
    
      
<tr><td class="num" id="LN89">89</td><td class="line"> <span class='comment'>* @param[in] nside Number of sides.</span></td></tr>


    		
  

  

    
      169
    
    
      
<tr><td class="num" id="LN90">90</td><td class="line"> <span class='comment'>* @param[in] order Pixel ordering ('RING' or 'NESTED').</span></td></tr>


    		
  

  

    
      170
    
    
      
<tr><td class="num" id="LN91">91</td><td class="line"> <span class='comment'>* @param[in] coords Coordinate system ('EQU' or 'GAL').</span></td></tr>


    		
  

  

    
      171
    
    
      
<tr><td class="num" id="LN92">92</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      172
    
    
      
<tr><td class="num" id="LN93">93</td><td class="line"> <span class='comment'>* @exception GException::wcs_hpx_bad_nside</span> </td></tr>


    		
  

  

    
      173
    
    
      
<tr><td class="num" id="LN94">94</td><td class="line"> <span class='comment'>*            Invalid nside parameter.</span></td></tr>


    		
  

  

    
      174
    
    
      
<tr><td class="num" id="LN95">95</td><td class="line"> <span class='comment'>* @exception GException::wcs_bad_coords</span> </td></tr>


    		
  

  

    
      175
    
    
      
<tr><td class="num" id="LN96">96</td><td class="line"> <span class='comment'>*            Invalid coordsys parameter.</span></td></tr>


    		
  

  

    
      176
    
    
      
<tr><td class="num" id="LN97">97</td><td class="line"> <span class='comment'>* @exception GException::wcs_hpx_bad_ordering</span> </td></tr>


    		
  

  

    
      177
    
    
      
<tr><td class="num" id="LN98">98</td><td class="line"> <span class='comment'>*            Invalid ordering parameter.</span></td></tr>


    		
  

  

    
      178
    
    
      
<tr><td class="num" id="LN99">99</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      179
    
    
      
<tr><td class="num" id="LN100">100</td><td class="line">GHealpix::GHealpix(<span class='keyword'>const</span> <span class='keyword'>int</span>&amp;         nside,</td></tr>


    		
  

  

    
      180
    
    
      
<tr><td class="num" id="LN101">101</td><td class="line">                   <span class='keyword'>const</span> std::string&amp; order,</td></tr>


    		
  

  

    
      181
    
    
      
<tr><td class="num" id="LN102">102</td><td class="line">                   <span class='keyword'>const</span> std::string&amp; coords) : GSkyProjection()</td></tr>


    		
  

  

    
      182
    
    
      
<tr><td class="num" id="LN103">103</td><td class="line">{</td></tr>


    		
  

  

    
      183
    
    
      
<tr><td class="num" id="LN104">104</td><td class="line">    <span class='comment'>// Initialise class members</span></td></tr>


    		
  

  

    
      184
    
    
      
<tr><td class="num" id="LN105">105</td><td class="line">    init_members();</td></tr>


    		
  

  

    
      185
    
    
      
<tr><td class="num" id="LN106">106</td><td class="line"> </td></tr>


    		
  

  

    
      186
    
    
      
<tr><td class="num" id="LN107">107</td><td class="line">    <span class='comment'>// Check nside parameter (power of 2 between 1 and 8192)</span></td></tr>


    		
  

  

    
      187
    
    
      
<tr><td class="num" id="LN108">108</td><td class="line">    <span class='keyword'>if</span> (nside2order(nside) == -1) {</td></tr>


    		
  

  

    
      188
    
    
      
<tr><td class="num" id="LN109">109</td><td class="line">        <span class='keyword'>throw</span> GException::wcs_hpx_bad_nside(<span class='macro'>G_CONSTRUCT<span class='expansion'>"GHealpix::GHealpix(int&amp; ,std::string&amp; ,std::string&amp;)"</span></span>, nside);</td></tr>


    		
  

  

    
      189
    
    
      
<tr><td class="num" id="LN110">110</td><td class="line">    }</td></tr>


    		
  

  

    
      190
    
    
      
<tr><td class="num" id="LN111">111</td><td class="line"> </td></tr>


    		
  

  

    
      191
    
    
      
<tr><td class="num" id="LN112">112</td><td class="line">    <span class='comment'>// Set coordinate system</span></td></tr>


    		
  

  

    
      192
    
    
      
<tr><td class="num" id="LN113">113</td><td class="line">    coordsys(coords);</td></tr>


    		
  

  

    
      193
    
    
      
<tr><td class="num" id="LN114">114</td><td class="line"> </td></tr>


    		
  

  

    
      194
    
    
      
<tr><td class="num" id="LN115">115</td><td class="line">    <span class='comment'>// Set pixel ordering</span></td></tr>


    		
  

  

    
      195
    
    
      
<tr><td class="num" id="LN116">116</td><td class="line">    ordering(order);</td></tr>


    		
  

  

    
      196
    
    
      
<tr><td class="num" id="LN117">117</td><td class="line"> </td></tr>


    		
  

  

    
      197
    
    
      
<tr><td class="num" id="LN118">118</td><td class="line">    <span class='comment'>// Set Healpix parameters</span></td></tr>


    		
  

  

    
      198
    
    
      
<tr><td class="num" id="LN119">119</td><td class="line">    m_nside      = nside;</td></tr>


    		
  

  

    
      199
    
    
      
<tr><td class="num" id="LN120">120</td><td class="line">    m_npface     = m_nside * m_nside;</td></tr>


    		
  

  

    
      200
    
    
      
<tr><td class="num" id="LN121">121</td><td class="line">    m_ncap       = 2 * (m_npface - m_nside);</td></tr>


    		
  

  

    
      201
    
    
      
<tr><td class="num" id="LN122">122</td><td class="line">    m_num_pixels = 12 * m_npface;</td></tr>


    		
  

  

    
      202
    
    
      
<tr><td class="num" id="LN123">123</td><td class="line">    m_fact2      = 4.0 / m_num_pixels;</td></tr>


    		
  

  

    
      203
    
    
      
<tr><td class="num" id="LN124">124</td><td class="line">    m_fact1      = 2 * m_nside * m_fact2;</td></tr>


    		
  

  

    
      204
    
    
      
<tr><td class="num" id="LN125">125</td><td class="line">    m_omega      = gammalib::fourpi / m_num_pixels;</td></tr>


    		
  

  

    
      205
    
    
      
<tr><td class="num" id="LN126">126</td><td class="line">    m_order      = nside2order(m_nside);</td></tr>


    		
  

  

    
      206
    
    
      
<tr><td class="num" id="LN127">127</td><td class="line"> </td></tr>


    		
  

  

    
      207
    
    
      
<tr><td class="num" id="LN128">128</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      208
    
    
      
<tr><td class="num" id="LN129">129</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      209
    
    
      
<tr><td class="num" id="LN130">130</td><td class="line">}</td></tr>


    		
  

  

    
      210
    
    
      
<tr><td class="num" id="LN131">131</td><td class="line"> </td></tr>


    		
  

  

    
      211
    
    
      
<tr><td class="num" id="LN132">132</td><td class="line"> </td></tr>


    		
  

  

    
      212
    
    
      
<tr><td class="num" id="LN133">133</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      213
    
    
      
<tr><td class="num" id="LN134">134</td><td class="line"> <span class='comment'>* @brief Constructor from FITS HDU table</span></td></tr>


    		
  

  

    
      214
    
    
      
<tr><td class="num" id="LN135">135</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      215
    
    
      
<tr><td class="num" id="LN136">136</td><td class="line"> <span class='comment'>* @param[in] hdu FITS HDU.</span></td></tr>


    		
  

  

    
      216
    
    
      
<tr><td class="num" id="LN137">137</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      217
    
    
      
<tr><td class="num" id="LN138">138</td><td class="line">GHealpix::GHealpix(<span class='keyword'>const</span> GFitsHDU&amp; hdu) : GSkyProjection()</td></tr>


    		
  

  

    
      218
    
    
      
<tr><td class="num" id="LN139">139</td><td class="line">{</td></tr>


    		
  

  

    
      219
    
    
      
<tr><td class="num" id="LN140">140</td><td class="line">    <span class='comment'>// Initialise class members</span></td></tr>


    		
  

  

    
      220
    
    
      
<tr><td class="num" id="LN141">141</td><td class="line">    init_members();</td></tr>


    		
  

  

    
      221
    
    
      
<tr><td class="num" id="LN142">142</td><td class="line"> </td></tr>


    		
  

  

    
      222
    
    
      
<tr><td class="num" id="LN143">143</td><td class="line">    <span class='comment'>// Read Healpix defintion from FITS table</span></td></tr>


    		
  

  

    
      223
    
    
      
<tr><td class="num" id="LN144">144</td><td class="line">    read(hdu);</td></tr>


    		
  

  

    
      224
    
    
      
<tr><td class="num" id="LN145">145</td><td class="line"> </td></tr>


    		
  

  

    
      225
    
    
      
<tr><td class="num" id="LN146">146</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      226
    
    
      
<tr><td class="num" id="LN147">147</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      227
    
    
      
<tr><td class="num" id="LN148">148</td><td class="line">}</td></tr>


    		
  

  

    
      228
    
    
      
<tr><td class="num" id="LN149">149</td><td class="line"> </td></tr>


    		
  

  

    
      229
    
    
      
<tr><td class="num" id="LN150">150</td><td class="line"> </td></tr>


    		
  

  

    
      230
    
    
      
<tr><td class="num" id="LN151">151</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      231
    
    
      
<tr><td class="num" id="LN152">152</td><td class="line"> <span class='comment'>* @brief Copy constructor</span></td></tr>


    		
  

  

    
      232
    
    
      
<tr><td class="num" id="LN153">153</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      233
    
    
      
<tr><td class="num" id="LN154">154</td><td class="line"> <span class='comment'>* @param[in] proj Healpix projection.</span></td></tr>


    		
  

  

    
      234
    
    
      
<tr><td class="num" id="LN155">155</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      235
    
    
      
<tr><td class="num" id="LN156">156</td><td class="line">GHealpix::GHealpix(<span class='keyword'>const</span> GHealpix&amp; proj) : GSkyProjection(proj)</td></tr>


    		
  

  

    
      236
    
    
      
<tr><td class="num" id="LN157">157</td><td class="line">{</td></tr>


    		
  

  

    
      237
    
    
      
<tr><td class="num" id="LN158">158</td><td class="line">    <span class='comment'>// Initialise class members for clean destruction</span></td></tr>


    		
  

  

    
      238
    
    
      
<tr><td class="num" id="LN159">159</td><td class="line">    init_members();</td></tr>


    		
  

  

    
      239
    
    
      
<tr><td class="num" id="LN160">160</td><td class="line"> </td></tr>


    		
  

  

    
      240
    
    
      
<tr><td class="num" id="LN161">161</td><td class="line">    <span class='comment'>// Copy members</span></td></tr>


    		
  

  

    
      241
    
    
      
<tr><td class="num" id="LN162">162</td><td class="line">    copy_members(proj);</td></tr>


    		
  

  

    
      242
    
    
      
<tr><td class="num" id="LN163">163</td><td class="line"> </td></tr>


    		
  

  

    
      243
    
    
      
<tr><td class="num" id="LN164">164</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      244
    
    
      
<tr><td class="num" id="LN165">165</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      245
    
    
      
<tr><td class="num" id="LN166">166</td><td class="line">}</td></tr>


    		
  

  

    
      246
    
    
      
<tr><td class="num" id="LN167">167</td><td class="line"> </td></tr>


    		
  

  

    
      247
    
    
      
<tr><td class="num" id="LN168">168</td><td class="line"> </td></tr>


    		
  

  

    
      248
    
    
      
<tr><td class="num" id="LN169">169</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      249
    
    
      
<tr><td class="num" id="LN170">170</td><td class="line"> <span class='comment'>* @brief Destructor</span></td></tr>


    		
  

  

    
      250
    
    
      
<tr><td class="num" id="LN171">171</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      251
    
    
      
<tr><td class="num" id="LN172">172</td><td class="line">GHealpix::~GHealpix()</td></tr>


    		
  

  

    
      252
    
    
      
<tr><td class="num" id="LN173">173</td><td class="line">{</td></tr>


    		
  

  

    
      253
    
    
      
<tr><td class="num" id="LN174">174</td><td class="line">    <span class='comment'>// Free members</span></td></tr>


    		
  

  

    
      254
    
    
      
<tr><td class="num" id="LN175">175</td><td class="line">    free_members();</td></tr>


    		
  

  

    
      255
    
    
      
<tr><td class="num" id="LN176">176</td><td class="line"> </td></tr>


    		
  

  

    
      256
    
    
      
<tr><td class="num" id="LN177">177</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      257
    
    
      
<tr><td class="num" id="LN178">178</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      258
    
    
      
<tr><td class="num" id="LN179">179</td><td class="line">}</td></tr>


    		
  

  

    
      259
    
    
      
<tr><td class="num" id="LN180">180</td><td class="line"> </td></tr>


    		
  

  

    
      260
    
    
      
<tr><td class="num" id="LN181">181</td><td class="line"> </td></tr>


    		
  

  

    
      261
    
    
      
<tr><td class="num" id="LN182">182</td><td class="line"><span class='comment'>/*==========================================================================</span></td></tr>


    		
  

  

    
      262
    
    
      
<tr><td class="num" id="LN183">183</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      263
    
    
      
<tr><td class="num" id="LN184">184</td><td class="line"> <span class='comment'>=                               Operators                                 =</span></td></tr>


    		
  

  

    
      264
    
    
      
<tr><td class="num" id="LN185">185</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      265
    
    
      
<tr><td class="num" id="LN186">186</td><td class="line"> <span class='comment'>==========================================================================*/</span></td></tr>


    		
  

  

    
      266
    
    
      
<tr><td class="num" id="LN187">187</td><td class="line"> </td></tr>


    		
  

  

    
      267
    
    
      
<tr><td class="num" id="LN188">188</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      268
    
    
      
<tr><td class="num" id="LN189">189</td><td class="line"> <span class='comment'>* @brief Assignment operator</span></td></tr>


    		
  

  

    
      269
    
    
      
<tr><td class="num" id="LN190">190</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      270
    
    
      
<tr><td class="num" id="LN191">191</td><td class="line"> <span class='comment'>* @param[in] proj Healpix projection.</span></td></tr>


    		
  

  

    
      271
    
    
      
<tr><td class="num" id="LN192">192</td><td class="line"> <span class='comment'>* @return Healpix projection.</span></td></tr>


    		
  

  

    
      272
    
    
      
<tr><td class="num" id="LN193">193</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      273
    
    
      
<tr><td class="num" id="LN194">194</td><td class="line">GHealpix&amp; GHealpix::<span class='keyword'>operator</span>=(<span class='keyword'>const</span> GHealpix&amp; proj)</td></tr>


    		
  

  

    
      274
    
    
      
<tr><td class="num" id="LN195">195</td><td class="line">{</td></tr>


    		
  

  

    
      275
    
    
      
<tr><td class="num" id="LN196">196</td><td class="line">    <span class='comment'>// Execute only if object is not identical</span></td></tr>


    		
  

  

    
      276
    
    
      
<tr><td class="num" id="LN197">197</td><td class="line">    <span class='keyword'>if</span> (<span class='keyword'>this</span> != &amp;proj) {</td></tr>


    		
  

  

    
      277
    
    
      
<tr><td class="num" id="LN198">198</td><td class="line"> </td></tr>


    		
  

  

    
      278
    
    
      
<tr><td class="num" id="LN199">199</td><td class="line">        <span class='comment'>// Copy base class members</span></td></tr>


    		
  

  

    
      279
    
    
      
<tr><td class="num" id="LN200">200</td><td class="line">        <span class='keyword'>this</span>-&gt;GSkyProjection::<span class='keyword'>operator</span>=(proj);</td></tr>


    		
  

  

    
      280
    
    
      
<tr><td class="num" id="LN201">201</td><td class="line"> </td></tr>


    		
  

  

    
      281
    
    
      
<tr><td class="num" id="LN202">202</td><td class="line">        <span class='comment'>// Free members</span></td></tr>


    		
  

  

    
      282
    
    
      
<tr><td class="num" id="LN203">203</td><td class="line">        free_members();</td></tr>


    		
  

  

    
      283
    
    
      
<tr><td class="num" id="LN204">204</td><td class="line"> </td></tr>


    		
  

  

    
      284
    
    
      
<tr><td class="num" id="LN205">205</td><td class="line">        <span class='comment'>// Initialise private members for clean destruction</span></td></tr>


    		
  

  

    
      285
    
    
      
<tr><td class="num" id="LN206">206</td><td class="line">        init_members();</td></tr>


    		
  

  

    
      286
    
    
      
<tr><td class="num" id="LN207">207</td><td class="line"> </td></tr>


    		
  

  

    
      287
    
    
      
<tr><td class="num" id="LN208">208</td><td class="line">        <span class='comment'>// Copy members</span></td></tr>


    		
  

  

    
      288
    
    
      
<tr><td class="num" id="LN209">209</td><td class="line">        copy_members(proj);</td></tr>


    		
  

  

    
      289
    
    
      
<tr><td class="num" id="LN210">210</td><td class="line"> </td></tr>


    		
  

  

    
      290
    
    
      
<tr><td class="num" id="LN211">211</td><td class="line">    } <span class='comment'>// endif: object was not identical</span></td></tr>


    		
  

  

    
      291
    
    
      
<tr><td class="num" id="LN212">212</td><td class="line"> </td></tr>


    		
  

  

    
      292
    
    
      
<tr><td class="num" id="LN213">213</td><td class="line">    <span class='comment'>// Return this object</span></td></tr>


    		
  

  

    
      293
    
    
      
<tr><td class="num" id="LN214">214</td><td class="line">    <span class='keyword'>return</span> *<span class='keyword'>this</span>;</td></tr>


    		
  

  

    
      294
    
    
      
<tr><td class="num" id="LN215">215</td><td class="line">}</td></tr>


    		
  

  

    
      295
    
    
      
<tr><td class="num" id="LN216">216</td><td class="line"> </td></tr>


    		
  

  

    
      296
    
    
      
<tr><td class="num" id="LN217">217</td><td class="line"> </td></tr>


    		
  

  

    
      297
    
    
      
<tr><td class="num" id="LN218">218</td><td class="line"><span class='comment'>/*==========================================================================</span></td></tr>


    		
  

  

    
      298
    
    
      
<tr><td class="num" id="LN219">219</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      299
    
    
      
<tr><td class="num" id="LN220">220</td><td class="line"> <span class='comment'>=                             Public methods                              =</span></td></tr>


    		
  

  

    
      300
    
    
      
<tr><td class="num" id="LN221">221</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      301
    
    
      
<tr><td class="num" id="LN222">222</td><td class="line"> <span class='comment'>==========================================================================*/</span></td></tr>


    		
  

  

    
      302
    
    
      
<tr><td class="num" id="LN223">223</td><td class="line"> </td></tr>


    		
  

  

    
      303
    
    
      
<tr><td class="num" id="LN224">224</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      304
    
    
      
<tr><td class="num" id="LN225">225</td><td class="line"> <span class='comment'>* @brief Clear object</span></td></tr>


    		
  

  

    
      305
    
    
      
<tr><td class="num" id="LN226">226</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      306
    
    
      
<tr><td class="num" id="LN227">227</td><td class="line"> <span class='comment'>* This method properly resets the object to an initial state.</span></td></tr>


    		
  

  

    
      307
    
    
      
<tr><td class="num" id="LN228">228</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      308
    
    
      
<tr><td class="num" id="LN229">229</td><td class="line"><span class='keyword'>void</span> GHealpix::clear(<span class='keyword'>void</span>)</td></tr>


    		
  

  

    
      309
    
    
      
<tr><td class="num" id="LN230">230</td><td class="line">{</td></tr>


    		
  

  

    
      310
    
    
      
<tr><td class="num" id="LN231">231</td><td class="line">    <span class='comment'>// Free class members (base and derived classes, derived class first)</span></td></tr>


    		
  

  

    
      311
    
    
      
<tr><td class="num" id="LN232">232</td><td class="line">    free_members();</td></tr>


    		
  

  

    
      312
    
    
      
<tr><td class="num" id="LN233">233</td><td class="line">    <span class='keyword'>this</span>-&gt;GSkyProjection::free_members();</td></tr>


    		
  

  

    
      313
    
    
      
<tr><td class="num" id="LN234">234</td><td class="line"> </td></tr>


    		
  

  

    
      314
    
    
      
<tr><td class="num" id="LN235">235</td><td class="line">    <span class='comment'>// Initialise members</span></td></tr>


    		
  

  

    
      315
    
    
      
<tr><td class="num" id="LN236">236</td><td class="line">    <span class='keyword'>this</span>-&gt;GSkyProjection::init_members();</td></tr>


    		
  

  

    
      316
    
    
      
<tr><td class="num" id="LN237">237</td><td class="line">    init_members();</td></tr>


    		
  

  

    
      317
    
    
      
<tr><td class="num" id="LN238">238</td><td class="line"> </td></tr>


    		
  

  

    
      318
    
    
      
<tr><td class="num" id="LN239">239</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      319
    
    
      
<tr><td class="num" id="LN240">240</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      320
    
    
      
<tr><td class="num" id="LN241">241</td><td class="line">}</td></tr>


    		
  

  

    
      321
    
    
      
<tr><td class="num" id="LN242">242</td><td class="line"> </td></tr>


    		
  

  

    
      322
    
    
      
<tr><td class="num" id="LN243">243</td><td class="line"> </td></tr>


    		
  

  

    
      323
    
    
      
<tr><td class="num" id="LN244">244</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      324
    
    
      
<tr><td class="num" id="LN245">245</td><td class="line"> <span class='comment'>* @brief Clone instance</span></td></tr>


    		
  

  

    
      325
    
    
      
<tr><td class="num" id="LN246">246</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      326
    
    
      
<tr><td class="num" id="LN247">247</td><td class="line">GHealpix* GHealpix::clone(<span class='keyword'>void</span>) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      327
    
    
      
<tr><td class="num" id="LN248">248</td><td class="line">{</td></tr>


    		
  

  

    
      328
    
    
      
<tr><td class="num" id="LN249">249</td><td class="line">    <span class='keyword'>return</span> <span class='keyword'>new</span> GHealpix(*<span class='keyword'>this</span>);</td></tr>


    		
  

  

    
      329
    
    
      
<tr><td class="num" id="LN250">250</td><td class="line">}</td></tr>


    		
  

  

    
      330
    
    
      
<tr><td class="num" id="LN251">251</td><td class="line"> </td></tr>


    		
  

  

    
      331
    
    
      
<tr><td class="num" id="LN252">252</td><td class="line"> </td></tr>


    		
  

  

    
      332
    
    
      
<tr><td class="num" id="LN253">253</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      333
    
    
      
<tr><td class="num" id="LN254">254</td><td class="line"> <span class='comment'>* @brief Read Healpix definiton from FITS header</span></td></tr>


    		
  

  

    
      334
    
    
      
<tr><td class="num" id="LN255">255</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      335
    
    
      
<tr><td class="num" id="LN256">256</td><td class="line"> <span class='comment'>* @param[in] hdu FITS HDU.</span></td></tr>


    		
  

  

    
      336
    
    
      
<tr><td class="num" id="LN257">257</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      337
    
    
      
<tr><td class="num" id="LN258">258</td><td class="line"> <span class='comment'>* @exception GException::wcs</span></td></tr>


    		
  

  

    
      338
    
    
      
<tr><td class="num" id="LN259">259</td><td class="line"> <span class='comment'>*            Unable to load Healpix definition from HDU.</span></td></tr>


    		
  

  

    
      339
    
    
      
<tr><td class="num" id="LN260">260</td><td class="line"> <span class='comment'>* @exception GException::wcs_bad_coords</span></td></tr>


    		
  

  

    
      340
    
    
      
<tr><td class="num" id="LN261">261</td><td class="line"> <span class='comment'>*            Invalid coordsys parameter.</span></td></tr>


    		
  

  

    
      341
    
    
      
<tr><td class="num" id="LN262">262</td><td class="line"> <span class='comment'>* @exception GException::wcs_hpx_bad_ordering</span></td></tr>


    		
  

  

    
      342
    
    
      
<tr><td class="num" id="LN263">263</td><td class="line"> <span class='comment'>*            Invalid ordering parameter.</span></td></tr>


    		
  

  

    
      343
    
    
      
<tr><td class="num" id="LN264">264</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      344
    
    
      
<tr><td class="num" id="LN265">265</td><td class="line"><span class='keyword'>void</span> GHealpix::read(<span class='keyword'>const</span> GFitsHDU&amp; hdu)</td></tr>


    		
  

  

    
      345
    
    
      
<tr><td class="num" id="LN266">266</td><td class="line">{</td></tr>


    		
  

  

    
      346
    
    
      
<tr><td class="num" id="LN267">267</td><td class="line">    <span class='comment'>// Clear object</span></td></tr>


    		
  

  

    
      347
    
    
      
<tr><td class="num" id="LN268">268</td><td class="line">    clear();</td></tr>


    		
  

  

    
      348
    
    
      
<tr><td class="num" id="LN269">269</td><td class="line"> </td></tr>


    		
  

  

    
      349
    
    
      
<tr><td class="num" id="LN270">270</td><td class="line">    <span class='comment'>// Check if we have a healpix representation</span></td></tr>


    		
  

  

    
      350
    
    
      
<tr><td class="num" id="LN271">271</td><td class="line">    <span class='keyword'>if</span> (hdu.string(<span class='string_literal'>"PIXTYPE"</span>) != <span class='string_literal'>"HEALPIX"</span>) {</td></tr>


    		
  

  

    
      351
    
    
      
<tr><td class="num" id="LN272">272</td><td class="line">        <span class='keyword'>throw</span> GException::wcs(<span class='macro'>G_READ<span class='expansion'>"GHealpix::read(GFitsHDU&amp;)"</span></span>, <span class='string_literal'>"HDU does not contain Healpix data"</span>);</td></tr>


    		
  

  

    
      352
    
    
      
<tr><td class="num" id="LN273">273</td><td class="line">    }</td></tr>


    		
  

  

    
      353
    
    
      
<tr><td class="num" id="LN274">274</td><td class="line"> </td></tr>


    		
  

  

    
      354
    
    
      
<tr><td class="num" id="LN275">275</td><td class="line">    <span class='comment'>// Get pixel ordering</span></td></tr>


    		
  

  

    
      355
    
    
      
<tr><td class="num" id="LN276">276</td><td class="line">    std::string order = hdu.string(<span class='string_literal'>"ORDERING"</span>);</td></tr>


    		
  

  

    
      356
    
    
      
<tr><td class="num" id="LN277">277</td><td class="line"> </td></tr>


    		
  

  

    
      357
    
    
      
<tr><td class="num" id="LN278">278</td><td class="line">    <span class='comment'>// Get coordinate system.</span></td></tr>


    		
  

  

    
      358
    
    
      
<tr><td class="num" id="LN279">279</td><td class="line">    <span class='comment'>// First search for HIER_CRD keyword (this has been used in older</span></td></tr>


    		
  

  

    
      359
    
    
      
<tr><td class="num" id="LN280">280</td><td class="line">    <span class='comment'>// versions of LAT exposure cubes). If not found then search for standard</span></td></tr>


    		
  

  

    
      360
    
    
      
<tr><td class="num" id="LN281">281</td><td class="line">    <span class='comment'>// COORDSYS keyword.</span></td></tr>


    		
  

  

    
      361
    
    
      
<tr><td class="num" id="LN282">282</td><td class="line">    std::string coords;</td></tr>


    		
  

  

    
      362
    
    
      
<tr><td class="num" id="LN283">283</td><td class="line">    <span class='keyword'>if</span> (hdu.hascard(<span class='string_literal'>"HIER_CRD"</span>)) {</td></tr>


    		
  

  

    
      363
    
    
      
<tr><td class="num" id="LN284">284</td><td class="line">        coords = hdu.string(<span class='string_literal'>"HIER_CRD"</span>);</td></tr>


    		
  

  

    
      364
    
    
      
<tr><td class="num" id="LN285">285</td><td class="line">    }</td></tr>


    		
  

  

    
      365
    
    
      
<tr><td class="num" id="LN286">286</td><td class="line">    <span class='keyword'>else</span> <span class='keyword'>if</span> (hdu.hascard(<span class='string_literal'>"COORDSYS"</span>)) {</td></tr>


    		
  

  

    
      366
    
    
      
<tr><td class="num" id="LN287">287</td><td class="line">        coords = hdu.string(<span class='string_literal'>"COORDSYS"</span>);</td></tr>


    		
  

  

    
      367
    
    
      
<tr><td class="num" id="LN288">288</td><td class="line">    }</td></tr>


    		
  

  

    
      368
    
    
      
<tr><td class="num" id="LN289">289</td><td class="line"> </td></tr>


    		
  

  

    
      369
    
    
      
<tr><td class="num" id="LN290">290</td><td class="line">    <span class='comment'>// Set coordinate system</span></td></tr>


    		
  

  

    
      370
    
    
      
<tr><td class="num" id="LN291">291</td><td class="line">    coordsys(coords);</td></tr>


    		
  

  

    
      371
    
    
      
<tr><td class="num" id="LN292">292</td><td class="line"> </td></tr>


    		
  

  

    
      372
    
    
      
<tr><td class="num" id="LN293">293</td><td class="line">    <span class='comment'>// Set pixel ordering</span></td></tr>


    		
  

  

    
      373
    
    
      
<tr><td class="num" id="LN294">294</td><td class="line">    ordering(order);</td></tr>


    		
  

  

    
      374
    
    
      
<tr><td class="num" id="LN295">295</td><td class="line"> </td></tr>


    		
  

  

    
      375
    
    
      
<tr><td class="num" id="LN296">296</td><td class="line">    <span class='comment'>// Get Healpix resolution and determine number of pixels and solid angle</span></td></tr>


    		
  

  

    
      376
    
    
      
<tr><td class="num" id="LN297">297</td><td class="line">    m_nside      = hdu.integer(<span class='string_literal'>"NSIDE"</span>);</td></tr>


    		
  

  

    
      377
    
    
      
<tr><td class="num" id="LN298">298</td><td class="line">    m_npface     = m_nside * m_nside;</td></tr>


    		
  

  

    
      378
    
    
      
<tr><td class="num" id="LN299">299</td><td class="line">    m_ncap       = 2 * (m_npface - m_nside);</td></tr>


    		
  

  

    
      379
    
    
      
<tr><td class="num" id="LN300">300</td><td class="line">    m_num_pixels = 12 * m_npface;</td></tr>


    		
  

  

    
      380
    
    
      
<tr><td class="num" id="LN301">301</td><td class="line">    m_fact2      = 4.0 / m_num_pixels;</td></tr>


    		
  

  

    
      381
    
    
      
<tr><td class="num" id="LN302">302</td><td class="line">    m_fact1      = 2 * m_nside * m_fact2;</td></tr>


    		
  

  

    
      382
    
    
      
<tr><td class="num" id="LN303">303</td><td class="line">    m_omega      = gammalib::fourpi / m_num_pixels;</td></tr>


    		
  

  

    
      383
    
    
      
<tr><td class="num" id="LN304">304</td><td class="line">    m_order      = nside2order(m_nside);</td></tr>


    		
  

  

    
      384
    
    
      
<tr><td class="num" id="LN305">305</td><td class="line"> </td></tr>


    		
  

  

    
      385
    
    
      
<tr><td class="num" id="LN306">306</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      386
    
    
      
<tr><td class="num" id="LN307">307</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      387
    
    
      
<tr><td class="num" id="LN308">308</td><td class="line">}</td></tr>


    		
  

  

    
      388
    
    
      
<tr><td class="num" id="LN309">309</td><td class="line"> </td></tr>


    		
  

  

    
      389
    
    
      
<tr><td class="num" id="LN310">310</td><td class="line"> </td></tr>


    		
  

  

    
      390
    
    
      
<tr><td class="num" id="LN311">311</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      391
    
    
      
<tr><td class="num" id="LN312">312</td><td class="line"> <span class='comment'>* @brief Write Healpix definiton into FITS HDU</span></td></tr>


    		
  

  

    
      392
    
    
      
<tr><td class="num" id="LN313">313</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      393
    
    
      
<tr><td class="num" id="LN314">314</td><td class="line"> <span class='comment'>* @param[in] hdu FITS HDU.</span></td></tr>


    		
  

  

    
      394
    
    
      
<tr><td class="num" id="LN315">315</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      395
    
    
      
<tr><td class="num" id="LN316">316</td><td class="line"> <span class='comment'>* Writes the following keywords in the FITS HDU:</span></td></tr>


    		
  

  

    
      396
    
    
      
<tr><td class="num" id="LN317">317</td><td class="line"> <span class='comment'>* EXTNAME  = HEALPIX</span></td></tr>


    		
  

  

    
      397
    
    
      
<tr><td class="num" id="LN318">318</td><td class="line"> <span class='comment'>* PIXTYPE  = HEALPIX</span></td></tr>


    		
  

  

    
      398
    
    
      
<tr><td class="num" id="LN319">319</td><td class="line"> <span class='comment'>* NSIDE    = nside() = m_nside</span></td></tr>


    		
  

  

    
      399
    
    
      
<tr><td class="num" id="LN320">320</td><td class="line"> <span class='comment'>* FIRSTPIX = 0</span></td></tr>


    		
  

  

    
      400
    
    
      
<tr><td class="num" id="LN321">321</td><td class="line"> <span class='comment'>* LASTPIX  = npix()-1 = m_num_pixels-1</span></td></tr>


    		
  

  

    
      401
    
    
      
<tr><td class="num" id="LN322">322</td><td class="line"> <span class='comment'>* ORDERING = ordering()</span></td></tr>


    		
  

  

    
      402
    
    
      
<tr><td class="num" id="LN323">323</td><td class="line"> <span class='comment'>* COORDSYS = coordsys()</span></td></tr>


    		
  

  

    
      403
    
    
      
<tr><td class="num" id="LN324">324</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      404
    
    
      
<tr><td class="num" id="LN325">325</td><td class="line"><span class='keyword'>void</span> GHealpix::write(GFitsHDU&amp; hdu) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      405
    
    
      
<tr><td class="num" id="LN326">326</td><td class="line">{</td></tr>


    		
  

  

    
      406
    
    
      
<tr><td class="num" id="LN327">327</td><td class="line">    <span class='comment'>// Set extension name</span></td></tr>


    		
  

  

    
      407
    
    
      
<tr><td class="num" id="LN328">328</td><td class="line">    hdu.extname(<span class='string_literal'>"HEALPIX"</span>);</td></tr>


    		
  

  

    
      408
    
    
      
<tr><td class="num" id="LN329">329</td><td class="line"> </td></tr>


    		
  

  

    
      409
    
    
      
<tr><td class="num" id="LN330">330</td><td class="line">    <span class='comment'>// Set keywords</span></td></tr>


    		
  

  

    
      410
    
    
      
<tr><td class="num" id="LN331">331</td><td class="line">    hdu.card(<span class='string_literal'>"PIXTYPE"</span>,  <span class='string_literal'>"HEALPIX"</span>,  <span class='string_literal'>"HEALPix pixelisation"</span>);</td></tr>


    		
  

  

    
      411
    
    
      
<tr><td class="num" id="LN332">332</td><td class="line">    hdu.card(<span class='string_literal'>"NSIDE"</span>,    nside(),    <span class='string_literal'>"HEALPix resolution parameter"</span>);</td></tr>


    		
  

  

    
      412
    
    
      
<tr><td class="num" id="LN333">333</td><td class="line">    hdu.card(<span class='string_literal'>"FIRSTPIX"</span>, 0,          <span class='string_literal'>"Index of first pixel"</span>);</td></tr>


    		
  

  

    
      413
    
    
      
<tr><td class="num" id="LN334">334</td><td class="line">    hdu.card(<span class='string_literal'>"LASTPIX"</span>,  npix()-1,   <span class='string_literal'>"Index of last pixel"</span>);</td></tr>


    		
  

  

    
      414
    
    
      
<tr><td class="num" id="LN335">335</td><td class="line">    hdu.card(<span class='string_literal'>"ORDERING"</span>, ordering(),</td></tr>


    		
  

  

    
      415
    
    
      
<tr><td class="num" id="LN336">336</td><td class="line">             <span class='string_literal'>"Pixel ordering scheme, either RING or NESTED"</span>);</td></tr>


    		
  

  

    
      416
    
    
      
<tr><td class="num" id="LN337">337</td><td class="line">    hdu.card(<span class='string_literal'>"COORDSYS"</span>, coordsys(),</td></tr>


    		
  

  

    
      417
    
    
      
<tr><td class="num" id="LN338">338</td><td class="line">             <span class='string_literal'>"Coordinate system, either EQU or GAL"</span>);</td></tr>


    		
  

  

    
      418
    
    
      
<tr><td class="num" id="LN339">339</td><td class="line"> </td></tr>


    		
  

  

    
      419
    
    
      
<tr><td class="num" id="LN340">340</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      420
    
    
      
<tr><td class="num" id="LN341">341</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      421
    
    
      
<tr><td class="num" id="LN342">342</td><td class="line">}</td></tr>


    		
  

  

    
      422
    
    
      
<tr><td class="num" id="LN343">343</td><td class="line"> </td></tr>


    		
  

  

    
      423
    
    
      
<tr><td class="num" id="LN344">344</td><td class="line"> </td></tr>


    		
  

  

    
      424
    
    
      
<tr><td class="num" id="LN345">345</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      425
    
    
      
<tr><td class="num" id="LN346">346</td><td class="line"> <span class='comment'>* @brief Returns sky direction of pixel</span></td></tr>


    		
  

  

    
      426
    
    
      
<tr><td class="num" id="LN347">347</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      427
    
    
      
<tr><td class="num" id="LN348">348</td><td class="line"> <span class='comment'>* @param[in] pixel Sky map pixel.</span></td></tr>


    		
  

  

    
      428
    
    
      
<tr><td class="num" id="LN349">349</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      429
    
    
      
<tr><td class="num" id="LN350">350</td><td class="line">GSkyDir GHealpix::pix2dir(<span class='keyword'>const</span> GSkyPixel&amp; pixel) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      430
    
    
      
<tr><td class="num" id="LN351">351</td><td class="line">{</td></tr>


    		
  

  

    
      431
    
    
      
<tr><td class="num" id="LN352">352</td><td class="line">    <span class='comment'>// Throw an exception if sky map pixel is not 1D</span></td></tr>


    		
  

  

    
      432
    
    
      
<tr><td class="num" id="LN353">353</td><td class="line">    <span class='keyword'>if</span> (!pixel.is1D()) {</td></tr>


    		
  

  

    
      433
    
    
      
<tr><td class="num" id="LN354">354</td><td class="line">        std::string msg = <span class='string_literal'>"Sky map pixel "</span>+pixel.print()+<span class='string_literal'>" is not"</span></td></tr>


    		
  

  

    
      434
    
    
      
<tr><td class="num" id="LN355">355</td><td class="line">                          <span class='string_literal'>" 1-dimensional.\n"</span></td></tr>


    		
  

  

    
      435
    
    
      
<tr><td class="num" id="LN356">356</td><td class="line">                          <span class='string_literal'>"Only 1-dimensional pixels are supported by the"</span></td></tr>


    		
  

  

    
      436
    
    
      
<tr><td class="num" id="LN357">357</td><td class="line">                          <span class='string_literal'>" Healpix projection."</span>;</td></tr>


    		
  

  

    
      437
    
    
      
<tr><td class="num" id="LN358">358</td><td class="line">        <span class='keyword'>throw</span> GException::invalid_argument(<span class='macro'>G_XY2DIR<span class='expansion'>"GHealpix::xy2dir(GSkyPixel&amp;)"</span></span>, msg);</td></tr>


    		
  

  

    
      438
    
    
      
<tr><td class="num" id="LN359">359</td><td class="line">    }</td></tr>


    		
  

  

    
      439
    
    
      
<tr><td class="num" id="LN360">360</td><td class="line"> </td></tr>


    		
  

  

    
      440
    
    
      
<tr><td class="num" id="LN361">361</td><td class="line">    <span class='comment'>// Perform ordering dependent conversion</span></td></tr>


    		
  

  

    
      441
    
    
      
<tr><td class="num" id="LN362">362</td><td class="line">    <span class='keyword'>double</span> theta = 0.0;</td></tr>


    		
  

  

    
      442
    
    
      
<tr><td class="num" id="LN363">363</td><td class="line">    <span class='keyword'>double</span> phi   = 0.0;</td></tr>


    		
  

  

    
      443
    
    
      
<tr><td class="num" id="LN364">364</td><td class="line">    <span class='keyword'>switch</span> (m_ordering) {</td></tr>


    		
  

  

    
      444
    
    
      
<tr><td class="num" id="LN365">365</td><td class="line">    <span class='keyword'>case</span> 0:</td></tr>


    		
  

  

    
      445
    
    
      
<tr><td class="num" id="LN366">366</td><td class="line">        pix2ang_ring(<span class='keyword'>int</span>(pixel), &amp;theta, &amp;phi);</td></tr>


    		
  

  

    
      446
    
    
      
<tr><td class="num" id="LN367">367</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      447
    
    
      
<tr><td class="num" id="LN368">368</td><td class="line">    <span class='keyword'>case</span> 1:</td></tr>


    		
  

  

    
      448
    
    
      
<tr><td class="num" id="LN369">369</td><td class="line">        pix2ang_nest(<span class='keyword'>int</span>(pixel), &amp;theta, &amp;phi);</td></tr>


    		
  

  

    
      449
    
    
      
<tr><td class="num" id="LN370">370</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      450
    
    
      
<tr><td class="num" id="LN371">371</td><td class="line">    <span class='keyword'>default</span>:</td></tr>


    		
  

  

    
      451
    
    
      
<tr><td class="num" id="LN372">372</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      452
    
    
      
<tr><td class="num" id="LN373">373</td><td class="line">    }</td></tr>


    		
  

  

    
      453
    
    
      
<tr><td class="num" id="LN374">374</td><td class="line"> </td></tr>


    		
  

  

    
      454
    
    
      
<tr><td class="num" id="LN375">375</td><td class="line">    <span class='comment'>// Store coordinate system dependent result</span></td></tr>


    		
  

  

    
      455
    
    
      
<tr><td class="num" id="LN376">376</td><td class="line">    GSkyDir result;</td></tr>


    		
  

  

    
      456
    
    
      
<tr><td class="num" id="LN377">377</td><td class="line">    <span class='keyword'>switch</span> (m_coordsys) {</td></tr>


    		
  

  

    
      457
    
    
      
<tr><td class="num" id="LN378">378</td><td class="line">    <span class='keyword'>case</span> 0:</td></tr>


    		
  

  

    
      458
    
    
      
<tr><td class="num" id="LN379">379</td><td class="line">        result.radec(phi, gammalib::pihalf - theta);</td></tr>


    		
  

  

    
      459
    
    
      
<tr><td class="num" id="LN380">380</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      460
    
    
      
<tr><td class="num" id="LN381">381</td><td class="line">    <span class='keyword'>case</span> 1:</td></tr>


    		
  

  

    
      461
    
    
      
<tr><td class="num" id="LN382">382</td><td class="line">        result.lb(phi, gammalib::pihalf - theta);</td></tr>


    		
  

  

    
      462
    
    
      
<tr><td class="num" id="LN383">383</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      463
    
    
      
<tr><td class="num" id="LN384">384</td><td class="line">    <span class='keyword'>default</span>:</td></tr>


    		
  

  

    
      464
    
    
      
<tr><td class="num" id="LN385">385</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      465
    
    
      
<tr><td class="num" id="LN386">386</td><td class="line">    }</td></tr>


    		
  

  

    
      466
    
    
      
<tr><td class="num" id="LN387">387</td><td class="line">    </td></tr>


    		
  

  

    
      467
    
    
      
<tr><td class="num" id="LN388">388</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      468
    
    
      
<tr><td class="num" id="LN389">389</td><td class="line">    <span class='keyword'>return</span> result;</td></tr>


    		
  

  

    
      469
    
    
      
<tr><td class="num" id="LN390">390</td><td class="line">}</td></tr>


    		
  

  

    
      470
    
    
      
<tr><td class="num" id="LN391">391</td><td class="line"> </td></tr>


    		
  

  

    
      471
    
    
      
<tr><td class="num" id="LN392">392</td><td class="line"> </td></tr>


    		
  

  

    
      472
    
    
      
<tr><td class="num" id="LN393">393</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      473
    
    
      
<tr><td class="num" id="LN394">394</td><td class="line"> <span class='comment'>* @brief Returns sky map pixel of sky coordinate</span></td></tr>


    		
  

  

    
      474
    
    
      
<tr><td class="num" id="LN395">395</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      475
    
    
      
<tr><td class="num" id="LN396">396</td><td class="line"> <span class='comment'>* @param[in] dir Sky coordinate.</span></td></tr>


    		
  

  

    
      476
    
    
      
<tr><td class="num" id="LN397">397</td><td class="line"> <span class='comment'>* @return Sky map pixel.</span></td></tr>


    		
  

  

    
      477
    
    
      
<tr><td class="num" id="LN398">398</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      478
    
    
      
<tr><td class="num" id="LN399">399</td><td class="line">GSkyPixel GHealpix::dir2pix(<span class='keyword'>const</span> GSkyDir&amp; dir) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      479
    
    
      
<tr><td class="num" id="LN400">400</td><td class="line">{</td></tr>


    		
  

  

    
      480
    
    
      
<tr><td class="num" id="LN401">401</td><td class="line">    <span class='comment'>// Compute coordinate system dependent (z,phi)</span></td></tr>


    		
  

  

    
      481
    
    
      
<tr><td class="num" id="LN402">402</td><td class="line">    <span class="mrange"><span class='keyword'>double</span> z</span>;</td></tr>


    		
  

  

    
      482
    
    
      
<tr><td class="num"></td><td class="line"><div id="Path1" class="msg msgEvent" style="margin-left:5ex"><table class="msgT"><tr><td valign="top"><div class="PathIndex PathIndexEvent">1</div></td><td>'z' declared without an initial value</td><td><div class="PathNav"><a href="#Path2" title="Next event (2)">&#x2192;</a></div></td></tr></table></div></td></tr>


    		
  

  

    
      483
    
    
      
<tr><td class="num" id="LN403">403</td><td class="line">    <span class='keyword'>double</span> phi;</td></tr>


    		
  

  

    
      484
    
    
      
<tr><td class="num" id="LN404">404</td><td class="line">    <span class='keyword'>switch</span> (m_coordsys) {</td></tr>


    		
  

  

    
      485
    
    
      
<tr><td class="num"></td><td class="line"><div id="Path2" class="msg msgControl" style="margin-left:5ex"><table class="msgT"><tr><td valign="top"><div class="PathIndex PathIndexControl">2</div></td><td><div class="PathNav"><a href="#Path1" title="Previous event (1)">&#x2190;</a></div></td></td><td>Control jumps to the 'default' case at line 413</td><td><div class="PathNav"><a href="#Path3" title="Next event (3)">&#x2192;</a></div></td></tr></table></div></td></tr>


    		
  

  

    
      486
    
    
      
<tr><td class="num" id="LN405">405</td><td class="line">    <span class='keyword'>case</span> 0:</td></tr>


    		
  

  

    
      487
    
    
      
<tr><td class="num" id="LN406">406</td><td class="line">        z   = cos(gammalib::pihalf - dir.dec());</td></tr>


    		
  

  

    
      488
    
    
      
<tr><td class="num" id="LN407">407</td><td class="line">        phi = dir.ra();</td></tr>


    		
  

  

    
      489
    
    
      
<tr><td class="num" id="LN408">408</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      490
    
    
      
<tr><td class="num" id="LN409">409</td><td class="line">    <span class='keyword'>case</span> 1:</td></tr>


    		
  

  

    
      491
    
    
      
<tr><td class="num" id="LN410">410</td><td class="line">        z   = cos(gammalib::pihalf - dir.b());</td></tr>


    		
  

  

    
      492
    
    
      
<tr><td class="num" id="LN411">411</td><td class="line">        phi = dir.l();</td></tr>


    		
  

  

    
      493
    
    
      
<tr><td class="num" id="LN412">412</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      494
    
    
      
<tr><td class="num" id="LN413">413</td><td class="line">    <span class='keyword'>default</span>:</td></tr>


    		
  

  

    
      495
    
    
      
<tr><td class="num" id="LN414">414</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      496
    
    
      
<tr><td class="num"></td><td class="line"><div id="Path3" class="msg msgControl" style="margin-left:9ex"><table class="msgT"><tr><td valign="top"><div class="PathIndex PathIndexControl">3</div></td><td><div class="PathNav"><a href="#Path2" title="Previous event (2)">&#x2190;</a></div></td></td><td> Execution continues on line 418</td><td><div class="PathNav"><a href="#Path4" title="Next event (4)">&#x2192;</a></div></td></tr></table></div></td></tr>


    		
  

  

    
      497
    
    
      
<tr><td class="num" id="LN415">415</td><td class="line">    }</td></tr>


    		
  

  

    
      498
    
    
      
<tr><td class="num" id="LN416">416</td><td class="line"> </td></tr>


    		
  

  

    
      499
    
    
      
<tr><td class="num" id="LN417">417</td><td class="line">    <span class='comment'>// Perform ordering dependent conversion</span></td></tr>


    		
  

  

    
      500
    
    
      
<tr><td class="num" id="LN418">418</td><td class="line">    <span class='keyword'>int</span> index;</td></tr>


    		
  

  

    
      501
    
    
      
<tr><td class="num" id="LN419">419</td><td class="line">    <span class='keyword'>switch</span> (m_ordering) {</td></tr>


    		
  

  

    
      502
    
    
      
<tr><td class="num"></td><td class="line"><div id="Path4" class="msg msgControl" style="margin-left:5ex"><table class="msgT"><tr><td valign="top"><div class="PathIndex PathIndexControl">4</div></td><td><div class="PathNav"><a href="#Path3" title="Previous event (3)">&#x2190;</a></div></td></td><td>Control jumps to 'case 0:'  at line 420</td><td><div class="PathNav"><a href="#EndPath" title="Next event (5)">&#x2192;</a></div></td></tr></table></div></td></tr>


    		
  

  

    
      503
    
    
      
<tr><td class="num" id="LN420">420</td><td class="line">    <span class='keyword'>case</span> 0:</td></tr>


    		
  

  

    
      504
    
    
      
<tr><td class="num" id="LN421">421</td><td class="line">        index = ang2pix_z_phi_ring(<span class="mrange">z</span>, phi);</td></tr>


    		
  

  

    
      505
    
    
      
<tr><td class="num"></td><td class="line"><div id="EndPath" class="msg msgEvent" style="margin-left:17ex"><table class="msgT"><tr><td valign="top"><div class="PathIndex PathIndexEvent">5</div></td><td><div class="PathNav"><a href="#Path4" title="Previous event (4)">&#x2190;</a></div></td></td><td>Function call argument is an uninitialized value</td></tr></table></div></td></tr>


    		
  

  

    
      506
    
    
      
<tr><td class="num" id="LN422">422</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      507
    
    
      
<tr><td class="num" id="LN423">423</td><td class="line">    <span class='keyword'>case</span> 1:</td></tr>


    		
  

  

    
      508
    
    
      
<tr><td class="num" id="LN424">424</td><td class="line">        index = ang2pix_z_phi_nest(z, phi);</td></tr>


    		
  

  

    
      509
    
    
      
<tr><td class="num" id="LN425">425</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      510
    
    
      
<tr><td class="num" id="LN426">426</td><td class="line">    <span class='keyword'>default</span>:</td></tr>


    		
  

  

    
      511
    
    
      
<tr><td class="num" id="LN427">427</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      512
    
    
      
<tr><td class="num" id="LN428">428</td><td class="line">    }</td></tr>


    		
  

  

    
      513
    
    
      
<tr><td class="num" id="LN429">429</td><td class="line"> </td></tr>


    		
  

  

    
      514
    
    
      
<tr><td class="num" id="LN430">430</td><td class="line">    <span class='comment'>// Return sky map pixel</span></td></tr>


    		
  

  

    
      515
    
    
      
<tr><td class="num" id="LN431">431</td><td class="line">    <span class='keyword'>return</span> (GSkyPixel(index));</td></tr>


    		
  

  

    
      516
    
    
      
<tr><td class="num" id="LN432">432</td><td class="line">}</td></tr>


    		
  

  

    
      517
    
    
      
<tr><td class="num" id="LN433">433</td><td class="line"> </td></tr>


    		
  

  

    
      518
    
    
      
<tr><td class="num" id="LN434">434</td><td class="line"> </td></tr>


    		
  

  

    
      519
    
    
      
<tr><td class="num" id="LN435">435</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      520
    
    
      
<tr><td class="num" id="LN436">436</td><td class="line"> <span class='comment'>* @brief Returns ordering parameter.</span></td></tr>


    		
  

  

    
      521
    
    
      
<tr><td class="num" id="LN437">437</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      522
    
    
      
<tr><td class="num" id="LN438">438</td><td class="line">std::string GHealpix::ordering(<span class='keyword'>void</span>) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      523
    
    
      
<tr><td class="num" id="LN439">439</td><td class="line">{</td></tr>


    		
  

  

    
      524
    
    
      
<tr><td class="num" id="LN440">440</td><td class="line">    <span class='comment'>// Set pixel ordering type</span></td></tr>


    		
  

  

    
      525
    
    
      
<tr><td class="num" id="LN441">441</td><td class="line">    std::string s_ordering;</td></tr>


    		
  

  

    
      526
    
    
      
<tr><td class="num" id="LN442">442</td><td class="line">    <span class='keyword'>switch</span> (m_ordering) {</td></tr>


    		
  

  

    
      527
    
    
      
<tr><td class="num" id="LN443">443</td><td class="line">    <span class='keyword'>case</span> 0:</td></tr>


    		
  

  

    
      528
    
    
      
<tr><td class="num" id="LN444">444</td><td class="line">        s_ordering = <span class='string_literal'>"RING"</span>;</td></tr>


    		
  

  

    
      529
    
    
      
<tr><td class="num" id="LN445">445</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      530
    
    
      
<tr><td class="num" id="LN446">446</td><td class="line">    <span class='keyword'>case</span> 1:</td></tr>


    		
  

  

    
      531
    
    
      
<tr><td class="num" id="LN447">447</td><td class="line">        s_ordering = <span class='string_literal'>"NESTED"</span>;</td></tr>


    		
  

  

    
      532
    
    
      
<tr><td class="num" id="LN448">448</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      533
    
    
      
<tr><td class="num" id="LN449">449</td><td class="line">    <span class='keyword'>default</span>:</td></tr>


    		
  

  

    
      534
    
    
      
<tr><td class="num" id="LN450">450</td><td class="line">        s_ordering = <span class='string_literal'>"UNKNOWN"</span>;</td></tr>


    		
  

  

    
      535
    
    
      
<tr><td class="num" id="LN451">451</td><td class="line">        <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      536
    
    
      
<tr><td class="num" id="LN452">452</td><td class="line">    }</td></tr>


    		
  

  

    
      537
    
    
      
<tr><td class="num" id="LN453">453</td><td class="line"> </td></tr>


    		
  

  

    
      538
    
    
      
<tr><td class="num" id="LN454">454</td><td class="line">    <span class='comment'>// Return ordering</span></td></tr>


    		
  

  

    
      539
    
    
      
<tr><td class="num" id="LN455">455</td><td class="line">    <span class='keyword'>return</span> s_ordering;</td></tr>


    		
  

  

    
      540
    
    
      
<tr><td class="num" id="LN456">456</td><td class="line">}</td></tr>


    		
  

  

    
      541
    
    
      
<tr><td class="num" id="LN457">457</td><td class="line"> </td></tr>


    		
  

  

    
      542
    
    
      
<tr><td class="num" id="LN458">458</td><td class="line"> </td></tr>


    		
  

  

    
      543
    
    
      
<tr><td class="num" id="LN459">459</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      544
    
    
      
<tr><td class="num" id="LN460">460</td><td class="line"> <span class='comment'>* @brief Set pixel ordering.</span></td></tr>


    		
  

  

    
      545
    
    
      
<tr><td class="num" id="LN461">461</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      546
    
    
      
<tr><td class="num" id="LN462">462</td><td class="line"> <span class='comment'>* @param[in] ordering Pixel ordering (RING or NEST/NESTED).</span></td></tr>


    		
  

  

    
      547
    
    
      
<tr><td class="num" id="LN463">463</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      548
    
    
      
<tr><td class="num" id="LN464">464</td><td class="line"> <span class='comment'>* @exception GException::wcs_hpx_bad_ordering</span></td></tr>


    		
  

  

    
      549
    
    
      
<tr><td class="num" id="LN465">465</td><td class="line"> <span class='comment'>*            Invalid ordering parameter.</span></td></tr>


    		
  

  

    
      550
    
    
      
<tr><td class="num" id="LN466">466</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      551
    
    
      
<tr><td class="num" id="LN467">467</td><td class="line"><span class='keyword'>void</span> GHealpix::ordering(<span class='keyword'>const</span> std::string&amp; ordering)</td></tr>


    		
  

  

    
      552
    
    
      
<tr><td class="num" id="LN468">468</td><td class="line">{</td></tr>


    		
  

  

    
      553
    
    
      
<tr><td class="num" id="LN469">469</td><td class="line">    <span class='comment'>// Convert argument to upper case</span></td></tr>


    		
  

  

    
      554
    
    
      
<tr><td class="num" id="LN470">470</td><td class="line">    std::string uordering = gammalib::toupper(ordering);</td></tr>


    		
  

  

    
      555
    
    
      
<tr><td class="num" id="LN471">471</td><td class="line"> </td></tr>


    		
  

  

    
      556
    
    
      
<tr><td class="num" id="LN472">472</td><td class="line">    <span class='comment'>// Set pixel ordering</span></td></tr>


    		
  

  

    
      557
    
    
      
<tr><td class="num" id="LN473">473</td><td class="line">    <span class='keyword'>if</span> (uordering == <span class='string_literal'>"RING"</span>) {</td></tr>


    		
  

  

    
      558
    
    
      
<tr><td class="num" id="LN474">474</td><td class="line">        m_ordering = 0;</td></tr>


    		
  

  

    
      559
    
    
      
<tr><td class="num" id="LN475">475</td><td class="line">    }</td></tr>


    		
  

  

    
      560
    
    
      
<tr><td class="num" id="LN476">476</td><td class="line">    <span class='keyword'>else</span> <span class='keyword'>if</span> (uordering == <span class='string_literal'>"NESTED"</span> || uordering == <span class='string_literal'>"NEST"</span>) {</td></tr>


    		
  

  

    
      561
    
    
      
<tr><td class="num" id="LN477">477</td><td class="line">        m_ordering = 1;</td></tr>


    		
  

  

    
      562
    
    
      
<tr><td class="num" id="LN478">478</td><td class="line">    }</td></tr>


    		
  

  

    
      563
    
    
      
<tr><td class="num" id="LN479">479</td><td class="line">    <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      564
    
    
      
<tr><td class="num" id="LN480">480</td><td class="line">        <span class='keyword'>throw</span> GException::wcs_hpx_bad_ordering(<span class='macro'>G_ORDERING_SET<span class='expansion'>"GHealpix::ordering(std::string&amp;)"</span></span>, ordering);</td></tr>


    		
  

  

    
      565
    
    
      
<tr><td class="num" id="LN481">481</td><td class="line">    }</td></tr>


    		
  

  

    
      566
    
    
      
<tr><td class="num" id="LN482">482</td><td class="line"> </td></tr>


    		
  

  

    
      567
    
    
      
<tr><td class="num" id="LN483">483</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      568
    
    
      
<tr><td class="num" id="LN484">484</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      569
    
    
      
<tr><td class="num" id="LN485">485</td><td class="line">}</td></tr>


    		
  

  

    
      570
    
    
      
<tr><td class="num" id="LN486">486</td><td class="line"> </td></tr>


    		
  

  

    
      571
    
    
      
<tr><td class="num" id="LN487">487</td><td class="line"> </td></tr>


    		
  

  

    
      572
    
    
      
<tr><td class="num" id="LN488">488</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      573
    
    
      
<tr><td class="num" id="LN489">489</td><td class="line"> <span class='comment'>* @brief Print WCS information</span></td></tr>


    		
  

  

    
      574
    
    
      
<tr><td class="num" id="LN490">490</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      575
    
    
      
<tr><td class="num" id="LN491">491</td><td class="line"> <span class='comment'>* @param[in] chatter Chattiness (defaults to NORMAL).</span></td></tr>


    		
  

  

    
      576
    
    
      
<tr><td class="num" id="LN492">492</td><td class="line"> <span class='comment'>* @return String containing WCS information.</span></td></tr>


    		
  

  

    
      577
    
    
      
<tr><td class="num" id="LN493">493</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      578
    
    
      
<tr><td class="num" id="LN494">494</td><td class="line">std::string GHealpix::print(<span class='keyword'>const</span> GChatter&amp; chatter) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      579
    
    
      
<tr><td class="num" id="LN495">495</td><td class="line">{</td></tr>


    		
  

  

    
      580
    
    
      
<tr><td class="num" id="LN496">496</td><td class="line">    <span class='comment'>// Initialise result string</span></td></tr>


    		
  

  

    
      581
    
    
      
<tr><td class="num" id="LN497">497</td><td class="line">    std::string result;</td></tr>


    		
  

  

    
      582
    
    
      
<tr><td class="num" id="LN498">498</td><td class="line"> </td></tr>


    		
  

  

    
      583
    
    
      
<tr><td class="num" id="LN499">499</td><td class="line">    <span class='comment'>// Continue only if chatter is not silent</span></td></tr>


    		
  

  

    
      584
    
    
      
<tr><td class="num" id="LN500">500</td><td class="line">    <span class='keyword'>if</span> (chatter != SILENT) {</td></tr>


    		
  

  

    
      585
    
    
      
<tr><td class="num" id="LN501">501</td><td class="line"> </td></tr>


    		
  

  

    
      586
    
    
      
<tr><td class="num" id="LN502">502</td><td class="line">        <span class='comment'>// Append header</span></td></tr>


    		
  

  

    
      587
    
    
      
<tr><td class="num" id="LN503">503</td><td class="line">        result.append(<span class='string_literal'>"=== GHealpix ==="</span>);</td></tr>


    		
  

  

    
      588
    
    
      
<tr><td class="num" id="LN504">504</td><td class="line"> </td></tr>


    		
  

  

    
      589
    
    
      
<tr><td class="num" id="LN505">505</td><td class="line">        <span class='comment'>// Append information</span></td></tr>


    		
  

  

    
      590
    
    
      
<tr><td class="num" id="LN506">506</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Coordinate system"</span>));</td></tr>


    		
  

  

    
      591
    
    
      
<tr><td class="num" id="LN507">507</td><td class="line">        result.append(coordsys());</td></tr>


    		
  

  

    
      592
    
    
      
<tr><td class="num" id="LN508">508</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Nside (# of divisions)"</span>));</td></tr>


    		
  

  

    
      593
    
    
      
<tr><td class="num" id="LN509">509</td><td class="line">        result.append(gammalib::str(m_nside));</td></tr>


    		
  

  

    
      594
    
    
      
<tr><td class="num" id="LN510">510</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Npface (pixels per face)"</span>));</td></tr>


    		
  

  

    
      595
    
    
      
<tr><td class="num" id="LN511">511</td><td class="line">        result.append(gammalib::str(m_npface));</td></tr>


    		
  

  

    
      596
    
    
      
<tr><td class="num" id="LN512">512</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Ncap (# of cap pixels)"</span>));</td></tr>


    		
  

  

    
      597
    
    
      
<tr><td class="num" id="LN513">513</td><td class="line">        result.append(gammalib::str(m_ncap));</td></tr>


    		
  

  

    
      598
    
    
      
<tr><td class="num" id="LN514">514</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Npix (# of pixels)"</span>));</td></tr>


    		
  

  

    
      599
    
    
      
<tr><td class="num" id="LN515">515</td><td class="line">        result.append(gammalib::str(m_num_pixels));</td></tr>


    		
  

  

    
      600
    
    
      
<tr><td class="num" id="LN516">516</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Order"</span>));</td></tr>


    		
  

  

    
      601
    
    
      
<tr><td class="num" id="LN517">517</td><td class="line">        result.append(gammalib::str(m_order));</td></tr>


    		
  

  

    
      602
    
    
      
<tr><td class="num" id="LN518">518</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Solid angle per pixel"</span>));</td></tr>


    		
  

  

    
      603
    
    
      
<tr><td class="num" id="LN519">519</td><td class="line">        result.append(gammalib::str(m_omega)+<span class='string_literal'>" sr"</span>);</td></tr>


    		
  

  

    
      604
    
    
      
<tr><td class="num" id="LN520">520</td><td class="line">        result.append(<span class='string_literal'>"\n"</span>+gammalib::parformat(<span class='string_literal'>"Ordering"</span>)+ordering());</td></tr>


    		
  

  

    
      605
    
    
      
<tr><td class="num" id="LN521">521</td><td class="line"> </td></tr>


    		
  

  

    
      606
    
    
      
<tr><td class="num" id="LN522">522</td><td class="line">    } <span class='comment'>// endif: chatter was not silent</span></td></tr>


    		
  

  

    
      607
    
    
      
<tr><td class="num" id="LN523">523</td><td class="line"> </td></tr>


    		
  

  

    
      608
    
    
      
<tr><td class="num" id="LN524">524</td><td class="line">    <span class='comment'>// Return result</span></td></tr>


    		
  

  

    
      609
    
    
      
<tr><td class="num" id="LN525">525</td><td class="line">    <span class='keyword'>return</span> result;</td></tr>


    		
  

  

    
      610
    
    
      
<tr><td class="num" id="LN526">526</td><td class="line">}</td></tr>


    		
  

  

    
      611
    
    
      
<tr><td class="num" id="LN527">527</td><td class="line"> </td></tr>


    		
  

  

    
      612
    
    
      
<tr><td class="num" id="LN528">528</td><td class="line"> </td></tr>


    		
  

  

    
      613
    
    
      
<tr><td class="num" id="LN529">529</td><td class="line"><span class='comment'>/*==========================================================================</span></td></tr>


    		
  

  

    
      614
    
    
      
<tr><td class="num" id="LN530">530</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      615
    
    
      
<tr><td class="num" id="LN531">531</td><td class="line"> <span class='comment'>=                             Private methods                             =</span></td></tr>


    		
  

  

    
      616
    
    
      
<tr><td class="num" id="LN532">532</td><td class="line"> <span class='comment'>=                                                                         =</span></td></tr>


    		
  

  

    
      617
    
    
      
<tr><td class="num" id="LN533">533</td><td class="line"> <span class='comment'>==========================================================================*/</span></td></tr>


    		
  

  

    
      618
    
    
      
<tr><td class="num" id="LN534">534</td><td class="line"> </td></tr>


    		
  

  

    
      619
    
    
      
<tr><td class="num" id="LN535">535</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      620
    
    
      
<tr><td class="num" id="LN536">536</td><td class="line"> <span class='comment'>* @brief Initialise class members</span></td></tr>


    		
  

  

    
      621
    
    
      
<tr><td class="num" id="LN537">537</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      622
    
    
      
<tr><td class="num" id="LN538">538</td><td class="line"><span class='keyword'>void</span> GHealpix::init_members(<span class='keyword'>void</span>)</td></tr>


    		
  

  

    
      623
    
    
      
<tr><td class="num" id="LN539">539</td><td class="line">{</td></tr>


    		
  

  

    
      624
    
    
      
<tr><td class="num" id="LN540">540</td><td class="line">    <span class='comment'>// Initialise members</span></td></tr>


    		
  

  

    
      625
    
    
      
<tr><td class="num" id="LN541">541</td><td class="line">    <span class='comment'>//m_type        = "HPX";</span></td></tr>


    		
  

  

    
      626
    
    
      
<tr><td class="num" id="LN542">542</td><td class="line">    m_nside       = 0;</td></tr>


    		
  

  

    
      627
    
    
      
<tr><td class="num" id="LN543">543</td><td class="line">    m_npface      = 0;</td></tr>


    		
  

  

    
      628
    
    
      
<tr><td class="num" id="LN544">544</td><td class="line">    m_ncap        = 0;</td></tr>


    		
  

  

    
      629
    
    
      
<tr><td class="num" id="LN545">545</td><td class="line">    m_order       = 0;</td></tr>


    		
  

  

    
      630
    
    
      
<tr><td class="num" id="LN546">546</td><td class="line">    m_ordering    = 0;</td></tr>


    		
  

  

    
      631
    
    
      
<tr><td class="num" id="LN547">547</td><td class="line">    m_num_pixels  = 0;</td></tr>


    		
  

  

    
      632
    
    
      
<tr><td class="num" id="LN548">548</td><td class="line">    m_fact1       = 0.0;</td></tr>


    		
  

  

    
      633
    
    
      
<tr><td class="num" id="LN549">549</td><td class="line">    m_fact2       = 0.0;</td></tr>


    		
  

  

    
      634
    
    
      
<tr><td class="num" id="LN550">550</td><td class="line">    m_omega       = 0.0;</td></tr>


    		
  

  

    
      635
    
    
      
<tr><td class="num" id="LN551">551</td><td class="line"> </td></tr>


    		
  

  

    
      636
    
    
      
<tr><td class="num" id="LN552">552</td><td class="line">    <span class='comment'>// Construct conversion arrays</span></td></tr>


    		
  

  

    
      637
    
    
      
<tr><td class="num" id="LN553">553</td><td class="line">    <span class='keyword'>for</span> (<span class='keyword'>int</span> m = 0; m &lt; 0x100; ++m) {</td></tr>


    		
  

  

    
      638
    
    
      
<tr><td class="num" id="LN554">554</td><td class="line">    ctab[m] =</td></tr>


    		
  

  

    
      639
    
    
      
<tr><td class="num" id="LN555">555</td><td class="line">         (m&amp;0x1 )       | ((m&amp;0x2 ) &lt;&lt; 7) | ((m&amp;0x4 ) &gt;&gt; 1) | ((m&amp;0x8 ) &lt;&lt; 6)</td></tr>


    		
  

  

    
      640
    
    
      
<tr><td class="num" id="LN556">556</td><td class="line">      | ((m&amp;0x10) &gt;&gt; 2) | ((m&amp;0x20) &lt;&lt; 5) | ((m&amp;0x40) &gt;&gt; 3) | ((m&amp;0x80) &lt;&lt; 4);</td></tr>


    		
  

  

    
      641
    
    
      
<tr><td class="num" id="LN557">557</td><td class="line">    utab[m] =</td></tr>


    		
  

  

    
      642
    
    
      
<tr><td class="num" id="LN558">558</td><td class="line">         (m&amp;0x1 )       | ((m&amp;0x2 ) &lt;&lt; 1) | ((m&amp;0x4 ) &lt;&lt; 2) | ((m&amp;0x8 ) &lt;&lt; 3)</td></tr>


    		
  

  

    
      643
    
    
      
<tr><td class="num" id="LN559">559</td><td class="line">      | ((m&amp;0x10) &lt;&lt; 4) | ((m&amp;0x20) &lt;&lt; 5) | ((m&amp;0x40) &lt;&lt; 6) | ((m&amp;0x80) &lt;&lt; 7);</td></tr>


    		
  

  

    
      644
    
    
      
<tr><td class="num" id="LN560">560</td><td class="line">    }</td></tr>


    		
  

  

    
      645
    
    
      
<tr><td class="num" id="LN561">561</td><td class="line"> </td></tr>


    		
  

  

    
      646
    
    
      
<tr><td class="num" id="LN562">562</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      647
    
    
      
<tr><td class="num" id="LN563">563</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      648
    
    
      
<tr><td class="num" id="LN564">564</td><td class="line">}</td></tr>


    		
  

  

    
      649
    
    
      
<tr><td class="num" id="LN565">565</td><td class="line"> </td></tr>


    		
  

  

    
      650
    
    
      
<tr><td class="num" id="LN566">566</td><td class="line"> </td></tr>


    		
  

  

    
      651
    
    
      
<tr><td class="num" id="LN567">567</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      652
    
    
      
<tr><td class="num" id="LN568">568</td><td class="line"> <span class='comment'>* @brief Copy class members</span></td></tr>


    		
  

  

    
      653
    
    
      
<tr><td class="num" id="LN569">569</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      654
    
    
      
<tr><td class="num" id="LN570">570</td><td class="line"> <span class='comment'>* @param[in] proj Healpix projection.</span></td></tr>


    		
  

  

    
      655
    
    
      
<tr><td class="num" id="LN571">571</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      656
    
    
      
<tr><td class="num" id="LN572">572</td><td class="line"><span class='keyword'>void</span> GHealpix::copy_members(<span class='keyword'>const</span> GHealpix&amp; proj)</td></tr>


    		
  

  

    
      657
    
    
      
<tr><td class="num" id="LN573">573</td><td class="line">{</td></tr>


    		
  

  

    
      658
    
    
      
<tr><td class="num" id="LN574">574</td><td class="line">    <span class='comment'>// Copy attributes</span></td></tr>


    		
  

  

    
      659
    
    
      
<tr><td class="num" id="LN575">575</td><td class="line">    m_nside      = proj.m_nside;</td></tr>


    		
  

  

    
      660
    
    
      
<tr><td class="num" id="LN576">576</td><td class="line">    m_npface     = proj.m_npface;</td></tr>


    		
  

  

    
      661
    
    
      
<tr><td class="num" id="LN577">577</td><td class="line">    m_ncap       = proj.m_ncap;</td></tr>


    		
  

  

    
      662
    
    
      
<tr><td class="num" id="LN578">578</td><td class="line">    m_order      = proj.m_order;</td></tr>


    		
  

  

    
      663
    
    
      
<tr><td class="num" id="LN579">579</td><td class="line">    m_ordering   = proj.m_ordering;</td></tr>


    		
  

  

    
      664
    
    
      
<tr><td class="num" id="LN580">580</td><td class="line">    m_num_pixels = proj.m_num_pixels;</td></tr>


    		
  

  

    
      665
    
    
      
<tr><td class="num" id="LN581">581</td><td class="line">    m_fact1      = proj.m_fact1;</td></tr>


    		
  

  

    
      666
    
    
      
<tr><td class="num" id="LN582">582</td><td class="line">    m_fact2      = proj.m_fact2;</td></tr>


    		
  

  

    
      667
    
    
      
<tr><td class="num" id="LN583">583</td><td class="line">    m_omega      = proj.m_omega;</td></tr>


    		
  

  

    
      668
    
    
      
<tr><td class="num" id="LN584">584</td><td class="line"> </td></tr>


    		
  

  

    
      669
    
    
      
<tr><td class="num" id="LN585">585</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      670
    
    
      
<tr><td class="num" id="LN586">586</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      671
    
    
      
<tr><td class="num" id="LN587">587</td><td class="line">}</td></tr>


    		
  

  

    
      672
    
    
      
<tr><td class="num" id="LN588">588</td><td class="line"> </td></tr>


    		
  

  

    
      673
    
    
      
<tr><td class="num" id="LN589">589</td><td class="line"> </td></tr>


    		
  

  

    
      674
    
    
      
<tr><td class="num" id="LN590">590</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      675
    
    
      
<tr><td class="num" id="LN591">591</td><td class="line"> <span class='comment'>* @brief Delete class members</span></td></tr>


    		
  

  

    
      676
    
    
      
<tr><td class="num" id="LN592">592</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      677
    
    
      
<tr><td class="num" id="LN593">593</td><td class="line"><span class='keyword'>void</span> GHealpix::free_members(<span class='keyword'>void</span>)</td></tr>


    		
  

  

    
      678
    
    
      
<tr><td class="num" id="LN594">594</td><td class="line">{</td></tr>


    		
  

  

    
      679
    
    
      
<tr><td class="num" id="LN595">595</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      680
    
    
      
<tr><td class="num" id="LN596">596</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      681
    
    
      
<tr><td class="num" id="LN597">597</td><td class="line">}</td></tr>


    		
  

  

    
      682
    
    
      
<tr><td class="num" id="LN598">598</td><td class="line"> </td></tr>


    		
  

  

    
      683
    
    
      
<tr><td class="num" id="LN599">599</td><td class="line"> </td></tr>


    		
  

  

    
      684
    
    
      
<tr><td class="num" id="LN600">600</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      685
    
    
      
<tr><td class="num" id="LN601">601</td><td class="line"> <span class='comment'>* @brief Returns true if argument is identical</span></td></tr>


    		
  

  

    
      686
    
    
      
<tr><td class="num" id="LN602">602</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      687
    
    
      
<tr><td class="num" id="LN603">603</td><td class="line"> <span class='comment'>* @param[in] proj Sky projection.</span></td></tr>


    		
  

  

    
      688
    
    
      
<tr><td class="num" id="LN604">604</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      689
    
    
      
<tr><td class="num" id="LN605">605</td><td class="line"> <span class='comment'>* This method is a helper for the sky projection friends.</span></td></tr>


    		
  

  

    
      690
    
    
      
<tr><td class="num" id="LN606">606</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      691
    
    
      
<tr><td class="num" id="LN607">607</td><td class="line"><span class='keyword'>bool</span> GHealpix::compare(<span class='keyword'>const</span> GSkyProjection&amp; proj) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      692
    
    
      
<tr><td class="num" id="LN608">608</td><td class="line">{</td></tr>


    		
  

  

    
      693
    
    
      
<tr><td class="num" id="LN609">609</td><td class="line">    <span class='comment'>// Initialise result</span></td></tr>


    		
  

  

    
      694
    
    
      
<tr><td class="num" id="LN610">610</td><td class="line">    <span class='keyword'>bool</span> result = <span class='keyword'>false</span>;</td></tr>


    		
  

  

    
      695
    
    
      
<tr><td class="num" id="LN611">611</td><td class="line">    </td></tr>


    		
  

  

    
      696
    
    
      
<tr><td class="num" id="LN612">612</td><td class="line">    <span class='comment'>// Continue only we compare to a GHealpix object</span></td></tr>


    		
  

  

    
      697
    
    
      
<tr><td class="num" id="LN613">613</td><td class="line">    <span class='keyword'>const</span> GHealpix* ptr = <span class='keyword'>dynamic_cast</span>&lt;<span class='keyword'>const</span> GHealpix*&gt;(&amp;proj);</td></tr>


    		
  

  

    
      698
    
    
      
<tr><td class="num" id="LN614">614</td><td class="line">    <span class='keyword'>if</span> (ptr != <span class='macro'>NULL<span class='expansion'>__null</span></span>) {</td></tr>


    		
  

  

    
      699
    
    
      
<tr><td class="num" id="LN615">615</td><td class="line">        result = ((m_coordsys   == ptr-&gt;m_coordsys) &amp;&amp;</td></tr>


    		
  

  

    
      700
    
    
      
<tr><td class="num" id="LN616">616</td><td class="line">                  (m_nside      == ptr-&gt;m_nside)    &amp;&amp;</td></tr>


    		
  

  

    
      701
    
    
      
<tr><td class="num" id="LN617">617</td><td class="line">                  (m_npface     == ptr-&gt;m_npface)   &amp;&amp;</td></tr>


    		
  

  

    
      702
    
    
      
<tr><td class="num" id="LN618">618</td><td class="line">                  (m_ncap       == ptr-&gt;m_ncap)     &amp;&amp;</td></tr>


    		
  

  

    
      703
    
    
      
<tr><td class="num" id="LN619">619</td><td class="line">                  (m_order      == ptr-&gt;m_order)    &amp;&amp;</td></tr>


    		
  

  

    
      704
    
    
      
<tr><td class="num" id="LN620">620</td><td class="line">                  (m_ordering   == ptr-&gt;m_ordering) &amp;&amp;</td></tr>


    		
  

  

    
      705
    
    
      
<tr><td class="num" id="LN621">621</td><td class="line">                  (m_num_pixels == ptr-&gt;m_num_pixels));</td></tr>


    		
  

  

    
      706
    
    
      
<tr><td class="num" id="LN622">622</td><td class="line">    }</td></tr>


    		
  

  

    
      707
    
    
      
<tr><td class="num" id="LN623">623</td><td class="line"> </td></tr>


    		
  

  

    
      708
    
    
      
<tr><td class="num" id="LN624">624</td><td class="line">    <span class='comment'>// Return result</span></td></tr>


    		
  

  

    
      709
    
    
      
<tr><td class="num" id="LN625">625</td><td class="line">    <span class='keyword'>return</span> result;</td></tr>


    		
  

  

    
      710
    
    
      
<tr><td class="num" id="LN626">626</td><td class="line">}</td></tr>


    		
  

  

    
      711
    
    
      
<tr><td class="num" id="LN627">627</td><td class="line"> </td></tr>


    		
  

  

    
      712
    
    
      
<tr><td class="num" id="LN628">628</td><td class="line"> </td></tr>


    		
  

  

    
      713
    
    
      
<tr><td class="num" id="LN629">629</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      714
    
    
      
<tr><td class="num" id="LN630">630</td><td class="line"> <span class='comment'>* @brief Convert nside to order</span></td></tr>


    		
  

  

    
      715
    
    
      
<tr><td class="num" id="LN631">631</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      716
    
    
      
<tr><td class="num" id="LN632">632</td><td class="line"> <span class='comment'>* @param[in] nside Number of sides.</span></td></tr>


    		
  

  

    
      717
    
    
      
<tr><td class="num" id="LN633">633</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      718
    
    
      
<tr><td class="num" id="LN634">634</td><td class="line"><span class='keyword'>int</span> GHealpix::nside2order(<span class='keyword'>int</span> nside)</td></tr>


    		
  

  

    
      719
    
    
      
<tr><td class="num" id="LN635">635</td><td class="line">{</td></tr>


    		
  

  

    
      720
    
    
      
<tr><td class="num" id="LN636">636</td><td class="line">    <span class='comment'>// Initialise order</span></td></tr>


    		
  

  

    
      721
    
    
      
<tr><td class="num" id="LN637">637</td><td class="line">    <span class='keyword'>int</span> order = -1;</td></tr>


    		
  

  

    
      722
    
    
      
<tr><td class="num" id="LN638">638</td><td class="line"> </td></tr>


    		
  

  

    
      723
    
    
      
<tr><td class="num" id="LN639">639</td><td class="line">    <span class='comment'>// Determine order</span></td></tr>


    		
  

  

    
      724
    
    
      
<tr><td class="num" id="LN640">640</td><td class="line">    <span class='keyword'>for</span> (<span class='keyword'>int</span> m = 0; m &lt;= order_max; ++m) {</td></tr>


    		
  

  

    
      725
    
    
      
<tr><td class="num" id="LN641">641</td><td class="line">        <span class='keyword'>int</span> nstest = 1 &lt;&lt; m;</td></tr>


    		
  

  

    
      726
    
    
      
<tr><td class="num" id="LN642">642</td><td class="line">        <span class='keyword'>if</span> (nside == nstest) {</td></tr>


    		
  

  

    
      727
    
    
      
<tr><td class="num" id="LN643">643</td><td class="line">            order = m;</td></tr>


    		
  

  

    
      728
    
    
      
<tr><td class="num" id="LN644">644</td><td class="line">            <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      729
    
    
      
<tr><td class="num" id="LN645">645</td><td class="line">        }</td></tr>


    		
  

  

    
      730
    
    
      
<tr><td class="num" id="LN646">646</td><td class="line">        <span class='keyword'>if</span> (nside &lt; nstest)</td></tr>


    		
  

  

    
      731
    
    
      
<tr><td class="num" id="LN647">647</td><td class="line">            <span class='keyword'>break</span>;</td></tr>


    		
  

  

    
      732
    
    
      
<tr><td class="num" id="LN648">648</td><td class="line">    }</td></tr>


    		
  

  

    
      733
    
    
      
<tr><td class="num" id="LN649">649</td><td class="line"> </td></tr>


    		
  

  

    
      734
    
    
      
<tr><td class="num" id="LN650">650</td><td class="line">    <span class='comment'>// Return order</span></td></tr>


    		
  

  

    
      735
    
    
      
<tr><td class="num" id="LN651">651</td><td class="line">    <span class='keyword'>return</span> order;</td></tr>


    		
  

  

    
      736
    
    
      
<tr><td class="num" id="LN652">652</td><td class="line">}</td></tr>


    		
  

  

    
      737
    
    
      
<tr><td class="num" id="LN653">653</td><td class="line"> </td></tr>


    		
  

  

    
      738
    
    
      
<tr><td class="num" id="LN654">654</td><td class="line"> </td></tr>


    		
  

  

    
      739
    
    
      
<tr><td class="num" id="LN655">655</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      740
    
    
      
<tr><td class="num" id="LN656">656</td><td class="line"> <span class='comment'>* @brief Convert pixel index to (x,y) coordinate</span></td></tr>


    		
  

  

    
      741
    
    
      
<tr><td class="num" id="LN657">657</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      742
    
    
      
<tr><td class="num" id="LN658">658</td><td class="line"> <span class='comment'>* @param[in] ipix Pixel index for which (x,y) are to be computed.</span></td></tr>


    		
  

  

    
      743
    
    
      
<tr><td class="num" id="LN659">659</td><td class="line"> <span class='comment'>* @param[out] x Pointer to x coordinate.</span></td></tr>


    		
  

  

    
      744
    
    
      
<tr><td class="num" id="LN660">660</td><td class="line"> <span class='comment'>* @param[out] y Pointer to y coordinate.</span></td></tr>


    		
  

  

    
      745
    
    
      
<tr><td class="num" id="LN661">661</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      746
    
    
      
<tr><td class="num" id="LN662">662</td><td class="line"><span class='keyword'>void</span> GHealpix::pix2xy(<span class='keyword'>const</span> <span class='keyword'>int</span>&amp; ipix, <span class='keyword'>int</span>* x, <span class='keyword'>int</span>* y) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      747
    
    
      
<tr><td class="num" id="LN663">663</td><td class="line">{</td></tr>


    		
  

  

    
      748
    
    
      
<tr><td class="num" id="LN664">664</td><td class="line">    <span class='comment'>// Set x coordinate</span></td></tr>


    		
  

  

    
      749
    
    
      
<tr><td class="num" id="LN665">665</td><td class="line">    <span class='keyword'>int</span> raw = (ipix &amp; 0x5555) | ((ipix &amp; 0x55550000) &gt;&gt; 15);</td></tr>


    		
  

  

    
      750
    
    
      
<tr><td class="num" id="LN666">666</td><td class="line">    *x      = ctab[raw &amp; 0xff] | (ctab[raw &gt;&gt; 8] &lt;&lt; 4);</td></tr>


    		
  

  

    
      751
    
    
      
<tr><td class="num" id="LN667">667</td><td class="line"> </td></tr>


    		
  

  

    
      752
    
    
      
<tr><td class="num" id="LN668">668</td><td class="line">    <span class='comment'>// Set y coordinate</span></td></tr>


    		
  

  

    
      753
    
    
      
<tr><td class="num" id="LN669">669</td><td class="line">    raw = ((ipix &amp; 0xaaaa) &gt;&gt; 1) | ((ipix &amp; 0xaaaa0000) &gt;&gt; 16);</td></tr>


    		
  

  

    
      754
    
    
      
<tr><td class="num" id="LN670">670</td><td class="line">    *y  = ctab[raw &amp; 0xff] | (ctab[raw &gt;&gt; 8] &lt;&lt; 4);</td></tr>


    		
  

  

    
      755
    
    
      
<tr><td class="num" id="LN671">671</td><td class="line"> </td></tr>


    		
  

  

    
      756
    
    
      
<tr><td class="num" id="LN672">672</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      757
    
    
      
<tr><td class="num" id="LN673">673</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      758
    
    
      
<tr><td class="num" id="LN674">674</td><td class="line">}</td></tr>


    		
  

  

    
      759
    
    
      
<tr><td class="num" id="LN675">675</td><td class="line"> </td></tr>


    		
  

  

    
      760
    
    
      
<tr><td class="num" id="LN676">676</td><td class="line"> </td></tr>


    		
  

  

    
      761
    
    
      
<tr><td class="num" id="LN677">677</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      762
    
    
      
<tr><td class="num" id="LN678">678</td><td class="line"> <span class='comment'>* @brief Convert (x,y) coordinate to pixel</span></td></tr>


    		
  

  

    
      763
    
    
      
<tr><td class="num" id="LN679">679</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      764
    
    
      
<tr><td class="num" id="LN680">680</td><td class="line"> <span class='comment'>* @param[in] x x coordinate.</span></td></tr>


    		
  

  

    
      765
    
    
      
<tr><td class="num" id="LN681">681</td><td class="line"> <span class='comment'>* @param[in] y y coordinate.</span></td></tr>


    		
  

  

    
      766
    
    
      
<tr><td class="num" id="LN682">682</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      767
    
    
      
<tr><td class="num" id="LN683">683</td><td class="line"><span class='keyword'>int</span> GHealpix::xy2pix(<span class='keyword'>int</span> x, <span class='keyword'>int</span> y) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      768
    
    
      
<tr><td class="num" id="LN684">684</td><td class="line">{</td></tr>


    		
  

  

    
      769
    
    
      
<tr><td class="num" id="LN685">685</td><td class="line">    <span class='comment'>// Return pixel</span></td></tr>


    		
  

  

    
      770
    
    
      
<tr><td class="num" id="LN686">686</td><td class="line">    <span class='keyword'>return</span> utab[x&amp;0xff] | (utab[x&gt;&gt;8]&lt;&lt;16) | (utab[y&amp;0xff]&lt;&lt;1) | (utab[y&gt;&gt;8]&lt;&lt;17);</td></tr>


    		
  

  

    
      771
    
    
      
<tr><td class="num" id="LN687">687</td><td class="line">}</td></tr>


    		
  

  

    
      772
    
    
      
<tr><td class="num" id="LN688">688</td><td class="line"> </td></tr>


    		
  

  

    
      773
    
    
      
<tr><td class="num" id="LN689">689</td><td class="line"> </td></tr>


    		
  

  

    
      774
    
    
      
<tr><td class="num" id="LN690">690</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      775
    
    
      
<tr><td class="num" id="LN691">691</td><td class="line"> <span class='comment'>* @brief Convert pixel index to (theta,phi) angles for ring ordering</span></td></tr>


    		
  

  

    
      776
    
    
      
<tr><td class="num" id="LN692">692</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      777
    
    
      
<tr><td class="num" id="LN693">693</td><td class="line"> <span class='comment'>* @param[in] ipix Pixel index for which (theta,phi) are to be computed.</span></td></tr>


    		
  

  

    
      778
    
    
      
<tr><td class="num" id="LN694">694</td><td class="line"> <span class='comment'>* @param[out] theta Pointer to result zenith angle in radians.</span></td></tr>


    		
  

  

    
      779
    
    
      
<tr><td class="num" id="LN695">695</td><td class="line"> <span class='comment'>* @param[out] phi Pointer to result azimuth angle in radians.</span></td></tr>


    		
  

  

    
      780
    
    
      
<tr><td class="num" id="LN696">696</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      781
    
    
      
<tr><td class="num" id="LN697">697</td><td class="line"> <span class='comment'>* @exception GException::out_of_range</span></td></tr>


    		
  

  

    
      782
    
    
      
<tr><td class="num" id="LN698">698</td><td class="line"> <span class='comment'>*            Pixel index is out of range.</span></td></tr>


    		
  

  

    
      783
    
    
      
<tr><td class="num" id="LN699">699</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      784
    
    
      
<tr><td class="num" id="LN700">700</td><td class="line"><span class='keyword'>void</span> GHealpix::pix2ang_ring(<span class='keyword'>int</span> ipix, <span class='keyword'>double</span>* theta, <span class='keyword'>double</span>* phi) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      785
    
    
      
<tr><td class="num" id="LN701">701</td><td class="line">{</td></tr>


    		
  

  

    
      786
    
    
      
<tr><td class="num" id="LN702">702</td><td class="line">    <span class='comment'>// Check if ipix is in range</span></td></tr>


    		
  

  

    
      787
    
    
      
<tr><td class="num" id="LN703">703</td><td class="line">    <span class='keyword'>if</span> (ipix &lt; 0 || ipix &gt;= m_num_pixels)</td></tr>


    		
  

  

    
      788
    
    
      
<tr><td class="num" id="LN704">704</td><td class="line">        <span class='keyword'>throw</span>  GException::out_of_range(<span class='macro'>G_PIX2ANG_RING<span class='expansion'>"GHealpix::pix2ang_ring(int, double*, double*)"</span></span>, ipix, 0, m_num_pixels-1);</td></tr>


    		
  

  

    
      789
    
    
      
<tr><td class="num" id="LN705">705</td><td class="line"> </td></tr>


    		
  

  

    
      790
    
    
      
<tr><td class="num" id="LN706">706</td><td class="line">    <span class='comment'>// Handle North Polar cap</span></td></tr>


    		
  

  

    
      791
    
    
      
<tr><td class="num" id="LN707">707</td><td class="line">    <span class='keyword'>if</span> (ipix &lt; m_ncap) {</td></tr>


    		
  

  

    
      792
    
    
      
<tr><td class="num" id="LN708">708</td><td class="line">        <span class='keyword'>int</span> iring = <span class='keyword'>int</span>(0.5*(1+isqrt(1+2*ipix))); <span class='comment'>// counted from North pole</span></td></tr>


    		
  

  

    
      793
    
    
      
<tr><td class="num" id="LN709">709</td><td class="line">        <span class='keyword'>int</span> iphi  = (ipix+1) - 2*iring*(iring-1);</td></tr>


    		
  

  

    
      794
    
    
      
<tr><td class="num" id="LN710">710</td><td class="line">        *theta    = acos(1.0 - (iring*iring) * m_fact2);</td></tr>


    		
  

  

    
      795
    
    
      
<tr><td class="num" id="LN711">711</td><td class="line">        *phi      = (iphi - 0.5) * gammalib::pi/(2.0*iring);</td></tr>


    		
  

  

    
      796
    
    
      
<tr><td class="num" id="LN712">712</td><td class="line">    }</td></tr>


    		
  

  

    
      797
    
    
      
<tr><td class="num" id="LN713">713</td><td class="line"> </td></tr>


    		
  

  

    
      798
    
    
      
<tr><td class="num" id="LN714">714</td><td class="line">    <span class='comment'>// Handle Equatorial region</span></td></tr>


    		
  

  

    
      799
    
    
      
<tr><td class="num" id="LN715">715</td><td class="line">    <span class='keyword'>else</span> <span class='keyword'>if</span> (ipix &lt; (m_num_pixels - m_ncap)) {</td></tr>


    		
  

  

    
      800
    
    
      
<tr><td class="num" id="LN716">716</td><td class="line">        <span class='keyword'>int</span>    ip    = ipix - m_ncap;</td></tr>


    		
  

  

    
      801
    
    
      
<tr><td class="num" id="LN717">717</td><td class="line">        <span class='keyword'>int</span>    iring = ip/(4*m_nside) + m_nside;   <span class='comment'>// counted from North pole</span></td></tr>


    		
  

  

    
      802
    
    
      
<tr><td class="num" id="LN718">718</td><td class="line">        <span class='keyword'>int</span>    iphi  = ip%(4*m_nside) + 1;</td></tr>


    		
  

  

    
      803
    
    
      
<tr><td class="num" id="LN719">719</td><td class="line">        <span class='keyword'>double</span> fodd  = ((iring+m_nside)&amp;1) ? 1 : 0.5;</td></tr>


    		
  

  

    
      804
    
    
      
<tr><td class="num" id="LN720">720</td><td class="line">        <span class='keyword'>int</span>    nl2   = 2*m_nside;</td></tr>


    		
  

  

    
      805
    
    
      
<tr><td class="num" id="LN721">721</td><td class="line">        *theta       = std::acos((nl2 - iring) * m_fact1);</td></tr>


    		
  

  

    
      806
    
    
      
<tr><td class="num" id="LN722">722</td><td class="line">        *phi         = (iphi - fodd) * gammalib::pi/nl2;</td></tr>


    		
  

  

    
      807
    
    
      
<tr><td class="num" id="LN723">723</td><td class="line">    }</td></tr>


    		
  

  

    
      808
    
    
      
<tr><td class="num" id="LN724">724</td><td class="line"> </td></tr>


    		
  

  

    
      809
    
    
      
<tr><td class="num" id="LN725">725</td><td class="line">    <span class='comment'>// Handle South Polar cap</span></td></tr>


    		
  

  

    
      810
    
    
      
<tr><td class="num" id="LN726">726</td><td class="line">    <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      811
    
    
      
<tr><td class="num" id="LN727">727</td><td class="line">        <span class='keyword'>int</span> ip    = m_num_pixels - ipix;</td></tr>


    		
  

  

    
      812
    
    
      
<tr><td class="num" id="LN728">728</td><td class="line">        <span class='keyword'>int</span> iring = <span class='keyword'>int</span>(0.5*(1+isqrt(2*ip-1)));    <span class='comment'>// Counted from South pole</span></td></tr>


    		
  

  

    
      813
    
    
      
<tr><td class="num" id="LN729">729</td><td class="line">        <span class='keyword'>int</span> iphi  = 4*iring + 1 - (ip - 2*iring*(iring-1));</td></tr>


    		
  

  

    
      814
    
    
      
<tr><td class="num" id="LN730">730</td><td class="line">        *theta    = std::acos(-1.0 + (iring*iring) * m_fact2);</td></tr>


    		
  

  

    
      815
    
    
      
<tr><td class="num" id="LN731">731</td><td class="line">        *phi      = (iphi - 0.5) * gammalib::pi/(2.*iring);</td></tr>


    		
  

  

    
      816
    
    
      
<tr><td class="num" id="LN732">732</td><td class="line">    }</td></tr>


    		
  

  

    
      817
    
    
      
<tr><td class="num" id="LN733">733</td><td class="line"> </td></tr>


    		
  

  

    
      818
    
    
      
<tr><td class="num" id="LN734">734</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      819
    
    
      
<tr><td class="num" id="LN735">735</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      820
    
    
      
<tr><td class="num" id="LN736">736</td><td class="line">}</td></tr>


    		
  

  

    
      821
    
    
      
<tr><td class="num" id="LN737">737</td><td class="line"> </td></tr>


    		
  

  

    
      822
    
    
      
<tr><td class="num" id="LN738">738</td><td class="line"> </td></tr>


    		
  

  

    
      823
    
    
      
<tr><td class="num" id="LN739">739</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      824
    
    
      
<tr><td class="num" id="LN740">740</td><td class="line"> <span class='comment'>* @brief Convert pixel index to (theta,phi) angles for nested ordering</span></td></tr>


    		
  

  

    
      825
    
    
      
<tr><td class="num" id="LN741">741</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      826
    
    
      
<tr><td class="num" id="LN742">742</td><td class="line"> <span class='comment'>* @param[in] ipix Pixel index for which (theta,phi) are to be computed.</span></td></tr>


    		
  

  

    
      827
    
    
      
<tr><td class="num" id="LN743">743</td><td class="line"> <span class='comment'>* @param[out] theta Pointer to result zenith angle in radians.</span></td></tr>


    		
  

  

    
      828
    
    
      
<tr><td class="num" id="LN744">744</td><td class="line"> <span class='comment'>* @param[out] phi Pointer to result azimuth angle in radians.</span></td></tr>


    		
  

  

    
      829
    
    
      
<tr><td class="num" id="LN745">745</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      830
    
    
      
<tr><td class="num" id="LN746">746</td><td class="line"> <span class='comment'>* @exception GException::out_of_range</span></td></tr>


    		
  

  

    
      831
    
    
      
<tr><td class="num" id="LN747">747</td><td class="line"> <span class='comment'>*            Pixel index is out of range.</span></td></tr>


    		
  

  

    
      832
    
    
      
<tr><td class="num" id="LN748">748</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      833
    
    
      
<tr><td class="num" id="LN749">749</td><td class="line"><span class='keyword'>void</span> GHealpix::pix2ang_nest(<span class='keyword'>int</span> ipix, <span class='keyword'>double</span>* theta, <span class='keyword'>double</span>* phi) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      834
    
    
      
<tr><td class="num" id="LN750">750</td><td class="line">{</td></tr>


    		
  

  

    
      835
    
    
      
<tr><td class="num" id="LN751">751</td><td class="line">    <span class='comment'>// Check if ipix is in range</span></td></tr>


    		
  

  

    
      836
    
    
      
<tr><td class="num" id="LN752">752</td><td class="line">    <span class='keyword'>if</span> (ipix &lt; 0 || ipix &gt;= m_num_pixels)</td></tr>


    		
  

  

    
      837
    
    
      
<tr><td class="num" id="LN753">753</td><td class="line">        <span class='keyword'>throw</span> GException::out_of_range(<span class='macro'>G_PIX2ANG_NEST<span class='expansion'>"GHealpix::pix2ang_nest(int, double*, double*)"</span></span>, ipix, 0, m_num_pixels-1);</td></tr>


    		
  

  

    
      838
    
    
      
<tr><td class="num" id="LN754">754</td><td class="line"> </td></tr>


    		
  

  

    
      839
    
    
      
<tr><td class="num" id="LN755">755</td><td class="line">    <span class='comment'>// Get face number and index in face</span></td></tr>


    		
  

  

    
      840
    
    
      
<tr><td class="num" id="LN756">756</td><td class="line">    <span class='keyword'>int</span> nl4      = 4 * m_nside;</td></tr>


    		
  

  

    
      841
    
    
      
<tr><td class="num" id="LN757">757</td><td class="line">    <span class='keyword'>int</span> face_num = ipix &gt;&gt; (2*m_order);      <span class='comment'>// Face number in {0,11}</span></td></tr>


    		
  

  

    
      842
    
    
      
<tr><td class="num" id="LN758">758</td><td class="line">    <span class='keyword'>int</span> ipf      = ipix &amp; (m_npface - 1);</td></tr>


    		
  

  

    
      843
    
    
      
<tr><td class="num" id="LN759">759</td><td class="line"> </td></tr>


    		
  

  

    
      844
    
    
      
<tr><td class="num" id="LN760">760</td><td class="line">    <span class='comment'>// Get pixel coordinates</span></td></tr>


    		
  

  

    
      845
    
    
      
<tr><td class="num" id="LN761">761</td><td class="line">    <span class='keyword'>int</span> ix;</td></tr>


    		
  

  

    
      846
    
    
      
<tr><td class="num" id="LN762">762</td><td class="line">    <span class='keyword'>int</span> iy;</td></tr>


    		
  

  

    
      847
    
    
      
<tr><td class="num" id="LN763">763</td><td class="line">    pix2xy(ipf, &amp;ix, &amp;iy);</td></tr>


    		
  

  

    
      848
    
    
      
<tr><td class="num" id="LN764">764</td><td class="line"> </td></tr>


    		
  

  

    
      849
    
    
      
<tr><td class="num" id="LN765">765</td><td class="line">    <span class='comment'>// Computes the z coordinate on the sphere</span></td></tr>


    		
  

  

    
      850
    
    
      
<tr><td class="num" id="LN766">766</td><td class="line">    <span class='keyword'>int</span> jr = (jrll[face_num] &lt;&lt; m_order) - ix - iy - 1;</td></tr>


    		
  

  

    
      851
    
    
      
<tr><td class="num" id="LN767">767</td><td class="line"> </td></tr>


    		
  

  

    
      852
    
    
      
<tr><td class="num" id="LN768">768</td><td class="line">    <span class='comment'>// Declare result variables</span></td></tr>


    		
  

  

    
      853
    
    
      
<tr><td class="num" id="LN769">769</td><td class="line">    <span class='keyword'>int</span>    nr;</td></tr>


    		
  

  

    
      854
    
    
      
<tr><td class="num" id="LN770">770</td><td class="line">    <span class='keyword'>double</span> z;</td></tr>


    		
  

  

    
      855
    
    
      
<tr><td class="num" id="LN771">771</td><td class="line">    <span class='keyword'>int</span>    kshift;</td></tr>


    		
  

  

    
      856
    
    
      
<tr><td class="num" id="LN772">772</td><td class="line"> </td></tr>


    		
  

  

    
      857
    
    
      
<tr><td class="num" id="LN773">773</td><td class="line">    <span class='comment'>// North pole region</span></td></tr>


    		
  

  

    
      858
    
    
      
<tr><td class="num" id="LN774">774</td><td class="line">    <span class='keyword'>if</span> (jr &lt; m_nside) {</td></tr>


    		
  

  

    
      859
    
    
      
<tr><td class="num" id="LN775">775</td><td class="line">        nr     = jr;</td></tr>


    		
  

  

    
      860
    
    
      
<tr><td class="num" id="LN776">776</td><td class="line">        z      = 1. - nr*nr*m_fact2;</td></tr>


    		
  

  

    
      861
    
    
      
<tr><td class="num" id="LN777">777</td><td class="line">        kshift = 0;</td></tr>


    		
  

  

    
      862
    
    
      
<tr><td class="num" id="LN778">778</td><td class="line">    }</td></tr>


    		
  

  

    
      863
    
    
      
<tr><td class="num" id="LN779">779</td><td class="line"> </td></tr>


    		
  

  

    
      864
    
    
      
<tr><td class="num" id="LN780">780</td><td class="line">    <span class='comment'>// South pole region</span></td></tr>


    		
  

  

    
      865
    
    
      
<tr><td class="num" id="LN781">781</td><td class="line">    <span class='keyword'>else</span> <span class='keyword'>if</span> (jr &gt; 3*m_nside) {</td></tr>


    		
  

  

    
      866
    
    
      
<tr><td class="num" id="LN782">782</td><td class="line">        nr     = nl4 - jr;</td></tr>


    		
  

  

    
      867
    
    
      
<tr><td class="num" id="LN783">783</td><td class="line">        z      = nr*nr*m_fact2 - 1;</td></tr>


    		
  

  

    
      868
    
    
      
<tr><td class="num" id="LN784">784</td><td class="line">        kshift = 0;</td></tr>


    		
  

  

    
      869
    
    
      
<tr><td class="num" id="LN785">785</td><td class="line">    }</td></tr>


    		
  

  

    
      870
    
    
      
<tr><td class="num" id="LN786">786</td><td class="line"> </td></tr>


    		
  

  

    
      871
    
    
      
<tr><td class="num" id="LN787">787</td><td class="line">    <span class='comment'>// Equatorial region</span></td></tr>


    		
  

  

    
      872
    
    
      
<tr><td class="num" id="LN788">788</td><td class="line">    <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      873
    
    
      
<tr><td class="num" id="LN789">789</td><td class="line">        nr     = m_nside;</td></tr>


    		
  

  

    
      874
    
    
      
<tr><td class="num" id="LN790">790</td><td class="line">        z      = (2*m_nside-jr) * m_fact1;</td></tr>


    		
  

  

    
      875
    
    
      
<tr><td class="num" id="LN791">791</td><td class="line">        kshift = (jr-m_nside) &amp; 1;</td></tr>


    		
  

  

    
      876
    
    
      
<tr><td class="num" id="LN792">792</td><td class="line">    }</td></tr>


    		
  

  

    
      877
    
    
      
<tr><td class="num" id="LN793">793</td><td class="line"> </td></tr>


    		
  

  

    
      878
    
    
      
<tr><td class="num" id="LN794">794</td><td class="line">    <span class='comment'>// Computes the phi coordinate on the sphere, in [0,2Pi]</span></td></tr>


    		
  

  

    
      879
    
    
      
<tr><td class="num" id="LN795">795</td><td class="line">    <span class='keyword'>int</span> jp = (jpll[face_num]*nr + ix - iy + 1 + kshift) / 2;</td></tr>


    		
  

  

    
      880
    
    
      
<tr><td class="num" id="LN796">796</td><td class="line">    <span class='keyword'>if</span> (jp &gt; nl4) jp -= nl4;</td></tr>


    		
  

  

    
      881
    
    
      
<tr><td class="num" id="LN797">797</td><td class="line">    <span class='keyword'>if</span> (jp &lt;   1) jp += nl4;</td></tr>


    		
  

  

    
      882
    
    
      
<tr><td class="num" id="LN798">798</td><td class="line"> </td></tr>


    		
  

  

    
      883
    
    
      
<tr><td class="num" id="LN799">799</td><td class="line">    <span class='comment'>// Computes Theta and Phi</span></td></tr>


    		
  

  

    
      884
    
    
      
<tr><td class="num" id="LN800">800</td><td class="line">    *theta = acos(z);</td></tr>


    		
  

  

    
      885
    
    
      
<tr><td class="num" id="LN801">801</td><td class="line">    *phi   = (jp - (kshift+1)*0.5) * (gammalib::pihalf / nr);</td></tr>


    		
  

  

    
      886
    
    
      
<tr><td class="num" id="LN802">802</td><td class="line"> </td></tr>


    		
  

  

    
      887
    
    
      
<tr><td class="num" id="LN803">803</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      888
    
    
      
<tr><td class="num" id="LN804">804</td><td class="line">    <span class='keyword'>return</span>;</td></tr>


    		
  

  

    
      889
    
    
      
<tr><td class="num" id="LN805">805</td><td class="line">}</td></tr>


    		
  

  

    
      890
    
    
      
<tr><td class="num" id="LN806">806</td><td class="line"> </td></tr>


    		
  

  

    
      891
    
    
      
<tr><td class="num" id="LN807">807</td><td class="line"> </td></tr>


    		
  

  

    
      892
    
    
      
<tr><td class="num" id="LN808">808</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      893
    
    
      
<tr><td class="num" id="LN809">809</td><td class="line"> <span class='comment'>* @brief Returns pixels which contains angular coordinates (z,phi)</span></td></tr>


    		
  

  

    
      894
    
    
      
<tr><td class="num" id="LN810">810</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      895
    
    
      
<tr><td class="num" id="LN811">811</td><td class="line"> <span class='comment'>* @param[in] z Cosine of zenith angle - cos(theta).</span></td></tr>


    		
  

  

    
      896
    
    
      
<tr><td class="num" id="LN812">812</td><td class="line"> <span class='comment'>* @param[in] phi Azimuth angle in radians.</span></td></tr>


    		
  

  

    
      897
    
    
      
<tr><td class="num" id="LN813">813</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      898
    
    
      
<tr><td class="num" id="LN814">814</td><td class="line"><span class='keyword'>int</span> GHealpix::ang2pix_z_phi_ring(<span class='keyword'>double</span> z, <span class='keyword'>double</span> phi) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      899
    
    
      
<tr><td class="num" id="LN815">815</td><td class="line">{</td></tr>


    		
  

  

    
      900
    
    
      
<tr><td class="num" id="LN816">816</td><td class="line">    <span class='comment'>// Initialise pixel</span></td></tr>


    		
  

  

    
      901
    
    
      
<tr><td class="num" id="LN817">817</td><td class="line">    <span class='keyword'>int</span> ipix = 0;</td></tr>


    		
  

  

    
      902
    
    
      
<tr><td class="num" id="LN818">818</td><td class="line"> </td></tr>


    		
  

  

    
      903
    
    
      
<tr><td class="num" id="LN819">819</td><td class="line">    <span class='comment'>// Setup</span></td></tr>


    		
  

  

    
      904
    
    
      
<tr><td class="num" id="LN820">820</td><td class="line">    <span class='keyword'>double</span> za = fabs(z);</td></tr>


    		
  

  

    
      905
    
    
      
<tr><td class="num" id="LN821">821</td><td class="line">    <span class='keyword'>double</span> tt = gammalib::modulo(phi, gammalib::twopi) * gammalib::inv_pihalf; <span class='comment'>// in [0,4)</span></td></tr>


    		
  

  

    
      906
    
    
      
<tr><td class="num" id="LN822">822</td><td class="line"> </td></tr>


    		
  

  

    
      907
    
    
      
<tr><td class="num" id="LN823">823</td><td class="line">    <span class='comment'>// Equatorial region</span></td></tr>


    		
  

  

    
      908
    
    
      
<tr><td class="num" id="LN824">824</td><td class="line">    <span class='keyword'>if</span> (za &lt;= gammalib::twothird) {</td></tr>


    		
  

  

    
      909
    
    
      
<tr><td class="num" id="LN825">825</td><td class="line">        <span class='keyword'>double</span> temp1  = m_nside*(0.5+tt);</td></tr>


    		
  

  

    
      910
    
    
      
<tr><td class="num" id="LN826">826</td><td class="line">        <span class='keyword'>double</span> temp2  = m_nside*z*0.75;</td></tr>


    		
  

  

    
      911
    
    
      
<tr><td class="num" id="LN827">827</td><td class="line">        <span class='keyword'>int</span>    jp     = <span class='keyword'>int</span>(temp1-temp2);           <span class='comment'>// index of ascending edge line</span></td></tr>


    		
  

  

    
      912
    
    
      
<tr><td class="num" id="LN828">828</td><td class="line">        <span class='keyword'>int</span>    jm     = <span class='keyword'>int</span>(temp1+temp2);           <span class='comment'>// index of descending edge line</span></td></tr>


    		
  

  

    
      913
    
    
      
<tr><td class="num" id="LN829">829</td><td class="line">        <span class='keyword'>int</span>    ir     = m_nside + 1 + jp - jm;      <span class='comment'>// in {1,2n+1}</span></td></tr>


    		
  

  

    
      914
    
    
      
<tr><td class="num" id="LN830">830</td><td class="line">        <span class='keyword'>int</span>    kshift = 1 - (ir &amp; 1);               <span class='comment'>// kshift=1 if ir even, 0 otherwise</span></td></tr>


    		
  

  

    
      915
    
    
      
<tr><td class="num" id="LN831">831</td><td class="line">        <span class='keyword'>int</span>    ip     = (jp+jm-m_nside+kshift+1)/2; <span class='comment'>// in {0,4n-1}</span></td></tr>


    		
  

  

    
      916
    
    
      
<tr><td class="num" id="LN832">832</td><td class="line">        ip            = <span class='keyword'>int</span>(gammalib::modulo(ip,4*m_nside));</td></tr>


    		
  

  

    
      917
    
    
      
<tr><td class="num" id="LN833">833</td><td class="line">        ipix          = m_ncap + (ir-1)*4*m_nside + ip;</td></tr>


    		
  

  

    
      918
    
    
      
<tr><td class="num" id="LN834">834</td><td class="line">    }</td></tr>


    		
  

  

    
      919
    
    
      
<tr><td class="num" id="LN835">835</td><td class="line"> </td></tr>


    		
  

  

    
      920
    
    
      
<tr><td class="num" id="LN836">836</td><td class="line">    <span class='comment'>// North &amp; South polar caps</span></td></tr>


    		
  

  

    
      921
    
    
      
<tr><td class="num" id="LN837">837</td><td class="line">    <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      922
    
    
      
<tr><td class="num" id="LN838">838</td><td class="line">        <span class='keyword'>double</span> tp  = tt - <span class='keyword'>int</span>(tt);</td></tr>


    		
  

  

    
      923
    
    
      
<tr><td class="num" id="LN839">839</td><td class="line">        <span class='keyword'>double</span> tmp = m_nside * std::sqrt(3*(1-za));</td></tr>


    		
  

  

    
      924
    
    
      
<tr><td class="num" id="LN840">840</td><td class="line">        <span class='keyword'>int</span>    jp  = <span class='keyword'>int</span>(tp*tmp);       <span class='comment'>// increasing edge line index</span></td></tr>


    		
  

  

    
      925
    
    
      
<tr><td class="num" id="LN841">841</td><td class="line">        <span class='keyword'>int</span>    jm  = <span class='keyword'>int</span>((1.0-tp)*tmp); <span class='comment'>// decreasing edge line index</span></td></tr>


    		
  

  

    
      926
    
    
      
<tr><td class="num" id="LN842">842</td><td class="line">        <span class='keyword'>int</span>    ir  = jp + jm + 1;       <span class='comment'>// ring number counted from the closest pole</span></td></tr>


    		
  

  

    
      927
    
    
      
<tr><td class="num" id="LN843">843</td><td class="line">        <span class='keyword'>int</span>    ip  = <span class='keyword'>int</span>(tt*ir);        <span class='comment'>// in {0,4*ir-1}</span></td></tr>


    		
  

  

    
      928
    
    
      
<tr><td class="num" id="LN844">844</td><td class="line">        ip = <span class='keyword'>int</span>(gammalib::modulo(ip,4*ir));</td></tr>


    		
  

  

    
      929
    
    
      
<tr><td class="num" id="LN845">845</td><td class="line">        <span class='keyword'>if</span> (z&gt;0)</td></tr>


    		
  

  

    
      930
    
    
      
<tr><td class="num" id="LN846">846</td><td class="line">            ipix = 2*ir*(ir-1) + ip;</td></tr>


    		
  

  

    
      931
    
    
      
<tr><td class="num" id="LN847">847</td><td class="line">        <span class='keyword'>else</span></td></tr>


    		
  

  

    
      932
    
    
      
<tr><td class="num" id="LN848">848</td><td class="line">            ipix = m_num_pixels - 2*ir*(ir+1) + ip;</td></tr>


    		
  

  

    
      933
    
    
      
<tr><td class="num" id="LN849">849</td><td class="line">    }</td></tr>


    		
  

  

    
      934
    
    
      
<tr><td class="num" id="LN850">850</td><td class="line"> </td></tr>


    		
  

  

    
      935
    
    
      
<tr><td class="num" id="LN851">851</td><td class="line">    <span class='comment'>// Return pixel</span></td></tr>


    		
  

  

    
      936
    
    
      
<tr><td class="num" id="LN852">852</td><td class="line">    <span class='keyword'>return</span> ipix;</td></tr>


    		
  

  

    
      937
    
    
      
<tr><td class="num" id="LN853">853</td><td class="line">}</td></tr>


    		
  

  

    
      938
    
    
      
<tr><td class="num" id="LN854">854</td><td class="line"> </td></tr>


    		
  

  

    
      939
    
    
      
<tr><td class="num" id="LN855">855</td><td class="line"> </td></tr>


    		
  

  

    
      940
    
    
      
<tr><td class="num" id="LN856">856</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      941
    
    
      
<tr><td class="num" id="LN857">857</td><td class="line"> <span class='comment'>* @brief Returns pixels which contains angular coordinates (z,phi)</span></td></tr>


    		
  

  

    
      942
    
    
      
<tr><td class="num" id="LN858">858</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      943
    
    
      
<tr><td class="num" id="LN859">859</td><td class="line"> <span class='comment'>* @param[in] z Cosine of zenith angle - cos(theta).</span></td></tr>


    		
  

  

    
      944
    
    
      
<tr><td class="num" id="LN860">860</td><td class="line"> <span class='comment'>* @param[in] phi Azimuth angle in radians.</span></td></tr>


    		
  

  

    
      945
    
    
      
<tr><td class="num" id="LN861">861</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      946
    
    
      
<tr><td class="num" id="LN862">862</td><td class="line"><span class='keyword'>int</span> GHealpix::ang2pix_z_phi_nest(<span class='keyword'>double</span> z, <span class='keyword'>double</span> phi) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      947
    
    
      
<tr><td class="num" id="LN863">863</td><td class="line">{</td></tr>


    		
  

  

    
      948
    
    
      
<tr><td class="num" id="LN864">864</td><td class="line">    <span class='comment'>// Initialise face and pixel numbers</span></td></tr>


    		
  

  

    
      949
    
    
      
<tr><td class="num" id="LN865">865</td><td class="line">    <span class='keyword'>int</span> face_num;</td></tr>


    		
  

  

    
      950
    
    
      
<tr><td class="num" id="LN866">866</td><td class="line">    <span class='keyword'>int</span> ix;</td></tr>


    		
  

  

    
      951
    
    
      
<tr><td class="num" id="LN867">867</td><td class="line">    <span class='keyword'>int</span> iy;</td></tr>


    		
  

  

    
      952
    
    
      
<tr><td class="num" id="LN868">868</td><td class="line"> </td></tr>


    		
  

  

    
      953
    
    
      
<tr><td class="num" id="LN869">869</td><td class="line">    <span class='comment'>// Setup</span></td></tr>


    		
  

  

    
      954
    
    
      
<tr><td class="num" id="LN870">870</td><td class="line">    <span class='keyword'>double</span> za = fabs(z);</td></tr>


    		
  

  

    
      955
    
    
      
<tr><td class="num" id="LN871">871</td><td class="line">    <span class='keyword'>double</span> tt = gammalib::modulo(phi, gammalib::twopi) * gammalib::inv_pihalf; <span class='comment'>// in [0,4)</span></td></tr>


    		
  

  

    
      956
    
    
      
<tr><td class="num" id="LN872">872</td><td class="line"> </td></tr>


    		
  

  

    
      957
    
    
      
<tr><td class="num" id="LN873">873</td><td class="line">    <span class='comment'>// Equatorial region</span></td></tr>


    		
  

  

    
      958
    
    
      
<tr><td class="num" id="LN874">874</td><td class="line">    <span class='keyword'>if</span> (za &lt;= gammalib::twothird) {</td></tr>


    		
  

  

    
      959
    
    
      
<tr><td class="num" id="LN875">875</td><td class="line">        <span class='keyword'>double</span> temp1 = ns_max*(0.5+tt);</td></tr>


    		
  

  

    
      960
    
    
      
<tr><td class="num" id="LN876">876</td><td class="line">        <span class='keyword'>double</span> temp2 = ns_max*z*0.75;</td></tr>


    		
  

  

    
      961
    
    
      
<tr><td class="num" id="LN877">877</td><td class="line">        <span class='keyword'>int</span>    jp    = <span class='keyword'>int</span>(temp1-temp2); <span class='comment'>// index of  ascending edge line</span></td></tr>


    		
  

  

    
      962
    
    
      
<tr><td class="num" id="LN878">878</td><td class="line">        <span class='keyword'>int</span>    jm    = <span class='keyword'>int</span>(temp1+temp2); <span class='comment'>// index of descending edge line</span></td></tr>


    		
  

  

    
      963
    
    
      
<tr><td class="num" id="LN879">879</td><td class="line">        <span class='keyword'>int</span>    ifp   = jp &gt;&gt; order_max;  <span class='comment'>// in {0,4}</span></td></tr>


    		
  

  

    
      964
    
    
      
<tr><td class="num" id="LN880">880</td><td class="line">        <span class='keyword'>int</span>    ifm   = jm &gt;&gt; order_max;</td></tr>


    		
  

  

    
      965
    
    
      
<tr><td class="num" id="LN881">881</td><td class="line">        <span class='keyword'>if</span> (ifp == ifm)                  <span class='comment'>// faces 4 to 7</span></td></tr>


    		
  

  

    
      966
    
    
      
<tr><td class="num" id="LN882">882</td><td class="line">            face_num = (ifp==4) ? 4: ifp+4;</td></tr>


    		
  

  

    
      967
    
    
      
<tr><td class="num" id="LN883">883</td><td class="line">        <span class='keyword'>else</span> <span class='keyword'>if</span> (ifp &lt; ifm)              <span class='comment'>// (half-)faces 0 to 3</span></td></tr>


    		
  

  

    
      968
    
    
      
<tr><td class="num" id="LN884">884</td><td class="line">            face_num = ifp;</td></tr>


    		
  

  

    
      969
    
    
      
<tr><td class="num" id="LN885">885</td><td class="line">        <span class='keyword'>else</span>                             <span class='comment'>// (half-)faces 8 to 11</span></td></tr>


    		
  

  

    
      970
    
    
      
<tr><td class="num" id="LN886">886</td><td class="line">            face_num = ifm + 8;</td></tr>


    		
  

  

    
      971
    
    
      
<tr><td class="num" id="LN887">887</td><td class="line">        ix = jm &amp; (ns_max-1);</td></tr>


    		
  

  

    
      972
    
    
      
<tr><td class="num" id="LN888">888</td><td class="line">        iy = ns_max - (jp &amp; (ns_max-1)) - 1;</td></tr>


    		
  

  

    
      973
    
    
      
<tr><td class="num" id="LN889">889</td><td class="line">    }</td></tr>


    		
  

  

    
      974
    
    
      
<tr><td class="num" id="LN890">890</td><td class="line"> </td></tr>


    		
  

  

    
      975
    
    
      
<tr><td class="num" id="LN891">891</td><td class="line">    <span class='comment'>// Polar region, za &gt; 2/3</span></td></tr>


    		
  

  

    
      976
    
    
      
<tr><td class="num" id="LN892">892</td><td class="line">    <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      977
    
    
      
<tr><td class="num" id="LN893">893</td><td class="line">        <span class='keyword'>int</span>    ntt = <span class='keyword'>int</span>(tt);</td></tr>


    		
  

  

    
      978
    
    
      
<tr><td class="num" id="LN894">894</td><td class="line">        <span class='keyword'>double</span> tp  = tt-ntt;</td></tr>


    		
  

  

    
      979
    
    
      
<tr><td class="num" id="LN895">895</td><td class="line">        <span class='keyword'>double</span> tmp = ns_max * std::sqrt(3*(1-za));</td></tr>


    		
  

  

    
      980
    
    
      
<tr><td class="num" id="LN896">896</td><td class="line">        <span class='keyword'>int</span>    jp  = <span class='keyword'>int</span>(tp*tmp);        <span class='comment'>// increasing edge line index</span></td></tr>


    		
  

  

    
      981
    
    
      
<tr><td class="num" id="LN897">897</td><td class="line">        <span class='keyword'>int</span>    jm  = <span class='keyword'>int</span>((1.0-tp)*tmp);  <span class='comment'>// decreasing edge line index</span></td></tr>


    		
  

  

    
      982
    
    
      
<tr><td class="num" id="LN898">898</td><td class="line">        <span class='keyword'>if</span> (jp &gt;= ns_max) jp = ns_max-1; <span class='comment'>// for points too close to the boundary</span></td></tr>


    		
  

  

    
      983
    
    
      
<tr><td class="num" id="LN899">899</td><td class="line">        <span class='keyword'>if</span> (jm &gt;= ns_max) jm = ns_max-1;</td></tr>


    		
  

  

    
      984
    
    
      
<tr><td class="num" id="LN900">900</td><td class="line">        <span class='keyword'>if</span> (z &gt;= 0) {</td></tr>


    		
  

  

    
      985
    
    
      
<tr><td class="num" id="LN901">901</td><td class="line">            face_num = ntt;              <span class='comment'>// in {0,3}</span></td></tr>


    		
  

  

    
      986
    
    
      
<tr><td class="num" id="LN902">902</td><td class="line">            ix       = ns_max - jm - 1;</td></tr>


    		
  

  

    
      987
    
    
      
<tr><td class="num" id="LN903">903</td><td class="line">            iy       = ns_max - jp - 1;</td></tr>


    		
  

  

    
      988
    
    
      
<tr><td class="num" id="LN904">904</td><td class="line">        }</td></tr>


    		
  

  

    
      989
    
    
      
<tr><td class="num" id="LN905">905</td><td class="line">        <span class='keyword'>else</span> {</td></tr>


    		
  

  

    
      990
    
    
      
<tr><td class="num" id="LN906">906</td><td class="line">            face_num = ntt + 8;          <span class='comment'>// in {8,11}</span></td></tr>


    		
  

  

    
      991
    
    
      
<tr><td class="num" id="LN907">907</td><td class="line">            ix       =  jp;</td></tr>


    		
  

  

    
      992
    
    
      
<tr><td class="num" id="LN908">908</td><td class="line">            iy       =  jm;</td></tr>


    		
  

  

    
      993
    
    
      
<tr><td class="num" id="LN909">909</td><td class="line">        }</td></tr>


    		
  

  

    
      994
    
    
      
<tr><td class="num" id="LN910">910</td><td class="line">    }</td></tr>


    		
  

  

    
      995
    
    
      
<tr><td class="num" id="LN911">911</td><td class="line"> </td></tr>


    		
  

  

    
      996
    
    
      
<tr><td class="num" id="LN912">912</td><td class="line">    <span class='comment'>// Get pixel</span></td></tr>


    		
  

  

    
      997
    
    
      
<tr><td class="num" id="LN913">913</td><td class="line">    <span class='keyword'>int</span> ipf = xy2pix(ix, iy);</td></tr>


    		
  

  

    
      998
    
    
      
<tr><td class="num" id="LN914">914</td><td class="line">    ipf &gt;&gt;= (2*(order_max - m_order));     <span class='comment'>// in {0, nside**2 - 1}</span></td></tr>


    		
  

  

    
      999
    
    
      
<tr><td class="num" id="LN915">915</td><td class="line">    <span class='keyword'>return</span> ipf + (face_num&lt;&lt;(2*m_order));  <span class='comment'>// in {0, 12*nside**2 - 1}</span></td></tr>


    		
  

  

    
      1000
    
    
      
<tr><td class="num" id="LN916">916</td><td class="line">}</td></tr>


    		
  

  

    
      1001
    
    
      
<tr><td class="num" id="LN917">917</td><td class="line"> </td></tr>


    		
  

  

    
      1002
    
    
      
<tr><td class="num" id="LN918">918</td><td class="line"> </td></tr>


    		
  

  

    
      1003
    
    
      
<tr><td class="num" id="LN919">919</td><td class="line"><span class='comment'>/***********************************************************************/</span><span class='comment'>/**</span></td></tr>


    		
  

  

    
      1004
    
    
      
<tr><td class="num" id="LN920">920</td><td class="line"> <span class='comment'>* @brief Integer n that fulfills n*n &lt;= arg &lt; (n+1)*(n+1)</span></td></tr>


    		
  

  

    
      1005
    
    
      
<tr><td class="num" id="LN921">921</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      1006
    
    
      
<tr><td class="num" id="LN922">922</td><td class="line"> <span class='comment'>* @param[in] arg Argument.</span></td></tr>


    		
  

  

    
      1007
    
    
      
<tr><td class="num" id="LN923">923</td><td class="line"> <span class='comment'>*</span></td></tr>


    		
  

  

    
      1008
    
    
      
<tr><td class="num" id="LN924">924</td><td class="line"> <span class='comment'>* Returns the integer \a n, which fulfills \a n*n &lt;= arg &lt; (n+1)*(n+1).</span></td></tr>


    		
  

  

    
      1009
    
    
      
<tr><td class="num" id="LN925">925</td><td class="line"> <span class='comment'>***************************************************************************/</span></td></tr>


    		
  

  

    
      1010
    
    
      
<tr><td class="num" id="LN926">926</td><td class="line"><span class='keyword'>unsigned</span> <span class='keyword'>int</span> GHealpix::isqrt(<span class='keyword'>unsigned</span> <span class='keyword'>int</span> arg) <span class='keyword'>const</span></td></tr>


    		
  

  

    
      1011
    
    
      
<tr><td class="num" id="LN927">927</td><td class="line">{</td></tr>


    		
  

  

    
      1012
    
    
      
<tr><td class="num" id="LN928">928</td><td class="line">    <span class='comment'>// Return</span></td></tr>


    		
  

  

    
      1013
    
    
      
<tr><td class="num" id="LN929">929</td><td class="line">    <span class='keyword'>return</span> <span class='keyword'>unsigned</span>(std::sqrt(arg+0.5));</td></tr>


    		
  

  

    
      1014
    
    
      
<tr><td class="num" id="LN930">930</td><td class="line">}</td></tr>


    		
  

  

    
      1015
    
    
      
</table></body></html>