itkFEMPArray.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFEMPArray.h,v $
  Language:  C++
  Date:      $Date: 2002/07/21 22:18:08 $
  Version:   $Revision: 1.10 $

  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 __itkFEMPArray_h
#define __itkFEMPArray_h

#include "itkFEMP.h"
#include "itkFEMException.h"
#include <vector>

namespace itk {
namespace fem {




template<class T>
class FEMPArray : public std::vector<FEMP<T> >
{
public:
  typedef FEMPArray Self;

  typedef std::vector<FEMP<T> > Superclass;

  typedef Self* Pointer;
  typedef const Self* ConstPointer;

  typedef T ClassType;
  typedef typename ClassType::Pointer ClassTypePointer;
  typedef typename ClassType::ConstPointer ClassTypeConstPointer;

  ClassTypePointer Find(int gn);
  ClassTypeConstPointer Find(int gn) const;

  ClassTypePointer operator() (int i) 
  {
    return &(*operator[](i));
  }

  ClassTypeConstPointer operator() (int i) const 
  {  
    return &(*operator[](i));
  }


  int Renumber();

};



template<class T>
typename FEMPArray<T>::ClassTypePointer FEMPArray<T>::Find(int gn)
{

  typename Superclass::iterator i;

  if( gn<0 ||
      gn>=(int)size() ||
      ( *( i=static_cast<typename Superclass::iterator>(&this->operator[](gn)) ) )->GN!=gn )
  {
    for(i=begin(); i!=end() && (*i)->GN!=gn; i++);

    if(i==end())
    {
      throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find()",typeid(T).name(),gn);
    }
  }

  return &(*(*i));

}




template<class T>
typename FEMPArray<T>::ClassTypeConstPointer FEMPArray<T>::Find(int gn) const
{

  typename Superclass::const_iterator i;

  if( gn<0 ||
      gn>=(int)size() ||
      ( *( i=static_cast<typename Superclass::const_iterator>(&this->operator[](gn)) ) )->GN!=gn )
  {
    for(i=begin(); i!=end() && (*i)->GN!=gn; i++);
    if(i==end())
    {
      throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find() const",typeid(T).name(),gn);
    }

  }

  return &(*(*i));

}




template<class T>
int FEMPArray<T>::Renumber() 
{

  typename Superclass::iterator i;
  int j=0;

  for(i=begin(); i!=end(); i++)
  {
    (*i)->GN=j;
    j++;
  }

  return j;

}




}} // end namespace itk::fem

#endif // #ifndef __itkFEMPArray_h

---