Point.h

/* -*- mode: c++ -*- */
/*
 * \@copyright This File is part of the Board library which is
 * licensed under the terms of the GNU Lesser General Public Licence.
 * See the LICENCE file for further details.
 */
#ifndef _BOARD_POINT_H_
#define _BOARD_POINT_H_

#include <cmath>
#include "DGtal/io/Color.h"

namespace LibBoard {


struct Point {

  double x;      
  double y;      
  Point():x(0.0),y(0.0) { } 

  Point( const Point & other ):x(other.x),y(other.y) { } 

  Point( double xc, double yc ):x(xc),y(yc) { } 
  
  inline Point & rotate( double angle );

  inline void get( double & xout, double & yout ) const;
  
  inline Point rotated( double angle ) const;

  inline Point & rotate( double angle, const Point & center );
  
  inline Point & rotated( double angle, const Point & center ) const;

  inline Point & operator+=( const Point & other );

  inline Point & operator-=( const Point & other );

  inline Point & operator*=( double s );
  
  inline Point & operator/=( double s );

  inline Point operator-();
  
  inline double norm() const;
  
};

inline void
Point::get( double & xout, double & yout ) const
{
  xout = x;
  yout = y;
}

inline Point
operator+( const Point & a, const Point & b )
{
  return Point( a.x + b.x, a.y + b.y );
}

inline Point
operator-( const Point & a, const Point & b )
{
  return Point( a.x - b.x, a.y - b.y );
}

inline double
operator*( const Point & a, const Point & b )
{
  return a.x * b.x + a.y * b.y;
}

inline Point
operator*( const Point & p, double s  )
{
  return Point( p.x * s, p.y * s );
}

inline Point
operator*( double s, const Point & p )
{
  return Point( s * p.x, s * p.y );
}

inline Point
operator/( const Point & p, double s )
{
  return Point( p.x / s, p.y / s );
}

inline Point &
Point::operator+=( const Point & other )
{
  x += other.x;
  y += other.y;
  return *this;
}

inline Point &
Point::operator-=( const Point & other )
{
  x -= other.x;
  y -= other.y;
  return *this;
}

inline Point &
Point::operator*=( double s )
{
  x *= s;
  y *= s;
  return *this;
}

inline Point &
Point::operator/=( double s )
{
  x /= s;
  y /= s;
  return *this;
}

inline bool
operator==( const Point & a, const Point & b )
{
  return ( a.x == b.x ) && ( a.y == b.y ) ;
}

inline bool
operator!=( const Point & a, const Point & b )
{
  return ( a.x != b.x ) || ( a.y != b.y ) ;
}

Point &
Point::rotate( double angle )
{
  double newx = cos( angle ) * Point::x - sin( angle ) * Point::y; 
  double newy = sin( angle ) * Point::x + cos( angle ) * Point::y;
  x = newx;
  y = newy;
  return *this;
}

Point
Point::rotated( double angle ) const
{
  return Point(*this).rotate( angle );
}

Point &
Point::rotate( double angle, const Point & center )
{
  (*this) -= center;
  (*this).rotate( angle );
  (*this) += center;
  return *this;
}

Point &
Point::rotated( double angle, const Point & center ) const
{
  return Point(*this).rotate( angle, center );
}

double
Point::norm() const
{
  return sqrt( x*x + y*y );
}

Point Point::operator-()
{
  return Point( -x, -y );
}

} // mamespace BoardLib

#endif // _POINT_H_