geometry.hpp

Go to the documentation of this file.

#ifndef H_GEO
#define H_GEO

#include <vector>
#include <iostream>
#include <gmpxx.h>
#include <boost/rational.hpp>
#include "integer_type.hpp"

#if USE_BIGINT

typedef mpz_class I;
typedef mpq_class rat;

#else 

typedef int I;
typedef boost::rational<int> rat;

#endif

class Point {
private:
  rat x; 
  rat y; 

public:
 
  Point() : x(), y() {};

  Point(const Point& p) : x(p.x), y(p.y) {};
  
  Point(const rat& xx, const rat& yy) : x(xx), y(yy) {};
  rat get_abscissa() const;
  rat get_ordinate() const;
  void assign(const rat&,const rat&);
};


enum slope_t {undef, rational, infty};

struct Slope {
  slope_t type; 
  rat value; 
  void make(const rat&,const rat&,const rat&,const rat&);
  bool operator<(const Slope&) const;
};

class Segment {
private:
  Point left; 
  Point right; 
  Slope slope; 

public:

  Segment() : left(), right() { slope.type = undef; };

  Segment(const Segment& seg) : left(seg.left), right(seg.right) {
    slope.type = seg.slope.type;
    if(slope.type == rational)
      slope.value = seg.slope.value;
  };

  Segment(const Point& a, const Point& b) : left(a), right(b) {
    slope.make(a.get_abscissa(), a.get_ordinate(), b.get_abscissa(), b.get_ordinate());
  };

  Segment(const rat& xa, const rat& ya, const rat& xb, const rat& yb)
    : left(xa, ya), right (xb, yb) {
    slope.make(xa,ya,xb,yb);
  };
  Point get_left() const;
  Point get_right() const;
  rat high(const Point&) const;
  bool less(const Segment&, const Point&) const;
  bool intersect(const Segment&, Point&) const;
  bool is_rend(const Point&) const;
  bool is_lend(const Point&) const;
  bool is_in(const Point&) const;
};

// Pretty-printing
std::ostream& operator<<(std::ostream&, const Segment&);
std::ostream& operator<<(std::ostream&, const Segment*&);
std::ostream& operator<<(std::ostream&, const Point&);
//std::ostream& operator<<(std::ostream&, const mpz_class&);
//std::ostream& operator<<(std::ostream&, const mpq_class&);

class Event {
private:
  Point point; 
public:
  Event() : point() {};

  Event(const Event& e) : point(e.point) {};

  Event(const Point& p) : point(p) {};
  Point get_point() const;
  bool operator<(const Event&) const;
};

#endif //H_GEO