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 &operator=( const Point & other ) = default;

  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_