402 lines
70 KiB
HTML
402 lines
70 KiB
HTML
<!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" lang="en-US">
|
||
<head>
|
||
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
|
||
<meta http-equiv="X-UA-Compatible" content="IE=11"/>
|
||
<meta name="generator" content="Doxygen 1.13.2"/>
|
||
<meta name="viewport" content="width=device-width, initial-scale=1"/>
|
||
<title>GridFire: src/lib/solver/solver.cpp Source File</title>
|
||
<link href="tabs.css" rel="stylesheet" type="text/css"/>
|
||
<script type="text/javascript" src="jquery.js"></script>
|
||
<script type="text/javascript" src="dynsections.js"></script>
|
||
<script type="text/javascript" src="clipboard.js"></script>
|
||
<link href="navtree.css" rel="stylesheet" type="text/css"/>
|
||
<script type="text/javascript" src="navtreedata.js"></script>
|
||
<script type="text/javascript" src="navtree.js"></script>
|
||
<script type="text/javascript" src="resize.js"></script>
|
||
<script type="text/javascript" src="cookie.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>
|
||
<link href="doxygen.css" rel="stylesheet" type="text/css" />
|
||
<link href="doxygen-awesome.css" rel="stylesheet" type="text/css"/>
|
||
<link href="doxygen-awesome-sidebar-only.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 id="projectrow">
|
||
<td id="projectalign">
|
||
<div id="projectname">GridFire<span id="projectnumber"> 0.0.1a</span>
|
||
</div>
|
||
<div id="projectbrief">General Purpose Nuclear Network</div>
|
||
</td>
|
||
</tr>
|
||
</tbody>
|
||
</table>
|
||
</div>
|
||
<!-- end header part -->
|
||
<!-- Generated by Doxygen 1.13.2 -->
|
||
<script type="text/javascript">
|
||
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
|
||
var searchBox = new SearchBox("searchBox", "search/",'.html');
|
||
/* @license-end */
|
||
</script>
|
||
<script type="text/javascript">
|
||
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
|
||
$(function() { codefold.init(0); });
|
||
/* @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:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
|
||
$(function() {
|
||
initMenu('',true,false,'search.php','Search',true);
|
||
$(function() { init_search(); });
|
||
});
|
||
/* @license-end */
|
||
</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">
|
||
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
|
||
$(function(){initNavTree('solver_8cpp_source.html',''); initResizable(true); });
|
||
/* @license-end */
|
||
</script>
|
||
<div id="doc-content">
|
||
<!-- window showing the filter options -->
|
||
<div id="MSearchSelectWindow"
|
||
onmouseover="return searchBox.OnSearchSelectShow()"
|
||
onmouseout="return searchBox.OnSearchSelectHide()"
|
||
onkeydown="return searchBox.OnSearchSelectKey(event)">
|
||
</div>
|
||
|
||
<!-- iframe showing the search results (closed by default) -->
|
||
<div id="MSearchResultsWindow">
|
||
<div id="MSearchResults">
|
||
<div class="SRPage">
|
||
<div id="SRIndex">
|
||
<div id="SRResults"></div>
|
||
<div class="SRStatus" id="Loading">Loading...</div>
|
||
<div class="SRStatus" id="Searching">Searching...</div>
|
||
<div class="SRStatus" id="NoMatches">No Matches</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div class="header">
|
||
<div class="headertitle"><div class="title">solver.cpp</div></div>
|
||
</div><!--header-->
|
||
<div class="contents">
|
||
<a href="solver_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a id="l00001" name="l00001"></a><span class="lineno"> 1</span><span class="preprocessor">#include "<a class="code" href="solver_8h.html">gridfire/solver/solver.h</a>"</span></div>
|
||
<div class="line"><a id="l00002" name="l00002"></a><span class="lineno"> 2</span><span class="preprocessor">#include "<a class="code" href="engine__graph_8h.html">gridfire/engine/engine_graph.h</a>"</span></div>
|
||
<div class="line"><a id="l00003" name="l00003"></a><span class="lineno"> 3</span><span class="preprocessor">#include "<a class="code" href="network_8h.html">gridfire/network.h</a>"</span></div>
|
||
<div class="line"><a id="l00004" name="l00004"></a><span class="lineno"> 4</span><span class="preprocessor">#include "<a class="code" href="error__engine_8h.html">gridfire/exceptions/error_engine.h</a>"</span></div>
|
||
<div class="line"><a id="l00005" name="l00005"></a><span class="lineno"> 5</span> </div>
|
||
<div class="line"><a id="l00006" name="l00006"></a><span class="lineno"> 6</span><span class="preprocessor">#include "fourdst/composition/atomicSpecies.h"</span></div>
|
||
<div class="line"><a id="l00007" name="l00007"></a><span class="lineno"> 7</span><span class="preprocessor">#include "fourdst/composition/composition.h"</span></div>
|
||
<div class="line"><a id="l00008" name="l00008"></a><span class="lineno"> 8</span><span class="preprocessor">#include "fourdst/config/config.h"</span></div>
|
||
<div class="line"><a id="l00009" name="l00009"></a><span class="lineno"> 9</span> </div>
|
||
<div class="line"><a id="l00010" name="l00010"></a><span class="lineno"> 10</span><span class="preprocessor">#include "unsupported/Eigen/NonLinearOptimization"</span></div>
|
||
<div class="line"><a id="l00011" name="l00011"></a><span class="lineno"> 11</span> </div>
|
||
<div class="line"><a id="l00012" name="l00012"></a><span class="lineno"> 12</span><span class="preprocessor">#include <boost/numeric/odeint.hpp></span></div>
|
||
<div class="line"><a id="l00013" name="l00013"></a><span class="lineno"> 13</span> </div>
|
||
<div class="line"><a id="l00014" name="l00014"></a><span class="lineno"> 14</span><span class="preprocessor">#include <vector></span></div>
|
||
<div class="line"><a id="l00015" name="l00015"></a><span class="lineno"> 15</span><span class="preprocessor">#include <string></span></div>
|
||
<div class="line"><a id="l00016" name="l00016"></a><span class="lineno"> 16</span><span class="preprocessor">#include <stdexcept></span></div>
|
||
<div class="line"><a id="l00017" name="l00017"></a><span class="lineno"> 17</span><span class="preprocessor">#include <iomanip></span></div>
|
||
<div class="line"><a id="l00018" name="l00018"></a><span class="lineno"> 18</span> </div>
|
||
<div class="line"><a id="l00019" name="l00019"></a><span class="lineno"> 19</span><span class="preprocessor">#include "quill/LogMacros.h"</span></div>
|
||
<div class="line"><a id="l00020" name="l00020"></a><span class="lineno"> 20</span> </div>
|
||
<div class="line"><a id="l00021" name="l00021"></a><span class="lineno"> 21</span><span class="keyword">namespace </span><a class="code hl_namespace" href="namespacegridfire_1_1solver.html">gridfire::solver</a> {</div>
|
||
<div class="foldopen" id="foldopen00022" data-start="{" data-end="}">
|
||
<div class="line"><a id="l00022" name="l00022"></a><span class="lineno"><a class="line" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a0e8a4b8ef656e0b084d11bea982e412a"> 22</a></span> <a class="code hl_struct" href="structgridfire_1_1_net_out.html">NetOut</a> <a class="code hl_function" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a0e8a4b8ef656e0b084d11bea982e412a">DirectNetworkSolver::evaluate</a>(<span class="keyword">const</span> <a class="code hl_struct" href="structgridfire_1_1_net_in.html">NetIn</a> &netIn) {</div>
|
||
<div class="line"><a id="l00023" name="l00023"></a><span class="lineno"> 23</span> <span class="keyword">namespace </span>ublas = boost::numeric::ublas;</div>
|
||
<div class="line"><a id="l00024" name="l00024"></a><span class="lineno"> 24</span> <span class="keyword">namespace </span>odeint = boost::numeric::odeint;</div>
|
||
<div class="line"><a id="l00025" name="l00025"></a><span class="lineno"> 25</span> <span class="keyword">using </span>fourdst::composition::Composition;</div>
|
||
<div class="line"><a id="l00026" name="l00026"></a><span class="lineno"> 26</span> </div>
|
||
<div class="line"><a id="l00027" name="l00027"></a><span class="lineno"> 27</span> </div>
|
||
<div class="line"><a id="l00028" name="l00028"></a><span class="lineno"> 28</span> <span class="keyword">const</span> <span class="keywordtype">double</span> T9 = netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a5be0f5195a5cd1dd177b9fc5ab83a7be">temperature</a> / 1e9; <span class="comment">// Convert temperature from Kelvin to T9 (T9 = T / 1e9)</span></div>
|
||
<div class="line"><a id="l00029" name="l00029"></a><span class="lineno"> 29</span> </div>
|
||
<div class="line"><a id="l00030" name="l00030"></a><span class="lineno"> 30</span> <span class="keyword">const</span> <span class="keyword">auto</span> absTol = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a2cc12e737a753a42b72a45be3fbfa8ab">m_config</a>.get<<span class="keywordtype">double</span>>(<span class="stringliteral">"gridfire:solver:DirectNetworkSolver:absTol"</span>, 1.0e-8);</div>
|
||
<div class="line"><a id="l00031" name="l00031"></a><span class="lineno"> 31</span> <span class="keyword">const</span> <span class="keyword">auto</span> relTol = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a2cc12e737a753a42b72a45be3fbfa8ab">m_config</a>.get<<span class="keywordtype">double</span>>(<span class="stringliteral">"gridfire:solver:DirectNetworkSolver:relTol"</span>, 1.0e-8);</div>
|
||
<div class="line"><a id="l00032" name="l00032"></a><span class="lineno"> 32</span> </div>
|
||
<div class="line"><a id="l00033" name="l00033"></a><span class="lineno"> 33</span> Composition equilibratedComposition = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.update(netIn);</div>
|
||
<div class="line"><a id="l00034" name="l00034"></a><span class="lineno"> 34</span> <span class="keywordtype">size_t</span> numSpecies = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies().size();</div>
|
||
<div class="line"><a id="l00035" name="l00035"></a><span class="lineno"> 35</span> ublas::vector<double> Y(numSpecies + 1);</div>
|
||
<div class="line"><a id="l00036" name="l00036"></a><span class="lineno"> 36</span> </div>
|
||
<div class="line"><a id="l00037" name="l00037"></a><span class="lineno"> 37</span> <a class="code hl_struct" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html">RHSManager</a> manager(<a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>, T9, netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a06f0dff9f8927b7cf2da3004c8fa1577">density</a>);</div>
|
||
<div class="line"><a id="l00038" name="l00038"></a><span class="lineno"> 38</span> <a class="code hl_struct" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html">JacobianFunctor</a> jacobianFunctor(<a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>, T9, netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a06f0dff9f8927b7cf2da3004c8fa1577">density</a>);</div>
|
||
<div class="line"><a id="l00039" name="l00039"></a><span class="lineno"> 39</span> </div>
|
||
<div class="line"><a id="l00040" name="l00040"></a><span class="lineno"> 40</span> <span class="keyword">auto</span> populateY = [&](<span class="keyword">const</span> Composition& comp) {</div>
|
||
<div class="line"><a id="l00041" name="l00041"></a><span class="lineno"> 41</span> <span class="keyword">const</span> <span class="keywordtype">size_t</span> numSpeciesInternal = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies().size();</div>
|
||
<div class="line"><a id="l00042" name="l00042"></a><span class="lineno"> 42</span> Y.resize(numSpeciesInternal + 1);</div>
|
||
<div class="line"><a id="l00043" name="l00043"></a><span class="lineno"> 43</span> <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i < numSpeciesInternal; i++) {</div>
|
||
<div class="line"><a id="l00044" name="l00044"></a><span class="lineno"> 44</span> <span class="keyword">const</span> <span class="keyword">auto</span>& species = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies()[i];</div>
|
||
<div class="line"><a id="l00045" name="l00045"></a><span class="lineno"> 45</span> <span class="keywordflow">if</span> (!comp.contains(species)) {</div>
|
||
<div class="line"><a id="l00046" name="l00046"></a><span class="lineno"> 46</span> <span class="keywordtype">double</span> lim = std::numeric_limits<double>::min();</div>
|
||
<div class="line"><a id="l00047" name="l00047"></a><span class="lineno"> 47</span> LOG_DEBUG(<a class="code hl_variable" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a093aa89fd23c2fe03266e286871c7079">m_logger</a>, <span class="stringliteral">"Species '{}' not found in composition. Setting abundance to {:0.3E}."</span>, species.name(), lim);</div>
|
||
<div class="line"><a id="l00048" name="l00048"></a><span class="lineno"> 48</span> Y(i) = lim; <span class="comment">// Species not in the composition, set to zero</span></div>
|
||
<div class="line"><a id="l00049" name="l00049"></a><span class="lineno"> 49</span> } <span class="keywordflow">else</span> {</div>
|
||
<div class="line"><a id="l00050" name="l00050"></a><span class="lineno"> 50</span> Y(i) = comp.getMolarAbundance(species);</div>
|
||
<div class="line"><a id="l00051" name="l00051"></a><span class="lineno"> 51</span> }</div>
|
||
<div class="line"><a id="l00052" name="l00052"></a><span class="lineno"> 52</span> }</div>
|
||
<div class="line"><a id="l00053" name="l00053"></a><span class="lineno"> 53</span> </div>
|
||
<div class="line"><a id="l00054" name="l00054"></a><span class="lineno"> 54</span> <span class="comment">// TODO: a good starting point to make the temperature, density, and energy self consistent would be to turn this into an accumulator</span></div>
|
||
<div class="line"><a id="l00055" name="l00055"></a><span class="lineno"> 55</span> Y(numSpeciesInternal) = 0.0; <span class="comment">// Specific energy rate, initialized to zero</span></div>
|
||
<div class="line"><a id="l00056" name="l00056"></a><span class="lineno"> 56</span> };</div>
|
||
<div class="line"><a id="l00057" name="l00057"></a><span class="lineno"> 57</span> </div>
|
||
<div class="line"><a id="l00058" name="l00058"></a><span class="lineno"> 58</span> <span class="comment">// This is a quick debug that can be turned on. For solar code input parameters (T~1.5e7K, ρ~1.5e3 g/cm^3) this should be near 8e-17</span></div>
|
||
<div class="line"><a id="l00059" name="l00059"></a><span class="lineno"> 59</span> <span class="comment">// std::cout << "D/H: " << equilibratedComposition.getMolarAbundance("H-2") / equilibratedComposition.getMolarAbundance("H-1") << std::endl;</span></div>
|
||
<div class="line"><a id="l00060" name="l00060"></a><span class="lineno"> 60</span> </div>
|
||
<div class="line"><a id="l00061" name="l00061"></a><span class="lineno"> 61</span> populateY(equilibratedComposition);</div>
|
||
<div class="line"><a id="l00062" name="l00062"></a><span class="lineno"> 62</span> <span class="keyword">const</span> <span class="keyword">auto</span> stepper = odeint::make_controlled<odeint::rosenbrock4<double>>(absTol, relTol);</div>
|
||
<div class="line"><a id="l00063" name="l00063"></a><span class="lineno"> 63</span> </div>
|
||
<div class="line"><a id="l00064" name="l00064"></a><span class="lineno"> 64</span> <span class="keywordtype">double</span> current_time = 0.0;</div>
|
||
<div class="line"><a id="l00065" name="l00065"></a><span class="lineno"> 65</span> <span class="keywordtype">double</span> current_initial_timestep = netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a4e556f7bb18f46654b3445476734076a">dt0</a>;</div>
|
||
<div class="line"><a id="l00066" name="l00066"></a><span class="lineno"> 66</span> <span class="keywordtype">double</span> accumulated_energy = 0.0;</div>
|
||
<div class="line"><a id="l00067" name="l00067"></a><span class="lineno"> 67</span> <span class="comment">// size_t total_update_stages_triggered = 0;</span></div>
|
||
<div class="line"><a id="l00068" name="l00068"></a><span class="lineno"> 68</span> </div>
|
||
<div class="line"><a id="l00069" name="l00069"></a><span class="lineno"> 69</span> <span class="keywordflow">while</span> (current_time < netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a0a8d820cfeaa92ee31f253795c57e0d1">tMax</a>) {</div>
|
||
<div class="line"><a id="l00070" name="l00070"></a><span class="lineno"> 70</span> <span class="keywordflow">try</span> {</div>
|
||
<div class="line"><a id="l00071" name="l00071"></a><span class="lineno"> 71</span> odeint::integrate_adaptive(</div>
|
||
<div class="line"><a id="l00072" name="l00072"></a><span class="lineno"> 72</span> stepper,</div>
|
||
<div class="line"><a id="l00073" name="l00073"></a><span class="lineno"> 73</span> std::make_pair(manager, jacobianFunctor),</div>
|
||
<div class="line"><a id="l00074" name="l00074"></a><span class="lineno"> 74</span> Y,</div>
|
||
<div class="line"><a id="l00075" name="l00075"></a><span class="lineno"> 75</span> current_time,</div>
|
||
<div class="line"><a id="l00076" name="l00076"></a><span class="lineno"> 76</span> netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a0a8d820cfeaa92ee31f253795c57e0d1">tMax</a>,</div>
|
||
<div class="line"><a id="l00077" name="l00077"></a><span class="lineno"> 77</span> current_initial_timestep,</div>
|
||
<div class="line"><a id="l00078" name="l00078"></a><span class="lineno"> 78</span> [&](<span class="keyword">const</span> <span class="keyword">auto</span>& state, <span class="keywordtype">double</span> t) {</div>
|
||
<div class="line"><a id="l00079" name="l00079"></a><span class="lineno"> 79</span> current_time = t;</div>
|
||
<div class="line"><a id="l00080" name="l00080"></a><span class="lineno"> 80</span> manager.observe(state, t);</div>
|
||
<div class="line"><a id="l00081" name="l00081"></a><span class="lineno"> 81</span> }</div>
|
||
<div class="line"><a id="l00082" name="l00082"></a><span class="lineno"> 82</span> );</div>
|
||
<div class="line"><a id="l00083" name="l00083"></a><span class="lineno"> 83</span> current_time = netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a0a8d820cfeaa92ee31f253795c57e0d1">tMax</a>;</div>
|
||
<div class="line"><a id="l00084" name="l00084"></a><span class="lineno"> 84</span> } <span class="keywordflow">catch</span> (<span class="keyword">const</span> <a class="code hl_class" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html">exceptions::StaleEngineTrigger</a> &e) {</div>
|
||
<div class="line"><a id="l00085" name="l00085"></a><span class="lineno"> 85</span> LOG_INFO(<a class="code hl_variable" href="classgridfire_1_1solver_1_1_direct_network_solver.html#a093aa89fd23c2fe03266e286871c7079">m_logger</a>, <span class="stringliteral">"Catching StaleEngineTrigger at t = {:0.3E} with T9 = {:0.3E}, rho = {:0.3E}. Triggering update stage (last stage took {} steps)."</span>, current_time, T9, netIn.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a06f0dff9f8927b7cf2da3004c8fa1577">density</a>, e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a0b7c627c3e69390808bef352b3875408">totalSteps</a>());</div>
|
||
<div class="line"><a id="l00086" name="l00086"></a><span class="lineno"> 86</span> <a class="code hl_struct" href="structgridfire_1_1exceptions_1_1_stale_engine_trigger_1_1state.html">exceptions::StaleEngineTrigger::state</a> staleState = e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a06eb9eb03750038c193c23c7f53668f5">getState</a>();</div>
|
||
<div class="line"><a id="l00087" name="l00087"></a><span class="lineno"> 87</span> accumulated_energy += e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#aeebfb529118f8dfcaf1422ae1768f2bf">energy</a>(); <span class="comment">// Add the specific energy rate to the accumulated energy</span></div>
|
||
<div class="line"><a id="l00088" name="l00088"></a><span class="lineno"> 88</span> <span class="comment">// total_update_stages_triggered++;</span></div>
|
||
<div class="line"><a id="l00089" name="l00089"></a><span class="lineno"> 89</span> </div>
|
||
<div class="line"><a id="l00090" name="l00090"></a><span class="lineno"> 90</span> Composition temp_comp;</div>
|
||
<div class="line"><a id="l00091" name="l00091"></a><span class="lineno"> 91</span> std::vector<double> mass_fractions;</div>
|
||
<div class="line"><a id="l00092" name="l00092"></a><span class="lineno"> 92</span> <span class="keywordtype">size_t</span> num_species_at_stop = e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a44ac2f7510ecf86cd5b556a842eee30c">numSpecies</a>();</div>
|
||
<div class="line"><a id="l00093" name="l00093"></a><span class="lineno"> 93</span> </div>
|
||
<div class="line"><a id="l00094" name="l00094"></a><span class="lineno"> 94</span> <span class="keywordflow">if</span> (num_species_at_stop != <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies().size()) {</div>
|
||
<div class="line"><a id="l00095" name="l00095"></a><span class="lineno"> 95</span> <span class="keywordflow">throw</span> std::runtime_error(</div>
|
||
<div class="line"><a id="l00096" name="l00096"></a><span class="lineno"> 96</span> <span class="stringliteral">"StaleEngineError state has a different number of species than the engine. This should not happen."</span></div>
|
||
<div class="line"><a id="l00097" name="l00097"></a><span class="lineno"> 97</span> );</div>
|
||
<div class="line"><a id="l00098" name="l00098"></a><span class="lineno"> 98</span> }</div>
|
||
<div class="line"><a id="l00099" name="l00099"></a><span class="lineno"> 99</span> mass_fractions.reserve(num_species_at_stop);</div>
|
||
<div class="line"><a id="l00100" name="l00100"></a><span class="lineno"> 100</span> </div>
|
||
<div class="line"><a id="l00101" name="l00101"></a><span class="lineno"> 101</span> <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i < num_species_at_stop; ++i) {</div>
|
||
<div class="line"><a id="l00102" name="l00102"></a><span class="lineno"> 102</span> <span class="keyword">const</span> <span class="keyword">auto</span>& species = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies()[i];</div>
|
||
<div class="line"><a id="l00103" name="l00103"></a><span class="lineno"> 103</span> temp_comp.registerSpecies(species);</div>
|
||
<div class="line"><a id="l00104" name="l00104"></a><span class="lineno"> 104</span> mass_fractions.push_back(e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a3183a0030b14253eb49d8304fb6665ea">getMolarAbundance</a>(i) * species.mass()); <span class="comment">// Convert from molar abundance to mass fraction</span></div>
|
||
<div class="line"><a id="l00105" name="l00105"></a><span class="lineno"> 105</span> }</div>
|
||
<div class="line"><a id="l00106" name="l00106"></a><span class="lineno"> 106</span> temp_comp.setMassFraction(<a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies(), mass_fractions);</div>
|
||
<div class="line"><a id="l00107" name="l00107"></a><span class="lineno"> 107</span> temp_comp.finalize(<span class="keyword">true</span>);</div>
|
||
<div class="line"><a id="l00108" name="l00108"></a><span class="lineno"> 108</span> </div>
|
||
<div class="line"><a id="l00109" name="l00109"></a><span class="lineno"> 109</span> <a class="code hl_struct" href="structgridfire_1_1_net_in.html">NetIn</a> netInTemp = netIn;</div>
|
||
<div class="line"><a id="l00110" name="l00110"></a><span class="lineno"> 110</span> netInTemp.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a5be0f5195a5cd1dd177b9fc5ab83a7be">temperature</a> = e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a2f5925b67562cebd08568fce76c739e9">temperature</a>();</div>
|
||
<div class="line"><a id="l00111" name="l00111"></a><span class="lineno"> 111</span> netInTemp.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a06f0dff9f8927b7cf2da3004c8fa1577">density</a> = e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#ae8156ed7e659cb629da24a5b6734e2dc">density</a>();</div>
|
||
<div class="line"><a id="l00112" name="l00112"></a><span class="lineno"> 112</span> netInTemp.<a class="code hl_variable" href="structgridfire_1_1_net_in.html#a13058f4929e72c1187abbebcddb8aed1">composition</a> = std::move(temp_comp);</div>
|
||
<div class="line"><a id="l00113" name="l00113"></a><span class="lineno"> 113</span> </div>
|
||
<div class="line"><a id="l00114" name="l00114"></a><span class="lineno"> 114</span> Composition currentComposition = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.update(netInTemp);</div>
|
||
<div class="line"><a id="l00115" name="l00115"></a><span class="lineno"> 115</span> populateY(currentComposition);</div>
|
||
<div class="line"><a id="l00116" name="l00116"></a><span class="lineno"> 116</span> Y(Y.size() - 1) = e.<a class="code hl_function" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#aeebfb529118f8dfcaf1422ae1768f2bf">energy</a>(); <span class="comment">// Set the specific energy rate from the stale state</span></div>
|
||
<div class="line"><a id="l00117" name="l00117"></a><span class="lineno"> 117</span> numSpecies = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies().size();</div>
|
||
<div class="line"><a id="l00118" name="l00118"></a><span class="lineno"> 118</span> </div>
|
||
<div class="line"><a id="l00119" name="l00119"></a><span class="lineno"> 119</span> <span class="comment">// current_initial_timestep = 0.001 * manager.m_last_step_time; // set the new timestep to the last successful timestep before repartitioning</span></div>
|
||
<div class="line"><a id="l00120" name="l00120"></a><span class="lineno"> 120</span> }</div>
|
||
<div class="line"><a id="l00121" name="l00121"></a><span class="lineno"> 121</span> }</div>
|
||
<div class="line"><a id="l00122" name="l00122"></a><span class="lineno"> 122</span> </div>
|
||
<div class="line"><a id="l00123" name="l00123"></a><span class="lineno"> 123</span> accumulated_energy += Y(Y.size() - 1); <span class="comment">// Add the specific energy rate to the accumulated energy</span></div>
|
||
<div class="line"><a id="l00124" name="l00124"></a><span class="lineno"> 124</span> </div>
|
||
<div class="line"><a id="l00125" name="l00125"></a><span class="lineno"> 125</span> std::vector<double> finalMassFractions(numSpecies);</div>
|
||
<div class="line"><a id="l00126" name="l00126"></a><span class="lineno"> 126</span> <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i < numSpecies; ++i) {</div>
|
||
<div class="line"><a id="l00127" name="l00127"></a><span class="lineno"> 127</span> <span class="keyword">const</span> <span class="keywordtype">double</span> molarMass = <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies()[i].mass();</div>
|
||
<div class="line"><a id="l00128" name="l00128"></a><span class="lineno"> 128</span> finalMassFractions[i] = Y(i) * molarMass; <span class="comment">// Convert from molar abundance to mass fraction</span></div>
|
||
<div class="line"><a id="l00129" name="l00129"></a><span class="lineno"> 129</span> <span class="keywordflow">if</span> (finalMassFractions[i] < <a class="code hl_variable" href="namespacegridfire.html#a96c062f94713921e5d7568ecedcdcb06">MIN_ABUNDANCE_THRESHOLD</a>) {</div>
|
||
<div class="line"><a id="l00130" name="l00130"></a><span class="lineno"> 130</span> finalMassFractions[i] = 0.0;</div>
|
||
<div class="line"><a id="l00131" name="l00131"></a><span class="lineno"> 131</span> }</div>
|
||
<div class="line"><a id="l00132" name="l00132"></a><span class="lineno"> 132</span> }</div>
|
||
<div class="line"><a id="l00133" name="l00133"></a><span class="lineno"> 133</span> </div>
|
||
<div class="line"><a id="l00134" name="l00134"></a><span class="lineno"> 134</span> std::vector<std::string> speciesNames;</div>
|
||
<div class="line"><a id="l00135" name="l00135"></a><span class="lineno"> 135</span> speciesNames.reserve(numSpecies);</div>
|
||
<div class="line"><a id="l00136" name="l00136"></a><span class="lineno"> 136</span> <span class="keywordflow">for</span> (<span class="keyword">const</span> <span class="keyword">auto</span>& species : <a class="code hl_variable" href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">m_engine</a>.getNetworkSpecies()) {</div>
|
||
<div class="line"><a id="l00137" name="l00137"></a><span class="lineno"> 137</span> speciesNames.push_back(std::string(species.name()));</div>
|
||
<div class="line"><a id="l00138" name="l00138"></a><span class="lineno"> 138</span> }</div>
|
||
<div class="line"><a id="l00139" name="l00139"></a><span class="lineno"> 139</span> </div>
|
||
<div class="line"><a id="l00140" name="l00140"></a><span class="lineno"> 140</span> Composition outputComposition(speciesNames);</div>
|
||
<div class="line"><a id="l00141" name="l00141"></a><span class="lineno"> 141</span> outputComposition.setMassFraction(speciesNames, finalMassFractions);</div>
|
||
<div class="line"><a id="l00142" name="l00142"></a><span class="lineno"> 142</span> outputComposition.finalize(<span class="keyword">true</span>);</div>
|
||
<div class="line"><a id="l00143" name="l00143"></a><span class="lineno"> 143</span> </div>
|
||
<div class="line"><a id="l00144" name="l00144"></a><span class="lineno"> 144</span> <a class="code hl_struct" href="structgridfire_1_1_net_out.html">NetOut</a> netOut;</div>
|
||
<div class="line"><a id="l00145" name="l00145"></a><span class="lineno"> 145</span> netOut.<a class="code hl_variable" href="structgridfire_1_1_net_out.html#a073529511ae0e52f868b47cce0e8ac0a">composition</a> = std::move(outputComposition);</div>
|
||
<div class="line"><a id="l00146" name="l00146"></a><span class="lineno"> 146</span> netOut.<a class="code hl_variable" href="structgridfire_1_1_net_out.html#a43d5a861708992c949f616aa2a035ec6">energy</a> = accumulated_energy; <span class="comment">// Specific energy rate</span></div>
|
||
<div class="line"><a id="l00147" name="l00147"></a><span class="lineno"> 147</span> netOut.<a class="code hl_variable" href="structgridfire_1_1_net_out.html#a51c16703132cf739ec2fd89eae7badd6">num_steps</a> = manager.<a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">m_num_steps</a>;</div>
|
||
<div class="line"><a id="l00148" name="l00148"></a><span class="lineno"> 148</span> </div>
|
||
<div class="line"><a id="l00149" name="l00149"></a><span class="lineno"> 149</span> <span class="keywordflow">return</span> netOut;</div>
|
||
<div class="line"><a id="l00150" name="l00150"></a><span class="lineno"> 150</span> }</div>
|
||
</div>
|
||
<div class="line"><a id="l00151" name="l00151"></a><span class="lineno"> 151</span> </div>
|
||
<div class="foldopen" id="foldopen00152" data-start="{" data-end="}">
|
||
<div class="line"><a id="l00152" name="l00152"></a><span class="lineno"><a class="line" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aec8c0a0b2fbb71cebb40c263f64385b3"> 152</a></span> <span class="keywordtype">void</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aec8c0a0b2fbb71cebb40c263f64385b3">DirectNetworkSolver::RHSManager::operator()</a>(</div>
|
||
<div class="line"><a id="l00153" name="l00153"></a><span class="lineno"> 153</span> <span class="keyword">const</span> boost::numeric::ublas::vector<double> &Y,</div>
|
||
<div class="line"><a id="l00154" name="l00154"></a><span class="lineno"> 154</span> boost::numeric::ublas::vector<double> &dYdt,</div>
|
||
<div class="line"><a id="l00155" name="l00155"></a><span class="lineno"> 155</span> <span class="keyword">const</span> <span class="keywordtype">double</span> t</div>
|
||
<div class="line"><a id="l00156" name="l00156"></a><span class="lineno"> 156</span> )<span class="keyword"> const </span>{</div>
|
||
<div class="line"><a id="l00157" name="l00157"></a><span class="lineno"> 157</span> <span class="keyword">const</span> <span class="keywordtype">size_t</span> numSpecies = <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a035962dfdfc13d255def98befefcccd9">m_engine</a>.getNetworkSpecies().size();</div>
|
||
<div class="line"><a id="l00158" name="l00158"></a><span class="lineno"> 158</span> <span class="keywordflow">if</span> (t != <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a17b83f2478395c934c4ec2c964e9d35e">m_cached_time</a> || !<a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#acfecb0ebb0429f112d503771764f27ec">m_cached_result</a>.has_value() || <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#acfecb0ebb0429f112d503771764f27ec">m_cached_result</a>.value().dydt.size() != numSpecies + 1) {</div>
|
||
<div class="line"><a id="l00159" name="l00159"></a><span class="lineno"> 159</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a595aa16333693ee2bbcac35aa85a1c2a">compute_and_cache</a>(Y, t);</div>
|
||
<div class="line"><a id="l00160" name="l00160"></a><span class="lineno"> 160</span> }</div>
|
||
<div class="line"><a id="l00161" name="l00161"></a><span class="lineno"> 161</span> <span class="keyword">const</span> <span class="keyword">auto</span>&[dydt, nuclearEnergyGenerationRate] = <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#acfecb0ebb0429f112d503771764f27ec">m_cached_result</a>.value();</div>
|
||
<div class="line"><a id="l00162" name="l00162"></a><span class="lineno"> 162</span> dYdt.resize(numSpecies + 1);</div>
|
||
<div class="line"><a id="l00163" name="l00163"></a><span class="lineno"> 163</span> std::ranges::copy(dydt, dYdt.begin());</div>
|
||
<div class="line"><a id="l00164" name="l00164"></a><span class="lineno"> 164</span> dYdt(numSpecies) = nuclearEnergyGenerationRate; <span class="comment">// Set the last element to the specific energy rate</span></div>
|
||
<div class="line"><a id="l00165" name="l00165"></a><span class="lineno"> 165</span> }</div>
|
||
</div>
|
||
<div class="line"><a id="l00166" name="l00166"></a><span class="lineno"> 166</span> </div>
|
||
<div class="foldopen" id="foldopen00167" data-start="{" data-end="}">
|
||
<div class="line"><a id="l00167" name="l00167"></a><span class="lineno"><a class="line" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a226b007bfc9960b5c0bb7b88b4f122da"> 167</a></span> <span class="keywordtype">void</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a226b007bfc9960b5c0bb7b88b4f122da">DirectNetworkSolver::RHSManager::observe</a>(</div>
|
||
<div class="line"><a id="l00168" name="l00168"></a><span class="lineno"> 168</span> <span class="keyword">const</span> boost::numeric::ublas::vector<double> &state,</div>
|
||
<div class="line"><a id="l00169" name="l00169"></a><span class="lineno"> 169</span> <span class="keyword">const</span> <span class="keywordtype">double</span> t</div>
|
||
<div class="line"><a id="l00170" name="l00170"></a><span class="lineno"> 170</span> )<span class="keyword"> const </span>{</div>
|
||
<div class="line"><a id="l00171" name="l00171"></a><span class="lineno"> 171</span> <span class="keywordtype">double</span> dt = t - <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a49268e65b89444c3caf1e69323ce545b">m_last_observed_time</a>;</div>
|
||
<div class="line"><a id="l00172" name="l00172"></a><span class="lineno"> 172</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a595aa16333693ee2bbcac35aa85a1c2a">compute_and_cache</a>(state, t);</div>
|
||
<div class="line"><a id="l00173" name="l00173"></a><span class="lineno"> 173</span> LOG_INFO(</div>
|
||
<div class="line"><a id="l00174" name="l00174"></a><span class="lineno"> 174</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a6cc605a83b5ac5ae048d1044be284ada">m_logger</a>,</div>
|
||
<div class="line"><a id="l00175" name="l00175"></a><span class="lineno"> 175</span> <span class="stringliteral">"(Step {}) Observed state at t = {:0.3E} (dt = {:0.3E})"</span>,</div>
|
||
<div class="line"><a id="l00176" name="l00176"></a><span class="lineno"> 176</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">m_num_steps</a>,</div>
|
||
<div class="line"><a id="l00177" name="l00177"></a><span class="lineno"> 177</span> t,</div>
|
||
<div class="line"><a id="l00178" name="l00178"></a><span class="lineno"> 178</span> dt</div>
|
||
<div class="line"><a id="l00179" name="l00179"></a><span class="lineno"> 179</span> );</div>
|
||
<div class="line"><a id="l00180" name="l00180"></a><span class="lineno"> 180</span> std::ostringstream oss;</div>
|
||
<div class="line"><a id="l00181" name="l00181"></a><span class="lineno"> 181</span> oss << std::scientific << std::setprecision(3);</div>
|
||
<div class="line"><a id="l00182" name="l00182"></a><span class="lineno"> 182</span> oss << <span class="stringliteral">"(Step: "</span> << std::setw(10) << <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">m_num_steps</a> << <span class="stringliteral">") t = "</span> << t << <span class="stringliteral">" (dt = "</span> << dt << <span class="stringliteral">", eps_nuc: "</span> << state(state.size() - 1) << <span class="stringliteral">" [erg])\n"</span>;</div>
|
||
<div class="line"><a id="l00183" name="l00183"></a><span class="lineno"> 183</span> std::cout << oss.str();</div>
|
||
<div class="line"><a id="l00184" name="l00184"></a><span class="lineno"> 184</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a49268e65b89444c3caf1e69323ce545b">m_last_observed_time</a> = t;</div>
|
||
<div class="line"><a id="l00185" name="l00185"></a><span class="lineno"> 185</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a69d773a1cfe4804876dbf23de1f212c9">m_last_step_time</a> = dt;</div>
|
||
<div class="line"><a id="l00186" name="l00186"></a><span class="lineno"> 186</span> </div>
|
||
<div class="line"><a id="l00187" name="l00187"></a><span class="lineno"> 187</span> }</div>
|
||
</div>
|
||
<div class="line"><a id="l00188" name="l00188"></a><span class="lineno"> 188</span> </div>
|
||
<div class="foldopen" id="foldopen00189" data-start="{" data-end="}">
|
||
<div class="line"><a id="l00189" name="l00189"></a><span class="lineno"><a class="line" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a595aa16333693ee2bbcac35aa85a1c2a"> 189</a></span> <span class="keywordtype">void</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a595aa16333693ee2bbcac35aa85a1c2a">DirectNetworkSolver::RHSManager::compute_and_cache</a>(</div>
|
||
<div class="line"><a id="l00190" name="l00190"></a><span class="lineno"> 190</span> <span class="keyword">const</span> boost::numeric::ublas::vector<double> &state,</div>
|
||
<div class="line"><a id="l00191" name="l00191"></a><span class="lineno"> 191</span> <span class="keywordtype">double</span> t</div>
|
||
<div class="line"><a id="l00192" name="l00192"></a><span class="lineno"> 192</span> )<span class="keyword"> const </span>{</div>
|
||
<div class="line"><a id="l00193" name="l00193"></a><span class="lineno"> 193</span> std::vector<double> y_vec(state.begin(), state.end() - 1);</div>
|
||
<div class="line"><a id="l00194" name="l00194"></a><span class="lineno"> 194</span> std::ranges::replace_if(</div>
|
||
<div class="line"><a id="l00195" name="l00195"></a><span class="lineno"> 195</span> y_vec,</div>
|
||
<div class="line"><a id="l00196" name="l00196"></a><span class="lineno"> 196</span> [](<span class="keyword">const</span> <span class="keywordtype">double</span> yi){</div>
|
||
<div class="line"><a id="l00197" name="l00197"></a><span class="lineno"> 197</span> <span class="keywordflow">return</span> yi < 0.0;</div>
|
||
<div class="line"><a id="l00198" name="l00198"></a><span class="lineno"> 198</span> },</div>
|
||
<div class="line"><a id="l00199" name="l00199"></a><span class="lineno"> 199</span> 0.0 <span class="comment">// Avoid negative abundances</span></div>
|
||
<div class="line"><a id="l00200" name="l00200"></a><span class="lineno"> 200</span> );</div>
|
||
<div class="line"><a id="l00201" name="l00201"></a><span class="lineno"> 201</span> </div>
|
||
<div class="line"><a id="l00202" name="l00202"></a><span class="lineno"> 202</span> <span class="keyword">const</span> <span class="keyword">auto</span> result = <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a035962dfdfc13d255def98befefcccd9">m_engine</a>.calculateRHSAndEnergy(y_vec, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a46e39ab9f9fd2f3822c72712173d7aef">m_T9</a>, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aa5d0316fa2fd7d817cc77303776ab446">m_rho</a>);</div>
|
||
<div class="line"><a id="l00203" name="l00203"></a><span class="lineno"> 203</span> <span class="keywordflow">if</span> (!result) {</div>
|
||
<div class="line"><a id="l00204" name="l00204"></a><span class="lineno"> 204</span> LOG_INFO(<a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a6cc605a83b5ac5ae048d1044be284ada">m_logger</a>,</div>
|
||
<div class="line"><a id="l00205" name="l00205"></a><span class="lineno"> 205</span> <span class="stringliteral">"Triggering update stage due to stale engine detected at t = {:0.3E} with T9 = {:0.3E}, rho = {:0.3E}, y_vec (size: {})"</span>,</div>
|
||
<div class="line"><a id="l00206" name="l00206"></a><span class="lineno"> 206</span> t, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a46e39ab9f9fd2f3822c72712173d7aef">m_T9</a>, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aa5d0316fa2fd7d817cc77303776ab446">m_rho</a>, y_vec.size());</div>
|
||
<div class="line"><a id="l00207" name="l00207"></a><span class="lineno"> 207</span> <span class="keywordflow">throw</span> <a class="code hl_class" href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html">exceptions::StaleEngineTrigger</a>({<a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a46e39ab9f9fd2f3822c72712173d7aef">m_T9</a>, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aa5d0316fa2fd7d817cc77303776ab446">m_rho</a>, y_vec, t, <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">m_num_steps</a>, state(state.size() - 1)});</div>
|
||
<div class="line"><a id="l00208" name="l00208"></a><span class="lineno"> 208</span> }</div>
|
||
<div class="line"><a id="l00209" name="l00209"></a><span class="lineno"> 209</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#acfecb0ebb0429f112d503771764f27ec">m_cached_result</a> = result.value();</div>
|
||
<div class="line"><a id="l00210" name="l00210"></a><span class="lineno"> 210</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a17b83f2478395c934c4ec2c964e9d35e">m_cached_time</a> = t;</div>
|
||
<div class="line"><a id="l00211" name="l00211"></a><span class="lineno"> 211</span> </div>
|
||
<div class="line"><a id="l00212" name="l00212"></a><span class="lineno"> 212</span> <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">m_num_steps</a>++;</div>
|
||
<div class="line"><a id="l00213" name="l00213"></a><span class="lineno"> 213</span> }</div>
|
||
</div>
|
||
<div class="line"><a id="l00214" name="l00214"></a><span class="lineno"> 214</span> </div>
|
||
<div class="foldopen" id="foldopen00215" data-start="{" data-end="}">
|
||
<div class="line"><a id="l00215" name="l00215"></a><span class="lineno"><a class="line" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#ad4fc0b406218aa2f567f5c61a4eda3c9"> 215</a></span> <span class="keywordtype">void</span> <a class="code hl_function" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#ad4fc0b406218aa2f567f5c61a4eda3c9">DirectNetworkSolver::JacobianFunctor::operator()</a>(</div>
|
||
<div class="line"><a id="l00216" name="l00216"></a><span class="lineno"> 216</span> <span class="keyword">const</span> boost::numeric::ublas::vector<double> &Y,</div>
|
||
<div class="line"><a id="l00217" name="l00217"></a><span class="lineno"> 217</span> boost::numeric::ublas::matrix<double> &J,</div>
|
||
<div class="line"><a id="l00218" name="l00218"></a><span class="lineno"> 218</span> <span class="keywordtype">double</span> t,</div>
|
||
<div class="line"><a id="l00219" name="l00219"></a><span class="lineno"> 219</span> boost::numeric::ublas::vector<double> &dfdt</div>
|
||
<div class="line"><a id="l00220" name="l00220"></a><span class="lineno"> 220</span> )<span class="keyword"> const </span>{</div>
|
||
<div class="line"><a id="l00221" name="l00221"></a><span class="lineno"> 221</span> <span class="keywordtype">size_t</span> numSpecies = <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#a56f8b2b222fb2a7dac190ead0babfdd0">m_engine</a>.getNetworkSpecies().size();</div>
|
||
<div class="line"><a id="l00222" name="l00222"></a><span class="lineno"> 222</span> J.resize(numSpecies+1, numSpecies+1);</div>
|
||
<div class="line"><a id="l00223" name="l00223"></a><span class="lineno"> 223</span> J.clear();</div>
|
||
<div class="line"><a id="l00224" name="l00224"></a><span class="lineno"> 224</span> <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i < numSpecies; ++i) {</div>
|
||
<div class="line"><a id="l00225" name="l00225"></a><span class="lineno"> 225</span> <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> j = 0; j < numSpecies; ++j) {</div>
|
||
<div class="line"><a id="l00226" name="l00226"></a><span class="lineno"> 226</span> J(i, j) = <a class="code hl_variable" href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#a56f8b2b222fb2a7dac190ead0babfdd0">m_engine</a>.getJacobianMatrixEntry(i, j);</div>
|
||
<div class="line"><a id="l00227" name="l00227"></a><span class="lineno"> 227</span> }</div>
|
||
<div class="line"><a id="l00228" name="l00228"></a><span class="lineno"> 228</span> }</div>
|
||
<div class="line"><a id="l00229" name="l00229"></a><span class="lineno"> 229</span> }</div>
|
||
</div>
|
||
<div class="line"><a id="l00230" name="l00230"></a><span class="lineno"> 230</span> </div>
|
||
<div class="line"><a id="l00231" name="l00231"></a><span class="lineno"> 231</span>}</div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html">gridfire::exceptions::StaleEngineTrigger</a></div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00010">error_engine.h:10</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_a06eb9eb03750038c193c23c7f53668f5"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a06eb9eb03750038c193c23c7f53668f5">gridfire::exceptions::StaleEngineTrigger::getState</a></div><div class="ttdeci">state getState() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00027">error_engine.h:27</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_a0b7c627c3e69390808bef352b3875408"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a0b7c627c3e69390808bef352b3875408">gridfire::exceptions::StaleEngineTrigger::totalSteps</a></div><div class="ttdeci">size_t totalSteps() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00035">error_engine.h:35</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_a2f5925b67562cebd08568fce76c739e9"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a2f5925b67562cebd08568fce76c739e9">gridfire::exceptions::StaleEngineTrigger::temperature</a></div><div class="ttdeci">double temperature() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00050">error_engine.h:50</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_a3183a0030b14253eb49d8304fb6665ea"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a3183a0030b14253eb49d8304fb6665ea">gridfire::exceptions::StaleEngineTrigger::getMolarAbundance</a></div><div class="ttdeci">double getMolarAbundance(const size_t index) const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00043">error_engine.h:43</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_a44ac2f7510ecf86cd5b556a842eee30c"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#a44ac2f7510ecf86cd5b556a842eee30c">gridfire::exceptions::StaleEngineTrigger::numSpecies</a></div><div class="ttdeci">size_t numSpecies() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00031">error_engine.h:31</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_ae8156ed7e659cb629da24a5b6734e2dc"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#ae8156ed7e659cb629da24a5b6734e2dc">gridfire::exceptions::StaleEngineTrigger::density</a></div><div class="ttdeci">double density() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00054">error_engine.h:54</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1exceptions_1_1_stale_engine_trigger_html_aeebfb529118f8dfcaf1422ae1768f2bf"><div class="ttname"><a href="classgridfire_1_1exceptions_1_1_stale_engine_trigger.html#aeebfb529118f8dfcaf1422ae1768f2bf">gridfire::exceptions::StaleEngineTrigger::energy</a></div><div class="ttdeci">double energy() const</div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00039">error_engine.h:39</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1solver_1_1_direct_network_solver_html_a093aa89fd23c2fe03266e286871c7079"><div class="ttname"><a href="classgridfire_1_1solver_1_1_direct_network_solver.html#a093aa89fd23c2fe03266e286871c7079">gridfire::solver::DirectNetworkSolver::m_logger</a></div><div class="ttdeci">quill::Logger * m_logger</div><div class="ttdoc">Logger instance.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00180">solver.h:180</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1solver_1_1_direct_network_solver_html_a0e8a4b8ef656e0b084d11bea982e412a"><div class="ttname"><a href="classgridfire_1_1solver_1_1_direct_network_solver.html#a0e8a4b8ef656e0b084d11bea982e412a">gridfire::solver::DirectNetworkSolver::evaluate</a></div><div class="ttdeci">NetOut evaluate(const NetIn &netIn) override</div><div class="ttdoc">Evaluates the network for a given timestep using direct integration.</div><div class="ttdef"><b>Definition</b> <a href="#l00022">solver.cpp:22</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1solver_1_1_direct_network_solver_html_a2cc12e737a753a42b72a45be3fbfa8ab"><div class="ttname"><a href="classgridfire_1_1solver_1_1_direct_network_solver.html#a2cc12e737a753a42b72a45be3fbfa8ab">gridfire::solver::DirectNetworkSolver::m_config</a></div><div class="ttdeci">Config & m_config</div><div class="ttdoc">Configuration instance.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00181">solver.h:181</a></div></div>
|
||
<div class="ttc" id="aclassgridfire_1_1solver_1_1_network_solver_strategy_html_a724924d94eaf82b67d9988a55c3261e8"><div class="ttname"><a href="classgridfire_1_1solver_1_1_network_solver_strategy.html#a724924d94eaf82b67d9988a55c3261e8">gridfire::solver::NetworkSolverStrategy< DynamicEngine >::m_engine</a></div><div class="ttdeci">DynamicEngine & m_engine</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00047">solver.h:47</a></div></div>
|
||
<div class="ttc" id="aengine__graph_8h_html"><div class="ttname"><a href="engine__graph_8h.html">engine_graph.h</a></div></div>
|
||
<div class="ttc" id="aerror__engine_8h_html"><div class="ttname"><a href="error__engine_8h.html">error_engine.h</a></div></div>
|
||
<div class="ttc" id="anamespacegridfire_1_1solver_html"><div class="ttname"><a href="namespacegridfire_1_1solver.html">gridfire::solver</a></div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00015">solver.h:15</a></div></div>
|
||
<div class="ttc" id="anamespacegridfire_html_a96c062f94713921e5d7568ecedcdcb06"><div class="ttname"><a href="namespacegridfire.html#a96c062f94713921e5d7568ecedcdcb06">gridfire::MIN_ABUNDANCE_THRESHOLD</a></div><div class="ttdeci">static constexpr double MIN_ABUNDANCE_THRESHOLD</div><div class="ttdoc">Minimum abundance threshold below which species are ignored.</div><div class="ttdef"><b>Definition</b> <a href="engine__graph_8h_source.html#l00066">engine_graph.h:66</a></div></div>
|
||
<div class="ttc" id="anetwork_8h_html"><div class="ttname"><a href="network_8h.html">network.h</a></div></div>
|
||
<div class="ttc" id="asolver_8h_html"><div class="ttname"><a href="solver_8h.html">solver.h</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html"><div class="ttname"><a href="structgridfire_1_1_net_in.html">gridfire::NetIn</a></div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00053">network.h:53</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html_a06f0dff9f8927b7cf2da3004c8fa1577"><div class="ttname"><a href="structgridfire_1_1_net_in.html#a06f0dff9f8927b7cf2da3004c8fa1577">gridfire::NetIn::density</a></div><div class="ttdeci">double density</div><div class="ttdoc">Density in g/cm^3.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00058">network.h:58</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html_a0a8d820cfeaa92ee31f253795c57e0d1"><div class="ttname"><a href="structgridfire_1_1_net_in.html#a0a8d820cfeaa92ee31f253795c57e0d1">gridfire::NetIn::tMax</a></div><div class="ttdeci">double tMax</div><div class="ttdoc">Maximum time.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00055">network.h:55</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html_a13058f4929e72c1187abbebcddb8aed1"><div class="ttname"><a href="structgridfire_1_1_net_in.html#a13058f4929e72c1187abbebcddb8aed1">gridfire::NetIn::composition</a></div><div class="ttdeci">fourdst::composition::Composition composition</div><div class="ttdoc">Composition of the network.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00054">network.h:54</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html_a4e556f7bb18f46654b3445476734076a"><div class="ttname"><a href="structgridfire_1_1_net_in.html#a4e556f7bb18f46654b3445476734076a">gridfire::NetIn::dt0</a></div><div class="ttdeci">double dt0</div><div class="ttdoc">Initial time step.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00056">network.h:56</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_in_html_a5be0f5195a5cd1dd177b9fc5ab83a7be"><div class="ttname"><a href="structgridfire_1_1_net_in.html#a5be0f5195a5cd1dd177b9fc5ab83a7be">gridfire::NetIn::temperature</a></div><div class="ttdeci">double temperature</div><div class="ttdoc">Temperature in Kelvin.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00057">network.h:57</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_out_html"><div class="ttname"><a href="structgridfire_1_1_net_out.html">gridfire::NetOut</a></div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00065">network.h:65</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_out_html_a073529511ae0e52f868b47cce0e8ac0a"><div class="ttname"><a href="structgridfire_1_1_net_out.html#a073529511ae0e52f868b47cce0e8ac0a">gridfire::NetOut::composition</a></div><div class="ttdeci">fourdst::composition::Composition composition</div><div class="ttdoc">Composition of the network after evaluation.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00066">network.h:66</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_out_html_a43d5a861708992c949f616aa2a035ec6"><div class="ttname"><a href="structgridfire_1_1_net_out.html#a43d5a861708992c949f616aa2a035ec6">gridfire::NetOut::energy</a></div><div class="ttdeci">double energy</div><div class="ttdoc">Energy in ergs after evaluation.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00068">network.h:68</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1_net_out_html_a51c16703132cf739ec2fd89eae7badd6"><div class="ttname"><a href="structgridfire_1_1_net_out.html#a51c16703132cf739ec2fd89eae7badd6">gridfire::NetOut::num_steps</a></div><div class="ttdeci">int num_steps</div><div class="ttdoc">Number of steps taken in the evaluation.</div><div class="ttdef"><b>Definition</b> <a href="network_8h_source.html#l00067">network.h:67</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1exceptions_1_1_stale_engine_trigger_1_1state_html"><div class="ttname"><a href="structgridfire_1_1exceptions_1_1_stale_engine_trigger_1_1state.html">gridfire::exceptions::StaleEngineTrigger::state</a></div><div class="ttdef"><b>Definition</b> <a href="error__engine_8h_source.html#l00012">error_engine.h:12</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor_html"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html">gridfire::solver::DirectNetworkSolver::JacobianFunctor</a></div><div class="ttdoc">Functor for calculating the Jacobian matrix.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00143">solver.h:143</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor_html_a56f8b2b222fb2a7dac190ead0babfdd0"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#a56f8b2b222fb2a7dac190ead0babfdd0">gridfire::solver::DirectNetworkSolver::JacobianFunctor::m_engine</a></div><div class="ttdeci">DynamicEngine & m_engine</div><div class="ttdoc">The engine used to evaluate the network.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00144">solver.h:144</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor_html_ad4fc0b406218aa2f567f5c61a4eda3c9"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_jacobian_functor.html#ad4fc0b406218aa2f567f5c61a4eda3c9">gridfire::solver::DirectNetworkSolver::JacobianFunctor::operator()</a></div><div class="ttdeci">void operator()(const boost::numeric::ublas::vector< double > &Y, boost::numeric::ublas::matrix< double > &J, double t, boost::numeric::ublas::vector< double > &dfdt) const</div><div class="ttdoc">Calculates the Jacobian matrix.</div><div class="ttdef"><b>Definition</b> <a href="#l00215">solver.cpp:215</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html">gridfire::solver::DirectNetworkSolver::RHSManager</a></div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00088">solver.h:88</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a035962dfdfc13d255def98befefcccd9"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a035962dfdfc13d255def98befefcccd9">gridfire::solver::DirectNetworkSolver::RHSManager::m_engine</a></div><div class="ttdeci">DynamicEngine & m_engine</div><div class="ttdoc">The engine used to evaluate the network.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00089">solver.h:89</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a17b83f2478395c934c4ec2c964e9d35e"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a17b83f2478395c934c4ec2c964e9d35e">gridfire::solver::DirectNetworkSolver::RHSManager::m_cached_time</a></div><div class="ttdeci">double m_cached_time</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00093">solver.h:93</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a226b007bfc9960b5c0bb7b88b4f122da"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a226b007bfc9960b5c0bb7b88b4f122da">gridfire::solver::DirectNetworkSolver::RHSManager::observe</a></div><div class="ttdeci">void observe(const boost::numeric::ublas::vector< double > &state, double t) const</div><div class="ttdef"><b>Definition</b> <a href="#l00167">solver.cpp:167</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a46e39ab9f9fd2f3822c72712173d7aef"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a46e39ab9f9fd2f3822c72712173d7aef">gridfire::solver::DirectNetworkSolver::RHSManager::m_T9</a></div><div class="ttdeci">const double m_T9</div><div class="ttdoc">Temperature in units of 10^9 K.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00090">solver.h:90</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a49268e65b89444c3caf1e69323ce545b"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a49268e65b89444c3caf1e69323ce545b">gridfire::solver::DirectNetworkSolver::RHSManager::m_last_observed_time</a></div><div class="ttdeci">double m_last_observed_time</div><div class="ttdoc">Last time the state was observed.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00096">solver.h:96</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a595aa16333693ee2bbcac35aa85a1c2a"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a595aa16333693ee2bbcac35aa85a1c2a">gridfire::solver::DirectNetworkSolver::RHSManager::compute_and_cache</a></div><div class="ttdeci">void compute_and_cache(const boost::numeric::ublas::vector< double > &state, double t) const</div><div class="ttdef"><b>Definition</b> <a href="#l00189">solver.cpp:189</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a69d773a1cfe4804876dbf23de1f212c9"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a69d773a1cfe4804876dbf23de1f212c9">gridfire::solver::DirectNetworkSolver::RHSManager::m_last_step_time</a></div><div class="ttdeci">double m_last_step_time</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00101">solver.h:101</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_a6cc605a83b5ac5ae048d1044be284ada"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#a6cc605a83b5ac5ae048d1044be284ada">gridfire::solver::DirectNetworkSolver::RHSManager::m_logger</a></div><div class="ttdeci">quill::Logger * m_logger</div><div class="ttdoc">Logger instance.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00099">solver.h:99</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_aa5d0316fa2fd7d817cc77303776ab446"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aa5d0316fa2fd7d817cc77303776ab446">gridfire::solver::DirectNetworkSolver::RHSManager::m_rho</a></div><div class="ttdeci">const double m_rho</div><div class="ttdoc">Density in g/cm^3.</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00091">solver.h:91</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_acfecb0ebb0429f112d503771764f27ec"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#acfecb0ebb0429f112d503771764f27ec">gridfire::solver::DirectNetworkSolver::RHSManager::m_cached_result</a></div><div class="ttdeci">std::optional< StepDerivatives< double > > m_cached_result</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00094">solver.h:94</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_ad9a07ff5cbe42a9455561903a0ae1708"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#ad9a07ff5cbe42a9455561903a0ae1708">gridfire::solver::DirectNetworkSolver::RHSManager::m_num_steps</a></div><div class="ttdeci">int m_num_steps</div><div class="ttdef"><b>Definition</b> <a href="solver_8h_source.html#l00100">solver.h:100</a></div></div>
|
||
<div class="ttc" id="astructgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager_html_aec8c0a0b2fbb71cebb40c263f64385b3"><div class="ttname"><a href="structgridfire_1_1solver_1_1_direct_network_solver_1_1_r_h_s_manager.html#aec8c0a0b2fbb71cebb40c263f64385b3">gridfire::solver::DirectNetworkSolver::RHSManager::operator()</a></div><div class="ttdeci">void operator()(const boost::numeric::ublas::vector< double > &Y, boost::numeric::ublas::vector< double > &dYdt, double t) const</div><div class="ttdoc">Calculates the time derivatives of the species abundances.</div><div class="ttdef"><b>Definition</b> <a href="#l00152">solver.cpp:152</a></div></div>
|
||
</div><!-- fragment --></div><!-- contents -->
|
||
</div><!-- doc-content -->
|
||
<!-- start footer part -->
|
||
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
|
||
<ul>
|
||
<li class="navelem"><a class="el" href="dir_68267d1309a1af8e8297ef4c3efbcdba.html">src</a></li><li class="navelem"><a class="el" href="dir_c85d3e3c5052e9ad9ce18c6863244a25.html">lib</a></li><li class="navelem"><a class="el" href="dir_4fd0dc9a50f7a53e22cb356c650f915e.html">solver</a></li><li class="navelem"><a class="el" href="solver_8cpp.html">solver.cpp</a></li>
|
||
<li class="footer">Generated by <a href="https://www.doxygen.org/index.html"><img class="footer" src="doxygen.svg" width="104" height="31" alt="doxygen"/></a> 1.13.2 </li>
|
||
</ul>
|
||
</div>
|
||
</body>
|
||
</html>
|