PageRenderTime 874ms CodeModel.GetById 48ms RepoModel.GetById 4ms app.codeStats 3ms

/doc/html/vabsSudo_8cpp_source.html

https://bitbucket.org/slawton/windturbinemdo-stevebitb
HTML | 3945 lines | 3932 code | 12 blank | 1 comment | 0 complexity | c0b71c54c03e62d2c468912abbc05e88 MD5 | raw file
Possible License(s): GPL-2.0, BSD-3-Clause

Large files files are truncated, but you can click here to view the full file

  1. <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
  2. <html xmlns="http://www.w3.org/1999/xhtml">
  3. <head>
  4. <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
  5. <title>WindTurbineMDO: structural/VABS_Calc/source/ver_2/vabsSudo.cpp Source File</title>
  6. <link href="tabs.css" rel="stylesheet" type="text/css"/>
  7. <link href="doxygen.css" rel="stylesheet" type="text/css" />
  8. <link href="navtree.css" rel="stylesheet" type="text/css"/>
  9. <script type="text/javascript" src="jquery.js"></script>
  10. <script type="text/javascript" src="resize.js"></script>
  11. <script type="text/javascript" src="navtree.js"></script>
  12. <script type="text/javascript">
  13. $(document).ready(initResizable);
  14. </script>
  15. <link href="search/search.css" rel="stylesheet" type="text/css"/>
  16. <script type="text/javascript" src="search/search.js"></script>
  17. <script type="text/javascript">
  18. $(document).ready(function() { searchBox.OnSelectItem(0); });
  19. </script>
  20. </head>
  21. <body>
  22. <div id="top"><!-- do not remove this div! -->
  23. <div id="titlearea">
  24. <table cellspacing="0" cellpadding="0">
  25. <tbody>
  26. <tr style="height: 56px;">
  27. <td style="padding-left: 0.5em;">
  28. <div id="projectname">WindTurbineMDO
  29. </div>
  30. </td>
  31. </tr>
  32. </tbody>
  33. </table>
  34. </div>
  35. <!-- Generated by Doxygen 1.7.5.1 -->
  36. <script type="text/javascript">
  37. var searchBox = new SearchBox("searchBox", "search",false,'Search');
  38. </script>
  39. <div id="navrow1" class="tabs">
  40. <ul class="tablist">
  41. <li><a href="index.html"><span>Main&#160;Page</span></a></li>
  42. <li><a href="namespaces.html"><span>Namespaces</span></a></li>
  43. <li><a href="annotated.html"><span>Classes</span></a></li>
  44. <li class="current"><a href="files.html"><span>Files</span></a></li>
  45. <li>
  46. <div id="MSearchBox" class="MSearchBoxInactive">
  47. <span class="left">
  48. <img id="MSearchSelect" src="search/mag_sel.png"
  49. onmouseover="return searchBox.OnSearchSelectShow()"
  50. onmouseout="return searchBox.OnSearchSelectHide()"
  51. alt=""/>
  52. <input type="text" id="MSearchField" value="Search" accesskey="S"
  53. onfocus="searchBox.OnSearchFieldFocus(true)"
  54. onblur="searchBox.OnSearchFieldFocus(false)"
  55. onkeyup="searchBox.OnSearchFieldChange(event)"/>
  56. </span><span class="right">
  57. <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
  58. </span>
  59. </div>
  60. </li>
  61. </ul>
  62. </div>
  63. <div id="navrow2" class="tabs2">
  64. <ul class="tablist">
  65. <li><a href="files.html"><span>File&#160;List</span></a></li>
  66. <li><a href="globals.html"><span>File&#160;Members</span></a></li>
  67. </ul>
  68. </div>
  69. </div>
  70. <div id="side-nav" class="ui-resizable side-nav-resizable">
  71. <div id="nav-tree">
  72. <div id="nav-tree-contents">
  73. </div>
  74. </div>
  75. <div id="splitbar" style="-moz-user-select:none;"
  76. class="ui-resizable-handle">
  77. </div>
  78. </div>
  79. <script type="text/javascript">
  80. initNavTree('vabsSudo_8cpp.html','');
  81. </script>
  82. <div id="doc-content">
  83. <div class="header">
  84. <div class="headertitle">
  85. <div class="title">structural/VABS_Calc/source/ver_2/vabsSudo.cpp</div> </div>
  86. </div>
  87. <div class="contents">
  88. <a href="vabsSudo_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
  89. <a name="l00002"></a>00002 <span class="comment">// #include &quot;matrix/matrixClass.h&quot;</span>
  90. <a name="l00003"></a>00003 <span class="comment">// #include &quot;/home/mike/work/UVic/NonlinearBeam/CAD_Calc/source/ver_1/meshSudo.cpp&quot;</span>
  91. <a name="l00004"></a>00004 <span class="comment">// #include &quot;/home/mike/work/UVic/NonlinearBeam/matrix_Solver/source/ver_1/matrixSolverSudo.cpp&quot;</span>
  92. <a name="l00005"></a>00005 <span class="comment">// #include &quot;/home/mike/work/UVic/NonlinearBeam/GEBT_Calc/source/ver_2/femSudoCode.cpp&quot;</span>
  93. <a name="l00006"></a>00006 <span class="comment">// #include &quot;material.cpp&quot;</span>
  94. <a name="l00007"></a>00007 <span class="preprocessor">#include &lt;algorithm&gt;</span>
  95. <a name="l00008"></a>00008 <span class="preprocessor">#include &lt;memory&gt;</span>
  96. <a name="l00009"></a>00009
  97. <a name="l00010"></a>00010 <span class="comment">// #define ADJOINT_VABS</span>
  98. <a name="l00011"></a>00011
  99. <a name="l00012"></a>00012 <span class="comment">// #define SVD_BASED_VABS_SOLUTION</span>
  100. <a name="l00013"></a>00013
  101. <a name="l00014"></a>00014 <span class="comment">// #define PRINT_VABS_PROGRESS</span>
  102. <a name="l00015"></a>00015
  103. <a name="l00016"></a>00016 <span class="comment">// #define TIME_ADJOINT_CALCULATION</span>
  104. <a name="l00017"></a>00017
  105. <a name="l00018"></a>00018 <span class="comment">// #include &lt;sys/time.h&gt;</span>
  106. <a name="l00019"></a>00019
  107. <a name="l00020"></a>00020 <span class="comment">// so lets talk about the VABS....</span>
  108. <a name="l00021"></a>00021
  109. <a name="l00022"></a>00022 <span class="comment">// ok the mesh program retrieves the minimal geometric information required to build a VABS model...</span>
  110. <a name="l00023"></a>00023
  111. <a name="l00024"></a>00024 <span class="comment">// A vabs element is only responsible for building a matrix</span>
  112. <a name="l00025"></a>00025
  113. <a name="l00026"></a>00026 <span class="comment">// Here is a structure</span>
  114. <a name="l00027"></a>00027 <span class="comment">// we have VABS elements that are like FEM elements with shape functions and matrix building abilities</span>
  115. <a name="l00028"></a>00028 <span class="comment">// we have VABS nodes that contain information from the mesh to build the matrices</span>
  116. <a name="l00029"></a>00029 <span class="comment">// we have a VABS assembly that takes a matrix solver and uses the assembly to get the matrices and solve the constituitive matrix</span>
  117. <a name="l00030"></a>00030
  118. <a name="l00031"></a>00031 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  119. <a name="l00032"></a><a class="code" href="classVABSNode.html">00032</a> <span class="keyword">class </span><a class="code" href="classVABSNode.html">VABSNode</a>
  120. <a name="l00033"></a>00033 {
  121. <a name="l00034"></a>00034 <span class="keyword">public</span>:
  122. <a name="l00035"></a>00035
  123. <a name="l00036"></a><a class="code" href="classVABSNode.html#afe63b8307d6bd3082ec5eddfd36ab28d">00036</a> <span class="keyword">virtual</span> <a class="code" href="classVABSNode.html#afe63b8307d6bd3082ec5eddfd36ab28d">~VABSNode</a>()
  124. <a name="l00037"></a>00037 {}
  125. <a name="l00038"></a>00038
  126. <a name="l00039"></a>00039 <span class="comment">// need to know position (2 scalars or 2-vect)</span>
  127. <a name="l00040"></a>00040 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode.html#a04a4be0511f905869894e34c0b2901d8">x</a>() <span class="keyword">const</span> =0;
  128. <a name="l00041"></a>00041 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode.html#ab8825545b32b55e848a2cdf394ebf8ed">y</a>() <span class="keyword">const</span> =0;
  129. <a name="l00042"></a>00042 <span class="comment">// need to know through thickness posisition (scalar)</span>
  130. <a name="l00043"></a>00043 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode.html#a47bb48aeeb59af55084e7743bd2963b5">h</a>() <span class="keyword">const</span> =0;
  131. <a name="l00044"></a>00044 <span class="comment">// need to know ply angle (scalar)</span>
  132. <a name="l00045"></a>00045 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode.html#acfd688e05d010009439d3459923ef7c0">plyAngle</a>() <span class="keyword">const</span> =0;
  133. <a name="l00046"></a>00046
  134. <a name="l00047"></a>00047 <span class="comment">// I am not sure if I need this as it could be stored in the material model as an h dependence</span>
  135. <a name="l00048"></a>00048 <span class="comment">// need to know fiber angle (scalar) (optional for now ...)</span>
  136. <a name="l00049"></a>00049
  137. <a name="l00050"></a>00050 <span class="comment">// this is where the results of the VABS calculate is stored</span>
  138. <a name="l00051"></a>00051 <span class="keyword">virtual</span> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;matrixStatic&lt;dataType,3,4&gt;</a> &gt;&amp; <a class="code" href="classVABSNode.html#a7bd73fd4e4208bb66e157c3882d01a29">warping</a>() <span class="keyword">const</span> =0;
  139. <a name="l00052"></a>00052
  140. <a name="l00053"></a>00053 };
  141. <a name="l00054"></a>00054
  142. <a name="l00055"></a>00055 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  143. <a name="l00056"></a><a class="code" href="classVABSNode__common.html">00056</a> <span class="keyword">class </span><a class="code" href="classVABSNode__common.html">VABSNode_common</a>
  144. <a name="l00057"></a>00057 {
  145. <a name="l00058"></a>00058 <span class="keyword">private</span>:
  146. <a name="l00059"></a>00059
  147. <a name="l00060"></a>00060 <a class="code" href="classdataSource__Data.html">dataSource_Data&lt;dataType&gt;</a> xSrc;
  148. <a name="l00061"></a>00061 <a class="code" href="classdataSource__Data.html">dataSource_Data&lt;dataType&gt;</a> ySrc;
  149. <a name="l00062"></a>00062 <a class="code" href="classdataSource__Data.html">dataSource_Data&lt;matrixStatic&lt;dataType,3,4&gt;</a>,<a class="code" href="classmatrixRdWtStatic__pData.html">matrixRdWtStatic_pData&lt;dataType,3,4&gt;</a> &gt; warpSrc;
  150. <a name="l00063"></a>00063
  151. <a name="l00064"></a>00064 <span class="keyword">public</span>:
  152. <a name="l00065"></a>00065
  153. <a name="l00066"></a><a class="code" href="classVABSNode__common.html#aa09906e00648ab7fb541a85432a53aff">00066</a> <a class="code" href="classVABSNode__common.html#aa09906e00648ab7fb541a85432a53aff">VABSNode_common</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> x_in=0.0,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> y_in=0.0) : xSrc(x_in) , ySrc(y_in)
  154. <a name="l00067"></a>00067 {}
  155. <a name="l00068"></a>00068
  156. <a name="l00069"></a><a class="code" href="classVABSNode__common.html#ac7174c77d18117ebb8304c9f325262d8">00069</a> <a class="code" href="classVABSNode__common.html">VABSNode_common</a>&amp; <a class="code" href="classVABSNode__common.html#ac7174c77d18117ebb8304c9f325262d8">setX</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> x_in)
  157. <a name="l00070"></a>00070 {
  158. <a name="l00071"></a>00071 xSrc.<a class="code" href="classdataSource__Data.html#aefa4c87ef5fb27413a5f2ce1356638d7">getData</a>()=x_in;
  159. <a name="l00072"></a>00072 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
  160. <a name="l00073"></a>00073 }
  161. <a name="l00074"></a>00074
  162. <a name="l00075"></a><a class="code" href="classVABSNode__common.html#ae9c8384acbcc1a44b7424cfee345df20">00075</a> <a class="code" href="classVABSNode__common.html">VABSNode_common</a>&amp; <a class="code" href="classVABSNode__common.html#ae9c8384acbcc1a44b7424cfee345df20">setY</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> y_in)
  163. <a name="l00076"></a>00076 {
  164. <a name="l00077"></a>00077 ySrc.<a class="code" href="classdataSource__Data.html#aefa4c87ef5fb27413a5f2ce1356638d7">getData</a>()=y_in;
  165. <a name="l00078"></a>00078 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
  166. <a name="l00079"></a>00079 }
  167. <a name="l00080"></a>00080
  168. <a name="l00081"></a><a class="code" href="classVABSNode__common.html#a08b644d8490442f08872ef2679d119f3">00081</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__common.html#a08b644d8490442f08872ef2679d119f3">x</a>()<span class="keyword"> const</span>
  169. <a name="l00082"></a>00082 <span class="keyword"> </span>{
  170. <a name="l00083"></a>00083 <span class="keywordflow">return</span> xSrc;
  171. <a name="l00084"></a>00084 }
  172. <a name="l00085"></a>00085
  173. <a name="l00086"></a><a class="code" href="classVABSNode__common.html#a091a6d771873b52c822d009c7bd8178e">00086</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__common.html#a091a6d771873b52c822d009c7bd8178e">y</a>()<span class="keyword"> const</span>
  174. <a name="l00087"></a>00087 <span class="keyword"> </span>{
  175. <a name="l00088"></a>00088 <span class="keywordflow">return</span> ySrc;
  176. <a name="l00089"></a>00089 }
  177. <a name="l00090"></a>00090
  178. <a name="l00091"></a><a class="code" href="classVABSNode__common.html#ab53b3d818bb3881738bed7b8a47c9eb2">00091</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;matrixStatic&lt;dataType,3,4&gt;</a> &gt;&amp; <a class="code" href="classVABSNode__common.html#ab53b3d818bb3881738bed7b8a47c9eb2">warping</a>()<span class="keyword"> const</span>
  179. <a name="l00092"></a>00092 <span class="keyword"> </span>{
  180. <a name="l00093"></a>00093 <span class="keywordflow">return</span> warpSrc;
  181. <a name="l00094"></a>00094 }
  182. <a name="l00095"></a>00095
  183. <a name="l00096"></a><a class="code" href="classVABSNode__common.html#af6524af5c58ff7a1c177aaf3aaba0cee">00096</a> <span class="keywordtype">void</span> <a class="code" href="classVABSNode__common.html#af6524af5c58ff7a1c177aaf3aaba0cee">setWarping</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>,<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> j,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a>* pntr)
  184. <a name="l00097"></a>00097 {
  185. <a name="l00098"></a>00098 warpSrc.<a class="code" href="classdataSource__Data.html#aefa4c87ef5fb27413a5f2ce1356638d7">getData</a>().setPointer(i,j,pntr,<span class="keyword">false</span>);
  186. <a name="l00099"></a>00099 }
  187. <a name="l00100"></a>00100 };
  188. <a name="l00101"></a>00101
  189. <a name="l00102"></a>00102 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  190. <a name="l00103"></a><a class="code" href="classVABSNode__element.html">00103</a> <span class="keyword">class </span><a class="code" href="classVABSNode__element.html">VABSNode_element</a> : <span class="keyword">public</span> <a class="code" href="classVABSNode.html">VABSNode</a>&lt;dataType&gt;
  191. <a name="l00104"></a>00104 {
  192. <a name="l00105"></a>00105 <span class="keyword">private</span>:
  193. <a name="l00106"></a>00106
  194. <a name="l00107"></a>00107 <a class="code" href="classdataSource__Data.html">dataSource_Data&lt;dataType&gt;</a> hSrc;
  195. <a name="l00108"></a>00108 <a class="code" href="classdataSource__Data.html">dataSource_Data&lt;dataType&gt;</a> angSrc;
  196. <a name="l00109"></a>00109 <a class="code" href="classVABSNode__common.html">VABSNode_common&lt;dataType&gt;</a>* common;
  197. <a name="l00110"></a>00110
  198. <a name="l00111"></a>00111 <span class="keyword">public</span>:
  199. <a name="l00112"></a>00112
  200. <a name="l00113"></a><a class="code" href="classVABSNode__element.html#a74b8f4f8455445ec22c902dd7190f69e">00113</a> <a class="code" href="classVABSNode__element.html#a74b8f4f8455445ec22c902dd7190f69e">VABSNode_element</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> h_in=0.5,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> ang_in=0.0,<a class="code" href="classVABSNode__common.html">VABSNode_common&lt;dataType&gt;</a>* common_in=NULL) : hSrc(h_in) , angSrc(ang_in) , common(common_in)
  201. <a name="l00114"></a>00114 {}
  202. <a name="l00115"></a>00115
  203. <a name="l00116"></a><a class="code" href="classVABSNode__element.html#a89fa980bf57dbff7ee2b322e4e7e37a1">00116</a> <a class="code" href="classVABSNode__element.html">VABSNode_element</a>&amp; <a class="code" href="classVABSNode__element.html#a89fa980bf57dbff7ee2b322e4e7e37a1">setH</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> h_in)
  204. <a name="l00117"></a>00117 {
  205. <a name="l00118"></a>00118 hSrc.<a class="code" href="classdataSource__Data.html#aefa4c87ef5fb27413a5f2ce1356638d7">getData</a>()=h_in;
  206. <a name="l00119"></a>00119 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
  207. <a name="l00120"></a>00120 }
  208. <a name="l00121"></a>00121
  209. <a name="l00122"></a><a class="code" href="classVABSNode__element.html#ab65defe78f48f9e7d9eb8f41f24a86dd">00122</a> <a class="code" href="classVABSNode__element.html">VABSNode_element</a>&amp; <a class="code" href="classVABSNode__element.html#ab65defe78f48f9e7d9eb8f41f24a86dd">setAng</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> ang_in)
  210. <a name="l00123"></a>00123 {
  211. <a name="l00124"></a>00124 angSrc.<a class="code" href="classdataSource__Data.html#aefa4c87ef5fb27413a5f2ce1356638d7">getData</a>()=ang_in;
  212. <a name="l00125"></a>00125 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
  213. <a name="l00126"></a>00126 }
  214. <a name="l00127"></a>00127
  215. <a name="l00128"></a><a class="code" href="classVABSNode__element.html#a86c4d761fb1b24efac2356d8527f5196">00128</a> <a class="code" href="classVABSNode__element.html">VABSNode_element</a>&amp; <a class="code" href="classVABSNode__element.html#a86c4d761fb1b24efac2356d8527f5196">setCommon</a>(<a class="code" href="classVABSNode__common.html">VABSNode_common&lt;dataType&gt;</a>* common_in)
  216. <a name="l00129"></a>00129 {
  217. <a name="l00130"></a>00130 common=common_in;
  218. <a name="l00131"></a>00131 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
  219. <a name="l00132"></a>00132 }
  220. <a name="l00133"></a>00133
  221. <a name="l00134"></a><a class="code" href="classVABSNode__element.html#a201c12a31d10e45f0036ecc3e7509fe4">00134</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__element.html#a201c12a31d10e45f0036ecc3e7509fe4">x</a>()<span class="keyword"> const</span>
  222. <a name="l00135"></a>00135 <span class="keyword"> </span>{
  223. <a name="l00136"></a>00136 <span class="keywordflow">if</span>(common==NULL)
  224. <a name="l00137"></a>00137 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;VABSNode_element::x(): The common node data has not been set&quot;</span>);
  225. <a name="l00138"></a>00138 <span class="keywordflow">return</span> common-&gt;x();
  226. <a name="l00139"></a>00139 }
  227. <a name="l00140"></a>00140
  228. <a name="l00141"></a><a class="code" href="classVABSNode__element.html#a5b4524069dfe5115920e9efe1aeebe13">00141</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__element.html#a5b4524069dfe5115920e9efe1aeebe13">y</a>()<span class="keyword"> const</span>
  229. <a name="l00142"></a>00142 <span class="keyword"> </span>{
  230. <a name="l00143"></a>00143 <span class="keywordflow">if</span>(common==NULL)
  231. <a name="l00144"></a>00144 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;VABSNode_element::y(): The common node data has not been set&quot;</span>);
  232. <a name="l00145"></a>00145 <span class="keywordflow">return</span> common-&gt;y();
  233. <a name="l00146"></a>00146 }
  234. <a name="l00147"></a>00147
  235. <a name="l00148"></a><a class="code" href="classVABSNode__element.html#ad2dbcf0c1d265e3a2edcaa542e8f22c2">00148</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;matrixStatic&lt;dataType,3,4&gt;</a> &gt;&amp; <a class="code" href="classVABSNode__element.html#ad2dbcf0c1d265e3a2edcaa542e8f22c2">warping</a>()<span class="keyword"> const</span>
  236. <a name="l00149"></a>00149 <span class="keyword"> </span>{
  237. <a name="l00150"></a>00150 <span class="keywordflow">if</span>(common==NULL)
  238. <a name="l00151"></a>00151 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;VABSNode_element::warping(): The common node data has not been set&quot;</span>);
  239. <a name="l00152"></a>00152 <span class="keywordflow">return</span> common-&gt;warping();
  240. <a name="l00153"></a>00153 }
  241. <a name="l00154"></a>00154
  242. <a name="l00155"></a><a class="code" href="classVABSNode__element.html#acfb36b6a9bdc5e5b39586adb35dc5514">00155</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__element.html#acfb36b6a9bdc5e5b39586adb35dc5514">h</a>()<span class="keyword"> const</span>
  243. <a name="l00156"></a>00156 <span class="keyword"> </span>{
  244. <a name="l00157"></a>00157 <span class="keywordflow">return</span> hSrc;
  245. <a name="l00158"></a>00158 }
  246. <a name="l00159"></a>00159
  247. <a name="l00160"></a><a class="code" href="classVABSNode__element.html#a71ba830ae4368a2d368624a48a4b0dad">00160</a> <span class="keyword">const</span> <a class="code" href="classdataSource.html">dataSource&lt;dataType&gt;</a>&amp; <a class="code" href="classVABSNode__element.html#a71ba830ae4368a2d368624a48a4b0dad">plyAngle</a>()<span class="keyword"> const</span>
  248. <a name="l00161"></a>00161 <span class="keyword"> </span>{
  249. <a name="l00162"></a>00162 <span class="keywordflow">return</span> angSrc;
  250. <a name="l00163"></a>00163 }
  251. <a name="l00164"></a>00164 };
  252. <a name="l00165"></a>00165
  253. <a name="l00166"></a>00166 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  254. <a name="l00167"></a><a class="code" href="classVABS1__MatrixSource.html">00167</a> <span class="keyword">class </span><a class="code" href="classVABS1__MatrixSource.html">VABS1_MatrixSource</a>
  255. <a name="l00168"></a>00168 {
  256. <a name="l00169"></a>00169 <span class="keyword">public</span>:
  257. <a name="l00170"></a>00170
  258. <a name="l00171"></a><a class="code" href="classVABS1__MatrixSource.html#a5a3abf044ac45bc7c7f19847bedc8db8">00171</a> <span class="keyword">virtual</span> <a class="code" href="classVABS1__MatrixSource.html#a5a3abf044ac45bc7c7f19847bedc8db8">~VABS1_MatrixSource</a>()
  259. <a name="l00172"></a>00172 {}
  260. <a name="l00173"></a>00173
  261. <a name="l00174"></a>00174 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classVABS1__MatrixSource.html#a161e7c3230c9372c9e70e68ebea02d34">addE</a>()=0;
  262. <a name="l00175"></a>00175 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classVABS1__MatrixSource.html#a2569e6b04d8c73c5615829db929066b9">addDee</a>()=0;
  263. <a name="l00176"></a>00176 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classVABS1__MatrixSource.html#a16927e196fd23ba8c8809097aee990cb">addDae</a>()=0;
  264. <a name="l00177"></a>00177 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classVABS1__MatrixSource.html#a73c217bf4f8c74b74223f724a7cad571">addH</a>()=0;
  265. <a name="l00178"></a>00178 <span class="comment">// virtual void addPsi()=0;</span>
  266. <a name="l00179"></a>00179 };
  267. <a name="l00180"></a>00180
  268. <a name="l00181"></a>00181 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  269. <a name="l00182"></a><a class="code" href="classfunction2D__Gamma__a.html">00182</a> <span class="keyword">class </span><a class="code" href="classfunction2D__Gamma__a.html">function2D_Gamma_a</a> : <span class="keyword">public</span> <a class="code" href="classfunction2D.html">function2D</a>&lt;matrixRdWtSz_Data&lt;dataType&gt;,dataType&gt;
  270. <a name="l00183"></a>00183 {
  271. <a name="l00184"></a>00184 <span class="keyword">private</span>:
  272. <a name="l00185"></a>00185
  273. <a name="l00186"></a>00186 std::vector&lt;isoGrad2D&lt;dataType,dataType&gt; &gt;* shpGrad;
  274. <a name="l00187"></a>00187
  275. <a name="l00188"></a>00188 <span class="keyword">public</span>:
  276. <a name="l00189"></a>00189
  277. <a name="l00190"></a><a class="code" href="classfunction2D__Gamma__a.html#a2fa3ec16927f9f8a113a5881c42abf25">00190</a> <a class="code" href="classfunction2D__Gamma__a.html#a2fa3ec16927f9f8a113a5881c42abf25">function2D_Gamma_a</a>(std::vector&lt;<a class="code" href="classisoGrad2D.html">isoGrad2D&lt;dataType,dataType&gt;</a> &gt;* shpGrad_in=NULL) : shpGrad(shpGrad_in)
  278. <a name="l00191"></a>00191 {}
  279. <a name="l00192"></a>00192
  280. <a name="l00193"></a><a class="code" href="classfunction2D__Gamma__a.html#ada070084c7e4862d43d00e273e866398">00193</a> <span class="keyword">virtual</span> <a class="code" href="classfunction2D__Gamma__a.html#ada070084c7e4862d43d00e273e866398">~function2D_Gamma_a</a>()
  281. <a name="l00194"></a>00194 {}
  282. <a name="l00195"></a>00195
  283. <a name="l00196"></a><a class="code" href="classfunction2D__Gamma__a.html#a87b84825cd3c347f8bc7681402a4b015">00196</a> <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> <a class="code" href="classfunction2D__Gamma__a.html#a87b84825cd3c347f8bc7681402a4b015">getValue2D</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at1,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at2)
  284. <a name="l00197"></a>00197 {
  285. <a name="l00198"></a>00198 <span class="keywordflow">if</span>(shpGrad==NULL)
  286. <a name="l00199"></a>00199 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;function2D_Gamma_a::getValue2D(): The shape grad has not been set&quot;</span>);
  287. <a name="l00200"></a>00200 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> colCnt=shpGrad-&gt;size();
  288. <a name="l00201"></a>00201 <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> retVal(6,colCnt*3);
  289. <a name="l00202"></a>00202 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>=0;<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>&lt;colCnt;<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>++)
  290. <a name="l00203"></a>00203 {
  291. <a name="l00204"></a>00204 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> j=0;j&lt;6;j++)
  292. <a name="l00205"></a>00205 {
  293. <a name="l00206"></a>00206 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> k=0;k&lt;3;k++)
  294. <a name="l00207"></a>00207 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(j,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+k)=0.0;
  295. <a name="l00208"></a>00208 }
  296. <a name="l00209"></a>00209 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(1,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+1)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad1(at1,at2); <span class="comment">// 22</span>
  297. <a name="l00210"></a>00210 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(2,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+2)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad2(at1,at2); <span class="comment">// 33</span>
  298. <a name="l00211"></a>00211 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(3,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad1(at1,at2); <span class="comment">// 12</span>
  299. <a name="l00212"></a>00212 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(4,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad2(at1,at2); <span class="comment">// 13</span>
  300. <a name="l00213"></a>00213 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(5,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+1)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad2(at1,at2); <span class="comment">// 23</span>
  301. <a name="l00214"></a>00214 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(5,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+2)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getGrad1(at1,at2); <span class="comment">// 23</span>
  302. <a name="l00215"></a>00215 }
  303. <a name="l00216"></a>00216 <span class="keywordflow">return</span> retVal;
  304. <a name="l00217"></a>00217 }
  305. <a name="l00218"></a>00218 };
  306. <a name="l00219"></a>00219
  307. <a name="l00220"></a>00220 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  308. <a name="l00221"></a><a class="code" href="classfunction2D__Gamma__e.html">00221</a> <span class="keyword">class </span><a class="code" href="classfunction2D__Gamma__e.html">function2D_Gamma_e</a> : <span class="keyword">public</span> <a class="code" href="classfunction2D.html">function2D</a>&lt;matrixRdWtStatic_Data&lt;dataType,6,4&gt;,dataType&gt;
  309. <a name="l00222"></a>00222 {
  310. <a name="l00223"></a>00223 <span class="keyword">private</span>:
  311. <a name="l00224"></a>00224
  312. <a name="l00225"></a>00225 <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* xLoc;
  313. <a name="l00226"></a>00226 <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* yLoc;
  314. <a name="l00227"></a>00227
  315. <a name="l00228"></a>00228 <span class="keyword">static</span> <span class="keywordtype">bool</span> isFirst;
  316. <a name="l00229"></a>00229
  317. <a name="l00230"></a>00230 <span class="keyword">public</span>:
  318. <a name="l00231"></a>00231
  319. <a name="l00232"></a><a class="code" href="classfunction2D__Gamma__e.html#aec36edf5f06e62dad57f89d47accc153">00232</a> <a class="code" href="classfunction2D__Gamma__e.html#aec36edf5f06e62dad57f89d47accc153">function2D_Gamma_e</a>(<a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* xLoc_in=NULL,<a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* yLoc_in=NULL) : xLoc(xLoc_in) , yLoc(yLoc_in)
  320. <a name="l00233"></a>00233 {}
  321. <a name="l00234"></a>00234
  322. <a name="l00235"></a><a class="code" href="classfunction2D__Gamma__e.html#a1fe0800dade51e4093b894ac7f30536f">00235</a> <a class="code" href="classmatrixRdWtStatic__Data.html">matrixRdWtStatic_Data&lt;dataType,6,4&gt;</a> <a class="code" href="classfunction2D__Gamma__e.html#a1fe0800dade51e4093b894ac7f30536f">getValue2D</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at1,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at2)
  323. <a name="l00236"></a>00236 {
  324. <a name="l00237"></a>00237 <span class="keywordflow">if</span>(xLoc==NULL || yLoc==NULL)
  325. <a name="l00238"></a>00238 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;function2D_Gamma_e::getValue2D(): The x and y locations have not been set&quot;</span>);
  326. <a name="l00239"></a>00239 <a class="code" href="classmatrixRdWtStatic__Data.html">matrixRdWtStatic_Data&lt;dataType,6,4&gt;</a> retVal;
  327. <a name="l00240"></a>00240 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> j=0;j&lt;6;j++)
  328. <a name="l00241"></a>00241 {
  329. <a name="l00242"></a>00242 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> k=0;k&lt;4;k++)
  330. <a name="l00243"></a>00243 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(j,k)=0.0;
  331. <a name="l00244"></a>00244 }
  332. <a name="l00245"></a>00245 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(0,0)=1.0; <span class="comment">// 11</span>
  333. <a name="l00246"></a>00246 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(0,2)=yLoc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2); <span class="comment">// 11</span>
  334. <a name="l00247"></a>00247 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(0,3)=-(xLoc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2)); <span class="comment">// 11</span>
  335. <a name="l00248"></a>00248 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(3,1)=-(yLoc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2)); <span class="comment">// 12</span>
  336. <a name="l00249"></a>00249 retVal.<a class="code" href="classmatrix.html#a81a746378de5ff16645b2cc51189e302">entry</a>(4,1)=xLoc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2); <span class="comment">// 13</span>
  337. <a name="l00250"></a>00250
  338. <a name="l00251"></a>00251 <span class="keywordflow">if</span>(isFirst)
  339. <a name="l00252"></a>00252 {
  340. <a name="l00253"></a>00253 std::cout &lt;&lt; <span class="stringliteral">&quot;Remember that you have to include the metric tensor into VABS for the case where there is initial curvature&quot;</span> &lt;&lt; std::endl;
  341. <a name="l00254"></a>00254 std::cout &lt;&lt; <span class="stringliteral">&quot;One other thing... the vabs blade does not shift the cross section to the desired origin. Just to the Xbar and Ybar&quot;</span> &lt;&lt; std::endl;
  342. <a name="l00255"></a>00255 isFirst=<span class="keyword">false</span>;
  343. <a name="l00256"></a>00256 }
  344. <a name="l00257"></a>00257
  345. <a name="l00258"></a>00258 <span class="keywordflow">return</span> retVal;
  346. <a name="l00259"></a>00259 }
  347. <a name="l00260"></a>00260 };
  348. <a name="l00261"></a>00261
  349. <a name="l00262"></a>00262 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  350. <a name="l00263"></a>00263 <span class="keywordtype">bool</span> <a class="code" href="classfunction2D__Gamma__e.html">function2D_Gamma_e&lt;dataType&gt;::isFirst</a>=<span class="keyword">true</span>;
  351. <a name="l00264"></a>00264
  352. <a name="l00265"></a>00265 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  353. <a name="l00266"></a><a class="code" href="classfunction2D__ShapeDiag.html">00266</a> <span class="keyword">class </span><a class="code" href="classfunction2D__ShapeDiag.html">function2D_ShapeDiag</a> : <span class="keyword">public</span> <a class="code" href="classfunction2D.html">function2D</a>&lt;matrixRdWtSz_Data&lt;dataType&gt;,dataType&gt;
  354. <a name="l00267"></a>00267 {
  355. <a name="l00268"></a>00268 <span class="keyword">private</span>:
  356. <a name="l00269"></a>00269
  357. <a name="l00270"></a>00270 std::vector&lt;isoGrad2D&lt;dataType,dataType&gt; &gt;* shpGrad;
  358. <a name="l00271"></a>00271
  359. <a name="l00272"></a>00272 <span class="keyword">public</span>:
  360. <a name="l00273"></a>00273
  361. <a name="l00274"></a><a class="code" href="classfunction2D__ShapeDiag.html#a4826cb643557ee732da15604ba877c4c">00274</a> <a class="code" href="classfunction2D__ShapeDiag.html#a4826cb643557ee732da15604ba877c4c">function2D_ShapeDiag</a>(std::vector&lt;<a class="code" href="classisoGrad2D.html">isoGrad2D&lt;dataType,dataType&gt;</a> &gt;* shpGrad_in=NULL) : shpGrad(shpGrad_in)
  362. <a name="l00275"></a>00275 {}
  363. <a name="l00276"></a>00276
  364. <a name="l00277"></a><a class="code" href="classfunction2D__ShapeDiag.html#a72675ea3d3e7fc00fc05fb34f3ff1d29">00277</a> <span class="keyword">virtual</span> <a class="code" href="classfunction2D__ShapeDiag.html#a72675ea3d3e7fc00fc05fb34f3ff1d29">~function2D_ShapeDiag</a>()
  365. <a name="l00278"></a>00278 {}
  366. <a name="l00279"></a>00279
  367. <a name="l00280"></a><a class="code" href="classfunction2D__ShapeDiag.html#a332be29feaa9b2b691036da939c2eb43">00280</a> <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> <a class="code" href="classfunction2D__ShapeDiag.html#a332be29feaa9b2b691036da939c2eb43">getValue2D</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at1,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at2)
  368. <a name="l00281"></a>00281 {
  369. <a name="l00282"></a>00282 <span class="keywordflow">if</span>(shpGrad==NULL)
  370. <a name="l00283"></a>00283 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;function2D_ShapeDiag::getValue2D(): The shape grad has not been set&quot;</span>);
  371. <a name="l00284"></a>00284 <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> colCnt=shpGrad-&gt;size();
  372. <a name="l00285"></a>00285 <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> retVal(3,colCnt*3);
  373. <a name="l00286"></a>00286 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>=0;<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>&lt;colCnt;<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>++)
  374. <a name="l00287"></a>00287 {
  375. <a name="l00288"></a>00288 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> j=0;j&lt;3;j++)
  376. <a name="l00289"></a>00289 {
  377. <a name="l00290"></a>00290 <span class="keywordflow">for</span>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> k=0;k&lt;3;k++)
  378. <a name="l00291"></a>00291 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(j,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+k)=0.0;
  379. <a name="l00292"></a>00292 }
  380. <a name="l00293"></a>00293 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(0,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getValue2D(at1,at2);
  381. <a name="l00294"></a>00294 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(1,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+1)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getValue2D(at1,at2);
  382. <a name="l00295"></a>00295 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(2,<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>*3+2)=(*shpGrad)[<a class="code" href="namespaceCPPL.html#aca67a31b860a54cc09b2519953aa2b39">i</a>].getValue2D(at1,at2);
  383. <a name="l00296"></a>00296 }
  384. <a name="l00297"></a>00297
  385. <a name="l00298"></a>00298 <span class="keywordflow">return</span> retVal;
  386. <a name="l00299"></a>00299 }
  387. <a name="l00300"></a>00300 };
  388. <a name="l00301"></a>00301
  389. <a name="l00302"></a>00302 <span class="keyword">template</span>&lt;<span class="keyword">class</span> dataType&gt;
  390. <a name="l00303"></a><a class="code" href="classfunction2D__vabsNull.html">00303</a> <span class="keyword">class </span><a class="code" href="classfunction2D__vabsNull.html">function2D_vabsNull</a> : <span class="keyword">public</span> <a class="code" href="classfunction2D.html">function2D</a>&lt;matrixRdWtSz_Data&lt;dataType&gt;,dataType&gt;
  391. <a name="l00304"></a>00304 {
  392. <a name="l00305"></a>00305 <span class="keyword">private</span>:
  393. <a name="l00306"></a>00306
  394. <a name="l00307"></a>00307 <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* xFunc;
  395. <a name="l00308"></a>00308 <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* yFunc;
  396. <a name="l00309"></a>00309 <span class="keywordtype">bool</span> deleteX;
  397. <a name="l00310"></a>00310 <span class="keywordtype">bool</span> deleteY;
  398. <a name="l00311"></a>00311
  399. <a name="l00312"></a>00312 <span class="keyword">public</span>:
  400. <a name="l00313"></a>00313
  401. <a name="l00314"></a><a class="code" href="classfunction2D__vabsNull.html#abed270cdd88dffa0f3520bc651b6e28e">00314</a> <a class="code" href="classfunction2D__vabsNull.html#abed270cdd88dffa0f3520bc651b6e28e">function2D_vabsNull</a>(
  402. <a name="l00315"></a>00315 <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* xFunc_in=NULL
  403. <a name="l00316"></a>00316 , <a class="code" href="classfunction2D.html">function2D&lt;dataType,dataType&gt;</a>* yFunc_in=NULL
  404. <a name="l00317"></a>00317 , <span class="keywordtype">bool</span> deleteX_in=<span class="keyword">false</span>
  405. <a name="l00318"></a>00318 , <span class="keywordtype">bool</span> deleteY_in=<span class="keyword">false</span>
  406. <a name="l00319"></a>00319 )
  407. <a name="l00320"></a>00320 : xFunc(xFunc_in)
  408. <a name="l00321"></a>00321 , yFunc(yFunc_in)
  409. <a name="l00322"></a>00322 , deleteX(deleteX_in)
  410. <a name="l00323"></a>00323 , deleteY(deleteY_in)
  411. <a name="l00324"></a>00324 {}
  412. <a name="l00325"></a>00325
  413. <a name="l00326"></a><a class="code" href="classfunction2D__vabsNull.html#a8895ac9819cb9b1044913567cb283325">00326</a> <span class="keyword">virtual</span> <a class="code" href="classfunction2D__vabsNull.html#a8895ac9819cb9b1044913567cb283325">~function2D_vabsNull</a>()
  414. <a name="l00327"></a>00327 {
  415. <a name="l00328"></a>00328 <span class="keywordflow">if</span>(deleteX &amp;&amp; xFunc!=NULL)
  416. <a name="l00329"></a>00329 <span class="keyword">delete</span> xFunc;
  417. <a name="l00330"></a>00330 <span class="keywordflow">if</span>(deleteY &amp;&amp; yFunc!=NULL)
  418. <a name="l00331"></a>00331 <span class="keyword">delete</span> yFunc;
  419. <a name="l00332"></a>00332 }
  420. <a name="l00333"></a>00333
  421. <a name="l00334"></a><a class="code" href="classfunction2D__vabsNull.html#a842b1f91ef7d5ed7c17eaadfb6de2400">00334</a> <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> <a class="code" href="classfunction2D__vabsNull.html#a842b1f91ef7d5ed7c17eaadfb6de2400">getValue2D</a>(<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at1,<a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> at2)
  422. <a name="l00335"></a>00335 {
  423. <a name="l00336"></a>00336 <span class="keywordflow">if</span>(xFunc==NULL)
  424. <a name="l00337"></a>00337 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;function2D_vabsNull::getValue2D(): The x function has not been specified&quot;</span>);
  425. <a name="l00338"></a>00338 <span class="keywordflow">if</span>(yFunc==NULL)
  426. <a name="l00339"></a>00339 <span class="keywordflow">throw</span> <a class="code" href="classeMsg.html">eMsg</a>(<span class="stringliteral">&quot;function2D_vabsNull::getValue2D(): The y function has not been specified&quot;</span>);
  427. <a name="l00340"></a>00340
  428. <a name="l00341"></a>00341 <a class="code" href="classmatrixRdWtSz__Data.html">matrixRdWtSz_Data&lt;dataType&gt;</a> retVal=<a class="code" href="classmatrixStatic__Zero.html">matrixStatic_Zero&lt;dataType,3,4&gt;</a>();
  429. <a name="l00342"></a>00342
  430. <a name="l00343"></a>00343 <a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> x=xFunc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2);
  431. <a name="l00344"></a>00344 <a class="code" href="timoshenko_8cpp.html#a926743e1c2bd9fc382c84014036246c2">dataType</a> y=yFunc-&gt;<a class="code" href="classfunction2D.html#ac48703a8b417ebe6ae34233ee4d81638">getValue2D</a>(at1,at2);
  432. <a name="l00345"></a>00345
  433. <a name="l00346"></a>00346 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(0,0)=1.0;
  434. <a name="l00347"></a>00347 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(1,1)=1.0;
  435. <a name="l00348"></a>00348 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(2,2)=1.0;
  436. <a name="l00349"></a>00349 retVal.<a class="code" href="classmatrixRdWt.html#a85aaa63333c688527b5e62dbce69cb61">entry</a>(1,3)=y;
  437. <a name="…

Large files files are truncated, but you can click here to view the full file