lorene-doc 0.0.0~cvs20161116+dfsg-1 / usr / share / doc / lorene-doc / refguide / scalar__import__symy_8C_source.html

This file is indexed.

/usr/share/doc/lorene-doc/refguide/scalar__import__symy_8C_source.html is in lorene-doc 0.0.0~cvs20161116+dfsg-1.

This file is owned by root:root, with mode 0o644.

The actual contents of the file can be viewed below. 

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.12"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>LORENE: scalar_import_symy.C Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
  $(document).ready(initResizable);
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  <td id="projectalign" style="padding-left: 0.5em;">
   <div id="projectname">LORENE
   </div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.12 -->
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
$(function() {
  initMenu('',false,false,'search.php','Search');
});
</script>
<div id="main-nav"></div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
      <div id="nav-sync" class="sync"></div>
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
$(document).ready(function(){initNavTree('scalar__import__symy_8C_source.html','');});
</script>
<div id="doc-content">
<div class="header">
  <div class="headertitle">
<div class="title">scalar_import_symy.C</div>  </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="comment">/*</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="comment"> * Member function of the Scalar class for initiating a Scalar from </span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment"> * a Scalar defined on another mapping.</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="comment"> * Case where both Scalar&#39;s are symmetric with respect to their y=0 plane.</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment"> */</span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;</div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="comment">/*</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="comment"> *   Copyright (c) 2003 Eric Gourgoulhon &amp; Jerome Novak</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="comment"> *   Copyright (c) 1999-2001 Eric Gourgoulhon (Cmp version)</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="comment"> *   This file is part of LORENE.</span></div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="comment"> *   LORENE is free software; you can redistribute it and/or modify</span></div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="comment"> *   it under the terms of the GNU General Public License as published by</span></div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="comment"> *   the Free Software Foundation; either version 2 of the License, or</span></div><div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="comment"> *   (at your option) any later version.</span></div><div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;<span class="comment"> *   LORENE is distributed in the hope that it will be useful,</span></div><div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="comment"> *   but WITHOUT ANY WARRANTY; without even the implied warranty of</span></div><div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;<span class="comment"> *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the</span></div><div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;<span class="comment"> *   GNU General Public License for more details.</span></div><div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;<span class="comment"> *   You should have received a copy of the GNU General Public License</span></div><div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;<span class="comment"> *   along with LORENE; if not, write to the Free Software</span></div><div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;<span class="comment"> *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA</span></div><div class="line"><a name="l00026"></a><span class="lineno">   26</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;<span class="comment"> */</span></div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;</div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;</div><div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;<span class="keywordtype">char</span> scalar_import_symy_C[] = <span class="stringliteral">&quot;$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_import_symy.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $&quot;</span> ;</div><div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;</div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;</div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;<span class="comment">/*</span></div><div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;<span class="comment"> * $Id: scalar_import_symy.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $</span></div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;<span class="comment"> * $Log: scalar_import_symy.C,v $</span></div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;<span class="comment"> * Revision 1.5  2014/10/13 08:53:46  j_novak</span></div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;<span class="comment"> * Lorene classes and functions now belong to the namespace Lorene.</span></div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;<span class="comment"> * Revision 1.4  2014/10/06 15:16:15  j_novak</span></div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;<span class="comment"> * Modified #include directives to use c++ syntax.</span></div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;<span class="comment"> * Revision 1.3  2003/10/10 15:57:29  j_novak</span></div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;<span class="comment"> * Added the state one (ETATUN) to the class Scalar</span></div><div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;<span class="comment"> * Revision 1.2  2003/10/01 13:04:44  e_gourgoulhon</span></div><div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;<span class="comment"> * The method Tensor::get_mp() returns now a reference (and not</span></div><div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;<span class="comment"> * a pointer) onto a mapping.</span></div><div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;<span class="comment"> * Revision 1.1  2003/09/25 09:07:05  j_novak</span></div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;<span class="comment"> * Added the functions for importing from another mapping (to be tested).</span></div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;<span class="comment"> * $Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_import_symy.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $</span></div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;<span class="comment"> *</span></div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;<span class="comment"> */</span></div><div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;</div><div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;<span class="comment">// Headers C</span></div><div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;<span class="preprocessor">#include &lt;cmath&gt;</span></div><div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;</div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;<span class="comment">// Headers Lorene</span></div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;<span class="preprocessor">#include &quot;tensor.h&quot;</span></div><div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;<span class="preprocessor">#include &quot;param.h&quot;</span></div><div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;<span class="preprocessor">#include &quot;nbr_spx.h&quot;</span></div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;</div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;            <span class="comment">//-------------------------------//</span></div><div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;            <span class="comment">//  Importation in all domains   //</span></div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;            <span class="comment">//-------------------------------//</span></div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;<span class="keyword">namespace </span><a class="code" href="namespaceLorene.html">Lorene</a> {</div><div class="line"><a name="l00072"></a><span class="lineno"><a class="line" href="classLorene_1_1Scalar.html#af17554dd17d8a77098b4f743ed86e38b">   72</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="classLorene_1_1Scalar.html#af17554dd17d8a77098b4f743ed86e38b">Scalar::import_symy</a>(<span class="keyword">const</span> <a class="code" href="classLorene_1_1Scalar.html">Scalar</a>&amp; ci) {</div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;    </div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;    <span class="keywordtype">int</span> nz = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4df8c6a943a85ef314ea636240ef31a7">get_mg</a>()-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a99fb0e7a3f1f0f3e405624855200055e">get_nzone</a>() ; </div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;    </div><div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#af17554dd17d8a77098b4f743ed86e38b">import_symy</a>(nz, ci) ; </div><div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;        </div><div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;}</div><div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;</div><div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;            <span class="comment">//--------------------------------------//</span></div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;            <span class="comment">//  Importation in inner domains only   //</span></div><div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;            <span class="comment">//--------------------------------------//</span></div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;</div><div class="line"><a name="l00084"></a><span class="lineno"><a class="line" href="classLorene_1_1Scalar.html#a8e44f80b2a24968ba3417a328828ebb1">   84</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="classLorene_1_1Scalar.html#af17554dd17d8a77098b4f743ed86e38b">Scalar::import_symy</a>(<span class="keywordtype">int</span> nzet, <span class="keyword">const</span> <a class="code" href="classLorene_1_1Scalar.html">Scalar</a>&amp; cm_d) {</div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;    </div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Map.html">Map</a>* mp_d = &amp;(cm_d.<a class="code" href="classLorene_1_1Tensor.html#a00ed03bf0b7364e5bec0cda699b07f26">get_mp</a>()) ; <span class="comment">// Departure mapping</span></div><div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;</div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;    <span class="comment">// Trivial case : mappings identical !</span></div><div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;    <span class="comment">// -----------------------------------</span></div><div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;    </div><div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;    <span class="keywordflow">if</span> (mp_d == <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>) {</div><div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;    *<span class="keyword">this</span> = cm_d ; </div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;    return ; </div><div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;    }</div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;    </div><div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;    <span class="comment">// Relative orientation of the two mappings</span></div><div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;    <span class="comment">// ----------------------------------------</span></div><div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160; </div><div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;    <span class="keywordtype">int</span> align_rel = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align()  </div><div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;            * (mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align() ;</div><div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;            </div><div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;    <span class="keywordflow">switch</span> (align_rel) {</div><div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160;</div><div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;    <span class="keywordflow">case</span> 1 : {  <span class="comment">// the two mappings have aligned Cartesian axis</span></div><div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;        <a class="code" href="classLorene_1_1Scalar.html#aa79d50444dc989294d4407354cc8c66b">import_align_symy</a>(nzet, cm_d) ; </div><div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;        break ; </div><div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;    }</div><div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;</div><div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;    <span class="keywordflow">case</span> -1 : {  <span class="comment">// the two mappings have anti-aligned Cartesian axis</span></div><div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;        <a class="code" href="classLorene_1_1Scalar.html#ae2425d68e625b5345ab518c058e5417b">import_anti_symy</a>(nzet, cm_d) ; </div><div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;        break ; </div><div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;    }</div><div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;</div><div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;    <span class="keywordflow">default</span> : {  </div><div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;        cout &lt;&lt; <span class="stringliteral">&quot;Scalar::import_symy : unexpected value of align_rel : &quot;</span> </div><div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;         &lt;&lt; align_rel &lt;&lt; endl ; </div><div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160;        abort() ; </div><div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;        break ; </div><div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;    }</div><div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;</div><div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;    }</div><div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;    </div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;}    </div><div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;</div><div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;</div><div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;        <span class="comment">//-----------------------------------------//</span></div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;        <span class="comment">//   Case of Cartesian axis anti-aligned   //</span></div><div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;        <span class="comment">//-----------------------------------------//</span></div><div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;</div><div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;</div><div class="line"><a name="l00131"></a><span class="lineno"><a class="line" href="classLorene_1_1Scalar.html#ae2425d68e625b5345ab518c058e5417b">  131</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="classLorene_1_1Scalar.html#ae2425d68e625b5345ab518c058e5417b">Scalar::import_anti_symy</a>(<span class="keywordtype">int</span> nzet, <span class="keyword">const</span> <a class="code" href="classLorene_1_1Scalar.html">Scalar</a>&amp; cm_d) {</div><div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;    </div><div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;    <span class="comment">// Trivial case : null Scalar</span></div><div class="line"><a name="l00134"></a><span class="lineno">  134</span>&#160;    <span class="comment">// ------------------------</span></div><div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160;    </div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATZERO) {</div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a68abd6fe37c8829980599db324fcd725">set_etat_zero</a>() ; </div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;    return ; </div><div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;    }</div><div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATUN) {</div><div class="line"><a name="l00141"></a><span class="lineno">  141</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a089209b1af70621abbf5ce370ecde0b7">set_etat_one</a>() ; </div><div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;    return ; </div><div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;    }</div><div class="line"><a name="l00144"></a><span class="lineno">  144</span>&#160;</div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Map.html">Map</a>* mp_d = &amp;(cm_d.<a class="code" href="classLorene_1_1Tensor.html#a00ed03bf0b7364e5bec0cda699b07f26">get_mp</a>()) ; <span class="comment">// Departure mapping</span></div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;</div><div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;    <span class="comment">// Protections</span></div><div class="line"><a name="l00148"></a><span class="lineno">  148</span>&#160;    <span class="comment">// -----------</span></div><div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;    <span class="keywordtype">int</span> align = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align() ;</div><div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;</div><div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;    assert( align  * (mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align() == -1 ) ; </div><div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;    </div><div class="line"><a name="l00153"></a><span class="lineno">  153</span>&#160;    assert(cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATQCQ) ;</div><div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;</div><div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a309bc82975514dbc9df435f7624e25fd">get_dzpuis</a>() != 0) {</div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;    cout &lt;&lt; </div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;    <span class="stringliteral">&quot;Scalar::import_anti_symy : the dzpuis of the Scalar to be imported&quot;</span></div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;    &lt;&lt; <span class="stringliteral">&quot;  must be zero !&quot;</span> &lt;&lt; endl ; </div><div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;    abort() ; </div><div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;    }</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;</div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;</div><div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mg3d.html">Mg3d</a>* mg_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4df8c6a943a85ef314ea636240ef31a7">get_mg</a>() ; </div><div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;    assert(mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a6fc60f226515b2c3ef3959ed9db0e1ce">get_type_p</a>() == NONSYM) ; </div><div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;</div><div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;    <span class="keywordtype">int</span> nz_a = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a99fb0e7a3f1f0f3e405624855200055e">get_nzone</a>() ; </div><div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;    assert(nzet &lt;= nz_a) ;     </div><div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;</div><div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Valeur.html">Valeur</a>&amp; va_d = cm_d.<a class="code" href="classLorene_1_1Scalar.html#aa8b9b16249b758df4fe8e1ad8689b25e">get_spectral_va</a>() ; </div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;    va_d.<a class="code" href="classLorene_1_1Valeur.html#aae49809b7bbd41b848e873ec076c025a">coef</a>() ;       <span class="comment">// The coefficients are required</span></div><div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;    </div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;</div><div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;    <span class="comment">// Preparations for storing the result in *this </span></div><div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;    <span class="comment">// --------------------------------------------</span></div><div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a908fc51b1306d30ba5525198e803f246">del_t</a>() ;   <span class="comment">// delete all previously computed derived quantities</span></div><div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;    </div><div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a39f70d51538ce68155c311f2eebeb82a">set_etat_qcq</a>() ;        <span class="comment">// Set the state to ETATQCQ</span></div><div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;    </div><div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a1e52288a6fd8f61dcf2a6c202756094e">set_etat_c_qcq</a>() ;   <span class="comment">// Allocates the memory for the Mtbl va.c</span></div><div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;                <span class="comment">//  if it does not exist already</span></div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#ab765d0b25416a6da5add6d2ea6b44acc">set_etat_qcq</a>() ;  <span class="comment">// Allocates the memory for the Tbl&#39;s in each</span></div><div class="line"><a name="l00182"></a><span class="lineno">  182</span>&#160;                <span class="comment">//  domain if they do not exist already</span></div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;</div><div class="line"><a name="l00184"></a><span class="lineno">  184</span>&#160;</div><div class="line"><a name="l00185"></a><span class="lineno">  185</span>&#160;    <span class="comment">// Departure (x,y,z) coordinates of the origin of the Arrival mapping :</span></div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;</div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;    <span class="keywordtype">double</span> xx_a, yy_a, zz_a ; </div><div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;    <span class="keywordflow">if</span> (align == 1) {</div><div class="line"><a name="l00189"></a><span class="lineno">  189</span>&#160;    xx_a = mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() - <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() ; </div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;    yy_a = mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() - <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() ; </div><div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;    }</div><div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;    xx_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() ; </div><div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;    yy_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() ; </div><div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;    }</div><div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;    zz_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a0159bdb0fe9f7fa3a78690ecce1658f1">get_ori_z</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a0159bdb0fe9f7fa3a78690ecce1658f1">get_ori_z</a>() ; </div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;</div><div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;    </div><div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;    <span class="comment">// r, theta, phi, x, y and z on the Arrival mapping </span></div><div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;    <span class="comment">//  update of the corresponding Coord&#39;s if necessary</span></div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;    </div><div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).fait() ; </div><div class="line"><a name="l00203"></a><span class="lineno">  203</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).fait() ; </div><div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).fait() ; </div><div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).fait() ; </div><div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).fait() ; </div><div class="line"><a name="l00207"></a><span class="lineno">  207</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).fait() ; </div><div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;</div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mr_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).c ; </div><div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mtet_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).c ; </div><div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mphi_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).c ; </div><div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mx_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).c ;</div><div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* my_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).c ;</div><div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mz_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).c ;</div><div class="line"><a name="l00215"></a><span class="lineno">  215</span>&#160;</div><div class="line"><a name="l00216"></a><span class="lineno">  216</span>&#160;    <a class="code" href="classLorene_1_1Param.html">Param</a> par_precis ;      <span class="comment">// Required precision in the method Map::val_lx</span></div><div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;    <span class="keywordtype">int</span> nitermax = 100 ;    <span class="comment">// Maximum number of iteration in the secant method</span></div><div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;    <span class="keywordtype">int</span> niter ; </div><div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;    <span class="keywordtype">double</span> precis = 1e-15 ; <span class="comment">// Absolute precision in the secant method</span></div><div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#aefd8de56123122b81d919e877a7a5cb7">add_int</a>(nitermax) ;  </div><div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#a8ea6e8bc0a7166c622764641bbc3744e">add_int_mod</a>(niter) ; </div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#a049f47c98bc2f7162a94b5a8c9215fde">add_double</a>(precis) ; </div><div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;</div><div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160;</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;    <span class="comment">// Loop of the Arrival domains where the computation is to be performed</span></div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;    <span class="comment">// --------------------------------------------------------------------</span></div><div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;    </div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> l=0; l &lt; nzet; l++) {</div><div class="line"><a name="l00229"></a><span class="lineno">  229</span>&#160;    </div><div class="line"><a name="l00230"></a><span class="lineno">  230</span>&#160;    <span class="keywordtype">int</span> nr = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#aad4b3f34fc58cc34de13540f3d37699a">get_nr</a>(l) ;</div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;    <span class="keywordtype">int</span> nt = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a7b8408df0cb45eb15d7ba838ee82e2ac">get_nt</a>(l) ;</div><div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;    <span class="keywordtype">int</span> np = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a6af09c50edf1dacdf9f9de24478746a4">get_np</a>(l) ;</div><div class="line"><a name="l00233"></a><span class="lineno">  233</span>&#160;    </div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;    </div><div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pr_a = mr_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of r</span></div><div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* ptet_a = mtet_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;  <span class="comment">// Pointer on the values of theta</span></div><div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pphi_a = mphi_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;  <span class="comment">// Pointer on the values of phi</span></div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* px_a = mx_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of X</span></div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* py_a = my_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of Y</span></div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pz_a = mz_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of Z</span></div><div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;    </div><div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;    (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;<a class="code" href="classLorene_1_1Scalar.html#a39f70d51538ce68155c311f2eebeb82a">set_etat_qcq</a>() ;       <span class="comment">// Allocates the array of double to</span></div><div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;                         <span class="comment">//  store the result </span></div><div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;    <span class="keywordtype">double</span>* ptx = (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;t ;      <span class="comment">// Pointer on the allocated array</span></div><div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;    </div><div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;    </div><div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;    <span class="comment">// Loop on half the grid points in the considered arrival domain</span></div><div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;    <span class="comment">// (the other half will be obtained by symmetry with respect to</span></div><div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;    <span class="comment">//  the y=0 plane). </span></div><div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;        </div><div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0 ; k&lt;np/2+1 ; k++) {    <span class="comment">// np/2+1 : half the grid</span></div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0 ; j&lt;nt ; j++) {</div><div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0 ; i&lt;nr ; i++) {</div><div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;</div><div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;            <span class="keywordtype">double</span> r = *pr_a ; </div><div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;            <span class="keywordtype">double</span> rd, tetd, phid ;</div><div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;            <span class="keywordflow">if</span> (r == __infinity) {</div><div class="line"><a name="l00258"></a><span class="lineno">  258</span>&#160;            rd = r ;</div><div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;            tetd = *ptet_a ;</div><div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;            phid = *pphi_a + M_PI ;</div><div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;            <span class="keywordflow">if</span> (phid &lt; 0) phid += 2*M_PI ;</div><div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;            }</div><div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;            <span class="keywordflow">else</span> {</div><div class="line"><a name="l00264"></a><span class="lineno">  264</span>&#160;</div><div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;            <span class="comment">// Cartesian coordinates on the Departure mapping</span></div><div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;            <span class="keywordtype">double</span> xd = - *px_a + xx_a ; </div><div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;            <span class="keywordtype">double</span> yd = - *py_a + yy_a ; </div><div class="line"><a name="l00268"></a><span class="lineno">  268</span>&#160;            <span class="keywordtype">double</span> zd = *pz_a + zz_a ; </div><div class="line"><a name="l00269"></a><span class="lineno">  269</span>&#160;</div><div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;            <span class="comment">// Spherical coordinates on the Departure mapping</span></div><div class="line"><a name="l00271"></a><span class="lineno">  271</span>&#160;            <span class="keywordtype">double</span> rhod2 = xd*xd + yd*yd ; </div><div class="line"><a name="l00272"></a><span class="lineno">  272</span>&#160;            <span class="keywordtype">double</span> rhod = <a class="code" href="classLorene_1_1Scalar.html#aae28013acc7f65d03bf107c0ded5f73e">sqrt</a>( rhod2 ) ;</div><div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;            rd = <a class="code" href="classLorene_1_1Scalar.html#aae28013acc7f65d03bf107c0ded5f73e">sqrt</a>(rhod2 + zd*zd) ;</div><div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;            tetd = atan2(rhod, zd) ;</div><div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;            phid = atan2(yd, xd) ; </div><div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;            <span class="keywordflow">if</span> (phid &lt; 0) phid += 2*M_PI ;</div><div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;            }           </div><div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;</div><div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;</div><div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;            <span class="comment">// NB: to increase the efficiency, the method Scalar::val_point</span></div><div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;            <span class="comment">//  is not invoked; the method Mtbl_cf::val_point is </span></div><div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;            <span class="comment">//  called directly instead. </span></div><div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;</div><div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;            <span class="comment">// Value of the grid coordinates (l,xi) corresponding to</span></div><div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;            <span class="comment">//  (rd,tetd,phid) :</span></div><div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;            </div><div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;            <span class="keywordtype">int</span> ld ;        <span class="comment">// domain index</span></div><div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;            <span class="keywordtype">double</span> xxd ;    <span class="comment">// radial coordinate xi in [0,1] or [-1,1]</span></div><div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;            mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a835fb3c3bfefe88135160b74ef7010b4">val_lx</a>(rd, tetd, phid, par_precis, ld, xxd) ;</div><div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;</div><div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;            <span class="comment">// Value of the Departure Scalar at the obtained point:</span></div><div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;            *ptx = va_d.<a class="code" href="classLorene_1_1Valeur.html#a9fb07214796f7f5c29b770c434c884c7">c_cf</a>-&gt;<a class="code" href="classLorene_1_1Mtbl__cf.html#a340f3b625c2961afa5ea7e6319f1d06b">val_point_symy</a>(ld, xxd, tetd, phid) ;</div><div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;</div><div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160;            <span class="comment">// Next point : </span></div><div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;            ptx++ ;   </div><div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;            pr_a++ ;  </div><div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;            ptet_a++ ;  </div><div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;            pphi_a++ ; </div><div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;            px_a++ ;  </div><div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;            py_a++ ;  </div><div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;            pz_a++ ;  </div><div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;</div><div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;        }</div><div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160;        }</div><div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;    }</div><div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;    </div><div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;    <span class="comment">// The remaining points are obtained by symmetry with rspect to the</span></div><div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;    <span class="comment">//  y=0 plane</span></div><div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;        </div><div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=np/2+1 ; k&lt;np ; k++) {</div><div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;</div><div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;        <span class="comment">// pointer on the value (already computed) at the point symmetric </span></div><div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;        <span class="comment">// with respect to the plane y=0</span></div><div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;        <span class="keywordtype">double</span>* ptx_symy = (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;t + (np-k)*nt*nr ;  </div><div class="line"><a name="l00315"></a><span class="lineno">  315</span>&#160;</div><div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;        <span class="comment">// copy : </span></div><div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0 ; j&lt;nt ; j++) {</div><div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0 ; i&lt;nr ; i++) {</div><div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;            *ptx = *ptx_symy ; </div><div class="line"><a name="l00320"></a><span class="lineno">  320</span>&#160;            ptx++ ;</div><div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;            ptx_symy++ ; </div><div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160;        }</div><div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;        }</div><div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;    }</div><div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;    </div><div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;        </div><div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;    }   <span class="comment">// End of the loop on the Arrival domains</span></div><div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;    </div><div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;    <span class="comment">// In the remaining domains, *this is set to zero:</span></div><div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;    <span class="comment">// ----------------------------------------------</span></div><div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;    </div><div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;    <span class="keywordflow">if</span> (nzet &lt; nz_a) {</div><div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a7d65bc7f080d1974a8429de810e0eff3">annule</a>(nzet, nz_a - 1) ; </div><div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;    }</div><div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;    </div><div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;    <span class="comment">// Treatment of dzpuis</span></div><div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;    <span class="comment">// -------------------</span></div><div class="line"><a name="l00338"></a><span class="lineno">  338</span>&#160;    </div><div class="line"><a name="l00339"></a><span class="lineno">  339</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a9c8f9055b9aa09e2f8c860e59c995f31">set_dzpuis</a>(0) ; </div><div class="line"><a name="l00340"></a><span class="lineno">  340</span>&#160;</div><div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;}</div><div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;</div><div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;</div><div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;        <span class="comment">//-------------------------------------//</span></div><div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;        <span class="comment">//   Case of aligned Cartesian axis    //</span></div><div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;        <span class="comment">//-------------------------------------//</span></div><div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;</div><div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;</div><div class="line"><a name="l00349"></a><span class="lineno"><a class="line" href="classLorene_1_1Scalar.html#aa79d50444dc989294d4407354cc8c66b">  349</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="classLorene_1_1Scalar.html#aa79d50444dc989294d4407354cc8c66b">Scalar::import_align_symy</a>(<span class="keywordtype">int</span> nzet, <span class="keyword">const</span> <a class="code" href="classLorene_1_1Scalar.html">Scalar</a>&amp; cm_d) {</div><div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;    </div><div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;    <span class="comment">// Trivial case : null Scalar</span></div><div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;    <span class="comment">// ------------------------</span></div><div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;    </div><div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATZERO) {</div><div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a68abd6fe37c8829980599db324fcd725">set_etat_zero</a>() ; </div><div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;    return ; </div><div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;    }</div><div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATUN) {</div><div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a089209b1af70621abbf5ce370ecde0b7">set_etat_one</a>() ; </div><div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;    return ; </div><div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;    }</div><div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;</div><div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Map.html">Map</a>* mp_d = &amp;(cm_d.<a class="code" href="classLorene_1_1Tensor.html#a00ed03bf0b7364e5bec0cda699b07f26">get_mp</a>()) ; <span class="comment">// Departure mapping</span></div><div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;</div><div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;    <span class="comment">// Protections</span></div><div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;    <span class="comment">// -----------</span></div><div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;    <span class="keywordtype">int</span> align = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align() ;</div><div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;</div><div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;    assert( align  * (mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">get_bvect_cart</a>()).get_align() == 1 ) ; </div><div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;    </div><div class="line"><a name="l00371"></a><span class="lineno">  371</span>&#160;    assert(cm_d.<a class="code" href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">get_etat</a>() == ETATQCQ) ;</div><div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;</div><div class="line"><a name="l00373"></a><span class="lineno">  373</span>&#160;    <span class="keywordflow">if</span> (cm_d.<a class="code" href="classLorene_1_1Scalar.html#a309bc82975514dbc9df435f7624e25fd">get_dzpuis</a>() != 0) {</div><div class="line"><a name="l00374"></a><span class="lineno">  374</span>&#160;    cout &lt;&lt; </div><div class="line"><a name="l00375"></a><span class="lineno">  375</span>&#160;    <span class="stringliteral">&quot;Scalar::import_align_symy : the dzpuis of the Scalar to be imported&quot;</span></div><div class="line"><a name="l00376"></a><span class="lineno">  376</span>&#160;    &lt;&lt; <span class="stringliteral">&quot; must be zero !&quot;</span> &lt;&lt; endl ; </div><div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160;    abort() ; </div><div class="line"><a name="l00378"></a><span class="lineno">  378</span>&#160;    }</div><div class="line"><a name="l00379"></a><span class="lineno">  379</span>&#160;</div><div class="line"><a name="l00380"></a><span class="lineno">  380</span>&#160;</div><div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mg3d.html">Mg3d</a>* mg_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4df8c6a943a85ef314ea636240ef31a7">get_mg</a>() ; </div><div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;    assert(mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a6fc60f226515b2c3ef3959ed9db0e1ce">get_type_p</a>() == NONSYM) ; </div><div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;</div><div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;    <span class="keywordtype">int</span> nz_a = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a99fb0e7a3f1f0f3e405624855200055e">get_nzone</a>() ; </div><div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;    assert(nzet &lt;= nz_a) ;     </div><div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;</div><div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Valeur.html">Valeur</a>&amp; va_d = cm_d.<a class="code" href="classLorene_1_1Scalar.html#aa8b9b16249b758df4fe8e1ad8689b25e">get_spectral_va</a>() ; </div><div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;    va_d.<a class="code" href="classLorene_1_1Valeur.html#aae49809b7bbd41b848e873ec076c025a">coef</a>() ;       <span class="comment">// The coefficients are required</span></div><div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;    </div><div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;</div><div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;    <span class="comment">// Preparations for storing the result in *this </span></div><div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160;    <span class="comment">// --------------------------------------------</span></div><div class="line"><a name="l00393"></a><span class="lineno">  393</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a908fc51b1306d30ba5525198e803f246">del_t</a>() ;   <span class="comment">// delete all previously computed derived quantities</span></div><div class="line"><a name="l00394"></a><span class="lineno">  394</span>&#160;    </div><div class="line"><a name="l00395"></a><span class="lineno">  395</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a39f70d51538ce68155c311f2eebeb82a">set_etat_qcq</a>() ;        <span class="comment">// Set the state to ETATQCQ</span></div><div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;    </div><div class="line"><a name="l00397"></a><span class="lineno">  397</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a1e52288a6fd8f61dcf2a6c202756094e">set_etat_c_qcq</a>() ;   <span class="comment">// Allocates the memory for the Mtbl va.c</span></div><div class="line"><a name="l00398"></a><span class="lineno">  398</span>&#160;                <span class="comment">//  if it does not exist already</span></div><div class="line"><a name="l00399"></a><span class="lineno">  399</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#ab765d0b25416a6da5add6d2ea6b44acc">set_etat_qcq</a>() ;  <span class="comment">// Allocates the memory for the Tbl&#39;s in each</span></div><div class="line"><a name="l00400"></a><span class="lineno">  400</span>&#160;                <span class="comment">//  domain if they do not exist already</span></div><div class="line"><a name="l00401"></a><span class="lineno">  401</span>&#160;</div><div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160;</div><div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;    <span class="comment">// Departure (x,y,z) coordinates of the origin of the Arrival mapping :</span></div><div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;</div><div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;    <span class="keywordtype">double</span> xx_a, yy_a, zz_a ; </div><div class="line"><a name="l00406"></a><span class="lineno">  406</span>&#160;    <span class="keywordflow">if</span> (align == 1) {</div><div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;    xx_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() ; </div><div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;    yy_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() ; </div><div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;    }</div><div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;    <span class="keywordflow">else</span> {</div><div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;    xx_a = mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() - <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">get_ori_x</a>() ; </div><div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;    yy_a = mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() - <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">get_ori_y</a>() ; </div><div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;    }</div><div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;    zz_a = <a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a0159bdb0fe9f7fa3a78690ecce1658f1">get_ori_z</a>() - mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a0159bdb0fe9f7fa3a78690ecce1658f1">get_ori_z</a>() ; </div><div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;</div><div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;    </div><div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;    <span class="comment">// r, theta, phi, x, y and z on the Arrival mapping </span></div><div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;    <span class="comment">//  update of the corresponding Coord&#39;s if necessary</span></div><div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;    </div><div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).fait() ; </div><div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).fait() ; </div><div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).fait() ; </div><div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).fait() ; </div><div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).fait() ; </div><div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;    <span class="keywordflow">if</span> ( (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).c == 0x0 ) (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).fait() ; </div><div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;</div><div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mr_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">r</a>).c ; </div><div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mtet_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">tet</a>).c ; </div><div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mphi_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">phi</a>).c ; </div><div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mx_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">x</a>).c ;</div><div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* my_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">y</a>).c ;</div><div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;    <span class="keyword">const</span> <a class="code" href="classLorene_1_1Mtbl.html">Mtbl</a>* mz_a = (<a class="code" href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">mp</a>-&gt;<a class="code" href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">z</a>).c ;</div><div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;</div><div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;    <a class="code" href="classLorene_1_1Param.html">Param</a> par_precis ;      <span class="comment">// Required precision in the method Map::val_lx</span></div><div class="line"><a name="l00435"></a><span class="lineno">  435</span>&#160;    <span class="keywordtype">int</span> nitermax = 100 ;    <span class="comment">// Maximum number of iteration in the secant method</span></div><div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;    <span class="keywordtype">int</span> niter ; </div><div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160;    <span class="keywordtype">double</span> precis = 1e-15 ; <span class="comment">// Absolute precision in the secant method</span></div><div class="line"><a name="l00438"></a><span class="lineno">  438</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#aefd8de56123122b81d919e877a7a5cb7">add_int</a>(nitermax) ;  </div><div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#a8ea6e8bc0a7166c622764641bbc3744e">add_int_mod</a>(niter) ; </div><div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;    par_precis.<a class="code" href="classLorene_1_1Param.html#a049f47c98bc2f7162a94b5a8c9215fde">add_double</a>(precis) ; </div><div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;</div><div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;</div><div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;    <span class="comment">// Loop of the Arrival domains where the computation is to be performed</span></div><div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160;    <span class="comment">// --------------------------------------------------------------------</span></div><div class="line"><a name="l00445"></a><span class="lineno">  445</span>&#160;    </div><div class="line"><a name="l00446"></a><span class="lineno">  446</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> l=0; l &lt; nzet; l++) {</div><div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;    </div><div class="line"><a name="l00448"></a><span class="lineno">  448</span>&#160;    <span class="keywordtype">int</span> nr = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#aad4b3f34fc58cc34de13540f3d37699a">get_nr</a>(l) ;</div><div class="line"><a name="l00449"></a><span class="lineno">  449</span>&#160;    <span class="keywordtype">int</span> nt = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a7b8408df0cb45eb15d7ba838ee82e2ac">get_nt</a>(l) ;</div><div class="line"><a name="l00450"></a><span class="lineno">  450</span>&#160;    <span class="keywordtype">int</span> np = mg_a-&gt;<a class="code" href="classLorene_1_1Mg3d.html#a6af09c50edf1dacdf9f9de24478746a4">get_np</a>(l) ;</div><div class="line"><a name="l00451"></a><span class="lineno">  451</span>&#160;    </div><div class="line"><a name="l00452"></a><span class="lineno">  452</span>&#160;    </div><div class="line"><a name="l00453"></a><span class="lineno">  453</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pr_a = mr_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of r</span></div><div class="line"><a name="l00454"></a><span class="lineno">  454</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* ptet_a = mtet_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;  <span class="comment">// Pointer on the values of theta</span></div><div class="line"><a name="l00455"></a><span class="lineno">  455</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pphi_a = mphi_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;  <span class="comment">// Pointer on the values of phi</span></div><div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* px_a = mx_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of X</span></div><div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* py_a = my_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of Y</span></div><div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;    <span class="keyword">const</span> <span class="keywordtype">double</span>* pz_a = mz_a-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l]-&gt;<a class="code" href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">t</a> ;      <span class="comment">// Pointer on the values of Z</span></div><div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;    </div><div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;    (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;<a class="code" href="classLorene_1_1Scalar.html#a39f70d51538ce68155c311f2eebeb82a">set_etat_qcq</a>() ;       <span class="comment">// Allocates the array of double to</span></div><div class="line"><a name="l00461"></a><span class="lineno">  461</span>&#160;                         <span class="comment">//  store the result </span></div><div class="line"><a name="l00462"></a><span class="lineno">  462</span>&#160;    <span class="keywordtype">double</span>* ptx = (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;t ;      <span class="comment">// Pointer on the allocated array</span></div><div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;    </div><div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;    </div><div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160;    </div><div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;    <span class="comment">// Loop on half the grid points in the considered arrival domain</span></div><div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;    <span class="comment">// (the other half will be obtained by symmetry with respect to</span></div><div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160;    <span class="comment">//  the y=0 plane). </span></div><div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;        </div><div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0 ; k&lt;np/2+1 ; k++) {    <span class="comment">// np/2+1 : half the grid</span></div><div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0 ; j&lt;nt ; j++) {</div><div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0 ; i&lt;nr ; i++) {</div><div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;</div><div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;            <span class="keywordtype">double</span> r = *pr_a ; </div><div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;            <span class="keywordtype">double</span> rd, tetd, phid ;</div><div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;            <span class="keywordflow">if</span> (r == __infinity) {</div><div class="line"><a name="l00477"></a><span class="lineno">  477</span>&#160;            rd = r ;</div><div class="line"><a name="l00478"></a><span class="lineno">  478</span>&#160;            tetd = *ptet_a ;</div><div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;            phid = *pphi_a ;</div><div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;            }</div><div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;            <span class="keywordflow">else</span> {</div><div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;</div><div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;            <span class="comment">// Cartesian coordinates on the Departure mapping</span></div><div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;            <span class="keywordtype">double</span> xd = *px_a + xx_a ; </div><div class="line"><a name="l00485"></a><span class="lineno">  485</span>&#160;            <span class="keywordtype">double</span> yd = *py_a + yy_a ; </div><div class="line"><a name="l00486"></a><span class="lineno">  486</span>&#160;            <span class="keywordtype">double</span> zd = *pz_a + zz_a ; </div><div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;</div><div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;            <span class="comment">// Spherical coordinates on the Departure mapping</span></div><div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;            <span class="keywordtype">double</span> rhod2 = xd*xd + yd*yd ; </div><div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160;            <span class="keywordtype">double</span> rhod = <a class="code" href="classLorene_1_1Scalar.html#aae28013acc7f65d03bf107c0ded5f73e">sqrt</a>( rhod2 ) ;</div><div class="line"><a name="l00491"></a><span class="lineno">  491</span>&#160;            rd = <a class="code" href="classLorene_1_1Scalar.html#aae28013acc7f65d03bf107c0ded5f73e">sqrt</a>(rhod2 + zd*zd) ;</div><div class="line"><a name="l00492"></a><span class="lineno">  492</span>&#160;            tetd = atan2(rhod, zd) ;</div><div class="line"><a name="l00493"></a><span class="lineno">  493</span>&#160;            phid = atan2(yd, xd) ; </div><div class="line"><a name="l00494"></a><span class="lineno">  494</span>&#160;            <span class="keywordflow">if</span> (phid &lt; 0) phid += 2*M_PI ;</div><div class="line"><a name="l00495"></a><span class="lineno">  495</span>&#160;            }           </div><div class="line"><a name="l00496"></a><span class="lineno">  496</span>&#160;</div><div class="line"><a name="l00497"></a><span class="lineno">  497</span>&#160;</div><div class="line"><a name="l00498"></a><span class="lineno">  498</span>&#160;            <span class="comment">// NB: to increase the efficiency, the method Scalar::val_point</span></div><div class="line"><a name="l00499"></a><span class="lineno">  499</span>&#160;            <span class="comment">//  is not invoked; the method Mtbl_cf::val_point is </span></div><div class="line"><a name="l00500"></a><span class="lineno">  500</span>&#160;            <span class="comment">//  called directly instead. </span></div><div class="line"><a name="l00501"></a><span class="lineno">  501</span>&#160;</div><div class="line"><a name="l00502"></a><span class="lineno">  502</span>&#160;            <span class="comment">// Value of the grid coordinates (l,xi) corresponding to</span></div><div class="line"><a name="l00503"></a><span class="lineno">  503</span>&#160;            <span class="comment">//  (rd,tetd,phid) :</span></div><div class="line"><a name="l00504"></a><span class="lineno">  504</span>&#160;            </div><div class="line"><a name="l00505"></a><span class="lineno">  505</span>&#160;            <span class="keywordtype">int</span> ld ;        <span class="comment">// domain index</span></div><div class="line"><a name="l00506"></a><span class="lineno">  506</span>&#160;            <span class="keywordtype">double</span> xxd ;    <span class="comment">// radial coordinate xi in [0,1] or [-1,1]</span></div><div class="line"><a name="l00507"></a><span class="lineno">  507</span>&#160;            mp_d-&gt;<a class="code" href="classLorene_1_1Map.html#a835fb3c3bfefe88135160b74ef7010b4">val_lx</a>(rd, tetd, phid, par_precis, ld, xxd) ;</div><div class="line"><a name="l00508"></a><span class="lineno">  508</span>&#160;</div><div class="line"><a name="l00509"></a><span class="lineno">  509</span>&#160;            <span class="comment">// Value of the Departure Scalar at the obtained point:</span></div><div class="line"><a name="l00510"></a><span class="lineno">  510</span>&#160;            *ptx = va_d.<a class="code" href="classLorene_1_1Valeur.html#a9fb07214796f7f5c29b770c434c884c7">c_cf</a>-&gt;<a class="code" href="classLorene_1_1Mtbl__cf.html#a340f3b625c2961afa5ea7e6319f1d06b">val_point_symy</a>(ld, xxd, tetd, phid) ;</div><div class="line"><a name="l00511"></a><span class="lineno">  511</span>&#160;</div><div class="line"><a name="l00512"></a><span class="lineno">  512</span>&#160;            <span class="comment">// Next point : </span></div><div class="line"><a name="l00513"></a><span class="lineno">  513</span>&#160;            ptx++ ;   </div><div class="line"><a name="l00514"></a><span class="lineno">  514</span>&#160;            pr_a++ ;  </div><div class="line"><a name="l00515"></a><span class="lineno">  515</span>&#160;            ptet_a++ ;  </div><div class="line"><a name="l00516"></a><span class="lineno">  516</span>&#160;            pphi_a++ ; </div><div class="line"><a name="l00517"></a><span class="lineno">  517</span>&#160;            px_a++ ;  </div><div class="line"><a name="l00518"></a><span class="lineno">  518</span>&#160;            py_a++ ;  </div><div class="line"><a name="l00519"></a><span class="lineno">  519</span>&#160;            pz_a++ ;  </div><div class="line"><a name="l00520"></a><span class="lineno">  520</span>&#160;</div><div class="line"><a name="l00521"></a><span class="lineno">  521</span>&#160;        }</div><div class="line"><a name="l00522"></a><span class="lineno">  522</span>&#160;        }</div><div class="line"><a name="l00523"></a><span class="lineno">  523</span>&#160;    }   </div><div class="line"><a name="l00524"></a><span class="lineno">  524</span>&#160;    </div><div class="line"><a name="l00525"></a><span class="lineno">  525</span>&#160;        </div><div class="line"><a name="l00526"></a><span class="lineno">  526</span>&#160;    <span class="comment">// The remaining points are obtained by symmetry with rspect to the</span></div><div class="line"><a name="l00527"></a><span class="lineno">  527</span>&#160;    <span class="comment">//  y=0 plane</span></div><div class="line"><a name="l00528"></a><span class="lineno">  528</span>&#160;        </div><div class="line"><a name="l00529"></a><span class="lineno">  529</span>&#160;    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=np/2+1 ; k&lt;np ; k++) {</div><div class="line"><a name="l00530"></a><span class="lineno">  530</span>&#160;</div><div class="line"><a name="l00531"></a><span class="lineno">  531</span>&#160;        <span class="comment">// pointer on the value (already computed) at the point symmetric </span></div><div class="line"><a name="l00532"></a><span class="lineno">  532</span>&#160;        <span class="comment">// with respect to the plane y=0</span></div><div class="line"><a name="l00533"></a><span class="lineno">  533</span>&#160;        <span class="keywordtype">double</span>* ptx_symy = (<a class="code" href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">va</a>.<a class="code" href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">c</a>-&gt;<a class="code" href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">t</a>[l])-&gt;t + (np-k)*nt*nr ;  </div><div class="line"><a name="l00534"></a><span class="lineno">  534</span>&#160;</div><div class="line"><a name="l00535"></a><span class="lineno">  535</span>&#160;        <span class="comment">// copy : </span></div><div class="line"><a name="l00536"></a><span class="lineno">  536</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0 ; j&lt;nt ; j++) {</div><div class="line"><a name="l00537"></a><span class="lineno">  537</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0 ; i&lt;nr ; i++) {</div><div class="line"><a name="l00538"></a><span class="lineno">  538</span>&#160;            *ptx = *ptx_symy ; </div><div class="line"><a name="l00539"></a><span class="lineno">  539</span>&#160;            ptx++ ;</div><div class="line"><a name="l00540"></a><span class="lineno">  540</span>&#160;            ptx_symy++ ; </div><div class="line"><a name="l00541"></a><span class="lineno">  541</span>&#160;        }</div><div class="line"><a name="l00542"></a><span class="lineno">  542</span>&#160;        }</div><div class="line"><a name="l00543"></a><span class="lineno">  543</span>&#160;    } </div><div class="line"><a name="l00544"></a><span class="lineno">  544</span>&#160;</div><div class="line"><a name="l00545"></a><span class="lineno">  545</span>&#160;    }   <span class="comment">// End of the loop on the Arrival domains</span></div><div class="line"><a name="l00546"></a><span class="lineno">  546</span>&#160;    </div><div class="line"><a name="l00547"></a><span class="lineno">  547</span>&#160;    <span class="comment">// In the remaining domains, *this is set to zero:</span></div><div class="line"><a name="l00548"></a><span class="lineno">  548</span>&#160;    <span class="comment">// ----------------------------------------------</span></div><div class="line"><a name="l00549"></a><span class="lineno">  549</span>&#160;    </div><div class="line"><a name="l00550"></a><span class="lineno">  550</span>&#160;    <span class="keywordflow">if</span> (nzet &lt; nz_a) {</div><div class="line"><a name="l00551"></a><span class="lineno">  551</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a7d65bc7f080d1974a8429de810e0eff3">annule</a>(nzet, nz_a - 1) ; </div><div class="line"><a name="l00552"></a><span class="lineno">  552</span>&#160;    }</div><div class="line"><a name="l00553"></a><span class="lineno">  553</span>&#160;    </div><div class="line"><a name="l00554"></a><span class="lineno">  554</span>&#160;    <span class="comment">// Treatment of dzpuis</span></div><div class="line"><a name="l00555"></a><span class="lineno">  555</span>&#160;    <span class="comment">// -------------------</span></div><div class="line"><a name="l00556"></a><span class="lineno">  556</span>&#160;    </div><div class="line"><a name="l00557"></a><span class="lineno">  557</span>&#160;    <a class="code" href="classLorene_1_1Scalar.html#a9c8f9055b9aa09e2f8c860e59c995f31">set_dzpuis</a>(0) ; </div><div class="line"><a name="l00558"></a><span class="lineno">  558</span>&#160;</div><div class="line"><a name="l00559"></a><span class="lineno">  559</span>&#160;}</div><div class="line"><a name="l00560"></a><span class="lineno">  560</span>&#160;}</div><div class="ttc" id="classLorene_1_1Mg3d_html_a6fc60f226515b2c3ef3959ed9db0e1ce"><div class="ttname"><a href="classLorene_1_1Mg3d.html#a6fc60f226515b2c3ef3959ed9db0e1ce">Lorene::Mg3d::get_type_p</a></div><div class="ttdeci">int get_type_p() const</div><div class="ttdoc">Returns the type of sampling in the  direction:   SYM :  : symmetry with respect to the transformatio...</div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00495">grilles.h:495</a></div></div>
<div class="ttc" id="classLorene_1_1Valeur_html_a9fb07214796f7f5c29b770c434c884c7"><div class="ttname"><a href="classLorene_1_1Valeur.html#a9fb07214796f7f5c29b770c434c884c7">Lorene::Valeur::c_cf</a></div><div class="ttdeci">Mtbl_cf * c_cf</div><div class="ttdoc">Coefficients of the spectral expansion of the function. </div><div class="ttdef"><b>Definition:</b> <a href="valeur_8h_source.html#l00302">valeur.h:302</a></div></div>
<div class="ttc" id="classLorene_1_1Param_html_aefd8de56123122b81d919e877a7a5cb7"><div class="ttname"><a href="classLorene_1_1Param.html#aefd8de56123122b81d919e877a7a5cb7">Lorene::Param::add_int</a></div><div class="ttdeci">void add_int(const int &amp;n, int position=0)</div><div class="ttdoc">Adds the address of a new int to the list. </div><div class="ttdef"><b>Definition:</b> <a href="param_8C_source.html#l00246">param.C:246</a></div></div>
<div class="ttc" id="classLorene_1_1Mg3d_html_a6af09c50edf1dacdf9f9de24478746a4"><div class="ttname"><a href="classLorene_1_1Mg3d.html#a6af09c50edf1dacdf9f9de24478746a4">Lorene::Mg3d::get_np</a></div><div class="ttdeci">int get_np(int l) const</div><div class="ttdoc">Returns the number of points in the azimuthal direction ( ) in domain no. l. </div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00462">grilles.h:462</a></div></div>
<div class="ttc" id="classLorene_1_1Valeur_html_aae49809b7bbd41b848e873ec076c025a"><div class="ttname"><a href="classLorene_1_1Valeur.html#aae49809b7bbd41b848e873ec076c025a">Lorene::Valeur::coef</a></div><div class="ttdeci">void coef() const</div><div class="ttdoc">Computes the coeffcients of *this. </div><div class="ttdef"><b>Definition:</b> <a href="valeur__coef_8C_source.html#l00148">valeur_coef.C:148</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a68abd6fe37c8829980599db324fcd725"><div class="ttname"><a href="classLorene_1_1Scalar.html#a68abd6fe37c8829980599db324fcd725">Lorene::Scalar::set_etat_zero</a></div><div class="ttdeci">virtual void set_etat_zero()</div><div class="ttdoc">Sets the logical state to ETATZERO (zero). </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00324">scalar.C:324</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a7d65bc7f080d1974a8429de810e0eff3"><div class="ttname"><a href="classLorene_1_1Scalar.html#a7d65bc7f080d1974a8429de810e0eff3">Lorene::Scalar::annule</a></div><div class="ttdeci">virtual void annule(int l_min, int l_max)</div><div class="ttdoc">Sets the Scalar to zero in several domains. </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00391">scalar.C:391</a></div></div>
<div class="ttc" id="classLorene_1_1Mtbl_html"><div class="ttname"><a href="classLorene_1_1Mtbl.html">Lorene::Mtbl</a></div><div class="ttdoc">Multi-domain array. </div><div class="ttdef"><b>Definition:</b> <a href="mtbl_8h_source.html#l00118">mtbl.h:118</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a751ddab3cb82b9bc3b9001ab65d7ccb8"><div class="ttname"><a href="classLorene_1_1Map.html#a751ddab3cb82b9bc3b9001ab65d7ccb8">Lorene::Map::get_ori_y</a></div><div class="ttdeci">double get_ori_y() const</div><div class="ttdoc">Returns the y coordinate of the origin. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00770">map.h:770</a></div></div>
<div class="ttc" id="namespaceLorene_html"><div class="ttname"><a href="namespaceLorene.html">Lorene</a></div><div class="ttdoc">Lorene prototypes. </div><div class="ttdef"><b>Definition:</b> <a href="app__hor_8h_source.html#l00064">app_hor.h:64</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a4df8c6a943a85ef314ea636240ef31a7"><div class="ttname"><a href="classLorene_1_1Map.html#a4df8c6a943a85ef314ea636240ef31a7">Lorene::Map::get_mg</a></div><div class="ttdeci">const Mg3d * get_mg() const</div><div class="ttdoc">Gives the Mg3d on which the mapping is defined. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00765">map.h:765</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html"><div class="ttname"><a href="classLorene_1_1Scalar.html">Lorene::Scalar</a></div><div class="ttdoc">Tensor field of valence 0 (or component of a tensorial field). </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8h_source.html#l00387">scalar.h:387</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html"><div class="ttname"><a href="classLorene_1_1Map.html">Lorene::Map</a></div><div class="ttdoc">Base class for coordinate mappings. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00670">map.h:670</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_aaed9fc290da48be795769157a4aacdb8"><div class="ttname"><a href="classLorene_1_1Map.html#aaed9fc290da48be795769157a4aacdb8">Lorene::Map::get_ori_x</a></div><div class="ttdeci">double get_ori_x() const</div><div class="ttdoc">Returns the x coordinate of the origin. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00768">map.h:768</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_ae2425d68e625b5345ab518c058e5417b"><div class="ttname"><a href="classLorene_1_1Scalar.html#ae2425d68e625b5345ab518c058e5417b">Lorene::Scalar::import_anti_symy</a></div><div class="ttdeci">void import_anti_symy(int nzet, const Scalar &amp;ci)</div><div class="ttdoc">Assignment to another Scalar defined on a different mapping, when the two mappings have anti-aligned ...</div><div class="ttdef"><b>Definition:</b> <a href="scalar__import__symy_8C_source.html#l00131">scalar_import_symy.C:131</a></div></div>
<div class="ttc" id="classLorene_1_1Valeur_html"><div class="ttname"><a href="classLorene_1_1Valeur.html">Lorene::Valeur</a></div><div class="ttdoc">Values and coefficients of a (real-value) function. </div><div class="ttdef"><b>Definition:</b> <a href="valeur_8h_source.html#l00287">valeur.h:287</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a957e9163ea5a022bef63bdd295cc473f"><div class="ttname"><a href="classLorene_1_1Scalar.html#a957e9163ea5a022bef63bdd295cc473f">Lorene::Scalar::get_etat</a></div><div class="ttdeci">int get_etat() const</div><div class="ttdoc">Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary). </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8h_source.html#l00554">scalar.h:554</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a39f70d51538ce68155c311f2eebeb82a"><div class="ttname"><a href="classLorene_1_1Scalar.html#a39f70d51538ce68155c311f2eebeb82a">Lorene::Scalar::set_etat_qcq</a></div><div class="ttdeci">virtual void set_etat_qcq()</div><div class="ttdoc">Sets the logical state to ETATQCQ (ordinary state). </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00353">scalar.C:353</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_aae28013acc7f65d03bf107c0ded5f73e"><div class="ttname"><a href="classLorene_1_1Scalar.html#aae28013acc7f65d03bf107c0ded5f73e">Lorene::Scalar::sqrt</a></div><div class="ttdeci">friend Scalar sqrt(const Scalar &amp;)</div><div class="ttdoc">Square root. </div><div class="ttdef"><b>Definition:</b> <a href="scalar__math_8C_source.html#l00263">scalar_math.C:263</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a801847722a3edf596667fd2603b07a29"><div class="ttname"><a href="classLorene_1_1Map.html#a801847722a3edf596667fd2603b07a29">Lorene::Map::tet</a></div><div class="ttdeci">Coord tet</div><div class="ttdoc"> coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00719">map.h:719</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a9c8f9055b9aa09e2f8c860e59c995f31"><div class="ttname"><a href="classLorene_1_1Scalar.html#a9c8f9055b9aa09e2f8c860e59c995f31">Lorene::Scalar::set_dzpuis</a></div><div class="ttdeci">void set_dzpuis(int)</div><div class="ttdoc">Modifies the dzpuis flag. </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00808">scalar.C:808</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_ac0430f88889bb9a4cd3e269a79d7fcdb"><div class="ttname"><a href="classLorene_1_1Map.html#ac0430f88889bb9a4cd3e269a79d7fcdb">Lorene::Map::phi</a></div><div class="ttdeci">Coord phi</div><div class="ttdoc"> coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00720">map.h:720</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a309bc82975514dbc9df435f7624e25fd"><div class="ttname"><a href="classLorene_1_1Scalar.html#a309bc82975514dbc9df435f7624e25fd">Lorene::Scalar::get_dzpuis</a></div><div class="ttdeci">int get_dzpuis() const</div><div class="ttdoc">Returns dzpuis. </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8h_source.html#l00557">scalar.h:557</a></div></div>
<div class="ttc" id="classLorene_1_1Tbl_html_ae35a88351a0826b159a34435c422c69a"><div class="ttname"><a href="classLorene_1_1Tbl.html#ae35a88351a0826b159a34435c422c69a">Lorene::Tbl::t</a></div><div class="ttdeci">double * t</div><div class="ttdoc">The array of double. </div><div class="ttdef"><b>Definition:</b> <a href="tbl_8h_source.html#l00173">tbl.h:173</a></div></div>
<div class="ttc" id="classLorene_1_1Valeur_html_a679447efc0853359043f254ae0a4e913"><div class="ttname"><a href="classLorene_1_1Valeur.html#a679447efc0853359043f254ae0a4e913">Lorene::Valeur::c</a></div><div class="ttdeci">Mtbl * c</div><div class="ttdoc">Values of the function at the points of the multi-grid. </div><div class="ttdef"><b>Definition:</b> <a href="valeur_8h_source.html#l00299">valeur.h:299</a></div></div>
<div class="ttc" id="classLorene_1_1Param_html"><div class="ttname"><a href="classLorene_1_1Param.html">Lorene::Param</a></div><div class="ttdoc">Parameter storage. </div><div class="ttdef"><b>Definition:</b> <a href="param_8h_source.html#l00125">param.h:125</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a089209b1af70621abbf5ce370ecde0b7"><div class="ttname"><a href="classLorene_1_1Scalar.html#a089209b1af70621abbf5ce370ecde0b7">Lorene::Scalar::set_etat_one</a></div><div class="ttdeci">void set_etat_one()</div><div class="ttdoc">Sets the logical state to ETATUN (one). </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00334">scalar.C:334</a></div></div>
<div class="ttc" id="classLorene_1_1Mg3d_html_a99fb0e7a3f1f0f3e405624855200055e"><div class="ttname"><a href="classLorene_1_1Mg3d.html#a99fb0e7a3f1f0f3e405624855200055e">Lorene::Mg3d::get_nzone</a></div><div class="ttdeci">int get_nzone() const</div><div class="ttdoc">Returns the number of domains. </div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00448">grilles.h:448</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_acb638e4a424f5a0ad2a839eff473f600"><div class="ttname"><a href="classLorene_1_1Scalar.html#acb638e4a424f5a0ad2a839eff473f600">Lorene::Scalar::va</a></div><div class="ttdeci">Valeur va</div><div class="ttdoc">The numerical value of the Scalar. </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8h_source.html#l00405">scalar.h:405</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_a908fc51b1306d30ba5525198e803f246"><div class="ttname"><a href="classLorene_1_1Scalar.html#a908fc51b1306d30ba5525198e803f246">Lorene::Scalar::del_t</a></div><div class="ttdeci">void del_t()</div><div class="ttdoc">Logical destructor. </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8C_source.html#l00279">scalar.C:279</a></div></div>
<div class="ttc" id="classLorene_1_1Mtbl_html_ab765d0b25416a6da5add6d2ea6b44acc"><div class="ttname"><a href="classLorene_1_1Mtbl.html#ab765d0b25416a6da5add6d2ea6b44acc">Lorene::Mtbl::set_etat_qcq</a></div><div class="ttdeci">void set_etat_qcq()</div><div class="ttdoc">Sets the logical state to ETATQCQ (ordinary state). </div><div class="ttdef"><b>Definition:</b> <a href="mtbl_8C_source.html#l00299">mtbl.C:299</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_aa79d50444dc989294d4407354cc8c66b"><div class="ttname"><a href="classLorene_1_1Scalar.html#aa79d50444dc989294d4407354cc8c66b">Lorene::Scalar::import_align_symy</a></div><div class="ttdeci">void import_align_symy(int nzet, const Scalar &amp;ci)</div><div class="ttdoc">Assignment to another Scalar defined on a different mapping, when the two mappings have aligned Carte...</div><div class="ttdef"><b>Definition:</b> <a href="scalar__import__symy_8C_source.html#l00349">scalar_import_symy.C:349</a></div></div>
<div class="ttc" id="classLorene_1_1Mg3d_html_aad4b3f34fc58cc34de13540f3d37699a"><div class="ttname"><a href="classLorene_1_1Mg3d.html#aad4b3f34fc58cc34de13540f3d37699a">Lorene::Mg3d::get_nr</a></div><div class="ttdeci">int get_nr(int l) const</div><div class="ttdoc">Returns the number of points in the radial direction ( ) in domain no. l. </div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00452">grilles.h:452</a></div></div>
<div class="ttc" id="classLorene_1_1Mg3d_html"><div class="ttname"><a href="classLorene_1_1Mg3d.html">Lorene::Mg3d</a></div><div class="ttdoc">Multi-domain grid. </div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00273">grilles.h:273</a></div></div>
<div class="ttc" id="classLorene_1_1Mtbl__cf_html_a340f3b625c2961afa5ea7e6319f1d06b"><div class="ttname"><a href="classLorene_1_1Mtbl__cf.html#a340f3b625c2961afa5ea7e6319f1d06b">Lorene::Mtbl_cf::val_point_symy</a></div><div class="ttdeci">double val_point_symy(int l, double x, double theta, double phi) const</div><div class="ttdoc">Computes the value of the field represented by *this at an arbitrary point, by means of the spectral ...</div><div class="ttdef"><b>Definition:</b> <a href="mtbl__cf__vp__symy_8C_source.html#l00069">mtbl_cf_vp_symy.C:69</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_adfc2b4aa1ef4f22a109bce235e7eb6b7"><div class="ttname"><a href="classLorene_1_1Map.html#adfc2b4aa1ef4f22a109bce235e7eb6b7">Lorene::Map::get_bvect_cart</a></div><div class="ttdeci">const Base_vect_cart &amp; get_bvect_cart() const</div><div class="ttdoc">Returns the Cartesian basis  associated with the coordinates (x,y,z) of the mapping, i.e. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00791">map.h:791</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a08ea71624c8555ab1e69adea5945b1db"><div class="ttname"><a href="classLorene_1_1Map.html#a08ea71624c8555ab1e69adea5945b1db">Lorene::Map::y</a></div><div class="ttdeci">Coord y</div><div class="ttdoc">y coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00727">map.h:727</a></div></div>
<div class="ttc" id="classLorene_1_1Param_html_a049f47c98bc2f7162a94b5a8c9215fde"><div class="ttname"><a href="classLorene_1_1Param.html#a049f47c98bc2f7162a94b5a8c9215fde">Lorene::Param::add_double</a></div><div class="ttdeci">void add_double(const double &amp;x, int position=0)</div><div class="ttdoc">Adds the the address of a new double to the list. </div><div class="ttdef"><b>Definition:</b> <a href="param_8C_source.html#l00315">param.C:315</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_abca4ab2ee51a2558ebd81e3c6b5bd04d"><div class="ttname"><a href="classLorene_1_1Map.html#abca4ab2ee51a2558ebd81e3c6b5bd04d">Lorene::Map::x</a></div><div class="ttdeci">Coord x</div><div class="ttdoc">x coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00726">map.h:726</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_af17554dd17d8a77098b4f743ed86e38b"><div class="ttname"><a href="classLorene_1_1Scalar.html#af17554dd17d8a77098b4f743ed86e38b">Lorene::Scalar::import_symy</a></div><div class="ttdeci">void import_symy(const Scalar &amp;ci)</div><div class="ttdoc">Assignment to another Scalar defined on a different mapping. </div><div class="ttdef"><b>Definition:</b> <a href="scalar__import__symy_8C_source.html#l00072">scalar_import_symy.C:72</a></div></div>
<div class="ttc" id="classLorene_1_1Valeur_html_a1e52288a6fd8f61dcf2a6c202756094e"><div class="ttname"><a href="classLorene_1_1Valeur.html#a1e52288a6fd8f61dcf2a6c202756094e">Lorene::Valeur::set_etat_c_qcq</a></div><div class="ttdeci">void set_etat_c_qcq()</div><div class="ttdoc">Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c )...</div><div class="ttdef"><b>Definition:</b> <a href="valeur_8C_source.html#l00701">valeur.C:701</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a0159bdb0fe9f7fa3a78690ecce1658f1"><div class="ttname"><a href="classLorene_1_1Map.html#a0159bdb0fe9f7fa3a78690ecce1658f1">Lorene::Map::get_ori_z</a></div><div class="ttdeci">double get_ori_z() const</div><div class="ttdoc">Returns the z coordinate of the origin. </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00772">map.h:772</a></div></div>
<div class="ttc" id="classLorene_1_1Mg3d_html_a7b8408df0cb45eb15d7ba838ee82e2ac"><div class="ttname"><a href="classLorene_1_1Mg3d.html#a7b8408df0cb45eb15d7ba838ee82e2ac">Lorene::Mg3d::get_nt</a></div><div class="ttdeci">int get_nt(int l) const</div><div class="ttdoc">Returns the number of points in the co-latitude direction ( ) in domain no. l. </div><div class="ttdef"><b>Definition:</b> <a href="grilles_8h_source.html#l00457">grilles.h:457</a></div></div>
<div class="ttc" id="classLorene_1_1Tensor_html_a1bb2cef406860a6ea7696d5101474ab1"><div class="ttname"><a href="classLorene_1_1Tensor.html#a1bb2cef406860a6ea7696d5101474ab1">Lorene::Tensor::mp</a></div><div class="ttdeci">const Map *const mp</div><div class="ttdoc">Mapping on which the numerical values at the grid points are defined. </div><div class="ttdef"><b>Definition:</b> <a href="tensor_8h_source.html#l00295">tensor.h:295</a></div></div>
<div class="ttc" id="classLorene_1_1Mtbl_html_aa294df3c399cb928068dea52a4632b9b"><div class="ttname"><a href="classLorene_1_1Mtbl.html#aa294df3c399cb928068dea52a4632b9b">Lorene::Mtbl::t</a></div><div class="ttdeci">Tbl ** t</div><div class="ttdoc">Array (size nzone ) of pointers on the Tbl &amp;#39;s. </div><div class="ttdef"><b>Definition:</b> <a href="mtbl_8h_source.html#l00132">mtbl.h:132</a></div></div>
<div class="ttc" id="classLorene_1_1Tensor_html_a00ed03bf0b7364e5bec0cda699b07f26"><div class="ttname"><a href="classLorene_1_1Tensor.html#a00ed03bf0b7364e5bec0cda699b07f26">Lorene::Tensor::get_mp</a></div><div class="ttdeci">const Map &amp; get_mp() const</div><div class="ttdoc">Returns the mapping. </div><div class="ttdef"><b>Definition:</b> <a href="tensor_8h_source.html#l00861">tensor.h:861</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_aaea595c163565f5878addd2959e3e901"><div class="ttname"><a href="classLorene_1_1Map.html#aaea595c163565f5878addd2959e3e901">Lorene::Map::z</a></div><div class="ttdeci">Coord z</div><div class="ttdoc">z coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00728">map.h:728</a></div></div>
<div class="ttc" id="classLorene_1_1Param_html_a8ea6e8bc0a7166c622764641bbc3744e"><div class="ttname"><a href="classLorene_1_1Param.html#a8ea6e8bc0a7166c622764641bbc3744e">Lorene::Param::add_int_mod</a></div><div class="ttdeci">void add_int_mod(int &amp;n, int position=0)</div><div class="ttdoc">Adds the address of a new modifiable int to the list. </div><div class="ttdef"><b>Definition:</b> <a href="param_8C_source.html#l00385">param.C:385</a></div></div>
<div class="ttc" id="classLorene_1_1Scalar_html_aa8b9b16249b758df4fe8e1ad8689b25e"><div class="ttname"><a href="classLorene_1_1Scalar.html#aa8b9b16249b758df4fe8e1ad8689b25e">Lorene::Scalar::get_spectral_va</a></div><div class="ttdeci">const Valeur &amp; get_spectral_va() const</div><div class="ttdoc">Returns va (read only version) </div><div class="ttdef"><b>Definition:</b> <a href="scalar_8h_source.html#l00601">scalar.h:601</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a4a82b11f6df0f0fa3a6a0f5a98e3b071"><div class="ttname"><a href="classLorene_1_1Map.html#a4a82b11f6df0f0fa3a6a0f5a98e3b071">Lorene::Map::r</a></div><div class="ttdeci">Coord r</div><div class="ttdoc">r coordinate centered on the grid </div><div class="ttdef"><b>Definition:</b> <a href="map_8h_source.html#l00718">map.h:718</a></div></div>
<div class="ttc" id="classLorene_1_1Map_html_a835fb3c3bfefe88135160b74ef7010b4"><div class="ttname"><a href="classLorene_1_1Map.html#a835fb3c3bfefe88135160b74ef7010b4">Lorene::Map::val_lx</a></div><div class="ttdeci">virtual void val_lx(double rr, double theta, double pphi, int &amp;l, double &amp;xi) const =0</div><div class="ttdoc">Computes the domain index l and the value of  corresponding to a point given by its physical coordina...</div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
  <ul>
    <li class="navelem"><a class="el" href="dir_2916f9fbcacdf2e5f0583aa159c707b6.html">C++</a></li><li class="navelem"><a class="el" href="dir_be81ae66a82603ee109e545ee9e7219c.html">Source</a></li><li class="navelem"><a class="el" href="dir_cf6b415a077e2ff00403c4fcfed5bf7f.html">Tensor</a></li><li class="navelem"><a class="el" href="dir_2863452ff65dfc8d646637482ced012e.html">Scalar</a></li><li class="navelem"><b>scalar_import_symy.C</b></li>
    <li class="footer">Generated by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.12 </li>
  </ul>
</div>
</body>
</html>

APT Browse - Built by Thomas Orozco.