<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <html><head><meta name="robots" content="noindex"> <meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1"> <title>ITK: itkFEMSolver.h Source File</title> <link href="DoxygenStyle.css" rel="stylesheet" type="text/css"> </head><body bgcolor="#ffffff"> <!-- Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 --> <center> <a href="index.html" class="qindex">Main Page</a> <a href="modules.html" class="qindex">Groups</a> <a href="namespaces.html" class="qindex">Namespace List</a> <a href="hierarchy.html" class="qindex">Class Hierarchy</a> <a href="classes.html" class="qindex">Alphabetical List</a> <a href="annotated.html" class="qindex">Compound List</a> <a href="files.html" class="qindex">File List</a> <a href="namespacemembers.html" class="qindex">Namespace Members</a> <a href="functions.html" class="qindex">Compound Members</a> <a href="globals.html" class="qindex">File Members</a> <a href="pages.html" class="qindex">Concepts</a></center> <!-- Generated by Doxygen 1.5.9 --> <div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a> » <a class="el" href="dir_2c33022d19e8077162952ea9ed57656d.html">Numerics</a> » <a class="el" href="dir_316e5ab5c8a9cb56a9dce1c9fc454a23.html">FEM</a> </div> <div class="contents"> <h1>itkFEMSolver.h</h1><a href="itkFEMSolver_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span> <a name="l00002"></a>00002 <span class="comment"></span> <a name="l00003"></a>00003 <span class="comment"> Program: Insight Segmentation & Registration Toolkit</span> <a name="l00004"></a>00004 <span class="comment"> Module: $RCSfile: itkFEMSolver.h,v $</span> <a name="l00005"></a>00005 <span class="comment"> Language: C++</span> <a name="l00006"></a>00006 <span class="comment"> Date: $Date: 2009-01-30 21:53:03 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.28 $</span> <a name="l00008"></a>00008 <span class="comment"></span> <a name="l00009"></a>00009 <span class="comment"> Copyright (c) Insight Software Consortium. All rights reserved.</span> <a name="l00010"></a>00010 <span class="comment"> See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span> <a name="l00011"></a>00011 <span class="comment"></span> <a name="l00012"></a>00012 <span class="comment"> This software is distributed WITHOUT ANY WARRANTY; without even </span> <a name="l00013"></a>00013 <span class="comment"> the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR </span> <a name="l00014"></a>00014 <span class="comment"> PURPOSE. See the above copyright notices for more information.</span> <a name="l00015"></a>00015 <span class="comment"></span> <a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span> <a name="l00017"></a>00017 <a name="l00018"></a>00018 <span class="preprocessor">#ifndef __itkFEMSolver_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span><span class="preprocessor">#define __itkFEMSolver_h</span> <a name="l00020"></a>00020 <span class="preprocessor"></span> <a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkFEMElementBase_8h.html">itkFEMElementBase.h</a>"</span> <a name="l00022"></a>00022 <span class="preprocessor">#include "<a class="code" href="itkFEMMaterialBase_8h.html">itkFEMMaterialBase.h</a>"</span> <a name="l00023"></a>00023 <span class="preprocessor">#include "<a class="code" href="itkFEMLoadBase_8h.html">itkFEMLoadBase.h</a>"</span> <a name="l00024"></a>00024 <a name="l00025"></a>00025 <span class="preprocessor">#include "<a class="code" href="itkFEMLinearSystemWrapper_8h.html">itkFEMLinearSystemWrapper.h</a>"</span> <a name="l00026"></a>00026 <span class="preprocessor">#include "<a class="code" href="itkFEMLinearSystemWrapperVNL_8h.html">itkFEMLinearSystemWrapperVNL.h</a>"</span> <a name="l00027"></a>00027 <a name="l00028"></a>00028 <span class="preprocessor">#include "<a class="code" href="itkImage_8h.html">itkImage.h</a>"</span> <a name="l00029"></a>00029 <a name="l00030"></a>00030 <span class="keyword">namespace </span>itk { <a name="l00031"></a>00031 <span class="keyword">namespace </span>fem { <a name="l00032"></a>00032 <a name="l00041"></a><a class="code" href="classitk_1_1fem_1_1Solver.html">00041</a> <span class="keyword">class </span><a class="code" href="classitk_1_1fem_1_1Solver.html" title="Main FEM solver.">Solver</a> <a name="l00042"></a>00042 { <a name="l00043"></a>00043 <span class="keyword">public</span>: <a name="l00044"></a>00044 <a name="l00048"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">00048</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1Element.html#c55f879c15daa302723012b7f9d21b94">Element::Float</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">Float</a>; <a name="l00049"></a>00049 <a name="l00054"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#2d5529ae60b8ce66e8540e49bbc23f2f">00054</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">Element::ArrayType</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">ElementArray</a>; <a name="l00055"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#6f2852ca227980b30bb1cc0ee3f5260f">00055</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">ElementArray</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#6f2852ca227980b30bb1cc0ee3f5260f">el</a>; <a name="l00056"></a>00056 <a name="l00060"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#e18f600e4f3b1782d74257b0f33c5636">00060</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">Node::ArrayType</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">NodeArray</a>; <a name="l00061"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#6762185b526b949d140efe0b36792f6b">00061</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">NodeArray</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#6762185b526b949d140efe0b36792f6b">node</a>; <a name="l00062"></a>00062 <a name="l00066"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#964128f515d7ba6b804c4cbee9bdc016">00066</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">Load::ArrayType</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">LoadArray</a>; <a name="l00067"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#7f6bb9739f97c41a075d25f02c66395b">00067</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">LoadArray</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#7f6bb9739f97c41a075d25f02c66395b">load</a>; <a name="l00068"></a>00068 <a name="l00072"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#ffae74f3d9d2de8d0453c25e7efb657f">00072</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">Material::ArrayType</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">MaterialArray</a>; <a name="l00073"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#75e35b0c24e7ac0be0b88ccba5cb23f6">00073</a> <a class="code" href="classitk_1_1fem_1_1FEMPArray.html" title="Array for FEMP objects.">MaterialArray</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#75e35b0c24e7ac0be0b88ccba5cb23f6">mat</a>; <a name="l00074"></a>00074 <a name="l00078"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#442958a53672d2899b2f81c8547777df">00078</a> <span class="keyword">typedef</span> <a class="code" href="classvnl__vector.html">Element::VectorType</a> <a class="code" href="classvnl__vector.html">VectorType</a>; <a name="l00079"></a>00079 <a name="l00090"></a>00090 <a class="code" href="classitk_1_1fem_1_1Solver.html#954e93629e4055eaff0c12b09b1ec712">itkStaticConstMacro</a>(MaxGridDimensions, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>, 3); <a name="l00091"></a>00091 <a name="l00095"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#1ed801173a178f4b73db7c1863cfed9f">00095</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Image.html" title="Templated n-dimensional image class.">itk::Image<Element::ConstPointer,MaxGridDimensions></a> <a class="code" href="classitk_1_1Image.html" title="Templated n-dimensional image class.">InterpolationGridType</a>; <a name="l00096"></a>00096 <a name="l00111"></a>00111 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#c648debf94f686eadb682c14ae8625ba">InitializeInterpolationGrid</a>(<span class="keyword">const</span> <a class="code" href="classvnl__vector.html">VectorType</a>& size, <span class="keyword">const</span> <a class="code" href="classvnl__vector.html">VectorType</a>& bb1, <span class="keyword">const</span> <a class="code" href="classvnl__vector.html">VectorType</a>& bb2); <a name="l00112"></a>00112 <a name="l00116"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#3e7dc50c626fbe771ab9dcc5d7f3daa3">00116</a> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#c648debf94f686eadb682c14ae8625ba">InitializeInterpolationGrid</a>(<span class="keyword">const</span> <a class="code" href="classvnl__vector.html">VectorType</a>& size) <a name="l00117"></a>00117 { <a name="l00118"></a>00118 <a class="code" href="classitk_1_1fem_1_1Solver.html#c648debf94f686eadb682c14ae8625ba">InitializeInterpolationGrid</a>(size, <a class="code" href="classitk_1_1fem_1_1Solver.html#442958a53672d2899b2f81c8547777df">VectorType</a>(size.size(),0.0), size-1.0); <a name="l00119"></a>00119 } <a name="l00121"></a>00121 <a name="l00132"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#b1cc210b38e6cc6193df4503c6813912">00132</a> <span class="keyword">const</span> <a class="code" href="classitk_1_1Image.html" title="Templated n-dimensional image class.">InterpolationGridType</a> * <a class="code" href="classitk_1_1fem_1_1Solver.html#b1cc210b38e6cc6193df4503c6813912">GetInterpolationGrid</a>(<span class="keywordtype">void</span>)<span class="keyword"> const</span> <a name="l00133"></a>00133 <span class="keyword"> </span>{ <span class="keywordflow">return</span> m_InterpolationGrid.<a class="code" href="classitk_1_1SmartPointer.html#0734d2308d93e0ca10199b79531db5de">GetPointer</a>(); } <a name="l00134"></a>00134 <a name="l00143"></a>00143 <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1Element.html" title="Abstract base element class.">Element</a> * <a class="code" href="classitk_1_1fem_1_1Solver.html#283bf84597aaa88dfc132d37d3059a79">GetElementAtPoint</a>(<span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#442958a53672d2899b2f81c8547777df">VectorType</a>& pt) <span class="keyword">const</span>; <a name="l00144"></a>00144 <a name="l00148"></a>00148 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#9292175bd5f735b99fd37a605f39a8aa">Read</a>( std::istream& f ); <a name="l00149"></a>00149 <a name="l00153"></a>00153 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#84bd13e3a4d716c9f8b5dd628b454d33">Write</a>( std::ostream& f ); <a name="l00154"></a>00154 <a name="l00158"></a>00158 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#7a255226c48b5bda03911d3f16d27c82">Clear</a>( <span class="keywordtype">void</span> ); <a name="l00159"></a>00159 <a name="l00160"></a>00160 <a name="l00169"></a>00169 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#e264fad911491bc0bd58efd7ebfbcbdf">GenerateGFN</a>( <span class="keywordtype">void</span> ); <a name="l00170"></a>00170 <a name="l00174"></a>00174 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#344ef3f19bfcf305af43d54597632782">AssembleK</a>( <span class="keywordtype">void</span> ); <a name="l00175"></a>00175 <a name="l00182"></a>00182 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#789cc6c7beebce0f37f5c3ceb290de7f">InitializeMatrixForAssembly</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> N); <a name="l00183"></a>00183 <a name="l00189"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#ddea79b5dbd2ea4d005e9a600a0ae56d">00189</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#ddea79b5dbd2ea4d005e9a600a0ae56d">FinalizeMatrixAfterAssembly</a>( <span class="keywordtype">void</span> ) <a name="l00190"></a>00190 { <a name="l00191"></a>00191 <span class="comment">// Apply the boundary conditions to the K matrix</span> <a name="l00192"></a>00192 this-><a class="code" href="classitk_1_1fem_1_1Solver.html#826f03c2fb925f35e59cdf5a9d4f241f">ApplyBC</a>(); <a name="l00193"></a>00193 } <a name="l00194"></a>00194 <a name="l00201"></a>00201 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#c7b511053301115cf987b7bbdba6ebc1">AssembleElementMatrix</a>(<a class="code" href="classitk_1_1fem_1_1Element.html" title="Abstract base element class.">Element::Pointer</a> e); <a name="l00202"></a>00202 <a name="l00210"></a>00210 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#4e4526c203e081c8f4b2f8b4adf7f5da">AssembleLandmarkContribution</a>(<a class="code" href="classitk_1_1fem_1_1Element.html" title="Abstract base element class.">Element::Pointer</a> e, <span class="keywordtype">float</span>); <a name="l00211"></a>00211 <a name="l00225"></a>00225 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#826f03c2fb925f35e59cdf5a9d4f241f">ApplyBC</a>(<span class="keywordtype">int</span> dim=0, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrix=0); <a name="l00226"></a>00226 <a name="l00234"></a>00234 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#53397f2e367f07721f566edf6ee7ed3f">AssembleF</a>(<span class="keywordtype">int</span> dim=0); <a name="l00235"></a>00235 <a name="l00239"></a>00239 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#b79d74274c60fde250264200bc47a76a">DecomposeK</a>( <span class="keywordtype">void</span> ); <a name="l00240"></a>00240 <a name="l00244"></a>00244 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#51a30382e69273da15b8f8967eb58f96">Solve</a>( <span class="keywordtype">void</span> ); <a name="l00245"></a>00245 <a name="l00250"></a>00250 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#ceddbdcb541ac5e99eb41fd8794b1275">UpdateDisplacements</a>( <span class="keywordtype">void</span> ); <a name="l00251"></a>00251 <a name="l00252"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#0f341208c8516aa18a2194133edbfc02">00252</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">Float</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#0f341208c8516aa18a2194133edbfc02">GetSolution</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i,<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> which=0) <a name="l00253"></a>00253 { <a name="l00254"></a>00254 <span class="keywordflow">return</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#22ce8e0f6c63c72542aa83a4055d76e4">m_ls</a>-><a class="code" href="classitk_1_1fem_1_1Solution.html#4d4ea59c159d3d809bd5f7cb83f44818">GetSolutionValue</a>(i,which); <a name="l00255"></a>00255 } <a name="l00256"></a>00256 <a name="l00257"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#b4976b11a31382a3e6e2f34a990a44fb">00257</a> <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#b4976b11a31382a3e6e2f34a990a44fb">GetNumberOfDegreesOfFreedom</a>( <span class="keywordtype">void</span> ) <a name="l00258"></a>00258 { <a name="l00259"></a>00259 <span class="keywordflow">return</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#0a7e098e5ae7324ae62afcdb62d391d2">NGFN</a>; <a name="l00260"></a>00260 } <a name="l00261"></a>00261 <a name="l00263"></a>00263 <a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">Float</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#fe57d9dd39f8ea6d748f53f3b8131c56">GetDeformationEnergy</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> SolutionIndex=0); <a name="l00264"></a>00264 <a name="l00265"></a>00265 <span class="keyword">public</span>: <a name="l00270"></a>00270 <a class="code" href="classitk_1_1fem_1_1Solver.html#7983575d156eebdf44afbb04844bf692">Solver</a>(); <a name="l00271"></a>00271 <a name="l00275"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#b6c84cc43b52be9eb0eafebc449bd35e">00275</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#b6c84cc43b52be9eb0eafebc449bd35e">~Solver</a>() {} <a name="l00276"></a>00276 <a name="l00291"></a>00291 <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#b6f9f5c887c1a4468a80d99006f5f694">SetLinearSystemWrapper</a>(<a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper::Pointer</a> ls); <a name="l00292"></a>00292 <a name="l00298"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#fe05ec0b4b1b5a52a5bfb712fc2a38ff">00298</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper::Pointer</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#fe05ec0b4b1b5a52a5bfb712fc2a38ff">GetLinearSystemWrapper</a>() { <span class="keywordflow">return</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#22ce8e0f6c63c72542aa83a4055d76e4">m_ls</a>; } <a name="l00299"></a>00299 <a name="l00304"></a>00304 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#78a8a446bf308b3ea8cc4c85ce049738">InitializeLinearSystemWrapper</a>(<span class="keywordtype">void</span>); <a name="l00305"></a>00305 <a name="l00309"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#75f96c6a8282cdc4728e481a57ca43cd">00309</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">Float</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#75f96c6a8282cdc4728e481a57ca43cd">GetTimeStep</a>( <span class="keywordtype">void</span> )<span class="keyword"> const </span>{ <span class="keywordflow">return</span> 0.0; } <a name="l00310"></a>00310 <a name="l00316"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#0b12daf400b92c77a18dba496e7df108">00316</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#0b12daf400b92c77a18dba496e7df108">SetTimeStep</a>(<a class="code" href="classitk_1_1fem_1_1Solver.html#335c6c3cfdae6e714302d4da3e9b3c37">Float</a>) {} <a name="l00317"></a>00317 <a name="l00318"></a>00318 <span class="keyword">protected</span>: <a name="l00319"></a>00319 <a name="l00323"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#0a7e098e5ae7324ae62afcdb62d391d2">00323</a> <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#0a7e098e5ae7324ae62afcdb62d391d2">NGFN</a>; <a name="l00324"></a>00324 <a name="l00329"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#7ad5f5085b838f93bf12e9a1736d9b44">00329</a> <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> <a class="code" href="classitk_1_1fem_1_1Solver.html#7ad5f5085b838f93bf12e9a1736d9b44">NMFC</a>; <a name="l00330"></a>00330 <a name="l00332"></a><a class="code" href="classitk_1_1fem_1_1Solver.html#22ce8e0f6c63c72542aa83a4055d76e4">00332</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper::Pointer</a> <a class="code" href="classitk_1_1fem_1_1Solver.html#22ce8e0f6c63c72542aa83a4055d76e4">m_ls</a>; <a name="l00333"></a>00333 <a name="l00334"></a>00334 <span class="keyword">private</span>: <a name="l00335"></a>00335 <a name="l00339"></a>00339 <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html" title="LinearSystemWrapper class that uses VNL numeric library functions to define a sparse...">LinearSystemWrapperVNL</a> m_lsVNL; <a name="l00340"></a>00340 <a name="l00346"></a>00346 <a class="code" href="classitk_1_1SmartPointer.html">InterpolationGridType::Pointer</a> m_InterpolationGrid; <a name="l00347"></a>00347 <a name="l00348"></a>00348 }; <a name="l00349"></a>00349 <a name="l00350"></a>00350 }} <span class="comment">// end namespace itk::fem</span> <a name="l00351"></a>00351 <a name="l00352"></a>00352 <span class="preprocessor">#endif // #ifndef __itkFEMSolver_h</span> <a name="l00353"></a>00353 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 22:42:07 2009 for ITK by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen" align="middle" border=0 width=110 height=53> </a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>, © 1997-2000</small></address> </body> </html>