geometry.hpp

Go to the documentation of this file.

#ifndef GEOMETRY_HPP
#define GEOMETRY_HPP 1
#include <vector>
#include <boost/array.hpp>
#ifdef RICH_MPI
#include "../misc/serializable.hpp"
#endif // RICH_MPI

class Vector2D
#ifdef RICH_MPI
  : public Serializable
#endif // RICH_MPI
{
public:

  Vector2D(void);

  Vector2D(double ix, double iy);

  Vector2D(const Vector2D& v);

  void Set(double ix, double iy);

  double x;

  double y;

  Vector2D& operator+=(Vector2D const& v);

  Vector2D& operator-=(Vector2D const& v);

  Vector2D& operator=(Vector2D const& v);

  Vector2D& operator*=(double s);

  void Rotate(double a);
  double distance(Vector2D const& v1) const;

#ifdef RICH_MPI

  template<class Archive>
  void serialize
  (Archive& ar, 
   const unsigned int /*version*/)
  {
    ar & x;
    ar & y;
  }

  vector<double> serialize(void) const;

  size_t getChunkSize(void) const;

  void unserialize
  (const vector<double>& data);
#endif // RICH_MPI
};

double abs(Vector2D const& v);

double fastabs(Vector2D const& v);

Vector2D operator+(Vector2D const& v1, Vector2D const& v2);

Vector2D operator-(Vector2D const& v1, Vector2D const& v2);

Vector2D operator*(double d, Vector2D const& v);

Vector2D operator*(Vector2D const& v, double d);

Vector2D operator/(Vector2D const& v, double d);

double ScalarProd(Vector2D const& v1, Vector2D const& v2);

double CalcAngle(Vector2D const& v1, Vector2D const& v2);

double Projection(Vector2D const& v1, Vector2D const& v2);

Vector2D Rotate(Vector2D const& v, double a);

Vector2D Reflect(Vector2D const& v, Vector2D const& axis);

double distance(Vector2D const& v1, Vector2D const& v2);

double CrossProduct(Vector2D const& v1, Vector2D const& v2);

Vector2D calc_mid_point(Vector2D const& v1, Vector2D const& v2);

Vector2D zcross(Vector2D const& v);

Vector2D pol2cart(double radius, double angle);

Vector2D normalize(const Vector2D& v);

double dist_sqr(const Vector2D& v);

#endif // GEOMETRY_HPP