-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathVecAsSymmMatrix_8h_source.html
333 lines (331 loc) · 35.3 KB
/
VecAsSymmMatrix_8h_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
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
<!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.9.8"/>
<title>G+Smo: /home/gitlab-runner/builds/t3_JWeu5_/0/gismo/gismo/src/gsMatrix/VecAsSymmMatrix.h Source File</title>
<link href="gismodoxy_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>
<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">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
$(document).ready(function() { init_search(); });
/* @license-end */
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
extensions: ["tex2jax.js"],
jax: ["input/TeX","output/HTML-CSS"],
});
</script>
<script type="text/javascript" async="async" src="https://people.ricam.oeaw.ac.at/gismo/mj/MathJax.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="gismodoxy.css" rel="stylesheet" type="text/css">
<!-- -->
<!-- <script type="text/javascript"> -->
<!-- </script> -->
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<!-- <a name="top"></a> -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectlogo"><img alt="Logo" src="gslogo-doxy.jpg"/></td>
<td style="padding-left: 0.5em;">
<div id="projectname"><a href="index.html"><font style="font-variant:small-caps;">G+S</font><font style="font-variant:small-caps;"
color="#000000">mo</font></a>
 <span id="projectnumber">25.01.0</span>
</div>
<div id="projectbrief">Geometry + Simulation Modules</div>
</td>
<td> <div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<span id="MSearchSelect" onmouseover="return searchBox.OnSearchSelectShow()" onmouseout="return searchBox.OnSearchSelectHide()"> </span>
<input type="text" id="MSearchField" value="" placeholder="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.svg" alt=""/></a>
</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.8 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
var searchBox = new SearchBox("searchBox", "search/",'.html');
/* @license-end */
</script>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
<div id="nav-tree">
<div id="nav-tree-contents">
<div id="nav-sync" class="sync"></div>
</div>
</div>
<div id="splitbar" style="-moz-user-select:none;"
class="ui-resizable-handle">
</div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
$(document).ready(function(){initNavTree('VecAsSymmMatrix_8h_source.html',''); initResizable(); });
/* @license-end */
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<div id="MSearchResults">
<div class="SRPage">
<div id="SRIndex">
<div id="SRResults"></div>
<div class="SRStatus" id="Loading">Loading...</div>
<div class="SRStatus" id="Searching">Searching...</div>
<div class="SRStatus" id="NoMatches">No Matches</div>
</div>
</div>
</div>
</div>
<div class="header">
<div class="headertitle"><div class="title">VecAsSymmMatrix.h</div></div>
</div><!--header-->
<div class="contents">
<a href="VecAsSymmMatrix_8h.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a id="l00001" name="l00001"></a><span class="lineno"> 1</span> </div>
<div class="line"><a id="l00015" name="l00015"></a><span class="lineno"> 15</span><span class="keyword">namespace </span><a class="code hl_namespaceRef" href="http://eigen.tuxfamily.org/dox/namespaceEigen.html">Eigen</a> { </div>
<div class="line"><a id="l00016" name="l00016"></a><span class="lineno"> 16</span> </div>
<div class="line"><a id="l00025" name="l00025"></a><span class="lineno"> 25</span><span class="keyword">namespace </span>internal {</div>
<div class="line"><a id="l00026" name="l00026"></a><span class="lineno"> 26</span><span class="keyword">template</span><<span class="keyword">typename</span> MatrixType,<span class="keywordtype">int</span> Dim></div>
<div class="line"><a id="l00027" name="l00027"></a><span class="lineno"> 27</span><span class="keyword">struct </span>traits<VecAsSymmMatrix<MatrixType,Dim> ></div>
<div class="line"><a id="l00028" name="l00028"></a><span class="lineno"> 28</span> : traits<MatrixType></div>
<div class="line"><a id="l00029" name="l00029"></a><span class="lineno"> 29</span>{</div>
<div class="line"><a id="l00030" name="l00030"></a><span class="lineno"> 30</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> MatrixType::Scalar Scalar;</div>
<div class="line"><a id="l00031" name="l00031"></a><span class="lineno"> 31</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> traits<MatrixType>::StorageKind StorageKind;</div>
<div class="line"><a id="l00032" name="l00032"></a><span class="lineno"> 32</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> traits<MatrixType>::XprKind XprKind;</div>
<div class="line"><a id="l00033" name="l00033"></a><span class="lineno"> 33</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> ref_selector<MatrixType>::type MatrixTypeNested;</div>
<div class="line"><a id="l00034" name="l00034"></a><span class="lineno"> 34</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> remove_reference<MatrixTypeNested>::type _MatrixTypeNested;</div>
<div class="line"><a id="l00035" name="l00035"></a><span class="lineno"> 35</span> <span class="keyword">enum</span> {</div>
<div class="line"><a id="l00036" name="l00036"></a><span class="lineno"> 36</span> RowsAtCompileTime = Dim==Dynamic || int(MatrixType::RowsAtCompileTime)==Dynamic</div>
<div class="line"><a id="l00037" name="l00037"></a><span class="lineno"> 37</span> ? Dynamic</div>
<div class="line"><a id="l00038" name="l00038"></a><span class="lineno"> 38</span> : Dim * MatrixType::RowsAtCompileTime,</div>
<div class="line"><a id="l00039" name="l00039"></a><span class="lineno"> 39</span> ColsAtCompileTime = Dim==Dynamic || int(MatrixType::ColsAtCompileTime)==Dynamic</div>
<div class="line"><a id="l00040" name="l00040"></a><span class="lineno"> 40</span> ? Dynamic</div>
<div class="line"><a id="l00041" name="l00041"></a><span class="lineno"> 41</span> : Dim * MatrixType::ColsAtCompileTime,</div>
<div class="line"><a id="l00042" name="l00042"></a><span class="lineno"> 42</span> <span class="comment">//FIXME we don't propagate the max sizes !!!</span></div>
<div class="line"><a id="l00043" name="l00043"></a><span class="lineno"> 43</span> MaxRowsAtCompileTime = RowsAtCompileTime,</div>
<div class="line"><a id="l00044" name="l00044"></a><span class="lineno"> 44</span> MaxColsAtCompileTime = ColsAtCompileTime,</div>
<div class="line"><a id="l00045" name="l00045"></a><span class="lineno"> 45</span> IsRowMajor = MaxRowsAtCompileTime==1 && MaxColsAtCompileTime!=1 ? 1</div>
<div class="line"><a id="l00046" name="l00046"></a><span class="lineno"> 46</span> : MaxColsAtCompileTime==1 && MaxRowsAtCompileTime!=1 ? 0</div>
<div class="line"><a id="l00047" name="l00047"></a><span class="lineno"> 47</span> : (MatrixType::Flags & RowMajorBit) ? 1 : 0,</div>
<div class="line"><a id="l00048" name="l00048"></a><span class="lineno"> 48</span> </div>
<div class="line"><a id="l00049" name="l00049"></a><span class="lineno"> 49</span> <span class="comment">// FIXME enable DirectAccess with negative strides?</span></div>
<div class="line"><a id="l00050" name="l00050"></a><span class="lineno"> 50</span> Flags = IsRowMajor ? RowMajorBit : 0</div>
<div class="line"><a id="l00051" name="l00051"></a><span class="lineno"> 51</span> };</div>
<div class="line"><a id="l00052" name="l00052"></a><span class="lineno"> 52</span>};</div>
<div class="line"><a id="l00053" name="l00053"></a><span class="lineno"> 53</span>}</div>
<div class="line"><a id="l00054" name="l00054"></a><span class="lineno"> 54</span> </div>
<div class="foldopen" id="foldopen00055" data-start="{" data-end="};">
<div class="line"><a id="l00055" name="l00055"></a><span class="lineno"><a class="line" href="classEigen_1_1VecAsSymmMatrix.html"> 55</a></span><span class="keyword">template</span><<span class="keyword">typename</span> MatrixType,<span class="keywordtype">int</span> Dim> <span class="keyword">class </span><a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix</a></div>
<div class="line"><a id="l00056" name="l00056"></a><span class="lineno"> 56</span> : <span class="keyword">public</span> internal::dense_xpr_base< VecAsSymmMatrix<MatrixType,Dim> >::type</div>
<div class="line"><a id="l00057" name="l00057"></a><span class="lineno"> 57</span>{</div>
<div class="line"><a id="l00058" name="l00058"></a><span class="lineno"> 58</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::traits<VecAsSymmMatrix>::MatrixTypeNested MatrixTypeNested;</div>
<div class="line"><a id="l00059" name="l00059"></a><span class="lineno"> 59</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::traits<VecAsSymmMatrix>::_MatrixTypeNested _MatrixTypeNested;</div>
<div class="line"><a id="l00060" name="l00060"></a><span class="lineno"> 60</span> <span class="keyword">public</span>:</div>
<div class="line"><a id="l00061" name="l00061"></a><span class="lineno"> 61</span> </div>
<div class="line"><a id="l00062" name="l00062"></a><span class="lineno"> 62</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::dense_xpr_base<VecAsSymmMatrix>::type Base;</div>
<div class="line"><a id="l00063" name="l00063"></a><span class="lineno"> 63</span> EIGEN_DENSE_PUBLIC_INTERFACE(<a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix</a>)</div>
<div class="line"><a id="l00064" name="l00064"></a><span class="lineno"> 64</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::remove_all<MatrixType>::type NestedExpression;</div>
<div class="line"><a id="l00065" name="l00065"></a><span class="lineno"> 65</span> </div>
<div class="line"><a id="l00066" name="l00066"></a><span class="lineno"> 66</span> <span class="keyword">template</span><<span class="keyword">typename</span> OriginalMatrixType></div>
<div class="line"><a id="l00067" name="l00067"></a><span class="lineno"> 67</span> EIGEN_DEVICE_FUNC</div>
<div class="line"><a id="l00068" name="l00068"></a><span class="lineno"> 68</span> <span class="keyword">inline</span> <span class="keyword">explicit</span> <a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix</a>(<span class="keyword">const</span> OriginalMatrixType& a_matrix)</div>
<div class="line"><a id="l00069" name="l00069"></a><span class="lineno"> 69</span> : m_matrix(a_matrix), m_mdim(Dim)</div>
<div class="line"><a id="l00070" name="l00070"></a><span class="lineno"> 70</span> {</div>
<div class="line"><a id="l00071" name="l00071"></a><span class="lineno"> 71</span> EIGEN_STATIC_ASSERT((internal::is_same<<span class="keyword">typename</span> internal::remove_const<MatrixType>::type,OriginalMatrixType>::value),</div>
<div class="line"><a id="l00072" name="l00072"></a><span class="lineno"> 72</span> THE_MATRIX_OR_EXPRESSION_THAT_YOU_PASSED_DOES_NOT_HAVE_THE_EXPECTED_TYPE)</div>
<div class="line"><a id="l00073" name="l00073"></a><span class="lineno"> 73</span> eigen_assert(Dim!=Dynamic);</div>
<div class="line"><a id="l00074" name="l00074"></a><span class="lineno"> 74</span> }</div>
<div class="line"><a id="l00075" name="l00075"></a><span class="lineno"> 75</span> </div>
<div class="line"><a id="l00076" name="l00076"></a><span class="lineno"> 76</span> <span class="keyword">template</span><<span class="keyword">typename</span> OriginalMatrixType></div>
<div class="line"><a id="l00077" name="l00077"></a><span class="lineno"> 77</span> EIGEN_DEVICE_FUNC</div>
<div class="line"><a id="l00078" name="l00078"></a><span class="lineno"> 78</span> <span class="keyword">inline</span> <a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix</a>(<span class="keyword">const</span> OriginalMatrixType& a_matrix, Index mdim)</div>
<div class="line"><a id="l00079" name="l00079"></a><span class="lineno"> 79</span> : m_matrix(a_matrix), m_mdim(mdim)</div>
<div class="line"><a id="l00080" name="l00080"></a><span class="lineno"> 80</span> {</div>
<div class="line"><a id="l00081" name="l00081"></a><span class="lineno"> 81</span> EIGEN_STATIC_ASSERT((internal::is_same<<span class="keyword">typename</span> internal::remove_const<MatrixType>::type,OriginalMatrixType>::value),</div>
<div class="line"><a id="l00082" name="l00082"></a><span class="lineno"> 82</span> THE_MATRIX_OR_EXPRESSION_THAT_YOU_PASSED_DOES_NOT_HAVE_THE_EXPECTED_TYPE)</div>
<div class="line"><a id="l00083" name="l00083"></a><span class="lineno"> 83</span> }</div>
<div class="line"><a id="l00084" name="l00084"></a><span class="lineno"> 84</span> </div>
<div class="line"><a id="l00085" name="l00085"></a><span class="lineno"> 85</span> EIGEN_DEVICE_FUNC</div>
<div class="line"><a id="l00086" name="l00086"></a><span class="lineno"> 86</span> <span class="keyword">inline</span> Index rows()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> m_mdim.value(); }</div>
<div class="line"><a id="l00087" name="l00087"></a><span class="lineno"> 87</span> </div>
<div class="line"><a id="l00088" name="l00088"></a><span class="lineno"> 88</span> EIGEN_DEVICE_FUNC</div>
<div class="line"><a id="l00089" name="l00089"></a><span class="lineno"> 89</span> <span class="keyword">inline</span> Index cols()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> m_mdim.value(); }</div>
<div class="line"><a id="l00090" name="l00090"></a><span class="lineno"> 90</span><span class="comment">/*</span></div>
<div class="line"><a id="l00091" name="l00091"></a><span class="lineno"> 91</span><span class="comment"> inline Scalar coeff(Index rowId, Index colId) const</span></div>
<div class="line"><a id="l00092" name="l00092"></a><span class="lineno"> 92</span><span class="comment"> {</span></div>
<div class="line"><a id="l00093" name="l00093"></a><span class="lineno"> 93</span><span class="comment"> if ( rowId / m_matrix.rows() != colId / m_matrix.cols() )</span></div>
<div class="line"><a id="l00094" name="l00094"></a><span class="lineno"> 94</span><span class="comment"> return Scalar(0);</span></div>
<div class="line"><a id="l00095" name="l00095"></a><span class="lineno"> 95</span><span class="comment"> // try to avoid using modulo; this is a pure optimization strategy</span></div>
<div class="line"><a id="l00096" name="l00096"></a><span class="lineno"> 96</span><span class="comment"> const Index actual_row = internal::traits<MatrixType>::RowsAtCompileTime==1 ? 0</span></div>
<div class="line"><a id="l00097" name="l00097"></a><span class="lineno"> 97</span><span class="comment"> : Dim==1 ? rowId</span></div>
<div class="line"><a id="l00098" name="l00098"></a><span class="lineno"> 98</span><span class="comment"> : rowId%m_matrix.rows();</span></div>
<div class="line"><a id="l00099" name="l00099"></a><span class="lineno"> 99</span><span class="comment"> const Index actual_col = internal::traits<MatrixType>::ColsAtCompileTime==1 ? 0</span></div>
<div class="line"><a id="l00100" name="l00100"></a><span class="lineno"> 100</span><span class="comment"> : Dim==1 ? colId</span></div>
<div class="line"><a id="l00101" name="l00101"></a><span class="lineno"> 101</span><span class="comment"> : colId%m_matrix.cols();</span></div>
<div class="line"><a id="l00102" name="l00102"></a><span class="lineno"> 102</span><span class="comment"></span> </div>
<div class="line"><a id="l00103" name="l00103"></a><span class="lineno"> 103</span><span class="comment"> return m_matrix.coeff(actual_row, actual_col);</span></div>
<div class="line"><a id="l00104" name="l00104"></a><span class="lineno"> 104</span><span class="comment"> }</span></div>
<div class="line"><a id="l00105" name="l00105"></a><span class="lineno"> 105</span><span class="comment"> template<int LoadMode></span></div>
<div class="line"><a id="l00106" name="l00106"></a><span class="lineno"> 106</span><span class="comment"> inline PacketScalar packet(Index rowId, Index colId) const</span></div>
<div class="line"><a id="l00107" name="l00107"></a><span class="lineno"> 107</span><span class="comment"> {</span></div>
<div class="line"><a id="l00108" name="l00108"></a><span class="lineno"> 108</span><span class="comment"> </span></div>
<div class="line"><a id="l00109" name="l00109"></a><span class="lineno"> 109</span><span class="comment"> if ( rowId / m_matrix.rows() != colId / m_matrix.cols() )</span></div>
<div class="line"><a id="l00110" name="l00110"></a><span class="lineno"> 110</span><span class="comment"> GISMO_ERROR("not implemented");</span></div>
<div class="line"><a id="l00111" name="l00111"></a><span class="lineno"> 111</span><span class="comment"> </span></div>
<div class="line"><a id="l00112" name="l00112"></a><span class="lineno"> 112</span><span class="comment"> const Index actual_row = internal::traits<MatrixType>::RowsAtCompileTime==1 ? 0</span></div>
<div class="line"><a id="l00113" name="l00113"></a><span class="lineno"> 113</span><span class="comment"> : Dim==1 ? rowId</span></div>
<div class="line"><a id="l00114" name="l00114"></a><span class="lineno"> 114</span><span class="comment"> : rowId%m_matrix.rows();</span></div>
<div class="line"><a id="l00115" name="l00115"></a><span class="lineno"> 115</span><span class="comment"> const Index actual_col = internal::traits<MatrixType>::ColsAtCompileTime==1 ? 0</span></div>
<div class="line"><a id="l00116" name="l00116"></a><span class="lineno"> 116</span><span class="comment"> : Dim==1 ? colId</span></div>
<div class="line"><a id="l00117" name="l00117"></a><span class="lineno"> 117</span><span class="comment"> : colId%m_matrix.cols();</span></div>
<div class="line"><a id="l00118" name="l00118"></a><span class="lineno"> 118</span><span class="comment"></span> </div>
<div class="line"><a id="l00119" name="l00119"></a><span class="lineno"> 119</span><span class="comment"> return m_matrix.template packet<LoadMode>(actual_row, actual_col);</span></div>
<div class="line"><a id="l00120" name="l00120"></a><span class="lineno"> 120</span><span class="comment"> }</span></div>
<div class="line"><a id="l00121" name="l00121"></a><span class="lineno"> 121</span><span class="comment">*/</span></div>
<div class="line"><a id="l00122" name="l00122"></a><span class="lineno"> 122</span> EIGEN_DEVICE_FUNC</div>
<div class="line"><a id="l00123" name="l00123"></a><span class="lineno"> 123</span> <span class="keyword">const</span> _MatrixTypeNested& nestedExpression()<span class="keyword"> const</span></div>
<div class="line"><a id="l00124" name="l00124"></a><span class="lineno"> 124</span><span class="keyword"> </span>{ </div>
<div class="line"><a id="l00125" name="l00125"></a><span class="lineno"> 125</span> <span class="keywordflow">return</span> m_matrix; </div>
<div class="line"><a id="l00126" name="l00126"></a><span class="lineno"> 126</span> }</div>
<div class="line"><a id="l00127" name="l00127"></a><span class="lineno"> 127</span> </div>
<div class="line"><a id="l00128" name="l00128"></a><span class="lineno"> 128</span> <span class="keyword">protected</span>:</div>
<div class="line"><a id="l00129" name="l00129"></a><span class="lineno"> 129</span> MatrixTypeNested m_matrix;</div>
<div class="line"><a id="l00130" name="l00130"></a><span class="lineno"> 130</span> <span class="keyword">const</span> internal::variable_if_dynamic<Index, Dim> m_mdim;</div>
<div class="line"><a id="l00131" name="l00131"></a><span class="lineno"> 131</span>};</div>
</div>
<div class="line"><a id="l00132" name="l00132"></a><span class="lineno"> 132</span> </div>
<div class="line"><a id="l00133" name="l00133"></a><span class="lineno"> 133</span><span class="keyword">namespace </span>internal {</div>
<div class="line"><a id="l00134" name="l00134"></a><span class="lineno"> 134</span><span class="keyword">template</span><<span class="keyword">typename</span> ArgType, <span class="keywordtype">int</span> Dim> </div>
<div class="line"><a id="l00135" name="l00135"></a><span class="lineno"> 135</span><span class="keyword">struct </span>unary_evaluator<<a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix</a><ArgType, Dim> ></div>
<div class="line"><a id="l00136" name="l00136"></a><span class="lineno"> 136</span> : evaluator_base<VecAsSymmMatrix<ArgType, Dim> ></div>
<div class="line"><a id="l00137" name="l00137"></a><span class="lineno"> 137</span>{</div>
<div class="line"><a id="l00138" name="l00138"></a><span class="lineno"> 138</span> <span class="keyword">typedef</span> <a class="code hl_class" href="classEigen_1_1VecAsSymmMatrix.html">VecAsSymmMatrix<ArgType, Dim></a> XprType;</div>
<div class="line"><a id="l00139" name="l00139"></a><span class="lineno"> 139</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> XprType::CoeffReturnType CoeffReturnType;</div>
<div class="line"><a id="l00140" name="l00140"></a><span class="lineno"> 140</span> <span class="keyword">enum</span> {</div>
<div class="line"><a id="l00141" name="l00141"></a><span class="lineno"> 141</span> Factor = (Dim==Dynamic) ? Dynamic : Dim*Dim</div>
<div class="line"><a id="l00142" name="l00142"></a><span class="lineno"> 142</span> };</div>
<div class="line"><a id="l00143" name="l00143"></a><span class="lineno"> 143</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::nested_eval<ArgType,Factor>::type ArgTypeNested;</div>
<div class="line"><a id="l00144" name="l00144"></a><span class="lineno"> 144</span> <span class="keyword">typedef</span> <span class="keyword">typename</span> internal::remove_all<ArgTypeNested>::type ArgTypeNestedCleaned;</div>
<div class="line"><a id="l00145" name="l00145"></a><span class="lineno"> 145</span> </div>
<div class="line"><a id="l00146" name="l00146"></a><span class="lineno"> 146</span> <span class="keyword">enum</span> {</div>
<div class="line"><a id="l00147" name="l00147"></a><span class="lineno"> 147</span> CoeffReadCost = evaluator<ArgTypeNestedCleaned>::CoeffReadCost,</div>
<div class="line"><a id="l00148" name="l00148"></a><span class="lineno"> 148</span> LinearAccessMask = XprType::IsVectorAtCompileTime ? LinearAccessBit : 0,</div>
<div class="line"><a id="l00149" name="l00149"></a><span class="lineno"> 149</span> Flags = (evaluator<ArgTypeNestedCleaned>::Flags & (HereditaryBits|LinearAccessMask) & ~RowMajorBit) | (traits<XprType>::Flags & RowMajorBit),</div>
<div class="line"><a id="l00150" name="l00150"></a><span class="lineno"> 150</span> </div>
<div class="line"><a id="l00151" name="l00151"></a><span class="lineno"> 151</span> Alignment = evaluator<ArgTypeNestedCleaned>::Alignment</div>
<div class="line"><a id="l00152" name="l00152"></a><span class="lineno"> 152</span> };</div>
<div class="line"><a id="l00153" name="l00153"></a><span class="lineno"> 153</span> </div>
<div class="line"><a id="l00154" name="l00154"></a><span class="lineno"> 154</span> EIGEN_DEVICE_FUNC <span class="keyword">explicit</span> unary_evaluator(<span class="keyword">const</span> XprType& xpr)</div>
<div class="line"><a id="l00155" name="l00155"></a><span class="lineno"> 155</span> : m_arg(xpr.nestedExpression()),</div>
<div class="line"><a id="l00156" name="l00156"></a><span class="lineno"> 156</span> m_argImpl(m_arg),</div>
<div class="line"><a id="l00157" name="l00157"></a><span class="lineno"> 157</span> m_rows(xpr.nestedExpression().rows()),</div>
<div class="line"><a id="l00158" name="l00158"></a><span class="lineno"> 158</span> m_cols(xpr.nestedExpression().cols())</div>
<div class="line"><a id="l00159" name="l00159"></a><span class="lineno"> 159</span> {}</div>
<div class="line"><a id="l00160" name="l00160"></a><span class="lineno"> 160</span> </div>
<div class="line"><a id="l00161" name="l00161"></a><span class="lineno"> 161</span> EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE</div>
<div class="line"><a id="l00162" name="l00162"></a><span class="lineno"> 162</span> CoeffReturnType coeff(Index rowId, Index colId)<span class="keyword"> const</span></div>
<div class="line"><a id="l00163" name="l00163"></a><span class="lineno"> 163</span><span class="keyword"> </span>{</div>
<div class="line"><a id="l00164" name="l00164"></a><span class="lineno"> 164</span> <span class="keywordflow">if</span> ( rowId / m_rows.value() != colId / m_cols.value() )</div>
<div class="line"><a id="l00165" name="l00165"></a><span class="lineno"> 165</span> <span class="keywordflow">return</span> CoeffReturnType(0);</div>
<div class="line"><a id="l00166" name="l00166"></a><span class="lineno"> 166</span> <span class="comment">// try to avoid using modulo; this is a pure optimization strategy</span></div>
<div class="line"><a id="l00167" name="l00167"></a><span class="lineno"> 167</span> <span class="keyword">const</span> Index actual_row = internal::traits<XprType>::RowsAtCompileTime==1 ? 0</div>
<div class="line"><a id="l00168" name="l00168"></a><span class="lineno"> 168</span> : Dim==1 ? rowId</div>
<div class="line"><a id="l00169" name="l00169"></a><span class="lineno"> 169</span> : rowId%m_rows.value();</div>
<div class="line"><a id="l00170" name="l00170"></a><span class="lineno"> 170</span> <span class="keyword">const</span> Index actual_col = internal::traits<XprType>::ColsAtCompileTime==1 ? 0</div>
<div class="line"><a id="l00171" name="l00171"></a><span class="lineno"> 171</span> : Dim==1 ? colId</div>
<div class="line"><a id="l00172" name="l00172"></a><span class="lineno"> 172</span> : colId%m_cols.value();</div>
<div class="line"><a id="l00173" name="l00173"></a><span class="lineno"> 173</span> </div>
<div class="line"><a id="l00174" name="l00174"></a><span class="lineno"> 174</span> <span class="keywordflow">return</span> m_argImpl.coeff(actual_row, actual_col);</div>
<div class="line"><a id="l00175" name="l00175"></a><span class="lineno"> 175</span> }</div>
<div class="line"><a id="l00176" name="l00176"></a><span class="lineno"> 176</span> </div>
<div class="line"><a id="l00177" name="l00177"></a><span class="lineno"> 177</span> <span class="keyword">template</span><<span class="keywordtype">int</span> LoadMode, <span class="keyword">typename</span> PacketType></div>
<div class="line"><a id="l00178" name="l00178"></a><span class="lineno"> 178</span> EIGEN_STRONG_INLINE</div>
<div class="line"><a id="l00179" name="l00179"></a><span class="lineno"> 179</span> PacketType packet(Index rowId, Index colId)<span class="keyword"> const</span></div>
<div class="line"><a id="l00180" name="l00180"></a><span class="lineno"> 180</span><span class="keyword"> </span>{</div>
<div class="line"><a id="l00181" name="l00181"></a><span class="lineno"> 181</span> assert( rowId / m_rows.value() != colId / m_cols.value() &&</div>
<div class="line"><a id="l00182" name="l00182"></a><span class="lineno"> 182</span> <span class="stringliteral">"Not implemented"</span>);</div>
<div class="line"><a id="l00183" name="l00183"></a><span class="lineno"> 183</span> </div>
<div class="line"><a id="l00184" name="l00184"></a><span class="lineno"> 184</span> <span class="keyword">const</span> Index actual_row = internal::traits<XprType>::RowsAtCompileTime==1 ? 0</div>
<div class="line"><a id="l00185" name="l00185"></a><span class="lineno"> 185</span> : Dim==1 ? rowId</div>
<div class="line"><a id="l00186" name="l00186"></a><span class="lineno"> 186</span> : rowId%m_rows.value();</div>
<div class="line"><a id="l00187" name="l00187"></a><span class="lineno"> 187</span> <span class="keyword">const</span> Index actual_col = internal::traits<XprType>::ColsAtCompileTime==1 ? 0</div>
<div class="line"><a id="l00188" name="l00188"></a><span class="lineno"> 188</span> : Dim==1 ? colId</div>
<div class="line"><a id="l00189" name="l00189"></a><span class="lineno"> 189</span> : colId%m_cols.value();</div>
<div class="line"><a id="l00190" name="l00190"></a><span class="lineno"> 190</span> </div>
<div class="line"><a id="l00191" name="l00191"></a><span class="lineno"> 191</span> <span class="keywordflow">return</span> m_argImpl.template packet<LoadMode>(actual_row, actual_col);</div>
<div class="line"><a id="l00192" name="l00192"></a><span class="lineno"> 192</span> }</div>
<div class="line"><a id="l00193" name="l00193"></a><span class="lineno"> 193</span> </div>
<div class="line"><a id="l00194" name="l00194"></a><span class="lineno"> 194</span><span class="keyword">protected</span>:</div>
<div class="line"><a id="l00195" name="l00195"></a><span class="lineno"> 195</span> ArgTypeNested m_arg;</div>
<div class="line"><a id="l00196" name="l00196"></a><span class="lineno"> 196</span> evaluator<ArgTypeNestedCleaned> m_argImpl;</div>
<div class="line"><a id="l00197" name="l00197"></a><span class="lineno"> 197</span> <span class="keyword">const</span> variable_if_dynamic<Index, ArgType::RowsAtCompileTime> m_rows;</div>
<div class="line"><a id="l00198" name="l00198"></a><span class="lineno"> 198</span> <span class="keyword">const</span> variable_if_dynamic<Index, ArgType::ColsAtCompileTime> m_cols;</div>
<div class="line"><a id="l00199" name="l00199"></a><span class="lineno"> 199</span>};</div>
<div class="line"><a id="l00200" name="l00200"></a><span class="lineno"> 200</span> </div>
<div class="line"><a id="l00201" name="l00201"></a><span class="lineno"> 201</span>} <span class="comment">// namespace internal</span></div>
<div class="line"><a id="l00202" name="l00202"></a><span class="lineno"> 202</span> </div>
<div class="line"><a id="l00211" name="l00211"></a><span class="lineno"> 211</span><span class="keyword">template</span><<span class="keyword">typename</span> Derived></div>
<div class="line"><a id="l00212" name="l00212"></a><span class="lineno"> 212</span><span class="keyword">const</span> <span class="keyword">typename</span> MatrixBase<Derived>::VecAsSymmMatrixReturnType</div>
<div class="line"><a id="l00213" name="l00213"></a><span class="lineno"> 213</span>MatrixBase<Derived>::blockDiag(Index mdim)<span class="keyword"> const</span></div>
<div class="line"><a id="l00214" name="l00214"></a><span class="lineno"> 214</span><span class="keyword"></span>{</div>
<div class="line"><a id="l00215" name="l00215"></a><span class="lineno"> 215</span> <span class="keywordflow">return</span> VecAsSymmMatrix<Derived,Dynamic>(derived(),mdim);</div>
<div class="line"><a id="l00216" name="l00216"></a><span class="lineno"> 216</span>}</div>
<div class="line"><a id="l00217" name="l00217"></a><span class="lineno"> 217</span> </div>
<div class="line"><a id="l00218" name="l00218"></a><span class="lineno"> 218</span>} <span class="comment">// namespace eigen</span></div>
<div class="ttc" id="aclassEigen_1_1VecAsSymmMatrix_html"><div class="ttname"><a href="classEigen_1_1VecAsSymmMatrix.html">Eigen::VecAsSymmMatrix</a></div><div class="ttdoc">Expression for presenting a vector as a DimxDim symmetric matrix.</div><div class="ttdef"><b>Definition</b> VecAsSymmMatrix.h:57</div></div>
<div class="ttc" id="anamespaceEigen_html"><div class="ttname"><a href="http://eigen.tuxfamily.org/dox/namespaceEigen.html">Eigen</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_a48fb26961da3afe39c877bacfaf3215.html">gismo</a></li><li class="navelem"><a class="el" href="dir_906eed49ab018d64c0970e7739089d77.html">src</a></li><li class="navelem"><a class="el" href="dir_4de348d47f42a77919db1107a94e3deb.html">gsMatrix</a></li><li class="navelem"><a class="el" href="VecAsSymmMatrix_8h.html">VecAsSymmMatrix.h</a></li>
<li class="footer">Generated on Sun Jan 12 2025 14:02:07 for G+Smo by Doxygen v1.9.8 </li>
</ul>
</div>
<!-- Piwik -->
<!--
<script type="text/javascript">
var pkBaseURL = (("https:" == document.location.protocol) ? "https://stats.sylphide-consulting.com/piwik/" : "http://stats.sylphide-consulting.com/piwik/");
document.write(unescape("%3Cscript src='" + pkBaseURL + "piwik.js' type='text/javascript'%3E%3C/script%3E"));
</script><script type="text/javascript">
try {
var piwikTracker = Piwik.getTracker(pkBaseURL + "piwik.php", 20);
piwikTracker.trackPageView();
piwikTracker.enableLinkTracking();
} catch( err ) {}
</script><noscript><p><img src="http://stats.sylphide-consulting.com/piwik/piwik.php?idsite=20" style="border:0" alt="" /></p></noscript>
-->
<!-- End Piwik Tracking Code -->
</body>
</html>