itkFEMPArray.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFEMPArray.h,v $
  Language:  C++
  Date:      $Date: 2003/08/06 17:15:47 $
  Version:   $Revision: 1.10.8.1 $

  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)
{

  if( gn < 0 || gn >= (int)this->size() )
    {
    throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find() const",typeid(T).name(),gn);
    }
      
  typedef typename Superclass::value_type     ValueType;
  const ValueType & value = this->operator[](gn);

  if( value->GN == gn )
    {
    return &(*(value));
    }
      
  typedef typename Superclass::iterator Iterator;
  Iterator it   = this->begin(); 
  Iterator iend = this->end();
  while( it != iend )
    {
    if( (*it)->GN == gn )
      {
      break;
      }
    it++;
    }  

  if( it == this->end() )
    {
    throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find() const",typeid(T).name(),gn);
    }

  return &(*(*it));

}




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

  if( gn < 0 || gn >= (int)this->size() )
    {
    throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find() const",typeid(T).name(),gn);
    }
      
  typedef typename Superclass::value_type     ValueType;
  const ValueType & value = this->operator[](gn);

  if( value->GN == gn )
    {
    return &(*(value));
    }
      
  typedef typename Superclass::const_iterator ConstIterator;
  ConstIterator it   = this->begin(); 
  ConstIterator iend = this->end();
  while( it != iend )
    {
    if( (*it)->GN == gn )
      {
      break;
      }
    it++;
    }  

  if( it == this->end() )
    {
    throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"FEMPArray::Find() const",typeid(T).name(),gn);
    }

  return &(*(*it));

}




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

---