itkFEMElementStd.h

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFEMElementStd.h,v $
  Language:  C++
  Date:      $Date: 2002/06/10 23:09:20 $
  Version:   $Revision: 1.15 $

  Copyright (c) 2002 Insight Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/

#ifndef __itkFEMElementStd_h
#define __itkFEMElementStd_h

#include "itkFEMElementBase.h"

namespace itk {
namespace fem {




template<unsigned int VNumberOfNodes, unsigned int VNumberOfSpatialDimensions, class TBaseClass=Element>
class ElementStd : public TBaseClass
{
FEM_ABSTRACT_CLASS(ElementStd,TBaseClass)
public:

// FIXME: Add concept cheking for TBaseClass, and TPointClass

  // Repeat typedefs and enums from parent class
  typedef typename Superclass::Float Float;
  typedef typename Superclass::MatrixType MatrixType;
  typedef typename Superclass::VectorType VectorType;
  typedef typename Superclass::LoadType LoadType;
  typedef typename Superclass::LoadPointer LoadPointer;
  typedef typename Superclass::NodeIDType NodeIDType;
  typedef typename Superclass::DegreeOfFreedomIDType DegreeOfFreedomIDType;
  typedef typename Superclass::Node Node;
  enum{ InvalidDegreeOfFreedomID = Superclass::InvalidDegreeOfFreedomID };

  enum { NumberOfNodes=VNumberOfNodes };

  enum { NumberOfSpatialDimensions=VNumberOfSpatialDimensions };




  ElementStd();



  /*
   * Methods that define the geometry of an element
   */
  virtual unsigned int GetNumberOfNodes( void ) const
  { return NumberOfNodes; }

  virtual NodeIDType GetNode(unsigned int n) const
  {
    if(n>=NumberOfNodes)
    {
      return 0;
    }
    return this->m_node[n];
  }

  virtual void SetNode(unsigned int n, NodeIDType node)
  {
    if(n>=NumberOfNodes) { return; }
    this->m_node[n]=node;
  }

  virtual const VectorType& GetNodeCoordinates( unsigned int n ) const
  {
    return m_node[n]->GetCoordinates();
  }

  virtual unsigned int GetNumberOfSpatialDimensions() const
  {
    return NumberOfSpatialDimensions;
  }




  /*
   * Methods related to I/O
   */

  virtual void Read( std::istream&, void* info );

  virtual void Write( std::ostream& f ) const;




protected:

  NodeIDType m_node[NumberOfNodes];

};




#ifdef _MSC_VER
// Declare a static dummy function to prevent a MSVC 6.0 SP5 from crashing.
// I have no idea why things don't work when this is not declared, but it
// looks like this declaration makes compiler forget about some of the
// troubles it has with templates.
static void Dummy( void );
#endif // #ifdef _MSC_VER

}} // end namespace itk::fem

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkFEMElementStd.txx"
#endif

#endif // #ifndef __itkFEMElementStd_h

---