itkVector.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkVector.h,v $
  Language:  C++
  Date:      $Date: 2002/11/13 03:51:29 $
  Version:   $Revision: 1.55 $

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

#include "itkFixedArray.h"
#include "vnl/vnl_vector_ref.h"
#include "itkIndent.h"

namespace itk
{

template<class T, unsigned int NVectorDimension=3>
class Vector : public FixedArray<T,NVectorDimension>
{
public:
  typedef Vector  Self;
  typedef FixedArray<T,NVectorDimension>  Superclass;
  
  typedef T ValueType;

  itkStaticConstMacro(Dimension, unsigned int, NVectorDimension);

  typedef Self VectorType;
  
  typedef FixedArray<T, NVectorDimension>                BaseArray;
    
  static unsigned int GetVectorDimension() 
    { return NVectorDimension; }  

  void Set_vnl_vector( const vnl_vector<T> & );

  vnl_vector_ref<T> Get_vnl_vector( void );

  vnl_vector<T> Get_vnl_vector( void ) const;

  Vector(): BaseArray() { }
  Vector(const ValueType& r);
  
  Vector(const Self& r): BaseArray(r) {}
  Vector(const ValueType r[Dimension]): BaseArray(r) {}  
    
  Vector& operator= (const Self& r);
  Vector& operator= (const ValueType r[Dimension]);
    
  const Self& operator*=(const ValueType &value);

  const Self& operator/=(const ValueType &value);

  const Self& operator+=(const Self &vec);

  const Self& operator-=(const Self &vec);

  Self operator-() const;

  Self operator+(const Self &vec) const;

  Self operator-(const Self &vec) const;

  ValueType operator*(const Self &vec) const;

  Self operator*(const ValueType& val) const;

  Self operator/(const ValueType& val) const;

  bool operator==(const Self& v) const
    { return Superclass::operator==(v); }
  bool operator!=(const Self& v) const
    { return !operator==(v); }
   
  ValueType GetNorm( void ) const;

  ValueType GetSquaredNorm( void ) const; 

  void Normalize(void);

  template < typename TCoordRepB >
  void CastFrom( const Vector<TCoordRepB,NVectorDimension> & pa )
  {
    for(unsigned int i=0; i<NVectorDimension; i++ )
      {
      (*this)[i] = static_cast<T>( pa[i] );
      }
  }

};

template< class T, unsigned int NVectorDimension >  
ITK_EXPORT std::ostream& operator<<(std::ostream& os, 
                                    const Vector<T,NVectorDimension> & v); 

template< class T, unsigned int NVectorDimension >  
ITK_EXPORT std::istream& operator>>(std::istream& is, 
                                    Vector<T,NVectorDimension> & v); 

ITK_EXPORT Vector<double,3> operator^( const Vector<double,3> &,
                                       const Vector<double,3> &  );

ITK_EXPORT Vector<float,3> operator^( const Vector<float,3> &,
                                      const Vector<float,3> &  );

ITK_EXPORT Vector<int,3> operator^( const Vector<int,3> &,
                                    const Vector<int,3> &  );

} // end namespace itk
  

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


#endif 

---