-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathsolvers_8h_source.html
More file actions
144 lines (142 loc) · 11 KB
/
solvers_8h_source.html
File metadata and controls
144 lines (142 loc) · 11 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://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.17"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Aether: include/solvers.h 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="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
extensions: ["tex2jax.js"],
jax: ["input/TeX","output/HTML-CSS"],
});
</script>
<script type="text/javascript" async="async" src="https://cdn.jsdelivr.net/npm/mathjax@2/MathJax.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectalign" style="padding-left: 0.5em;">
<div id="projectname">Aether
 <span id="projectnumber">0.0</span>
</div>
<div id="projectbrief">Ionosphere-Thermosphere model</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.17 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "search",false,'Search');
/* @license-end */
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
$(function() {
initMenu('',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
/* @license-end */</script>
<div id="main-nav"></div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_d44c64559bbebec7f509842c48db8b23.html">include</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">solvers.h</div> </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">// Copyright 2020, the Aether Development Team (see doc/dev_team.md for members)</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment">// Full license can be found in License.md</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span>  </div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="preprocessor">#ifndef INCLUDE_SOLVERS_H_</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="preprocessor">#define INCLUDE_SOLVERS_H_</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span>  </div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment">// The armadillo library is to allow the use of 3d cubes and other</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment">// array types, with array math built in. This eliminates loops!</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="preprocessor">#include <armadillo></span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="keyword">using namespace </span>arma;</div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span>  </div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> arma_vec solver_conduction(arma_vec value,</div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span>  arma_vec lambda,</div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span>  arma_vec front,</div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  precision_t dt,</div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  arma_vec dx);</div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span>  </div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  </div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> arma_cube solver_chemistry(arma_cube density,</div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  arma_cube source,</div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  arma_cube loss,</div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  precision_t dt);</div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  </div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="keyword">const</span> <span class="keywordtype">int</span> iPrevious_ = 1;</div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="keyword">const</span> <span class="keywordtype">int</span> iNext_ = 2;</div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="keyword">const</span> <span class="keywordtype">int</span> iClosest_ = 3; </div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="keyword">const</span> <span class="keywordtype">int</span> iInterp_ = 4;</div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  </div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="keywordtype">double</span> interpolate_1d_get_index_doubles(<span class="keywordtype">double</span> intime,</div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span>  std::vector<double> times);</div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  </div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment">// Overloading the interpolation function:</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="keywordtype">double</span> interpolate_1d_w_index(std::vector<double> values,</div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  <span class="keywordtype">double</span> interpolation_index,</div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="keywordtype">int</span> interpolation_type);</div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="keywordtype">double</span> interpolate_1d_w_index(std::vector<float> values,</div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  <span class="keywordtype">double</span> interpolation_index,</div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="keywordtype">int</span> interpolation_type);</div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="keywordtype">double</span> interpolate_1d_w_index(std::vector<float> values,</div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  <span class="keywordtype">float</span> interpolation_index,</div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordtype">int</span> interpolation_type);</div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="keywordtype">double</span> interpolate_1d_w_index(arma_vec values,</div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <span class="keywordtype">double</span> interpolation_index,</div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="keywordtype">int</span> interpolation_type);</div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> fmat interpolate_1d_w_index(std::vector<fmat> values,</div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="keywordtype">double</span> interpolation_index,</div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="keywordtype">int</span> interpolation_type);</div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  </div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span> arma_cube calc_gradient_lon(arma_cube value, <a class="code" href="classGrid.html">Grid</a> grid);</div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> arma_cube calc_gradient_lat(arma_cube value, <a class="code" href="classGrid.html">Grid</a> grid);</div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> arma_cube calc_gradient_alt(arma_cube value, <a class="code" href="classGrid.html">Grid</a> grid);</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> std::vector<arma_cube> calc_gradient_vector(arma_cube value_scgc, <a class="code" href="classGrid.html">Grid</a> grid);</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  </div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="preprocessor">#endif // INCLUDE_SOLVERS_H_</span></div>
</div><!-- fragment --></div><!-- contents -->
<div class="ttc" id="aclassGrid_html"><div class="ttname"><a href="classGrid.html">Grid</a></div><div class="ttdef"><b>Definition:</b> arm_vars.h:11</div></div>
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated by  <a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.17
</small></address>
</body>
</html>