<!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: itkFEMElement2DC0QuadraticTriangular.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>itkFEMElement2DC0QuadraticTriangular.h</h1><a href="itkFEMElement2DC0QuadraticTriangular_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: itkFEMElement2DC0QuadraticTriangular.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-28 21:44:35 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.7 $</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 __itkFEMElement2DC0QuadraticTriangular_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span><span class="preprocessor">#define __itkFEMElement2DC0QuadraticTriangular_h</span> <a name="l00020"></a>00020 <span class="preprocessor"></span> <a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkFEMElementStd_8h.html">itkFEMElementStd.h</a>"</span> <a name="l00022"></a>00022 <a name="l00023"></a>00023 <span class="keyword">namespace </span>itk { <a name="l00024"></a>00024 <span class="keyword">namespace </span>fem { <a name="l00025"></a>00025 <a name="l00030"></a><a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html">00030</a> <span class="keyword">class </span><a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html" title="3-noded, quadratic, C0 continuous finite element in 2D space.">Element2DC0QuadraticTriangular</a> : <span class="keyword">public</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html" title="Implements standard node management in the element classes.">ElementStd</a><6,2> <a name="l00031"></a>00031 { <a name="l00032"></a>00032 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html">ElementStd<3,2></a> <a class="code" href="classitk_1_1fem_1_1ElementStd.html">TemplatedParentClass</a>; <a name="l00033"></a><a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#f9e040fbdd21edd9b63dfd4880030f46">00033</a> <a class="code" href="itkFEMMacro_8h.html#fa7d760fae69e20af6c94904fb3366fa" title="If defined, FEM classes will use smart pointers.">FEM_ABSTRACT_CLASS</a>( <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html" title="3-noded, quadratic, C0 continuous finite element in 2D space.">Element2DC0QuadraticTriangular</a>, <a class="code" href="classitk_1_1fem_1_1ElementStd.html">TemplatedParentClass</a> ) <a name="l00034"></a>00034 public: <a name="l00035"></a>00035 <a name="l00037"></a>00037 <a name="l00041"></a><a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#6a2703c162c928b852b37ba6fdfef3e64c62b993e73f132ccf412777ca07b141">00041</a> enum { <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#6a2703c162c928b852b37ba6fdfef3e64c62b993e73f132ccf412777ca07b141">DefaultIntegrationOrder</a> = 2 }; <a name="l00042"></a>00042 <a name="l00043"></a>00043 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#0f1e6a9410d35cfc1cbe08c00e663bd8">GetIntegrationPointAndWeight</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a>& pt, <a class="code" href="classitk_1_1fem_1_1ElementStd.html#48b84bd9306df346dc091abf8607cc92">Float</a>& w, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> order) <span class="keyword">const</span>; <a name="l00044"></a>00044 <a name="l00045"></a>00045 <span class="keyword">virtual</span> <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#97264d1596e5276d030feb6211baf92d">GetNumberOfIntegrationPoints</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> order) <span class="keyword">const</span>; <a name="l00046"></a>00046 <a name="l00048"></a>00048 <a name="l00052"></a>00052 <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#c078cab1e5ac523f225cba115eb077cf">ShapeFunctions</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a>& pt ) <span class="keyword">const</span>; <a name="l00053"></a>00053 <a name="l00054"></a>00054 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#3f5dc6c513867f5d9360e57318f32413">ShapeFunctionDerivatives</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a>& pt, <a class="code" href="classitk_1_1fem_1_1ElementStd.html#c0c59eed49f0d89a7245da721c4c905d">MatrixType</a>& shapeD ) <span class="keyword">const</span>; <a name="l00055"></a>00055 <a name="l00056"></a>00056 <span class="comment">// FIXME: Write a proper implementation</span> <a name="l00057"></a><a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#f866d5da1b6229692564afff0bd9736b">00057</a> <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#f866d5da1b6229692564afff0bd9736b">GetLocalFromGlobalCoordinates</a>( <span class="keyword">const</span> <a class="code" href="classvnl__vector.html">VectorType</a>&, <a class="code" href="classvnl__vector.html">VectorType</a>& )<span class="keyword"> const</span> <a name="l00058"></a>00058 <span class="keyword"> </span>{ <a name="l00059"></a>00059 <span class="keywordflow">throw</span>; <a name="l00060"></a>00060 <span class="keywordflow">return</span> <span class="keyword">false</span>; <a name="l00061"></a>00061 } <a name="l00062"></a>00062 <a name="l00063"></a>00063 <span class="comment">// Since the Jacobian is not quadratic, we need to provide our</span> <a name="l00064"></a>00064 <span class="comment">// own implementation of calculating the determinant and inverse.</span> <a name="l00065"></a>00065 <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#48b84bd9306df346dc091abf8607cc92">Float</a> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#bef222daa438ebdf47d0abe23aa3cf4f">JacobianDeterminant</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a>& pt, <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#c0c59eed49f0d89a7245da721c4c905d">MatrixType</a>* pJ = 0 ) <span class="keyword">const</span>; <a name="l00066"></a>00066 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1Element2DC0QuadraticTriangular.html#27546b8f15cf083b4168f901ffea3a98">JacobianInverse</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#237fdb10334dab686c6f3da0d872ac31">VectorType</a>& pt, <a class="code" href="classitk_1_1fem_1_1ElementStd.html#c0c59eed49f0d89a7245da721c4c905d">MatrixType</a>& invJ, <span class="keyword">const</span> <a class="code" href="classitk_1_1fem_1_1ElementStd.html#c0c59eed49f0d89a7245da721c4c905d">MatrixType</a>* pJ = 0 ) <span class="keyword">const</span>; <a name="l00067"></a>00067 <a name="l00071"></a>00071 <span class="preprocessor">#ifdef FEM_BUILD_VISUALIZATION</span> <a name="l00072"></a>00072 <span class="preprocessor"></span> <span class="keywordtype">void</span> Draw(CDC* pDC, <a class="code" href="classitk_1_1fem_1_1Solution.html" title="Provides functions to access the values of the solution vector.">Solution::ConstPointer</a> sol) <span class="keyword">const</span>; <a name="l00073"></a>00073 <span class="preprocessor">#endif</span> <a name="l00074"></a>00074 <span class="preprocessor"></span> <a name="l00075"></a>00075 }; <a name="l00076"></a>00076 <a name="l00077"></a>00077 }} <span class="comment">// end namespace itk::fem</span> <a name="l00078"></a>00078 <a name="l00079"></a>00079 <span class="preprocessor">#endif // #ifndef __itkFEMElement2DC0QuadraticTriangular_h</span> <a name="l00080"></a>00080 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 22:41:04 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>