itkFEMElementStd.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFEMElementStd.h,v $
  Language:  C++
  Date:      $Date: 2009-01-29 21:28:16 $
  Version:   $Revision: 1.17 $

  Copyright (c) Insight Software 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();


  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;
    }
 

  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