-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathassembly__base_8cpp_source.html
196 lines (194 loc) · 122 KB
/
assembly__base_8cpp_source.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
<!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.11"/>
<title>MAST: /home/travis/build/MASTmultiphysics/mast-multiphysics/src/base/assembly_base.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="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);
$(window).load(resizeHeight);
</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/javascript">
$(document).ready(function() { init_search(); });
</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">MAST
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.11 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
</div><!-- 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('assembly__base_8cpp_source.html','');});
</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">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div class="header">
<div class="headertitle">
<div class="title">assembly_base.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="assembly__base_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/*</span></div><div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit</span></div><div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment"> * Copyright (C) 2013-2019 Manav Bhatia</span></div><div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> *</span></div><div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> * This library is free software; you can redistribute it and/or</span></div><div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> * modify it under the terms of the GNU Lesser General Public</span></div><div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> * License as published by the Free Software Foundation; either</span></div><div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> * version 2.1 of the License, or (at your option) any later version.</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> *</span></div><div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment"> * This library is distributed in the hope that it will be useful,</span></div><div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> * but WITHOUT ANY WARRANTY; without even the implied warranty of</span></div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment"> * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU</span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"> * Lesser General Public License for more details.</span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"> *</span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> * You should have received a copy of the GNU Lesser General Public</span></div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment"> * License along with this library; if not, write to the Free Software</span></div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"> * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA</span></div><div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment"> */</span></div><div class="line"><a name="l00019"></a><span class="lineno"> 19</span> </div><div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment">// MAST includes</span></div><div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="preprocessor">#include "<a class="code" href="assembly__base_8h.html">base/assembly_base.h</a>"</span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor">#include "<a class="code" href="system__initialization_8h.html">base/system_initialization.h</a>"</span></div><div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="preprocessor">#include "<a class="code" href="mesh__field__function_8h.html">base/mesh_field_function.h</a>"</span></div><div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="preprocessor">#include "<a class="code" href="elem__base_8h.html">base/elem_base.h</a>"</span></div><div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#include "<a class="code" href="physics__discipline__base_8h.html">base/physics_discipline_base.h</a>"</span></div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="preprocessor">#include "<a class="code" href="nonlinear__system_8h.html">base/nonlinear_system.h</a>"</span></div><div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="preprocessor">#include "<a class="code" href="assembly__elem__operation_8h.html">base/assembly_elem_operation.h</a>"</span></div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor">#include "<a class="code" href="output__assembly__elem__operations_8h.html">base/output_assembly_elem_operations.h</a>"</span></div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="preprocessor">#include "<a class="code" href="fe__base_8h.html">mesh/fe_base.h</a>"</span></div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor">#include "<a class="code" href="geom__elem_8h.html">mesh/geom_elem.h</a>"</span></div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="preprocessor">#include "<a class="code" href="utility_8h.html">numerics/utility.h</a>"</span></div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> </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">// libMesh includes</span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="preprocessor">#include "libmesh/numeric_vector.h"</span></div><div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="preprocessor">#include "libmesh/dof_map.h"</span></div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span> </div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span> </div><div class="line"><a name="l00039"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a9a4d4c88402ce1d98278c7b1aac25d1b"> 39</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a9a4d4c88402ce1d98278c7b1aac25d1b">MAST::AssemblyBase::AssemblyBase</a>():</div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span> close_matrix (true),</div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span> _elem_ops (nullptr),</div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span> _discipline (nullptr),</div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span> _system (nullptr),</div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span> _sol_function (nullptr),</div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span> _solver_monitor (nullptr),</div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span> _param_dependence (nullptr) {</div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  </div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span> }</div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span> </div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span> </div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span> </div><div class="line"><a name="l00052"></a><span class="lineno"> 52</span> </div><div class="line"><a name="l00053"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a2034db25d6074eef21899ab93c66a833"> 53</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a2034db25d6074eef21899ab93c66a833">MAST::AssemblyBase::~AssemblyBase</a>() {</div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  </div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a9f7f65a2be3dc8f6edb9000020aa8ccf">clear_discipline_and_system</a>();</div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad607ffca1c704155f56f503cbfe6350a">clear_elem_operation_object</a>();</div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span> }</div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span> </div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span> </div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span> </div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_physics_discipline_base.html">MAST::PhysicsDisciplineBase</a>&</div><div class="line"><a name="l00062"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462"> 62</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462">MAST::AssemblyBase::discipline</a>()<span class="keyword"> const </span>{</div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  </div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>,</div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <span class="stringliteral">"Error: Discipline not yet attached to Assembly."</span>);</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <span class="keywordflow">return</span> *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>;</div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span> }</div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span> </div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span> </div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span> </div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span> <a class="code" href="class_m_a_s_t_1_1_physics_discipline_base.html">MAST::PhysicsDisciplineBase</a>&</div><div class="line"><a name="l00072"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#ad8d7a29701f9a3b51f878178b87591c2"> 72</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462">MAST::AssemblyBase::discipline</a>() {</div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  </div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>,</div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="stringliteral">"Error: Discipline not yet attached to Assembly."</span>);</div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <span class="keywordflow">return</span> *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>;</div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span> }</div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span> </div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span> </div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span> </div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span> </div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span> <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>&</div><div class="line"><a name="l00083"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3"> 83</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">MAST::AssemblyBase::system</a>()<span class="keyword"> const </span>{</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  </div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>,</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="stringliteral">"Error: System not yet attached to Assembly."</span>);</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <span class="keywordflow">return</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span> }</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span> </div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span> </div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span> <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>&</div><div class="line"><a name="l00092"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a0e73ea60b9606d702e9f9af93f4fd423"> 92</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">MAST::AssemblyBase::system</a>() {</div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  </div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>,</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="stringliteral">"Error: System not yet attached to Assembly."</span>);</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordflow">return</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span> }</div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span> </div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span> </div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span> <a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html">MAST::AssemblyElemOperations</a>&</div><div class="line"><a name="l00101"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a0b4114f497996001f1e54f452a926786"> 101</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0b4114f497996001f1e54f452a926786">MAST::AssemblyBase::get_elem_ops</a>() {</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  </div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>,</div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  <span class="stringliteral">"Error: Not yet initialized."</span>);</div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  <span class="keywordflow">return</span> *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>;</div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span> }</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span> </div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span> </div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span> <a class="code" href="class_m_a_s_t_1_1_system_initialization.html">MAST::SystemInitialization</a>&</div><div class="line"><a name="l00110"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a12fa42f74de999be67da555082a32c5d"> 110</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a12fa42f74de999be67da555082a32c5d">MAST::AssemblyBase::system_init</a>() {</div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  </div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  libmesh_assert_msg(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>,</div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="stringliteral">"Error: System not yet attached to Assembly."</span>);</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  <span class="keywordflow">return</span> *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>;</div><div class="line"><a name="l00115"></a><span class="lineno"> 115</span> }</div><div class="line"><a name="l00116"></a><span class="lineno"> 116</span> </div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span> </div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00119"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#aff8ece1161fffa7fc8afc31a8f261455"> 119</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aff8ece1161fffa7fc8afc31a8f261455">MAST::AssemblyBase::set_solver_monitor</a>(<a class="code" href="class_m_a_s_t_1_1_assembly_base_1_1_solver_monitor.html">MAST::AssemblyBase::SolverMonitor</a>& monitor) {</div><div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  </div><div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  libmesh_assert(!<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abb647f77c49b15b60250bf1a5081dbe4">_solver_monitor</a>);</div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abb647f77c49b15b60250bf1a5081dbe4">_solver_monitor</a> = &monitor;</div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span> }</div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div><div class="line"><a name="l00125"></a><span class="lineno"> 125</span> </div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span> </div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7a236e900be0ba6dbfe630ebb51e3152">MAST::AssemblyBase::attach_elem_parameter_dependence_object</a></div><div class="line"><a name="l00129"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a7a236e900be0ba6dbfe630ebb51e3152"> 129</a></span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html">MAST::AssemblyBase::ElemParameterDependence</a>& dep) {</div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span> </div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  libmesh_assert(!<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a>);</div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  </div><div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a> = &dep;</div><div class="line"><a name="l00134"></a><span class="lineno"> 134</span> }</div><div class="line"><a name="l00135"></a><span class="lineno"> 135</span> </div><div class="line"><a name="l00136"></a><span class="lineno"> 136</span> </div><div class="line"><a name="l00137"></a><span class="lineno"> 137</span> </div><div class="line"><a name="l00138"></a><span class="lineno"> 138</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00139"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#af8281c5eece58746b64293f8ccea49f3"> 139</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#af8281c5eece58746b64293f8ccea49f3">MAST::AssemblyBase::clear_elem_parameter_dependence_object</a>() {</div><div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  </div><div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00142"></a><span class="lineno"> 142</span> }</div><div class="line"><a name="l00143"></a><span class="lineno"> 143</span> </div><div class="line"><a name="l00144"></a><span class="lineno"> 144</span> </div><div class="line"><a name="l00145"></a><span class="lineno"> 145</span> </div><div class="line"><a name="l00146"></a><span class="lineno"> 146</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base_1_1_solver_monitor.html">MAST::AssemblyBase::SolverMonitor</a>*</div><div class="line"><a name="l00147"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a652097c4673958f03ff6cca8cf037847"> 147</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a652097c4673958f03ff6cca8cf037847">MAST::AssemblyBase::get_solver_monitor</a>() {</div><div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  </div><div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <span class="keywordflow">return</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abb647f77c49b15b60250bf1a5081dbe4">_solver_monitor</a>;</div><div class="line"><a name="l00150"></a><span class="lineno"> 150</span> }</div><div class="line"><a name="l00151"></a><span class="lineno"> 151</span> </div><div class="line"><a name="l00152"></a><span class="lineno"> 152</span> </div><div class="line"><a name="l00153"></a><span class="lineno"> 153</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00154"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a5eb48a95d52bb9454c9efac225c73ba0"> 154</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a5eb48a95d52bb9454c9efac225c73ba0">MAST::AssemblyBase::clear_solver_monitor</a>() {</div><div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  </div><div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abb647f77c49b15b60250bf1a5081dbe4">_solver_monitor</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00157"></a><span class="lineno"> 157</span> }</div><div class="line"><a name="l00158"></a><span class="lineno"> 158</span> </div><div class="line"><a name="l00159"></a><span class="lineno"> 159</span> </div><div class="line"><a name="l00160"></a><span class="lineno"> 160</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00161"></a><span class="lineno"> 161</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0644e12cd53adc64d9b7e027becc08b5">MAST::AssemblyBase::</a></div><div class="line"><a name="l00162"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a0644e12cd53adc64d9b7e027becc08b5"> 162</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0644e12cd53adc64d9b7e027becc08b5">set_discipline_and_system</a>(<a class="code" href="class_m_a_s_t_1_1_physics_discipline_base.html">MAST::PhysicsDisciplineBase</a> &<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462">discipline</a>,</div><div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <a class="code" href="class_m_a_s_t_1_1_system_initialization.html">MAST::SystemInitialization</a> &<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</a>) {</div><div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  </div><div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  libmesh_assert_msg(!<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a> && !<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>,</div><div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  <span class="stringliteral">"Error: Assembly should be cleared before attaching System."</span>);</div><div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  </div><div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a> = &<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462">discipline</a>;</div><div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a> = &<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</a>;</div><div class="line"><a name="l00170"></a><span class="lineno"> 170</span> }</div><div class="line"><a name="l00171"></a><span class="lineno"> 171</span> </div><div class="line"><a name="l00172"></a><span class="lineno"> 172</span> </div><div class="line"><a name="l00173"></a><span class="lineno"> 173</span> </div><div class="line"><a name="l00174"></a><span class="lineno"> 174</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00175"></a><span class="lineno"> 175</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a9f7f65a2be3dc8f6edb9000020aa8ccf">MAST::AssemblyBase::</a></div><div class="line"><a name="l00176"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a9f7f65a2be3dc8f6edb9000020aa8ccf"> 176</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a9f7f65a2be3dc8f6edb9000020aa8ccf">clear_discipline_and_system</a>() {</div><div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  </div><div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a13f08f11353ef98dda92f3ccacc92902">close_matrix</a> = <span class="keyword">true</span>;</div><div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00182"></a><span class="lineno"> 182</span> }</div><div class="line"><a name="l00183"></a><span class="lineno"> 183</span> </div><div class="line"><a name="l00184"></a><span class="lineno"> 184</span> </div><div class="line"><a name="l00185"></a><span class="lineno"> 185</span> </div><div class="line"><a name="l00186"></a><span class="lineno"> 186</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00187"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#adcc03f1662807d86d866ee3a4e2b60df"> 187</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#adcc03f1662807d86d866ee3a4e2b60df">MAST::AssemblyBase::set_elem_operation_object</a>(<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html">MAST::AssemblyElemOperations</a>& elem_ops) {</div><div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  </div><div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  libmesh_assert_msg(!<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>,</div><div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  <span class="stringliteral">"Error: Assembly should be cleared before attaching Elem."</span>);</div><div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  </div><div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a> = &elem_ops;</div><div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#ac8c536e86f20680c1e6e41317bd12d2b">set_assembly</a>(*<span class="keyword">this</span>);</div><div class="line"><a name="l00194"></a><span class="lineno"> 194</span> }</div><div class="line"><a name="l00195"></a><span class="lineno"> 195</span> </div><div class="line"><a name="l00196"></a><span class="lineno"> 196</span> </div><div class="line"><a name="l00197"></a><span class="lineno"> 197</span> </div><div class="line"><a name="l00198"></a><span class="lineno"> 198</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00199"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#ad607ffca1c704155f56f503cbfe6350a"> 199</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad607ffca1c704155f56f503cbfe6350a">MAST::AssemblyBase::clear_elem_operation_object</a>() {</div><div class="line"><a name="l00200"></a><span class="lineno"> 200</span> </div><div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>) {</div><div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a4e6d01c9391e881f98a6caeb9a14503b">clear_assembly</a>();</div><div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">_elem_ops</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  }</div><div class="line"><a name="l00205"></a><span class="lineno"> 205</span> }</div><div class="line"><a name="l00206"></a><span class="lineno"> 206</span> </div><div class="line"><a name="l00207"></a><span class="lineno"> 207</span> </div><div class="line"><a name="l00208"></a><span class="lineno"> 208</span> </div><div class="line"><a name="l00209"></a><span class="lineno"> 209</span> std::unique_ptr<libMesh::NumericVector<Real> ></div><div class="line"><a name="l00210"></a><span class="lineno"> 210</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">MAST::AssemblyBase::</a></div><div class="line"><a name="l00211"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60"> 211</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">build_localized_vector</a>(<span class="keyword">const</span> libMesh::System& sys,</div><div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="keyword">const</span> libMesh::NumericVector<Real>& global)<span class="keyword"> const </span>{</div><div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  </div><div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  libMesh::NumericVector<Real>* local =</div><div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  libMesh::NumericVector<Real>::build(sys.comm()).release();</div><div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  </div><div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  <span class="keyword">const</span> std::vector<libMesh::dof_id_type>& send_list =</div><div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  sys.get_dof_map().get_send_list();</div><div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  </div><div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  local->init(sys.n_dofs(),</div><div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  sys.n_local_dofs(),</div><div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  send_list,</div><div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  <span class="keyword">false</span>,</div><div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  libMesh::GHOSTED);</div><div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  global.localize(*local, send_list);</div><div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  </div><div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  <span class="keywordflow">return</span> std::unique_ptr<libMesh::NumericVector<Real> >(local);</div><div class="line"><a name="l00228"></a><span class="lineno"> 228</span> }</div><div class="line"><a name="l00229"></a><span class="lineno"> 229</span> </div><div class="line"><a name="l00230"></a><span class="lineno"> 230</span> </div><div class="line"><a name="l00231"></a><span class="lineno"> 231</span> </div><div class="line"><a name="l00232"></a><span class="lineno"> 232</span> </div><div class="line"><a name="l00233"></a><span class="lineno"> 233</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00234"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#abb38e8d6b62b5332974755ec4587c1e8"> 234</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abb38e8d6b62b5332974755ec4587c1e8">MAST::AssemblyBase::attach_solution_function</a>(<a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html">MAST::MeshFieldFunction</a>& f){</div><div class="line"><a name="l00235"></a><span class="lineno"> 235</span>  </div><div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  <span class="comment">// make sure that no prior association is specified</span></div><div class="line"><a name="l00237"></a><span class="lineno"> 237</span>  libmesh_assert(!<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>);</div><div class="line"><a name="l00238"></a><span class="lineno"> 238</span>  </div><div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a> = &f;</div><div class="line"><a name="l00240"></a><span class="lineno"> 240</span> }</div><div class="line"><a name="l00241"></a><span class="lineno"> 241</span> </div><div class="line"><a name="l00242"></a><span class="lineno"> 242</span> </div><div class="line"><a name="l00243"></a><span class="lineno"> 243</span> </div><div class="line"><a name="l00244"></a><span class="lineno"> 244</span> </div><div class="line"><a name="l00245"></a><span class="lineno"> 245</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00246"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a97d3049ca6034e65f34cacb9cf5dd75b"> 246</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a97d3049ca6034e65f34cacb9cf5dd75b">MAST::AssemblyBase::detach_solution_function</a>() {</div><div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a> = <span class="keyword">nullptr</span>;</div><div class="line"><a name="l00248"></a><span class="lineno"> 248</span> }</div><div class="line"><a name="l00249"></a><span class="lineno"> 249</span> </div><div class="line"><a name="l00250"></a><span class="lineno"> 250</span> </div><div class="line"><a name="l00251"></a><span class="lineno"> 251</span> </div><div class="line"><a name="l00252"></a><span class="lineno"> 252</span> </div><div class="line"><a name="l00253"></a><span class="lineno"> 253</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00254"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#abbbc4fd5c13c458b2ef6c27636ee7510"> 254</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abbbc4fd5c13c458b2ef6c27636ee7510">MAST::AssemblyBase::calculate_output</a>(<span class="keyword">const</span> libMesh::NumericVector<Real>& X,</div><div class="line"><a name="l00255"></a><span class="lineno"> 255</span>  <a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html">MAST::OutputAssemblyElemOperations</a>& output) {</div><div class="line"><a name="l00256"></a><span class="lineno"> 256</span>  </div><div class="line"><a name="l00257"></a><span class="lineno"> 257</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>);</div><div class="line"><a name="l00258"></a><span class="lineno"> 258</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00259"></a><span class="lineno"> 259</span>  </div><div class="line"><a name="l00260"></a><span class="lineno"> 260</span>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>& nonlin_sys = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00261"></a><span class="lineno"> 261</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#ac8c536e86f20680c1e6e41317bd12d2b">set_assembly</a>(*<span class="keyword">this</span>);</div><div class="line"><a name="l00262"></a><span class="lineno"> 262</span>  </div><div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a6a81c128526ae2dff490652a39d60b11">zero_for_analysis</a>();</div><div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  </div><div class="line"><a name="l00265"></a><span class="lineno"> 265</span>  <span class="comment">// iterate over each element, initialize it and get the relevant</span></div><div class="line"><a name="l00266"></a><span class="lineno"> 266</span>  <span class="comment">// analysis quantities</span></div><div class="line"><a name="l00267"></a><span class="lineno"> 267</span>  <a class="code" href="mast__data__types_8h.html#abc8ad2bc4221d7400e5306159e28e588">RealVectorX</a> sol;</div><div class="line"><a name="l00268"></a><span class="lineno"> 268</span>  </div><div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  std::vector<libMesh::dof_id_type> dof_indices;</div><div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  <span class="keyword">const</span> libMesh::DofMap& dof_map = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>().get_dof_map();</div><div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  </div><div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  std::unique_ptr<libMesh::NumericVector<Real> > localized_solution;</div><div class="line"><a name="l00273"></a><span class="lineno"> 273</span>  localized_solution.reset(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">build_localized_vector</a>(nonlin_sys,</div><div class="line"><a name="l00274"></a><span class="lineno"> 274</span>  X).release());</div><div class="line"><a name="l00275"></a><span class="lineno"> 275</span>  </div><div class="line"><a name="l00276"></a><span class="lineno"> 276</span>  </div><div class="line"><a name="l00277"></a><span class="lineno"> 277</span>  <span class="comment">// if a solution function is attached, initialize it</span></div><div class="line"><a name="l00278"></a><span class="lineno"> 278</span>  <span class="comment">//if (_sol_function)</span></div><div class="line"><a name="l00279"></a><span class="lineno"> 279</span>  <span class="comment">// _sol_function->init( X);</span></div><div class="line"><a name="l00280"></a><span class="lineno"> 280</span>  </div><div class="line"><a name="l00281"></a><span class="lineno"> 281</span>  libMesh::MeshBase::const_element_iterator el =</div><div class="line"><a name="l00282"></a><span class="lineno"> 282</span>  nonlin_sys.get_mesh().active_local_elements_begin();</div><div class="line"><a name="l00283"></a><span class="lineno"> 283</span>  <span class="keyword">const</span> libMesh::MeshBase::const_element_iterator end_el =</div><div class="line"><a name="l00284"></a><span class="lineno"> 284</span>  nonlin_sys.get_mesh().active_local_elements_end();</div><div class="line"><a name="l00285"></a><span class="lineno"> 285</span>  </div><div class="line"><a name="l00286"></a><span class="lineno"> 286</span> </div><div class="line"><a name="l00287"></a><span class="lineno"> 287</span>  <span class="keywordflow">for</span> ( ; el != end_el; ++el) {</div><div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  </div><div class="line"><a name="l00289"></a><span class="lineno"> 289</span>  <span class="keyword">const</span> libMesh::Elem* elem = *el;</div><div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  </div><div class="line"><a name="l00291"></a><span class="lineno"> 291</span>  dof_map.dof_indices (elem, dof_indices);</div><div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  </div><div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  <span class="comment">// get the solution</span></div><div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ndofs = (<span class="keywordtype">unsigned</span> int)dof_indices.size();</div><div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  sol.setZero(ndofs);</div><div class="line"><a name="l00296"></a><span class="lineno"> 296</span>  </div><div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i=0; i<dof_indices.size(); i++)</div><div class="line"><a name="l00298"></a><span class="lineno"> 298</span>  sol(i) = (*localized_solution)(dof_indices[i]);</div><div class="line"><a name="l00299"></a><span class="lineno"> 299</span>  </div><div class="line"><a name="l00300"></a><span class="lineno"> 300</span>  </div><div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  <span class="comment">//if (_sol_function)</span></div><div class="line"><a name="l00302"></a><span class="lineno"> 302</span>  <span class="comment">// physics_elem->attach_active_solution_function(*_sol_function);</span></div><div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  </div><div class="line"><a name="l00304"></a><span class="lineno"> 304</span>  <a class="code" href="class_m_a_s_t_1_1_geom_elem.html">MAST::GeomElem</a> geom_elem;</div><div class="line"><a name="l00305"></a><span class="lineno"> 305</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a97d7148b7421efdfc8eadcd4e821ee1b">set_elem_data</a>(elem->dim(), *elem, geom_elem);</div><div class="line"><a name="l00306"></a><span class="lineno"> 306</span>  geom_elem.init(*elem, *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00307"></a><span class="lineno"> 307</span>  </div><div class="line"><a name="l00308"></a><span class="lineno"> 308</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a73b35bea5577b08a6a8ee70b19dac59c">init</a>(geom_elem);</div><div class="line"><a name="l00309"></a><span class="lineno"> 309</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a8cdd4717b7d4fab6fd7f864d629432e5">set_elem_solution</a>(sol);</div><div class="line"><a name="l00310"></a><span class="lineno"> 310</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a9d75826c8f6d322ddb4b40a33e6829cf">evaluate</a>();</div><div class="line"><a name="l00311"></a><span class="lineno"> 311</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#aad8893093e07905b7e5071468b23514c">clear_elem</a>();</div><div class="line"><a name="l00312"></a><span class="lineno"> 312</span> </div><div class="line"><a name="l00313"></a><span class="lineno"> 313</span>  <span class="comment">//physics_elem->detach_active_solution_function();</span></div><div class="line"><a name="l00314"></a><span class="lineno"> 314</span>  }</div><div class="line"><a name="l00315"></a><span class="lineno"> 315</span>  </div><div class="line"><a name="l00316"></a><span class="lineno"> 316</span>  <span class="comment">// if a solution function is attached, clear it</span></div><div class="line"><a name="l00317"></a><span class="lineno"> 317</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>)</div><div class="line"><a name="l00318"></a><span class="lineno"> 318</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>-><a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html#a45a0f5907f61898b83d649050617357d">clear</a>();</div><div class="line"><a name="l00319"></a><span class="lineno"> 319</span>  </div><div class="line"><a name="l00320"></a><span class="lineno"> 320</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a4e6d01c9391e881f98a6caeb9a14503b">clear_assembly</a>();</div><div class="line"><a name="l00321"></a><span class="lineno"> 321</span> }</div><div class="line"><a name="l00322"></a><span class="lineno"> 322</span> </div><div class="line"><a name="l00323"></a><span class="lineno"> 323</span> </div><div class="line"><a name="l00324"></a><span class="lineno"> 324</span> </div><div class="line"><a name="l00325"></a><span class="lineno"> 325</span> </div><div class="line"><a name="l00326"></a><span class="lineno"> 326</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00327"></a><span class="lineno"> 327</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#acb7019e60e22b2f28799de154a5571da">MAST::AssemblyBase::</a></div><div class="line"><a name="l00328"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#acb7019e60e22b2f28799de154a5571da"> 328</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#acb7019e60e22b2f28799de154a5571da">calculate_output_derivative</a>(<span class="keyword">const</span> libMesh::NumericVector<Real>& X,</div><div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  <a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html">MAST::OutputAssemblyElemOperations</a>& output,</div><div class="line"><a name="l00330"></a><span class="lineno"> 330</span>  libMesh::NumericVector<Real>& dq_dX) {</div><div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  </div><div class="line"><a name="l00332"></a><span class="lineno"> 332</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>);</div><div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00334"></a><span class="lineno"> 334</span> </div><div class="line"><a name="l00335"></a><span class="lineno"> 335</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a8deb480f8b65fd855bca1742f19829be">zero_for_sensitivity</a>();</div><div class="line"><a name="l00336"></a><span class="lineno"> 336</span>  </div><div class="line"><a name="l00337"></a><span class="lineno"> 337</span>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>& nonlin_sys = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00338"></a><span class="lineno"> 338</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#ac8c536e86f20680c1e6e41317bd12d2b">set_assembly</a>(*<span class="keyword">this</span>);</div><div class="line"><a name="l00339"></a><span class="lineno"> 339</span>  </div><div class="line"><a name="l00340"></a><span class="lineno"> 340</span>  dq_dX.zero();</div><div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  </div><div class="line"><a name="l00342"></a><span class="lineno"> 342</span>  <span class="comment">// iterate over each element, initialize it and get the relevant</span></div><div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <span class="comment">// analysis quantities</span></div><div class="line"><a name="l00344"></a><span class="lineno"> 344</span>  <a class="code" href="mast__data__types_8h.html#abc8ad2bc4221d7400e5306159e28e588">RealVectorX</a> vec, sol;</div><div class="line"><a name="l00345"></a><span class="lineno"> 345</span>  </div><div class="line"><a name="l00346"></a><span class="lineno"> 346</span>  std::vector<libMesh::dof_id_type> dof_indices;</div><div class="line"><a name="l00347"></a><span class="lineno"> 347</span>  <span class="keyword">const</span> libMesh::DofMap& dof_map = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>().get_dof_map();</div><div class="line"><a name="l00348"></a><span class="lineno"> 348</span>  </div><div class="line"><a name="l00349"></a><span class="lineno"> 349</span>  </div><div class="line"><a name="l00350"></a><span class="lineno"> 350</span>  std::unique_ptr<libMesh::NumericVector<Real> > localized_solution;</div><div class="line"><a name="l00351"></a><span class="lineno"> 351</span>  localized_solution.reset(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">build_localized_vector</a>(nonlin_sys,</div><div class="line"><a name="l00352"></a><span class="lineno"> 352</span>  X).release());</div><div class="line"><a name="l00353"></a><span class="lineno"> 353</span>  </div><div class="line"><a name="l00354"></a><span class="lineno"> 354</span>  </div><div class="line"><a name="l00355"></a><span class="lineno"> 355</span>  <span class="comment">// if a solution function is attached, initialize it</span></div><div class="line"><a name="l00356"></a><span class="lineno"> 356</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>)</div><div class="line"><a name="l00357"></a><span class="lineno"> 357</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>-><a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html#ad05c1a6c795fbee828c4a3daae01fecd">init</a>( X);</div><div class="line"><a name="l00358"></a><span class="lineno"> 358</span>  </div><div class="line"><a name="l00359"></a><span class="lineno"> 359</span>  </div><div class="line"><a name="l00360"></a><span class="lineno"> 360</span>  libMesh::MeshBase::const_element_iterator el =</div><div class="line"><a name="l00361"></a><span class="lineno"> 361</span>  nonlin_sys.get_mesh().active_local_elements_begin();</div><div class="line"><a name="l00362"></a><span class="lineno"> 362</span>  <span class="keyword">const</span> libMesh::MeshBase::const_element_iterator end_el =</div><div class="line"><a name="l00363"></a><span class="lineno"> 363</span>  nonlin_sys.get_mesh().active_local_elements_end();</div><div class="line"><a name="l00364"></a><span class="lineno"> 364</span>  </div><div class="line"><a name="l00365"></a><span class="lineno"> 365</span>  </div><div class="line"><a name="l00366"></a><span class="lineno"> 366</span>  <span class="keywordflow">for</span> ( ; el != end_el; ++el) {</div><div class="line"><a name="l00367"></a><span class="lineno"> 367</span>  </div><div class="line"><a name="l00368"></a><span class="lineno"> 368</span>  <span class="keyword">const</span> libMesh::Elem* elem = *el;</div><div class="line"><a name="l00369"></a><span class="lineno"> 369</span>  </div><div class="line"><a name="l00370"></a><span class="lineno"> 370</span>  dof_map.dof_indices (elem, dof_indices);</div><div class="line"><a name="l00371"></a><span class="lineno"> 371</span>  </div><div class="line"><a name="l00372"></a><span class="lineno"> 372</span>  <span class="comment">// get the solution</span></div><div class="line"><a name="l00373"></a><span class="lineno"> 373</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ndofs = (<span class="keywordtype">unsigned</span> int)dof_indices.size();</div><div class="line"><a name="l00374"></a><span class="lineno"> 374</span>  sol.setZero(ndofs);</div><div class="line"><a name="l00375"></a><span class="lineno"> 375</span>  vec.setZero(ndofs);</div><div class="line"><a name="l00376"></a><span class="lineno"> 376</span>  </div><div class="line"><a name="l00377"></a><span class="lineno"> 377</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i=0; i<dof_indices.size(); i++)</div><div class="line"><a name="l00378"></a><span class="lineno"> 378</span>  sol(i) = (*localized_solution)(dof_indices[i]);</div><div class="line"><a name="l00379"></a><span class="lineno"> 379</span>  </div><div class="line"><a name="l00380"></a><span class="lineno"> 380</span> <span class="comment">// if (_sol_function)</span></div><div class="line"><a name="l00381"></a><span class="lineno"> 381</span> <span class="comment">// physics_elem->attach_active_solution_function(*_sol_function);</span></div><div class="line"><a name="l00382"></a><span class="lineno"> 382</span>  </div><div class="line"><a name="l00383"></a><span class="lineno"> 383</span>  <a class="code" href="class_m_a_s_t_1_1_geom_elem.html">MAST::GeomElem</a> geom_elem;</div><div class="line"><a name="l00384"></a><span class="lineno"> 384</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a97d7148b7421efdfc8eadcd4e821ee1b">set_elem_data</a>(elem->dim(), *elem, geom_elem);</div><div class="line"><a name="l00385"></a><span class="lineno"> 385</span>  geom_elem.init(*elem, *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00386"></a><span class="lineno"> 386</span>  </div><div class="line"><a name="l00387"></a><span class="lineno"> 387</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a73b35bea5577b08a6a8ee70b19dac59c">init</a>(geom_elem);</div><div class="line"><a name="l00388"></a><span class="lineno"> 388</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a8cdd4717b7d4fab6fd7f864d629432e5">set_elem_solution</a>(sol);</div><div class="line"><a name="l00389"></a><span class="lineno"> 389</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a16a069336ddd4ccc3b71bf737b555483">output_derivative_for_elem</a>(vec);</div><div class="line"><a name="l00390"></a><span class="lineno"> 390</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#aad8893093e07905b7e5071468b23514c">clear_elem</a>();</div><div class="line"><a name="l00391"></a><span class="lineno"> 391</span>  </div><div class="line"><a name="l00392"></a><span class="lineno"> 392</span>  <a class="code" href="mast__data__types_8h.html#a579078c9edfe62906476a5e85dea6361">DenseRealVector</a> v;</div><div class="line"><a name="l00393"></a><span class="lineno"> 393</span>  <a class="code" href="namespace_m_a_s_t.html#a659caf861e2d27b164e12aef83e5ac1e">MAST::copy</a>(v, vec);</div><div class="line"><a name="l00394"></a><span class="lineno"> 394</span>  dof_map.constrain_element_vector(v, dof_indices);</div><div class="line"><a name="l00395"></a><span class="lineno"> 395</span>  dq_dX.add_vector(v, dof_indices);</div><div class="line"><a name="l00396"></a><span class="lineno"> 396</span>  dof_indices.clear();</div><div class="line"><a name="l00397"></a><span class="lineno"> 397</span>  }</div><div class="line"><a name="l00398"></a><span class="lineno"> 398</span>  </div><div class="line"><a name="l00399"></a><span class="lineno"> 399</span>  <span class="comment">// if a solution function is attached, clear it</span></div><div class="line"><a name="l00400"></a><span class="lineno"> 400</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>)</div><div class="line"><a name="l00401"></a><span class="lineno"> 401</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>-><a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html#a45a0f5907f61898b83d649050617357d">clear</a>();</div><div class="line"><a name="l00402"></a><span class="lineno"> 402</span>  </div><div class="line"><a name="l00403"></a><span class="lineno"> 403</span>  dq_dX.close();</div><div class="line"><a name="l00404"></a><span class="lineno"> 404</span>  </div><div class="line"><a name="l00405"></a><span class="lineno"> 405</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a4e6d01c9391e881f98a6caeb9a14503b">clear_assembly</a>();</div><div class="line"><a name="l00406"></a><span class="lineno"> 406</span> }</div><div class="line"><a name="l00407"></a><span class="lineno"> 407</span> </div><div class="line"><a name="l00408"></a><span class="lineno"> 408</span> </div><div class="line"><a name="l00409"></a><span class="lineno"> 409</span> </div><div class="line"><a name="l00410"></a><span class="lineno"> 410</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00411"></a><span class="lineno"> 411</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0c8b29a632eb1b975005798b3862ad60">MAST::AssemblyBase::</a></div><div class="line"><a name="l00412"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#a0c8b29a632eb1b975005798b3862ad60"> 412</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0c8b29a632eb1b975005798b3862ad60">calculate_output_direct_sensitivity</a>(<span class="keyword">const</span> libMesh::NumericVector<Real>& X,</div><div class="line"><a name="l00413"></a><span class="lineno"> 413</span>  <span class="keyword">const</span> libMesh::NumericVector<Real>* dXdp,</div><div class="line"><a name="l00414"></a><span class="lineno"> 414</span>  <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_function_base.html">MAST::FunctionBase</a>& p,</div><div class="line"><a name="l00415"></a><span class="lineno"> 415</span>  <a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html">MAST::OutputAssemblyElemOperations</a>& output) {</div><div class="line"><a name="l00416"></a><span class="lineno"> 416</span> </div><div class="line"><a name="l00417"></a><span class="lineno"> 417</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>);</div><div class="line"><a name="l00418"></a><span class="lineno"> 418</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00419"></a><span class="lineno"> 419</span> </div><div class="line"><a name="l00420"></a><span class="lineno"> 420</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a8deb480f8b65fd855bca1742f19829be">zero_for_sensitivity</a>();</div><div class="line"><a name="l00421"></a><span class="lineno"> 421</span> </div><div class="line"><a name="l00422"></a><span class="lineno"> 422</span>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>& nonlin_sys = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00423"></a><span class="lineno"> 423</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#ac8c536e86f20680c1e6e41317bd12d2b">set_assembly</a>(*<span class="keyword">this</span>);</div><div class="line"><a name="l00424"></a><span class="lineno"> 424</span>  </div><div class="line"><a name="l00425"></a><span class="lineno"> 425</span>  <span class="comment">// iterate over each element, initialize it and get the relevant</span></div><div class="line"><a name="l00426"></a><span class="lineno"> 426</span>  <span class="comment">// analysis quantities</span></div><div class="line"><a name="l00427"></a><span class="lineno"> 427</span>  <a class="code" href="mast__data__types_8h.html#abc8ad2bc4221d7400e5306159e28e588">RealVectorX</a></div><div class="line"><a name="l00428"></a><span class="lineno"> 428</span>  sol,</div><div class="line"><a name="l00429"></a><span class="lineno"> 429</span>  dsol;</div><div class="line"><a name="l00430"></a><span class="lineno"> 430</span>  </div><div class="line"><a name="l00431"></a><span class="lineno"> 431</span>  std::vector<libMesh::dof_id_type> dof_indices;</div><div class="line"><a name="l00432"></a><span class="lineno"> 432</span>  <span class="keyword">const</span> libMesh::DofMap& dof_map = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>().get_dof_map();</div><div class="line"><a name="l00433"></a><span class="lineno"> 433</span>  </div><div class="line"><a name="l00434"></a><span class="lineno"> 434</span>  </div><div class="line"><a name="l00435"></a><span class="lineno"> 435</span>  std::unique_ptr<libMesh::NumericVector<Real> ></div><div class="line"><a name="l00436"></a><span class="lineno"> 436</span>  localized_solution,</div><div class="line"><a name="l00437"></a><span class="lineno"> 437</span>  localized_solution_sens;</div><div class="line"><a name="l00438"></a><span class="lineno"> 438</span>  localized_solution.reset(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">build_localized_vector</a>(nonlin_sys,</div><div class="line"><a name="l00439"></a><span class="lineno"> 439</span>  X).release());</div><div class="line"><a name="l00440"></a><span class="lineno"> 440</span>  <span class="keywordflow">if</span> (dXdp)</div><div class="line"><a name="l00441"></a><span class="lineno"> 441</span>  localized_solution_sens.reset(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">build_localized_vector</a>(nonlin_sys,</div><div class="line"><a name="l00442"></a><span class="lineno"> 442</span>  *dXdp).release());</div><div class="line"><a name="l00443"></a><span class="lineno"> 443</span> </div><div class="line"><a name="l00444"></a><span class="lineno"> 444</span>  </div><div class="line"><a name="l00445"></a><span class="lineno"> 445</span>  <span class="comment">// if a solution function is attached, initialize it</span></div><div class="line"><a name="l00446"></a><span class="lineno"> 446</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>)</div><div class="line"><a name="l00447"></a><span class="lineno"> 447</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>-><a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html#ad05c1a6c795fbee828c4a3daae01fecd">init</a>( X);</div><div class="line"><a name="l00448"></a><span class="lineno"> 448</span>  </div><div class="line"><a name="l00449"></a><span class="lineno"> 449</span>  </div><div class="line"><a name="l00450"></a><span class="lineno"> 450</span>  libMesh::MeshBase::const_element_iterator el =</div><div class="line"><a name="l00451"></a><span class="lineno"> 451</span>  nonlin_sys.get_mesh().active_local_elements_begin();</div><div class="line"><a name="l00452"></a><span class="lineno"> 452</span>  <span class="keyword">const</span> libMesh::MeshBase::const_element_iterator end_el =</div><div class="line"><a name="l00453"></a><span class="lineno"> 453</span>  nonlin_sys.get_mesh().active_local_elements_end();</div><div class="line"><a name="l00454"></a><span class="lineno"> 454</span>  </div><div class="line"><a name="l00455"></a><span class="lineno"> 455</span>  </div><div class="line"><a name="l00456"></a><span class="lineno"> 456</span>  <span class="keywordflow">for</span> ( ; el != end_el; ++el) {</div><div class="line"><a name="l00457"></a><span class="lineno"> 457</span>  </div><div class="line"><a name="l00458"></a><span class="lineno"> 458</span>  <span class="keyword">const</span> libMesh::Elem* elem = *el;</div><div class="line"><a name="l00459"></a><span class="lineno"> 459</span>  </div><div class="line"><a name="l00460"></a><span class="lineno"> 460</span>  <span class="comment">// no sensitivity computation assembly is neeed in these cases</span></div><div class="line"><a name="l00461"></a><span class="lineno"> 461</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a> &&</div><div class="line"><a name="l00462"></a><span class="lineno"> 462</span>  <span class="comment">// if object is specified and elem does not depend on it</span></div><div class="line"><a name="l00463"></a><span class="lineno"> 463</span>  !<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html#a96e7f90a492e4799f1acccca0a7680a9">if_elem_depends_on_parameter</a>(*elem, p) &&</div><div class="line"><a name="l00464"></a><span class="lineno"> 464</span>  <span class="comment">// and if no sol_sens is given</span></div><div class="line"><a name="l00465"></a><span class="lineno"> 465</span>  (!dXdp ||</div><div class="line"><a name="l00466"></a><span class="lineno"> 466</span>  <span class="comment">// or if it can be ignored for elem</span></div><div class="line"><a name="l00467"></a><span class="lineno"> 467</span>  (dXdp && <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">_param_dependence</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html#af2dbdf12da9e6e2e885106de152fc953">override_flag</a>)))</div><div class="line"><a name="l00468"></a><span class="lineno"> 468</span>  <span class="keywordflow">continue</span>;</div><div class="line"><a name="l00469"></a><span class="lineno"> 469</span>  </div><div class="line"><a name="l00470"></a><span class="lineno"> 470</span>  dof_map.dof_indices (elem, dof_indices);</div><div class="line"><a name="l00471"></a><span class="lineno"> 471</span>  </div><div class="line"><a name="l00472"></a><span class="lineno"> 472</span>  <span class="comment">// get the solution</span></div><div class="line"><a name="l00473"></a><span class="lineno"> 473</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ndofs = (<span class="keywordtype">unsigned</span> int)dof_indices.size();</div><div class="line"><a name="l00474"></a><span class="lineno"> 474</span>  sol.setZero(ndofs);</div><div class="line"><a name="l00475"></a><span class="lineno"> 475</span>  dsol.setZero(ndofs);</div><div class="line"><a name="l00476"></a><span class="lineno"> 476</span>  </div><div class="line"><a name="l00477"></a><span class="lineno"> 477</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i=0; i<dof_indices.size(); i++) {</div><div class="line"><a name="l00478"></a><span class="lineno"> 478</span>  sol(i) = (*localized_solution)(dof_indices[i]);</div><div class="line"><a name="l00479"></a><span class="lineno"> 479</span>  <span class="keywordflow">if</span> (dXdp)</div><div class="line"><a name="l00480"></a><span class="lineno"> 480</span>  dsol(i) = (*localized_solution_sens)(dof_indices[i]);</div><div class="line"><a name="l00481"></a><span class="lineno"> 481</span>  }</div><div class="line"><a name="l00482"></a><span class="lineno"> 482</span>  </div><div class="line"><a name="l00483"></a><span class="lineno"> 483</span>  <span class="comment">// if (_sol_function)</span></div><div class="line"><a name="l00484"></a><span class="lineno"> 484</span>  <span class="comment">// physics_elem->attach_active_solution_function(*_sol_function);</span></div><div class="line"><a name="l00485"></a><span class="lineno"> 485</span>  </div><div class="line"><a name="l00486"></a><span class="lineno"> 486</span>  <a class="code" href="class_m_a_s_t_1_1_geom_elem.html">MAST::GeomElem</a> geom_elem;</div><div class="line"><a name="l00487"></a><span class="lineno"> 487</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a97d7148b7421efdfc8eadcd4e821ee1b">set_elem_data</a>(elem->dim(), *elem, geom_elem);</div><div class="line"><a name="l00488"></a><span class="lineno"> 488</span>  geom_elem.init(*elem, *<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00489"></a><span class="lineno"> 489</span>  </div><div class="line"><a name="l00490"></a><span class="lineno"> 490</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a73b35bea5577b08a6a8ee70b19dac59c">init</a>(geom_elem);</div><div class="line"><a name="l00491"></a><span class="lineno"> 491</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a8cdd4717b7d4fab6fd7f864d629432e5">set_elem_solution</a>(sol);</div><div class="line"><a name="l00492"></a><span class="lineno"> 492</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#af939a3e0e70c537ddb2cd2aa16b81935">set_elem_solution_sensitivity</a>(dsol);</div><div class="line"><a name="l00493"></a><span class="lineno"> 493</span>  output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#ae895cc8daf0b4d118ea0db9a7b036004">evaluate_sensitivity</a>(p);</div><div class="line"><a name="l00494"></a><span class="lineno"> 494</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#aad8893093e07905b7e5071468b23514c">clear_elem</a>();</div><div class="line"><a name="l00495"></a><span class="lineno"> 495</span>  </div><div class="line"><a name="l00496"></a><span class="lineno"> 496</span>  <span class="comment">// physics_elem->detach_active_solution_function();</span></div><div class="line"><a name="l00497"></a><span class="lineno"> 497</span>  }</div><div class="line"><a name="l00498"></a><span class="lineno"> 498</span>  </div><div class="line"><a name="l00499"></a><span class="lineno"> 499</span>  <span class="comment">// if a solution function is attached, clear it</span></div><div class="line"><a name="l00500"></a><span class="lineno"> 500</span>  <span class="keywordflow">if</span> (<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>)</div><div class="line"><a name="l00501"></a><span class="lineno"> 501</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">_sol_function</a>-><a class="code" href="class_m_a_s_t_1_1_mesh_field_function.html#a45a0f5907f61898b83d649050617357d">clear</a>();</div><div class="line"><a name="l00502"></a><span class="lineno"> 502</span>  </div><div class="line"><a name="l00503"></a><span class="lineno"> 503</span>  output.<a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html#a4e6d01c9391e881f98a6caeb9a14503b">clear_assembly</a>();</div><div class="line"><a name="l00504"></a><span class="lineno"> 504</span> }</div><div class="line"><a name="l00505"></a><span class="lineno"> 505</span> </div><div class="line"><a name="l00506"></a><span class="lineno"> 506</span> </div><div class="line"><a name="l00507"></a><span class="lineno"> 507</span> </div><div class="line"><a name="l00508"></a><span class="lineno"> 508</span> </div><div class="line"><a name="l00509"></a><span class="lineno"> 509</span> <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00510"></a><span class="lineno"> 510</span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ab6d7c383e893dca4de61657ee7f042cf">MAST::AssemblyBase::</a></div><div class="line"><a name="l00511"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_assembly_base.html#ab6d7c383e893dca4de61657ee7f042cf"> 511</a></span> <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ab6d7c383e893dca4de61657ee7f042cf">calculate_output_adjoint_sensitivity</a>(<span class="keyword">const</span> libMesh::NumericVector<Real>& X,</div><div class="line"><a name="l00512"></a><span class="lineno"> 512</span>  <span class="keyword">const</span> libMesh::NumericVector<Real>& dq_dX,</div><div class="line"><a name="l00513"></a><span class="lineno"> 513</span>  <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_function_base.html">MAST::FunctionBase</a>& p,</div><div class="line"><a name="l00514"></a><span class="lineno"> 514</span>  <a class="code" href="class_m_a_s_t_1_1_assembly_elem_operations.html">MAST::AssemblyElemOperations</a>& elem_ops,</div><div class="line"><a name="l00515"></a><span class="lineno"> 515</span>  <a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html">MAST::OutputAssemblyElemOperations</a>& output,</div><div class="line"><a name="l00516"></a><span class="lineno"> 516</span>  <span class="keyword">const</span> <span class="keywordtype">bool</span> include_partial_sens) {</div><div class="line"><a name="l00517"></a><span class="lineno"> 517</span> </div><div class="line"><a name="l00518"></a><span class="lineno"> 518</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">_discipline</a>);</div><div class="line"><a name="l00519"></a><span class="lineno"> 519</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>);</div><div class="line"><a name="l00520"></a><span class="lineno"> 520</span>  </div><div class="line"><a name="l00521"></a><span class="lineno"> 521</span>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>& nonlin_sys = <a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">_system</a>-><a class="code" href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">system</a>();</div><div class="line"><a name="l00522"></a><span class="lineno"> 522</span> </div><div class="line"><a name="l00523"></a><span class="lineno"> 523</span>  libMesh::NumericVector<Real></div><div class="line"><a name="l00524"></a><span class="lineno"> 524</span>  &dres_dp = nonlin_sys.add_sensitivity_rhs();</div><div class="line"><a name="l00525"></a><span class="lineno"> 525</span>  </div><div class="line"><a name="l00526"></a><span class="lineno"> 526</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#adcc03f1662807d86d866ee3a4e2b60df">set_elem_operation_object</a>(elem_ops);</div><div class="line"><a name="l00527"></a><span class="lineno"> 527</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a4eb0ce2595ad5f7b9279825fba621d79">sensitivity_assemble</a>(p, dres_dp);</div><div class="line"><a name="l00528"></a><span class="lineno"> 528</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#ad607ffca1c704155f56f503cbfe6350a">clear_elem_operation_object</a>();</div><div class="line"><a name="l00529"></a><span class="lineno"> 529</span> </div><div class="line"><a name="l00530"></a><span class="lineno"> 530</span>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00531"></a><span class="lineno"> 531</span>  dq_dp = dq_dX.dot(dres_dp);</div><div class="line"><a name="l00532"></a><span class="lineno"> 532</span> </div><div class="line"><a name="l00533"></a><span class="lineno"> 533</span>  <span class="keywordflow">if</span> (include_partial_sens) {</div><div class="line"><a name="l00534"></a><span class="lineno"> 534</span> </div><div class="line"><a name="l00535"></a><span class="lineno"> 535</span>  <span class="comment">// calculate the partial sensitivity of the output, which is done</span></div><div class="line"><a name="l00536"></a><span class="lineno"> 536</span>  <span class="comment">// with zero solution vector</span></div><div class="line"><a name="l00537"></a><span class="lineno"> 537</span>  this-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a0c8b29a632eb1b975005798b3862ad60">calculate_output_direct_sensitivity</a>(X, <span class="keyword">nullptr</span>, p, output);</div><div class="line"><a name="l00538"></a><span class="lineno"> 538</span> </div><div class="line"><a name="l00539"></a><span class="lineno"> 539</span>  dq_dp += output.<a class="code" href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a918765a6fc86ba98c02f84bc8ce8367b">output_sensitivity_total</a>(p);</div><div class="line"><a name="l00540"></a><span class="lineno"> 540</span>  }</div><div class="line"><a name="l00541"></a><span class="lineno"> 541</span> </div><div class="line"><a name="l00542"></a><span class="lineno"> 542</span>  <span class="keywordflow">return</span> dq_dp;</div><div class="line"><a name="l00543"></a><span class="lineno"> 543</span> }</div><div class="line"><a name="l00544"></a><span class="lineno"> 544</span> </div><div class="ttc" id="class_m_a_s_t_1_1_function_base_html"><div class="ttname"><a href="class_m_a_s_t_1_1_function_base.html">MAST::FunctionBase</a></div><div class="ttdef"><b>Definition:</b> <a href="function__base_8h_source.html#l00035">function_base.h:35</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_ad47d566c625e170029729b87c50a83a7"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#ad47d566c625e170029729b87c50a83a7">MAST::AssemblyBase::_elem_ops</a></div><div class="ttdeci">MAST::AssemblyElemOperations * _elem_ops</div><div class="ttdoc">provides assembly elem operations for use by this class </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00306">assembly_base.h:306</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_system_initialization_html_a40ed0a5367703c7bac16cc8d0917fc29"><div class="ttname"><a href="class_m_a_s_t_1_1_system_initialization.html#a40ed0a5367703c7bac16cc8d0917fc29">MAST::SystemInitialization::system</a></div><div class="ttdeci">MAST::NonlinearSystem & system()</div><div class="ttdef"><b>Definition:</b> <a href="system__initialization_8h_source.html#l00073">system_initialization.h:73</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence_html_a96e7f90a492e4799f1acccca0a7680a9"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html#a96e7f90a492e4799f1acccca0a7680a9">MAST::AssemblyBase::ElemParameterDependence::if_elem_depends_on_parameter</a></div><div class="ttdeci">virtual bool if_elem_depends_on_parameter(const libMesh::Elem &e, const MAST::FunctionBase &p) const =0</div></div>
<div class="ttc" id="assembly__elem__operation_8h_html"><div class="ttname"><a href="assembly__elem__operation_8h.html">assembly_elem_operation.h</a></div></div>
<div class="ttc" id="fe__base_8h_html"><div class="ttname"><a href="fe__base_8h.html">fe_base.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a9a4d4c88402ce1d98278c7b1aac25d1b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a9a4d4c88402ce1d98278c7b1aac25d1b">MAST::AssemblyBase::AssemblyBase</a></div><div class="ttdeci">AssemblyBase()</div><div class="ttdoc">constructor takes a reference to the discipline that provides the boundary conditions, volume loads, properties, etc. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00039">assembly_base.cpp:39</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a5eb48a95d52bb9454c9efac225c73ba0"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a5eb48a95d52bb9454c9efac225c73ba0">MAST::AssemblyBase::clear_solver_monitor</a></div><div class="ttdeci">void clear_solver_monitor()</div><div class="ttdoc">clears the monitor object </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00154">assembly_base.cpp:154</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a></div><div class="ttdoc">This class implements a system for solution of nonlinear systems. </div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00056">nonlinear_system.h:56</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence_html_af2dbdf12da9e6e2e885106de152fc953"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html#af2dbdf12da9e6e2e885106de152fc953">MAST::AssemblyBase::ElemParameterDependence::override_flag</a></div><div class="ttdeci">bool override_flag</div><div class="ttdoc">if true, assume zero solution sensitivity when elem does not dependent on parameter. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00097">assembly_base.h:97</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_ad607ffca1c704155f56f503cbfe6350a"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#ad607ffca1c704155f56f503cbfe6350a">MAST::AssemblyBase::clear_elem_operation_object</a></div><div class="ttdeci">virtual void clear_elem_operation_object()</div><div class="ttdoc">clears the association of this object with the assembly element operation object. ...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00199">assembly_base.cpp:199</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a4eb0ce2595ad5f7b9279825fba621d79"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a4eb0ce2595ad5f7b9279825fba621d79">MAST::AssemblyBase::sensitivity_assemble</a></div><div class="ttdeci">virtual bool sensitivity_assemble(const MAST::FunctionBase &f, libMesh::NumericVector< Real > &sensitivity_rhs)</div><div class="ttdoc">Assembly function. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00239">assembly_base.h:239</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_1_1_solver_monitor_html"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base_1_1_solver_monitor.html">MAST::AssemblyBase::SolverMonitor</a></div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00067">assembly_base.h:67</a></div></div>
<div class="ttc" id="physics__discipline__base_8h_html"><div class="ttname"><a href="physics__discipline__base_8h.html">physics_discipline_base.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_mesh_field_function_html"><div class="ttname"><a href="class_m_a_s_t_1_1_mesh_field_function.html">MAST::MeshFieldFunction</a></div><div class="ttdoc">This provides a wrapper FieldFunction compatible class that interpolates the solution using libMesh&#39;s...</div><div class="ttdef"><b>Definition:</b> <a href="mesh__field__function_8h_source.html#l00044">mesh_field_function.h:44</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html">MAST::AssemblyElemOperations</a></div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8h_source.html#l00041">assembly_elem_operation.h:41</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_abb38e8d6b62b5332974755ec4587c1e8"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#abb38e8d6b62b5332974755ec4587c1e8">MAST::AssemblyBase::attach_solution_function</a></div><div class="ttdeci">void attach_solution_function(MAST::MeshFieldFunction &f)</div><div class="ttdoc">tells the assembly object that this function is will need to be initialized before each residual eval...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00234">assembly_base.cpp:234</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_acb7019e60e22b2f28799de154a5571da"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#acb7019e60e22b2f28799de154a5571da">MAST::AssemblyBase::calculate_output_derivative</a></div><div class="ttdeci">virtual void calculate_output_derivative(const libMesh::NumericVector< Real > &X, MAST::OutputAssemblyElemOperations &output, libMesh::NumericVector< Real > &dq_dX)</div><div class="ttdoc">calculates </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00328">assembly_base.cpp:328</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a7a236e900be0ba6dbfe630ebb51e3152"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a7a236e900be0ba6dbfe630ebb51e3152">MAST::AssemblyBase::attach_elem_parameter_dependence_object</a></div><div class="ttdeci">void attach_elem_parameter_dependence_object(MAST::AssemblyBase::ElemParameterDependence &dep)</div><div class="ttdoc">This object, if provided by user, will be used to reduce unnecessary computations in sensitivity anal...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00129">assembly_base.cpp:129</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_ab6d7c383e893dca4de61657ee7f042cf"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#ab6d7c383e893dca4de61657ee7f042cf">MAST::AssemblyBase::calculate_output_adjoint_sensitivity</a></div><div class="ttdeci">virtual Real calculate_output_adjoint_sensitivity(const libMesh::NumericVector< Real > &X, const libMesh::NumericVector< Real > &dq_dX, const MAST::FunctionBase &p, MAST::AssemblyElemOperations &elem_ops, MAST::OutputAssemblyElemOperations &output, const bool include_partial_sens=true)</div><div class="ttdoc">Evaluates the total sensitivity of output wrt p using the adjoint solution provided in dq_dX for a li...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00511">assembly_base.cpp:511</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_af939a3e0e70c537ddb2cd2aa16b81935"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#af939a3e0e70c537ddb2cd2aa16b81935">MAST::AssemblyElemOperations::set_elem_solution_sensitivity</a></div><div class="ttdeci">virtual void set_elem_solution_sensitivity(const RealVectorX &sol)</div><div class="ttdoc">sets the element solution sensitivity </div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8cpp_source.html#l00117">assembly_elem_operation.cpp:117</a></div></div>
<div class="ttc" id="output__assembly__elem__operations_8h_html"><div class="ttname"><a href="output__assembly__elem__operations_8h.html">output_assembly_elem_operations.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_aff8ece1161fffa7fc8afc31a8f261455"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#aff8ece1161fffa7fc8afc31a8f261455">MAST::AssemblyBase::set_solver_monitor</a></div><div class="ttdeci">void set_solver_monitor(MAST::AssemblyBase::SolverMonitor &monitor)</div><div class="ttdoc">attaches the solver monitor, which is a user provided routine that is called each time ...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00119">assembly_base.cpp:119</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_ac8c536e86f20680c1e6e41317bd12d2b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#ac8c536e86f20680c1e6e41317bd12d2b">MAST::AssemblyElemOperations::set_assembly</a></div><div class="ttdeci">virtual void set_assembly(MAST::AssemblyBase &assembly)</div><div class="ttdoc">sets the assembly object </div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8cpp_source.html#l00070">assembly_elem_operation.cpp:70</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html">MAST::OutputAssemblyElemOperations</a></div><div class="ttdoc">This provides the base class for definitin of element level contribution of output quantity in an ana...</div><div class="ttdef"><b>Definition:</b> <a href="output__assembly__elem__operations_8h_source.html#l00054">output_assembly_elem_operations.h:54</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_a8deb480f8b65fd855bca1742f19829be"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a8deb480f8b65fd855bca1742f19829be">MAST::OutputAssemblyElemOperations::zero_for_sensitivity</a></div><div class="ttdeci">virtual void zero_for_sensitivity()=0</div><div class="ttdoc">zeroes the output quantity values stored inside this object so that assembly process can begin...</div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a13f08f11353ef98dda92f3ccacc92902"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a13f08f11353ef98dda92f3ccacc92902">MAST::AssemblyBase::close_matrix</a></div><div class="ttdeci">bool close_matrix</div><div class="ttdoc">flag to control the closing fo the Jacobian after assembly </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00103">assembly_base.h:103</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_system_initialization_html"><div class="ttname"><a href="class_m_a_s_t_1_1_system_initialization.html">MAST::SystemInitialization</a></div><div class="ttdef"><b>Definition:</b> <a href="system__initialization_8h_source.html#l00041">system_initialization.h:41</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_adcc03f1662807d86d866ee3a4e2b60df"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#adcc03f1662807d86d866ee3a4e2b60df">MAST::AssemblyBase::set_elem_operation_object</a></div><div class="ttdeci">virtual void set_elem_operation_object(MAST::AssemblyElemOperations &elem_ops)</div><div class="ttdoc">attaches a element operation to this object, and associated this with the element operation object...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00187">assembly_base.cpp:187</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_a8cdd4717b7d4fab6fd7f864d629432e5"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#a8cdd4717b7d4fab6fd7f864d629432e5">MAST::AssemblyElemOperations::set_elem_solution</a></div><div class="ttdeci">virtual void set_elem_solution(const RealVectorX &sol)</div><div class="ttdoc">sets the element solution </div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8cpp_source.html#l00108">assembly_elem_operation.cpp:108</a></div></div>
<div class="ttc" id="elem__base_8h_html"><div class="ttname"><a href="elem__base_8h.html">elem_base.h</a></div></div>
<div class="ttc" id="mast__data__types_8h_html_a0fc21c6ff882da286a44d42c85f5378f"><div class="ttname"><a href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="ttdeci">libMesh::Real Real</div><div class="ttdef"><b>Definition:</b> <a href="mast__data__types_8h_source.html#l00032">mast_data_types.h:32</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_a9d75826c8f6d322ddb4b40a33e6829cf"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a9d75826c8f6d322ddb4b40a33e6829cf">MAST::OutputAssemblyElemOperations::evaluate</a></div><div class="ttdeci">virtual void evaluate()=0</div><div class="ttdoc">this is the abstract interface to be implemented by derived classes. </div></div>
<div class="ttc" id="nonlinear__system_8h_html"><div class="ttname"><a href="nonlinear__system_8h.html">nonlinear_system.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a42583481abcb696e0436541b244dc51d"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a42583481abcb696e0436541b244dc51d">MAST::AssemblyBase::_system</a></div><div class="ttdeci">MAST::SystemInitialization * _system</div><div class="ttdoc">System for which this assembly is performed. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00318">assembly_base.h:318</a></div></div>
<div class="ttc" id="system__initialization_8h_html"><div class="ttname"><a href="system__initialization_8h.html">system_initialization.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a97d3049ca6034e65f34cacb9cf5dd75b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a97d3049ca6034e65f34cacb9cf5dd75b">MAST::AssemblyBase::detach_solution_function</a></div><div class="ttdeci">void detach_solution_function()</div><div class="ttdoc">removes the attachment of the solution function </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00246">assembly_base.cpp:246</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a7c42a8d69afb9d5c63dd9a0d234defae"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a7c42a8d69afb9d5c63dd9a0d234defae">MAST::AssemblyBase::_discipline</a></div><div class="ttdeci">MAST::PhysicsDisciplineBase * _discipline</div><div class="ttdoc">PhysicsDisciplineBase object for which this class is assembling. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00312">assembly_base.h:312</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_physics_discipline_base_html"><div class="ttname"><a href="class_m_a_s_t_1_1_physics_discipline_base.html">MAST::PhysicsDisciplineBase</a></div><div class="ttdef"><b>Definition:</b> <a href="physics__discipline__base_8h_source.html#l00056">physics_discipline_base.h:56</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_a97d7148b7421efdfc8eadcd4e821ee1b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#a97d7148b7421efdfc8eadcd4e821ee1b">MAST::AssemblyElemOperations::set_elem_data</a></div><div class="ttdeci">virtual void set_elem_data(unsigned int dim, const libMesh::Elem &ref_elem, MAST::GeomElem &elem) const =0</div><div class="ttdoc">some analyses may want to set additional element data before initialization of the GeomElem...</div></div>
<div class="ttc" id="geom__elem_8h_html"><div class="ttname"><a href="geom__elem_8h.html">geom_elem.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a9f7f65a2be3dc8f6edb9000020aa8ccf"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a9f7f65a2be3dc8f6edb9000020aa8ccf">MAST::AssemblyBase::clear_discipline_and_system</a></div><div class="ttdeci">virtual void clear_discipline_and_system()</div><div class="ttdoc">clears association with a system to this discipline </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00176">assembly_base.cpp:176</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_abbbc4fd5c13c458b2ef6c27636ee7510"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#abbbc4fd5c13c458b2ef6c27636ee7510">MAST::AssemblyBase::calculate_output</a></div><div class="ttdeci">virtual void calculate_output(const libMesh::NumericVector< Real > &X, MAST::OutputAssemblyElemOperations &output)</div><div class="ttdoc">calculates the value of quantity . </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00254">assembly_base.cpp:254</a></div></div>
<div class="ttc" id="mast__data__types_8h_html_a579078c9edfe62906476a5e85dea6361"><div class="ttname"><a href="mast__data__types_8h.html#a579078c9edfe62906476a5e85dea6361">DenseRealVector</a></div><div class="ttdeci">libMesh::DenseVector< Real > DenseRealVector</div><div class="ttdef"><b>Definition:</b> <a href="mast__data__types_8h_source.html#l00056">mast_data_types.h:56</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_a4e6d01c9391e881f98a6caeb9a14503b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#a4e6d01c9391e881f98a6caeb9a14503b">MAST::AssemblyElemOperations::clear_assembly</a></div><div class="ttdeci">virtual void clear_assembly()</div><div class="ttdoc">clears the assembly object </div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8cpp_source.html#l00078">assembly_elem_operation.cpp:78</a></div></div>
<div class="ttc" id="namespace_m_a_s_t_html_a659caf861e2d27b164e12aef83e5ac1e"><div class="ttname"><a href="namespace_m_a_s_t.html#a659caf861e2d27b164e12aef83e5ac1e">MAST::copy</a></div><div class="ttdeci">void copy(DenseRealMatrix &m1, const RealMatrixX &m2)</div><div class="ttdef"><b>Definition:</b> <a href="utility_8h_source.html#l00167">utility.h:167</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_mesh_field_function_html_a45a0f5907f61898b83d649050617357d"><div class="ttname"><a href="class_m_a_s_t_1_1_mesh_field_function.html#a45a0f5907f61898b83d649050617357d">MAST::MeshFieldFunction::clear</a></div><div class="ttdeci">void clear()</div><div class="ttdoc">clear the solution </div><div class="ttdef"><b>Definition:</b> <a href="mesh__field__function_8cpp_source.html#l00247">mesh_field_function.cpp:247</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a0c8b29a632eb1b975005798b3862ad60"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a0c8b29a632eb1b975005798b3862ad60">MAST::AssemblyBase::calculate_output_direct_sensitivity</a></div><div class="ttdeci">virtual void calculate_output_direct_sensitivity(const libMesh::NumericVector< Real > &X, const libMesh::NumericVector< Real > *dXdp, const MAST::FunctionBase &p, MAST::OutputAssemblyElemOperations &output)</div><div class="ttdoc">evaluates the sensitivity of the outputs in the attached discipline with respect to the parametrs in ...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00412">assembly_base.cpp:412</a></div></div>
<div class="ttc" id="mast__data__types_8h_html_abc8ad2bc4221d7400e5306159e28e588"><div class="ttname"><a href="mast__data__types_8h.html#abc8ad2bc4221d7400e5306159e28e588">RealVectorX</a></div><div class="ttdeci">Matrix< Real, Dynamic, 1 > RealVectorX</div><div class="ttdef"><b>Definition:</b> <a href="mast__data__types_8h_source.html#l00035">mast_data_types.h:35</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_a918765a6fc86ba98c02f84bc8ce8367b"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a918765a6fc86ba98c02f84bc8ce8367b">MAST::OutputAssemblyElemOperations::output_sensitivity_total</a></div><div class="ttdeci">virtual Real output_sensitivity_total(const MAST::FunctionBase &p)=0</div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_a73b35bea5577b08a6a8ee70b19dac59c"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#a73b35bea5577b08a6a8ee70b19dac59c">MAST::AssemblyElemOperations::init</a></div><div class="ttdeci">virtual void init(const MAST::GeomElem &elem)=0</div><div class="ttdoc">initializes the object for calculation of element quantities for the specified elem. </div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_a16a069336ddd4ccc3b71bf737b555483"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a16a069336ddd4ccc3b71bf737b555483">MAST::OutputAssemblyElemOperations::output_derivative_for_elem</a></div><div class="ttdeci">virtual void output_derivative_for_elem(RealVectorX &dq_dX)=0</div><div class="ttdoc">returns the output quantity derivative with respect to state vector in dq_dX. </div></div>
<div class="ttc" id="assembly__base_8h_html"><div class="ttname"><a href="assembly__base_8h.html">assembly_base.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_af8281c5eece58746b64293f8ccea49f3"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#af8281c5eece58746b64293f8ccea49f3">MAST::AssemblyBase::clear_elem_parameter_dependence_object</a></div><div class="ttdeci">void clear_elem_parameter_dependence_object()</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00139">assembly_base.cpp:139</a></div></div>
<div class="ttc" id="mesh__field__function_8h_html"><div class="ttname"><a href="mesh__field__function_8h.html">mesh_field_function.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_ad37944186c61ab97dfc315d53f873e60"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#ad37944186c61ab97dfc315d53f873e60">MAST::AssemblyBase::build_localized_vector</a></div><div class="ttdeci">std::unique_ptr< libMesh::NumericVector< Real > > build_localized_vector(const libMesh::System &sys, const libMesh::NumericVector< Real > &global) const </div><div class="ttdoc">localizes the parallel vector so that the local copy stores all values necessary for calculation of t...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00211">assembly_base.cpp:211</a></div></div>
<div class="ttc" id="utility_8h_html"><div class="ttname"><a href="utility_8h.html">utility.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a652097c4673958f03ff6cca8cf037847"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a652097c4673958f03ff6cca8cf037847">MAST::AssemblyBase::get_solver_monitor</a></div><div class="ttdeci">MAST::AssemblyBase::SolverMonitor * get_solver_monitor()</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00147">assembly_base.cpp:147</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_geom_elem_html"><div class="ttname"><a href="class_m_a_s_t_1_1_geom_elem.html">MAST::GeomElem</a></div><div class="ttdoc">This class acts as a wrapper around libMesh::Elem for the purpose of providing a uniform interface fo...</div><div class="ttdef"><b>Definition:</b> <a href="geom__elem_8h_source.html#l00059">geom_elem.h:59</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence_html"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base_1_1_elem_parameter_dependence.html">MAST::AssemblyBase::ElemParameterDependence</a></div><div class="ttdoc">Inherited objects from this class can be provided by the user provide assessment of whether or not an...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00086">assembly_base.h:86</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a12fa42f74de999be67da555082a32c5d"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a12fa42f74de999be67da555082a32c5d">MAST::AssemblyBase::system_init</a></div><div class="ttdeci">MAST::SystemInitialization & system_init()</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00110">assembly_base.cpp:110</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_a6a81c128526ae2dff490652a39d60b11"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#a6a81c128526ae2dff490652a39d60b11">MAST::OutputAssemblyElemOperations::zero_for_analysis</a></div><div class="ttdeci">virtual void zero_for_analysis()=0</div><div class="ttdoc">zeroes the output quantity values stored inside this object so that assembly process can begin...</div></div>
<div class="ttc" id="class_m_a_s_t_1_1_mesh_field_function_html_ad05c1a6c795fbee828c4a3daae01fecd"><div class="ttname"><a href="class_m_a_s_t_1_1_mesh_field_function.html#ad05c1a6c795fbee828c4a3daae01fecd">MAST::MeshFieldFunction::init</a></div><div class="ttdeci">void init(const libMesh::NumericVector< Real > &sol, const libMesh::NumericVector< Real > *dsol=nullptr)</div><div class="ttdoc">initializes the data structures to perform the interpolation function of sol. </div><div class="ttdef"><b>Definition:</b> <a href="mesh__field__function_8cpp_source.html#l00188">mesh_field_function.cpp:188</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_aca93947fd87f364b33d7944deb98b99b"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#aca93947fd87f364b33d7944deb98b99b">MAST::AssemblyBase::_param_dependence</a></div><div class="ttdeci">MAST::AssemblyBase::ElemParameterDependence * _param_dependence</div><div class="ttdoc">If provided by user, this object is used by sensitiivty analysis to check for whether or the current ...</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00336">assembly_base.h:336</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a0644e12cd53adc64d9b7e027becc08b5"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a0644e12cd53adc64d9b7e027becc08b5">MAST::AssemblyBase::set_discipline_and_system</a></div><div class="ttdeci">virtual void set_discipline_and_system(MAST::PhysicsDisciplineBase &discipline, MAST::SystemInitialization &system)</div><div class="ttdoc">attaches a system to this discipline </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00162">assembly_base.cpp:162</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_elem_operations_html_aad8893093e07905b7e5071468b23514c"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_elem_operations.html#aad8893093e07905b7e5071468b23514c">MAST::AssemblyElemOperations::clear_elem</a></div><div class="ttdeci">virtual void clear_elem()</div><div class="ttdoc">clears the element initialization </div><div class="ttdef"><b>Definition:</b> <a href="assembly__elem__operation_8cpp_source.html#l00191">assembly_elem_operation.cpp:191</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a2034db25d6074eef21899ab93c66a833"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a2034db25d6074eef21899ab93c66a833">MAST::AssemblyBase::~AssemblyBase</a></div><div class="ttdeci">virtual ~AssemblyBase()</div><div class="ttdoc">virtual destructor </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00053">assembly_base.cpp:53</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a0b4114f497996001f1e54f452a926786"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a0b4114f497996001f1e54f452a926786">MAST::AssemblyBase::get_elem_ops</a></div><div class="ttdeci">MAST::AssemblyElemOperations & get_elem_ops()</div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00101">assembly_base.cpp:101</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_abb647f77c49b15b60250bf1a5081dbe4"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#abb647f77c49b15b60250bf1a5081dbe4">MAST::AssemblyBase::_solver_monitor</a></div><div class="ttdeci">MAST::AssemblyBase::SolverMonitor * _solver_monitor</div><div class="ttdoc">User provided solver monitor is attached to the linear nonlinear solvers, if provided. </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00329">assembly_base.h:329</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_abdff8292b4929a3aa7276c639a8d05e3"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">MAST::AssemblyBase::system</a></div><div class="ttdeci">const MAST::NonlinearSystem & system() const </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00083">assembly_base.cpp:83</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_output_assembly_elem_operations_html_ae895cc8daf0b4d118ea0db9a7b036004"><div class="ttname"><a href="class_m_a_s_t_1_1_output_assembly_elem_operations.html#ae895cc8daf0b4d118ea0db9a7b036004">MAST::OutputAssemblyElemOperations::evaluate_sensitivity</a></div><div class="ttdeci">virtual void evaluate_sensitivity(const MAST::FunctionBase &f)=0</div><div class="ttdoc">this evaluates all relevant sensitivity components on the element. </div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a8f878daffaae4a0eb77db2f3de1d45dc"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a8f878daffaae4a0eb77db2f3de1d45dc">MAST::AssemblyBase::_sol_function</a></div><div class="ttdeci">MAST::MeshFieldFunction * _sol_function</div><div class="ttdoc">system solution that will be initialized before each solution </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8h_source.html#l00323">assembly_base.h:323</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_assembly_base_html_a447c284ed86a61935a8498a5fee23462"><div class="ttname"><a href="class_m_a_s_t_1_1_assembly_base.html#a447c284ed86a61935a8498a5fee23462">MAST::AssemblyBase::discipline</a></div><div class="ttdeci">const MAST::PhysicsDisciplineBase & discipline() const </div><div class="ttdef"><b>Definition:</b> <a href="assembly__base_8cpp_source.html#l00062">assembly_base.cpp:62</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_0af1587c8378955de40f48b4bd1869f0.html">base</a></li><li class="navelem"><a class="el" href="assembly__base_8cpp.html">assembly_base.cpp</a></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.11 </li>
</ul>
</div>
</body>
</html>