-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy patharclength__continuation__solver_8cpp_source.html
161 lines (159 loc) · 60.2 KB
/
arclength__continuation__solver_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
<!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/solver/arclength_continuation_solver.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<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('arclength__continuation__solver_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">arclength_continuation_solver.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="arclength__continuation__solver_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> </div><div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment">// MAST includes</span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor">#include "<a class="code" href="arclength__continuation__solver_8h.html">solver/arclength_continuation_solver.h</a>"</span></div><div class="line"><a name="l00023"></a><span class="lineno"> 23</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="l00024"></a><span class="lineno"> 24</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="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#include "<a class="code" href="parameter_8h.html">base/parameter.h</a>"</span></div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> </div><div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment">// libMesh includes</span></div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor">#include "libmesh/linear_solver.h"</span></div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="preprocessor">#include "libmesh/dof_map.h"</span></div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> </div><div class="line"><a name="l00031"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#af1d214c7ee6365b93d69a13ec77f0479"> 31</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#af1d214c7ee6365b93d69a13ec77f0479">MAST::ArclengthContinuationSolver::ArclengthContinuationSolver</a>():</div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <a class="code" href="namespace_m_a_s_t.html">MAST</a>::<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html">ContinuationSolverBase</a>(),</div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span> _dpds_sign (1.) {</div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  </div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> }</div><div class="line"><a name="l00036"></a><span class="lineno"> 36</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"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#adbf6af0ff37d94cace5d33968d886877"> 38</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#adbf6af0ff37d94cace5d33968d886877">MAST::ArclengthContinuationSolver::~ArclengthContinuationSolver</a>() {</div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  </div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span> }</div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span> </div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span> </div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00044"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a1a9b89d1c18d2ead43401d1af97daf3b"> 44</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a1a9b89d1c18d2ead43401d1af97daf3b">MAST::ArclengthContinuationSolver::initialize</a>(<a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a> dp) {</div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  </div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  libmesh_assert(!<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a261972057673fdb59cd8e2cdec85bccd">_initialized</a>);</div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a69af9ff1784113826021d9e6a47fa4e8">_elem_ops</a>);</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>  (*_p)() += dp;</div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="keywordflow">if</span> (dp < 0.) <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a257bfb32ae3cb87ef8994d836ce74247">_dpds_sign</a> = -1.;</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>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a>&</div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  system = <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</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>  std::unique_ptr<libMesh::NumericVector<Real>></div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  dX(system.solution->clone().release());</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>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</a>().<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#aaaed1f33d18591f0f380516e1e2fbb8c">solve</a>(*<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a69af9ff1784113826021d9e6a47fa4e8">_elem_ops</a>, *<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>);</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>  dX->add(-1., *system.solution);</div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  dX->scale(-1.);</div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  dX->close();</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>  <span class="comment">// initialize scaling factors</span></div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a> = 1./dX->l2_norm();</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a> = 1./std::fabs(dp);</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>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ad85b6542eb3b6d922e4e32ec771578f2">arc_length</a> = std::sqrt(std::pow(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a>, 2) * dp*dp +</div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  std::pow(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a>, 2) * std::pow(dX->l2_norm(), 2));</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>  libmesh_assert_greater(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ad85b6542eb3b6d922e4e32ec771578f2">arc_length</a>, 0.);</div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  </div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a261972057673fdb59cd8e2cdec85bccd">_initialized</a> = <span class="keyword">true</span>;</div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span> }</div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span> </div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span> </div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a13c3ace24f2b9838235eb58705abf5de">MAST::ArclengthContinuationSolver::</a></div><div class="line"><a name="l00079"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a13c3ace24f2b9838235eb58705abf5de"> 79</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a13c3ace24f2b9838235eb58705abf5de">_solve_NR_iterate</a>(libMesh::NumericVector<Real> &X,</div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <a class="code" href="class_m_a_s_t_1_1_parameter.html">MAST::Parameter</a> &p) {</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>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a261972057673fdb59cd8e2cdec85bccd">_initialized</a>);</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span> </div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  std::unique_ptr<libMesh::NumericVector<Real>></div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  f(X.zero_clone().release()),</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  dgdX(X.zero_clone().release()),</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  dfdp(X.zero_clone().release()),</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  dXdp(X.zero_clone().release()),</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  dX(X.zero_clone().release());</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>  libMesh::SparseMatrix<Real></div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  &jac = *<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</a>().matrix;</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>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  g = 0.,</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  dgdp = 0.,</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  dp = 0.;</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>  <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42">_g</a>(X, p, *dfdp, *dXdp, g, dgdp, dgdX.get());</div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  </div><div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="keywordflow">if</span> (!<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a81321bdc20b5e491976904f669f00ee8">schur_factorization</a>)</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ac8c15f4d9ee862883e0c2d2b3e32ab79">_solve</a>(X, p,</div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  *f, <span class="keyword">true</span>, <span class="comment">// update f</span></div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  *dfdp, <span class="keyword">false</span>, <span class="comment">// update dfdp</span></div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  *dgdX, dgdp, g,</div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  *dX, dp);</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  <span class="keywordflow">else</span></div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a0fcf588b7f70edca8efeb5ca0461b0f9">_solve_schur_factorization</a>(X, p,</div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  jac, <span class="keyword">true</span>, <span class="comment">// update jac</span></div><div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  *f, <span class="keyword">true</span>, <span class="comment">// update f</span></div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  *dfdp, <span class="keyword">false</span>, <span class="comment">// update dfdp</span></div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  *dXdp, <span class="keyword">false</span>, <span class="comment">// update dXdp</span></div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  *dgdX, dgdp, g,</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  *dX, dp);</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>  <span class="comment">// update the solution and load parameter</span></div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  p() += dp;</div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  X.add(1., *dX);</div><div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  X.close();</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> </div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span> </div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ace5a4e6bd5814c11b3da802c50e0cad0">MAST::ArclengthContinuationSolver::</a></div><div class="line"><a name="l00125"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ace5a4e6bd5814c11b3da802c50e0cad0"> 125</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ace5a4e6bd5814c11b3da802c50e0cad0">_dXdp</a>(<span class="keyword">const</span> libMesh::NumericVector<Real> &X,</div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_parameter.html">MAST::Parameter</a> &p,</div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  libMesh::NumericVector<Real> &dfdp,</div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  libMesh::NumericVector<Real> &dXdp) {</div><div class="line"><a name="l00129"></a><span class="lineno"> 129</span> </div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a69af9ff1784113826021d9e6a47fa4e8">_elem_ops</a>);</div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  </div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html">MAST::NonlinearSystem</a></div><div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  &system = <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#abdff8292b4929a3aa7276c639a8d05e3">system</a>();</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>  libmesh_assert(system.<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a0cc6d230801fe6ca4d257d39e4a62fe4">operation</a>() == <a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a8efc4c9cce173ac724bd1a8b1dd1bf9fa9e013a2b2c6719553de9e894a237f5df">MAST::NonlinearSystem::NONE</a>);</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>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#adcc03f1662807d86d866ee3a4e2b60df">set_elem_operation_object</a>(*<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a69af9ff1784113826021d9e6a47fa4e8">_elem_ops</a>);</div><div class="line"><a name="l00138"></a><span class="lineno"> 138</span> </div><div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  system.<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a0dc9f61dc43a4f983eb9a73627e1832a">set_operation</a>(<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a8efc4c9cce173ac724bd1a8b1dd1bf9fab8a38b72e840ee8536f72f3f63eaa910">MAST::NonlinearSystem::FORWARD_SENSITIVITY_SOLVE</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_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><a class="code" href="class_m_a_s_t_1_1_assembly_base.html#a4eb0ce2595ad5f7b9279825fba621d79">sensitivity_assemble</a>(p, dfdp);</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>  libMesh::SparseMatrix<Real></div><div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  *pc = system.request_matrix(<span class="stringliteral">"Preconditioner"</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>  std::pair<unsigned int, Real></div><div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  solver_params = system.<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a29f8ecb7fe86ca26e7aa4bdfbafd400b">get_linear_solve_parameters</a>(),</div><div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  rval;</div><div class="line"><a name="l00149"></a><span class="lineno"> 149</span> </div><div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  dXdp.zero();</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>  rval = system.<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a4f4f650179ea9c8de5fc16edc3509290">linear_solver</a>->solve (*system.matrix, pc,</div><div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  dXdp,</div><div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  dfdp,</div><div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  solver_params.second,</div><div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  solver_params.first);</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>  dXdp.scale(-1.);</div><div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  dXdp.close();</div><div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  </div><div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  <span class="comment">// The linear solver may not have fit our constraints exactly</span></div><div class="line"><a name="l00162"></a><span class="lineno"> 162</span> <span class="preprocessor">#ifdef LIBMESH_ENABLE_CONSTRAINTS</span></div><div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  system.get_dof_map().enforce_constraints_exactly (system,</div><div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  &dXdp,</div><div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="comment">/* homogeneous = */</span> <span class="keyword">true</span>);</div><div class="line"><a name="l00166"></a><span class="lineno"> 166</span> <span class="preprocessor">#endif</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_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">_assembly</a>-><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="l00169"></a><span class="lineno"> 169</span>  system.<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a0dc9f61dc43a4f983eb9a73627e1832a">set_operation</a>(<a class="code" href="class_m_a_s_t_1_1_nonlinear_system.html#a8efc4c9cce173ac724bd1a8b1dd1bf9fa9e013a2b2c6719553de9e894a237f5df">MAST::NonlinearSystem::NONE</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> </div><div class="line"><a name="l00175"></a><span class="lineno"> 175</span> <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00176"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42"> 176</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42">MAST::ArclengthContinuationSolver::_g</a>(<span class="keyword">const</span> libMesh::NumericVector<Real> &X,</div><div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_parameter.html">MAST::Parameter</a> &p) {</div><div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  </div><div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  std::unique_ptr<libMesh::NumericVector<Real>></div><div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  dfdp(X.zero_clone().release()),</div><div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  dXdp(X.zero_clone().release());</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>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  g = 0.,</div><div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  dgdp = 0.;</div><div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  </div><div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42">_g</a>(X, p, *dfdp, *dXdp, g, dgdp, <span class="keyword">nullptr</span>);</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>  <span class="keywordflow">return</span> g;</div><div class="line"><a name="l00190"></a><span class="lineno"> 190</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> </div><div class="line"><a name="l00193"></a><span class="lineno"> 193</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00194"></a><span class="lineno"><a class="line" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ab1ccea1201cd8e98a83e5b7b996a3004"> 194</a></span> <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42">MAST::ArclengthContinuationSolver::_g</a>(<span class="keyword">const</span> libMesh::NumericVector<Real> &X,</div><div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="keyword">const</span> <a class="code" href="class_m_a_s_t_1_1_parameter.html">MAST::Parameter</a> &p,</div><div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  libMesh::NumericVector<Real> &dfdp,</div><div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  libMesh::NumericVector<Real> &dXdp,</div><div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a> &g,</div><div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a> &dgdp,</div><div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  libMesh::NumericVector<Real> *dgdX) {</div><div class="line"><a name="l00201"></a><span class="lineno"> 201</span> </div><div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  libmesh_assert(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a261972057673fdb59cd8e2cdec85bccd">_initialized</a>);</div><div class="line"><a name="l00203"></a><span class="lineno"> 203</span> </div><div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="comment">// update the constraint data</span></div><div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ace5a4e6bd5814c11b3da802c50e0cad0">_dXdp</a>(X, p, dfdp, dXdp);</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>  <span class="comment">// this includes scaling of X and p</span></div><div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <a class="code" href="mast__data__types_8h.html#a0fc21c6ff882da286a44d42c85f5378f">Real</a></div><div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  dpds = <a class="code" href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a257bfb32ae3cb87ef8994d836ce74247">_dpds_sign</a> * std::sqrt(1./ ( std::pow(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a>/<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a>,2) * dXdp.dot(dXdp) + 1.));</div><div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  </div><div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  std::unique_ptr<libMesh::NumericVector<Real>></div><div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  dX(X.clone().release());</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>  dX->add(-1., *<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a64931a46675460f49be09939b715861d">_X0</a>);</div><div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  dX->close();</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="comment">// (dX/ds)_scaled = (dX/dp)_scaled * (dp/ds)_scaled</span></div><div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  <span class="comment">// = (dX_scaled/dX) * dX/dp * (dp/dp_scaled) * (dp/ds)_scaled</span></div><div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  <span class="comment">// = X_scale/p_scale * dX/dp * (dp/ds)_scaled</span></div><div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  g = (<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a> * (dX->dot(dXdp) * dpds * <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a>/<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a>) +</div><div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a> * (p() - <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a1b887e0b6d743fc4f0f29b4e6789fd55">_p0</a>) * dpds) - <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ad85b6542eb3b6d922e4e32ec771578f2">arc_length</a>;</div><div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  dgdp = <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a> * dpds;</div><div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  </div><div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  <span class="keywordflow">if</span> (dgdX) {</div><div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  </div><div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  dgdX->zero();</div><div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  dgdX->add(<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a> * (dpds * <a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">_X_scale</a>/<a class="code" href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">_p_scale</a>), dXdp);</div><div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  dgdX->close();</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="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_af1d214c7ee6365b93d69a13ec77f0479"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#af1d214c7ee6365b93d69a13ec77f0479">MAST::ArclengthContinuationSolver::ArclengthContinuationSolver</a></div><div class="ttdeci">ArclengthContinuationSolver()</div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00031">arclength_continuation_solver.cpp:31</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_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_continuation_solver_base_html_a1b887e0b6d743fc4f0f29b4e6789fd55"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a1b887e0b6d743fc4f0f29b4e6789fd55">MAST::ContinuationSolverBase::_p0</a></div><div class="ttdeci">Real _p0</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00228">continuation_solver_base.h:228</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_a13c3ace24f2b9838235eb58705abf5de"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a13c3ace24f2b9838235eb58705abf5de">MAST::ArclengthContinuationSolver::_solve_NR_iterate</a></div><div class="ttdeci">virtual void _solve_NR_iterate(libMesh::NumericVector< Real > &X, MAST::Parameter &p)</div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00079">arclength_continuation_solver.cpp:79</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_parameter_html"><div class="ttname"><a href="class_m_a_s_t_1_1_parameter.html">MAST::Parameter</a></div><div class="ttdoc">This is a scalar function whose value can be changed and one that can be used as a design variable in...</div><div class="ttdef"><b>Definition:</b> <a href="parameter_8h_source.html#l00035">parameter.h:35</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="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="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_continuation_solver_base_html_ab2cc153ec1c0d9c4c725acf1a4a3f7bd"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#ab2cc153ec1c0d9c4c725acf1a4a3f7bd">MAST::ContinuationSolverBase::_p_scale</a></div><div class="ttdeci">Real _p_scale</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00228">continuation_solver_base.h:228</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_a0fcf588b7f70edca8efeb5ca0461b0f9"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a0fcf588b7f70edca8efeb5ca0461b0f9">MAST::ContinuationSolverBase::_solve_schur_factorization</a></div><div class="ttdeci">void _solve_schur_factorization(const libMesh::NumericVector< Real > &X, const MAST::Parameter &p, libMesh::SparseMatrix< Real > &jac, bool update_jac, libMesh::NumericVector< Real > &f, bool update_f, libMesh::NumericVector< Real > &dfdp, bool update_dfdp, libMesh::NumericVector< Real > &dXdp, bool update_dXdp, const libMesh::NumericVector< Real > &dgdX, const Real dgdp, const Real g, libMesh::NumericVector< Real > &dX, Real &dp)</div><div class="ttdoc">solves for the linear system of equation using Schur factorization. </div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8cpp_source.html#l00421">continuation_solver_base.cpp:421</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a8efc4c9cce173ac724bd1a8b1dd1bf9fab8a38b72e840ee8536f72f3f63eaa910"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a8efc4c9cce173ac724bd1a8b1dd1bf9fab8a38b72e840ee8536f72f3f63eaa910">MAST::NonlinearSystem::FORWARD_SENSITIVITY_SOLVE</a></div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00075">nonlinear_system.h:75</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_adbf6af0ff37d94cace5d33968d886877"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#adbf6af0ff37d94cace5d33968d886877">MAST::ArclengthContinuationSolver::~ArclengthContinuationSolver</a></div><div class="ttdeci">virtual ~ArclengthContinuationSolver()</div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00038">arclength_continuation_solver.cpp:38</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_ac8c15f4d9ee862883e0c2d2b3e32ab79"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#ac8c15f4d9ee862883e0c2d2b3e32ab79">MAST::ContinuationSolverBase::_solve</a></div><div class="ttdeci">void _solve(const libMesh::NumericVector< Real > &X, const MAST::Parameter &p, libMesh::NumericVector< Real > &f, bool update_f, libMesh::NumericVector< Real > &dfdp, bool update_dfdp, const libMesh::NumericVector< Real > &dgdX, const Real dgdp, const Real g, libMesh::NumericVector< Real > &dX, Real &dp)</div><div class="ttdoc">solves for the linear system of equation as a monolithic system dX and dp are returned from the solu...</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8cpp_source.html#l00184">continuation_solver_base.cpp:184</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html">MAST::ContinuationSolverBase</a></div><div class="ttdoc">the equation set is: the N-R updates are calculated such that This equation is solved using Schur-f...</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00053">continuation_solver_base.h:53</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a29f8ecb7fe86ca26e7aa4bdfbafd400b"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a29f8ecb7fe86ca26e7aa4bdfbafd400b">MAST::NonlinearSystem::get_linear_solve_parameters</a></div><div class="ttdeci">virtual std::pair< unsigned int, Real > get_linear_solve_parameters()</div><div class="ttdoc">calls NonlinearImplicitSystem::set_solver_parameters() before accessing the values. </div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8cpp_source.html#l00201">nonlinear_system.cpp:201</a></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_continuation_solver_base_html_ad85b6542eb3b6d922e4e32ec771578f2"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#ad85b6542eb3b6d922e4e32ec771578f2">MAST::ContinuationSolverBase::arc_length</a></div><div class="ttdeci">Real arc_length</div><div class="ttdoc">arc length that the solver is required to satisfy for the update. </div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00109">continuation_solver_base.h:109</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_a257bfb32ae3cb87ef8994d836ce74247"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a257bfb32ae3cb87ef8994d836ce74247">MAST::ArclengthContinuationSolver::_dpds_sign</a></div><div class="ttdeci">Real _dpds_sign</div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8h_source.html#l00104">arclength_continuation_solver.h:104</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_aa426e3c559573ff4fd95eb37592a63ca"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#aa426e3c559573ff4fd95eb37592a63ca">MAST::ContinuationSolverBase::_assembly</a></div><div class="ttdeci">MAST::AssemblyBase * _assembly</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00224">continuation_solver_base.h:224</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_a7191c6cd76fa7972b7227f2de6849ac6"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a7191c6cd76fa7972b7227f2de6849ac6">MAST::ContinuationSolverBase::_X_scale</a></div><div class="ttdeci">Real _X_scale</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00228">continuation_solver_base.h:228</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_a64931a46675460f49be09939b715861d"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a64931a46675460f49be09939b715861d">MAST::ContinuationSolverBase::_X0</a></div><div class="ttdeci">std::unique_ptr< libMesh::NumericVector< Real > > _X0</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00233">continuation_solver_base.h:233</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a0dc9f61dc43a4f983eb9a73627e1832a"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a0dc9f61dc43a4f983eb9a73627e1832a">MAST::NonlinearSystem::set_operation</a></div><div class="ttdeci">void set_operation(MAST::NonlinearSystem::Operation op)</div><div class="ttdoc">sets the current operation of the system </div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00093">nonlinear_system.h:93</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a4f4f650179ea9c8de5fc16edc3509290"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a4f4f650179ea9c8de5fc16edc3509290">MAST::NonlinearSystem::linear_solver</a></div><div class="ttdeci">std::unique_ptr< libMesh::LinearSolver< Real > > linear_solver</div><div class="ttdoc">The LinearSolver for solution of the linear equations. </div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00280">nonlinear_system.h:280</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_ae7881a0930358642ce2c3e47cc8aaf42"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ae7881a0930358642ce2c3e47cc8aaf42">MAST::ArclengthContinuationSolver::_g</a></div><div class="ttdeci">virtual Real _g(const libMesh::NumericVector< Real > &X, const MAST::Parameter &p)</div><div class="ttdoc"> </div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00176">arclength_continuation_solver.cpp:176</a></div></div>
<div class="ttc" id="arclength__continuation__solver_8h_html"><div class="ttname"><a href="arclength__continuation__solver_8h.html">arclength_continuation_solver.h</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a0cc6d230801fe6ca4d257d39e4a62fe4"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a0cc6d230801fe6ca4d257d39e4a62fe4">MAST::NonlinearSystem::operation</a></div><div class="ttdeci">MAST::NonlinearSystem::Operation operation()</div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00084">nonlinear_system.h:84</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_a69af9ff1784113826021d9e6a47fa4e8"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a69af9ff1784113826021d9e6a47fa4e8">MAST::ContinuationSolverBase::_elem_ops</a></div><div class="ttdeci">MAST::AssemblyElemOperations * _elem_ops</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00223">continuation_solver_base.h:223</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_continuation_solver_base_html_a261972057673fdb59cd8e2cdec85bccd"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a261972057673fdb59cd8e2cdec85bccd">MAST::ContinuationSolverBase::_initialized</a></div><div class="ttdeci">bool _initialized</div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00221">continuation_solver_base.h:221</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_a1a9b89d1c18d2ead43401d1af97daf3b"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#a1a9b89d1c18d2ead43401d1af97daf3b">MAST::ArclengthContinuationSolver::initialize</a></div><div class="ttdeci">virtual void initialize(Real dp)</div><div class="ttdoc">sets the arc length using a nonlinear solution using a step dp. </div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00044">arclength_continuation_solver.cpp:44</a></div></div>
<div class="ttc" id="parameter_8h_html"><div class="ttname"><a href="parameter_8h.html">parameter.h</a></div></div>
<div class="ttc" id="namespace_m_a_s_t_html"><div class="ttname"><a href="namespace_m_a_s_t.html">MAST</a></div><div class="ttdef"><b>Definition:</b> <a href="flutter__root__base_8h_source.html#l00027">flutter_root_base.h:27</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_continuation_solver_base_html_a81321bdc20b5e491976904f669f00ee8"><div class="ttname"><a href="class_m_a_s_t_1_1_continuation_solver_base.html#a81321bdc20b5e491976904f669f00ee8">MAST::ContinuationSolverBase::schur_factorization</a></div><div class="ttdeci">bool schur_factorization</div><div class="ttdoc">flag to use Schur-factorizaiton (default) or monolithic solver </div><div class="ttdef"><b>Definition:</b> <a href="continuation__solver__base_8h_source.html#l00141">continuation_solver_base.h:141</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_arclength_continuation_solver_html_ace5a4e6bd5814c11b3da802c50e0cad0"><div class="ttname"><a href="class_m_a_s_t_1_1_arclength_continuation_solver.html#ace5a4e6bd5814c11b3da802c50e0cad0">MAST::ArclengthContinuationSolver::_dXdp</a></div><div class="ttdeci">virtual void _dXdp(const libMesh::NumericVector< Real > &X, const MAST::Parameter &p, libMesh::NumericVector< Real > &dfdp, libMesh::NumericVector< Real > &dXdp)</div><div class="ttdef"><b>Definition:</b> <a href="arclength__continuation__solver_8cpp_source.html#l00125">arclength_continuation_solver.cpp:125</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_a8efc4c9cce173ac724bd1a8b1dd1bf9fa9e013a2b2c6719553de9e894a237f5df"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#a8efc4c9cce173ac724bd1a8b1dd1bf9fa9e013a2b2c6719553de9e894a237f5df">MAST::NonlinearSystem::NONE</a></div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8h_source.html#l00077">nonlinear_system.h:77</a></div></div>
<div class="ttc" id="class_m_a_s_t_1_1_nonlinear_system_html_aaaed1f33d18591f0f380516e1e2fbb8c"><div class="ttname"><a href="class_m_a_s_t_1_1_nonlinear_system.html#aaaed1f33d18591f0f380516e1e2fbb8c">MAST::NonlinearSystem::solve</a></div><div class="ttdeci">virtual void solve(MAST::AssemblyElemOperations &elem_ops, MAST::AssemblyBase &assembly)</div><div class="ttdoc">solves the nonlinear problem with the specified assembly operation object </div><div class="ttdef"><b>Definition:</b> <a href="nonlinear__system_8cpp_source.html#l00209">nonlinear_system.cpp:209</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_041db6ad6a3d44280e40e5af409429cb.html">solver</a></li><li class="navelem"><a class="el" href="arclength__continuation__solver_8cpp.html">arclength_continuation_solver.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>