lemon-main: comparison lemon/dim2.h

equal deleted inserted replaced

-:0ec0a16fabe7
+:d87a1b88ab1a
-/* -*- C++ -*-
+/* -*- mode: C++; indent-tabs-mode: nil; -*-
 *
-* This file is a part of LEMON, a generic C++ optimization library
+* This file is a part of LEMON, a generic C++ optimization library.
 *
 * Copyright (C) 2003-2008
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
 *
 #include <iostream>
 #include <lemon/bits/utility.h>
 ///\ingroup misc
 ///\file
 ///\brief A simple two dimensional vector and a bounding box implementation
 ///
 /// The class \ref lemon::dim2::Point "dim2::Point" implements
 /// a two dimensional vector with the usual operations.
 ///
 /// The class \ref lemon::dim2::BoundingBox "dim2::BoundingBox"
 typedef T Value;
 ///First coordinate
 T x;
 ///Second coordinate
 T y;
 ///Default constructor
 Point() {}
 ///Construct an instance from coordinates
 Point(T a, T b) : x(a), y(b) { }
 ///Returns the dimension of the vector (i.e. returns 2).
 ///The dimension of the vector.
 ///This function always returns 2.
 int size() const { return 2; }
 ///Subscripting operator
 ///\c p[0] is \c p.x and \c p[1] is \c p.y
 ///Give back the square of the norm of the vector
 T normSquare() const {
 return x*x+y*y;
 }
 ///Increment the left hand side by \c u
 Point<T>& operator +=(const Point<T>& u) {
 x += u.x;
 y += u.y;
 return *this;
 }
 ///Decrement the left hand side by \c u
 Point<T>& operator -=(const Point<T>& u) {
 x -= u.x;
 y -= u.y;
 return *this;
 Point<T>& operator /=(const T &u) {
 x /= u;
 y /= u;
 return *this;
 }
 ///Return the scalar product of two vectors
 T operator *(const Point<T>& u) const {
 return x*u.x+y*u.y;
 }
 ///Return the sum of two vectors
 Point<T> operator+(const Point<T> &u) const {
 Point<T> b=*this;
 return b+=u;
 }
 return  (x!=u.x) || (y!=u.y);
 }
 };
 ///Return a Point
 ///Return a Point.
 ///\relates Point
 template <typename T>
 inline Point<T> makePoint(const T& x, const T& y) {
 inline Point<T> rot270(const Point<T> &z)
 {
 return Point<T>(z.y,-z.x);
 }
 /// A class to calculate or store the bounding box of plainvectors.
 /// A class to calculate or store the bounding box of plainvectors.
 ///
 template<typename T>
 class BoundingBox {
 Point<T> bottom_left, top_right;
 bool _empty;
 public:
 ///Default constructor: creates an empty bounding box
 BoundingBox() { _empty = true; }
 ///Construct an instance from one point
 BoundingBox(Point<T> a) { bottom_left=top_right=a; _empty = false; }
 ///Construct an instance from two points
 ///Construct an instance from two points.
 ///\param a The bottom left corner.
 ///\param b The top right corner.
 ///\warning The coordinates of the bottom left corner must be no more
 ///than those of the top right one.
 BoundingBox(Point<T> a,Point<T> b)
 {
-	bottom_left=a;
+bottom_left=a;
-	top_right=b;
+top_right=b;
-	_empty = false;
+_empty = false;
 }
 ///Construct an instance from four numbers
 ///Construct an instance from four numbers.
 ///\param l The left side of the box.
 ///\param b The bottom of the box.
 ///\param r The right side of the box.
 ///\param t The top of the box.
 ///\warning The left side must be no more than the right side and
 ///bottom must be no more than the top.
 BoundingBox(T l,T b,T r,T t)
 {
-	bottom_left=Point<T>(l,b);
+bottom_left=Point<T>(l,b);
-	top_right=Point<T>(r,t);
+top_right=Point<T>(r,t);
-	_empty = false;
+_empty = false;
 }
 ///Return \c true if the bounding box is empty.
 ///Return \c true if the bounding box is empty (i.e. return \c false
 ///if at least one point was added to the box or the coordinates of
 ///the box were set).
 ///
 ///The coordinates of an empty bounding box are not defined.
 bool empty() const {
 return _empty;
 }
 ///Make the BoundingBox empty
 void clear() {
 _empty=1;
 }
 ///Set the bottom left corner of the box
 ///Set the bottom left corner of the box.
 ///It should only be used for non-empty box.
 void bottomLeft(Point<T> p) {
-	bottom_left = p;
+bottom_left = p;
 }
 ///Give back the top right corner of the box
 ///Give back the top right corner of the box.
 ///Set the top right corner of the box
 ///Set the top right corner of the box.
 ///It should only be used for non-empty box.
 void topRight(Point<T> p) {
-	top_right = p;
+top_right = p;
 }
 ///Give back the bottom right corner of the box
 ///Give back the bottom right corner of the box.
 ///Set the bottom right corner of the box
 ///Set the bottom right corner of the box.
 ///It should only be used for non-empty box.
 void bottomRight(Point<T> p) {
-	top_right.x = p.x;
+top_right.x = p.x;
-	bottom_left.y = p.y;
+bottom_left.y = p.y;
 }
 ///Give back the top left corner of the box
 ///Give back the top left corner of the box.
 ///If the bounding box is empty, then the return value is not defined.
 Point<T> topLeft() const {
 ///Set the top left corner of the box
 ///Set the top left corner of the box.
 ///It should only be used for non-empty box.
 void topLeft(Point<T> p) {
-	top_right.y = p.y;
+top_right.y = p.y;
-	bottom_left.x = p.x;
+bottom_left.x = p.x;
 }
 ///Give back the bottom of the box
 ///Give back the bottom of the box.
 ///Set the bottom of the box
 ///Set the bottom of the box.
 ///It should only be used for non-empty box.
 void bottom(T t) {
-	bottom_left.y = t;
+bottom_left.y = t;
 }
 ///Give back the top of the box
 ///Give back the top of the box.
 ///Set the top of the box
 ///Set the top of the box.
 ///It should only be used for non-empty box.
 void top(T t) {
-	top_right.y = t;
+top_right.y = t;
 }
 ///Give back the left side of the box
 ///Give back the left side of the box.
 ///If the bounding box is empty, then the return value is not defined.
 T left() const {
 return bottom_left.x;
 }
 ///Set the left side of the box
 ///Set the left side of the box.
 ///It should only be used for non-empty box.
 void left(T t) {
-	bottom_left.x = t;
+bottom_left.x = t;
 }
 /// Give back the right side of the box
 /// Give back the right side of the box.
 ///Set the right side of the box
 ///Set the right side of the box.
 ///It should only be used for non-empty box.
 void right(T t) {
-	top_right.x = t;
+top_right.x = t;
 }
 ///Give back the height of the box
 ///Give back the height of the box.
 else{
 return ((u.x-bottom_left.x)*(top_right.x-u.x) >= 0 &&
 (u.y-bottom_left.y)*(top_right.y-u.y) >= 0 );
 }
 }
 ///Increments a bounding box with a point
 ///Increments a bounding box with a point.
 ///
 BoundingBox& add(const Point<T>& u){
 if (top_right.x < u.x) top_right.x = u.x;
 if (top_right.y < u.y) top_right.y = u.y;
 }
 return *this;
 }
 ///Increments a bounding box to contain another bounding box
 ///Increments a bounding box to contain another bounding box.
 ///
 BoundingBox& add(const BoundingBox &u){
 if ( !u.empty() ){
 this->add(u.bottomLeft());
-	  this->add(u.topRight());
+this->add(u.topRight());
 }
 return *this;
 }
 ///Intersection of two bounding boxes
 ///Intersection of two bounding boxes.
 ///
 BoundingBox operator&(const BoundingBox& u) const {
 BoundingBox b;
 if (this->_empty || u._empty) {
-	  b._empty = true;
+b._empty = true;
-	} else {
+} else {
-	  b.bottom_left.x = std::max(this->bottom_left.x,u.bottom_left.x);
+b.bottom_left.x = std::max(this->bottom_left.x,u.bottom_left.x);
-	  b.bottom_left.y = std::max(this->bottom_left.y,u.bottom_left.y);
+b.bottom_left.y = std::max(this->bottom_left.y,u.bottom_left.y);
-	  b.top_right.x = std::min(this->top_right.x,u.top_right.x);
+b.top_right.x = std::min(this->top_right.x,u.top_right.x);
-	  b.top_right.y = std::min(this->top_right.y,u.top_right.y);
+b.top_right.y = std::min(this->top_right.y,u.top_right.y);
-	  b._empty = b.bottom_left.x > b.top_right.x ||
+b._empty = b.bottom_left.x > b.top_right.x ||
-	             b.bottom_left.y > b.top_right.y;
+b.bottom_left.y > b.top_right.y;
-	}
+}
 return b;
 }
 };//class Boundingbox
 ///\ingroup maps
 ///Map of x-coordinates of a \ref Point "Point"-map.
 ///
 template<class M>
 class XMap
 {
 M& _map;
 public:
 typedef typename M::Value::Value Value;
 ///\e
 XMap(M& map) : _map(map) {}
 Value operator[](Key k) const {return _map[k].x;}
 void set(Key k,Value v) {_map.set(k,typename M::Value(v,_map[k].y));}
 };
 ///Returns an \ref XMap class
 ///This function just returns an \ref XMap class.
 ///
 ///\ingroup maps
 ///\relates XMap
 template<class M>
 inline XMap<M> xMap(M &m)
 {
 return XMap<M>(m);
 }
 template<class M>
 inline XMap<M> xMap(const M &m)
 {
 return XMap<M>(m);
 }
 ///Constant (read only) version of \ref XMap
 ///\ingroup maps
 ///Constant (read only) version of \ref XMap
 ///
 template<class M>
 class ConstXMap
 {
 const M& _map;
 public:
 typedef typename M::Value::Value Value;
 typedef typename M::Key Key;
 ///\e
 ConstXMap(const M &map) : _map(map) {}
 Value operator[](Key k) const {return _map[k].x;}
 };
 ///Returns a \ref ConstXMap class
 ///This function just returns a \ref ConstXMap class.
 ///
 ///\ingroup maps
 ///\relates ConstXMap
 template<class M>
 inline ConstXMap<M> xMap(const M &m)
 {
 return ConstXMap<M>(m);
 }
 ///Map of y-coordinates of a \ref Point "Point"-map
 ///\ingroup maps
 ///Map of y-coordinates of a \ref Point "Point"-map.
 ///
 template<class M>
 class YMap
 {
 M& _map;
 public:
 typedef typename M::Value::Value Value;
 ///This function just returns a \ref YMap class.
 ///
 ///\ingroup maps
 ///\relates YMap
 template<class M>
 inline YMap<M> yMap(M &m)
 {
 return YMap<M>(m);
 }
 template<class M>
 inline YMap<M> yMap(const M &m)
 {
 return YMap<M>(m);
 }
 ///Constant (read only) version of \ref YMap
 ///\ingroup maps
 ///Constant (read only) version of \ref YMap
 ///
 template<class M>
 class ConstYMap
 {
 const M& _map;
 public:
 typedef typename M::Value::Value Value;
 typedef typename M::Key Key;
 ///\e
 ConstYMap(const M &map) : _map(map) {}
 Value operator[](Key k) const {return _map[k].y;}
 };
 ///Returns a \ref ConstYMap class
 ///This function just returns a \ref ConstYMap class.
 ///
 ///\ingroup maps
 ///\relates ConstYMap
 template<class M>
 inline ConstYMap<M> yMap(const M &m)
 {
 return ConstYMap<M>(m);
 }
 ///
 ///Map of the \ref Point::normSquare() "normSquare()"
 ///of a \ref Point "Point"-map.
 ///\ingroup maps
 template<class M>
 class NormSquareMap
 {
 const M& _map;
 public:
 typedef typename M::Value::Value Value;
 typedef typename M::Key Key;
 ///\e
 NormSquareMap(const M &map) : _map(map) {}
 Value operator[](Key k) const {return _map[k].normSquare();}
 };
 ///Returns a \ref NormSquareMap class
 ///This function just returns a \ref NormSquareMap class.
 ///
 ///\ingroup maps
 ///\relates NormSquareMap
 template<class M>
 inline NormSquareMap<M> normSquareMap(const M &m)
 {
 return NormSquareMap<M>(m);
 }
 /// @}
 } //namespce dim2
 } //namespace lemon
 #endif //LEMON_DIM2_H

changeset 225	c5a40fc54f1a
parent 49	9a556af88710
child 220	a5d8c039f218