point3d.h

Go to the documentation of this file.

/***************************************************************************
                          point3d.h  -  description
                             -------------------
    begin                : Tue Mar 15 2005
    copyright            : (C) 2005-2006 by Ben Swerts
    email                : bswerts@users.sourceforge.net
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/





#ifndef POINT3D_H
#define POINT3D_H

#include <cmath>


template<class T> class Point3D
{
  public:
    Point3D();                          // default constructor
    Point3D(const T x, const T y, const T z);     // constructs a point with coordinates x, y, z
    Point3D(const Point3D& p);          // copy constructor
    ~Point3D();                         // destructor

    T x() const;                        // returns the x-value
    T y() const;                        // returns the y-value
    T z() const;                        // returns the z-value
    unsigned int id() const;            // returns the point's ID (used in IsoSurface and PlotMapBase)
    bool operator==(const Point3D<double>& p) const         
    {
      if(xCoord == p.xCoord && yCoord == p.yCoord && zCoord == p.zCoord)
        return true;
      return fabs(xCoord - p.xCoord) < TOLERANCE && fabs(yCoord - p.yCoord) < TOLERANCE && fabs(zCoord - p.zCoord) < TOLERANCE;
    }
    bool operator==(const Point3D<float>& p) const          
    {
      if(xCoord == p.xCoord && yCoord == p.yCoord && zCoord == p.zCoord)
        return true;
      return fabs(xCoord - p.xCoord) < TOLERANCE && fabs(yCoord - p.yCoord) < TOLERANCE && fabs(zCoord - p.zCoord) < TOLERANCE;
    }
    bool operator==(const Point3D<unsigned int>& p) const   
    {  
      return xCoord == p.xCoord && yCoord == p.yCoord && zCoord == p.zCoord;
    }

    void setValues(const T x, const T y, const T z);           // sets the soordinates (cartesian by default)
    void setCartesian(const T x, const T y, const T z);        // sets the cartesian coordinates directly
    void setPolar(const T theta, const T phi, const T r);      // sets the coordinates from polar coordinates
    void setID(const unsigned int id);  // sets the point's ID
    void add(const Point3D& p);         // adds the point p (maybe do the same for adding a Vector3D to it)

    static const T PI;                  
    static const T DEGTORAD;            
    static const T RADTODEG;            
    static const T TOLERANCE;           

  private:
    T xCoord, yCoord, zCoord;           
    unsigned int ID;                    

};


template <class T> Point3D<T>::Point3D()
{
  setValues(0, 0, 0);
}

template <class T> Point3D<T>::Point3D(const T x, const T y, const T z)
{
  setValues(x, y, z);
}

template <class T> Point3D<T>::Point3D(const Point3D& p)
{
  xCoord = p.xCoord;
  yCoord = p.yCoord;
  zCoord = p.zCoord;
  ID = p.ID;
}

template <class T> Point3D<T>::~Point3D()
{

}

template <class T> T Point3D<T>::x() const
{
  return xCoord;
}

template <class T> T Point3D<T>::y() const
{
  return yCoord;
}

template <class T> T Point3D<T>::z() const
{
  return zCoord;
}

template <class T> unsigned int Point3D<T>::id() const
{
  return ID;
}

template <class T> void Point3D<T>::setValues(const T x, const T y, const T z)
{
  xCoord = x;
  yCoord = y;
  zCoord = z;
}

template <class T> void Point3D<T>::setCartesian(const T x, const T y, const T z)
{
  setValues(x, y, z);
}

template <class T> void Point3D<T>::setPolar(const T theta, const T phi, const T r)
{
  setValues(r*sin(theta*DEGTORAD)*cos(phi*DEGTORAD), r*sin(theta*DEGTORAD)*sin(phi*DEGTORAD), r*cos(theta*DEGTORAD));
}

template <class T> void Point3D<T>::setID(const unsigned int id)
{
  ID = id;
}

template <class T> void Point3D<T>::add(const Point3D& p)
{
  xCoord += p.xCoord;
  yCoord += p.yCoord;
  zCoord += p.zCoord;
}

#endif

---