Angle.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 1999 by Russell Edwards
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/Util/genutil/Angle.h
 
// redwards 17 Mar 99 -- Time for a definitive C++ suite of
// angle and sky coordinate functions
 
#ifndef ANGLE_H
#define ANGLE_H
 
#include <string>
#include <iostream>
#include <limits>
#include <math.h>
 
static const double MilliSecin12Hours = 4.32e7;
 
#define ANGLE_STRLEN 128
 
class Angle
{
 
  friend class AnglePair;
 
 protected:
  double radians;  // angle in radians
  double wrap_point;
  virtual void wrap(); // make sure angle is +- pi radians
  void init() {radians=0.0;wrap_point=M_PI;}
 public:
 
  static bool verbose;
 
  enum Type { Degrees, Radians, Turns };
 
  static Type default_type;
 
  Angle (const Angle& a)  
    {init(); radians=a.radians;wrap_point=a.wrap_point;}
  Angle (const double& rad = 0);
 
  virtual ~Angle () {}
 
  void setWrapPoint(double wp){wrap_point=wp; wrap();}
  double getWrapPoint() {return wrap_point;}
 
  // WvS - this not only changes the wrap-point but also changes the scale
  // void makeArctangle() { setWrapPoint(0.5*M_PI); }
 
  void setHMS (int hours, int minutes, double seconds);
  void getHMS (int& hours, int& minutes, double& seconds) const;
 
  int     setHMS (const char *);
 
 
  char*   getHMS (char* str, int places=3) const;
 
  std::string  getHMS (int places = 3) const;
 
  void setDMS (int degrees, int minutes, double seconds);
  void getDMS (int& degrees, int& minutes, double& seconds) const;
 
  int     setDMS (const char *);
 
 
  char*   getDMS (char* str, int places=3) const;
 
  std::string  getDMS (int places = 3) const;
 
  // ms is given in milliseconds of an hour
  void  setRadMS(long int ms)
    { radians = double(ms) * M_PI / MilliSecin12Hours; };
  long  getRadMS() const
    { return long(radians * MilliSecin12Hours / M_PI); };
 
  void setDegrees(double deg)
    { radians = deg * M_PI/180.0; wrap();};
 
  double getDegrees() const
    { return radians * 180.0/M_PI; };
 
  void setTurns(double turns)
    { radians = turns * 2.0*M_PI; wrap();};
 
  double getTurns() const
    { return radians / (2.0*M_PI); };
 
  void setRadians(double rad)
    { radians = rad; wrap(); };
 
  double getRadians() const
    { return radians; };
  
  void setDegMS (double deg_mmss);
  double getDegMS () const;
 
  void setHourMS (double hour_mmss);
  double getHourMS () const;
 
  void setTurnMS (double turn_mmss);
  double getTurnMS () const;
 
  Angle & operator= (const Angle & a);
  Angle & operator= (const double &val);
 
  Angle & operator+= (const Angle & a);
  Angle & operator-= (const Angle & a);
  Angle & operator+= (const double & d);
  Angle & operator-= (const double & d);
  Angle & operator*= (const double & d);
  Angle & operator/= (const double & d);
  friend const Angle operator + (const Angle &, const Angle &);
  friend const Angle operator - (const Angle &, const Angle &);
  friend const Angle operator + (const Angle &, double);
  friend const Angle operator - (const Angle &, double);
  friend const Angle operator / (const Angle &, double);
  friend const Angle operator * (const Angle &, double);
  friend const Angle operator - (const Angle &);
 
  friend double operator * (const Angle &, const Angle &);
 
  friend int operator > (const Angle &, const Angle &);
  friend int operator < (const Angle &, const Angle &);
  friend int operator >= (const Angle &, const Angle &);
  friend int operator <= (const Angle &, const Angle &);
  friend int operator == (const Angle &, const Angle &);
  friend int operator != (const Angle &, const Angle &);
 
  inline friend double cast_double(const Angle &a) {return a.radians;}
};
 
std::ostream& operator << (std::ostream& os, const Angle& angle);
std::istream& operator >> (std::istream& is, Angle& angle);
 
namespace std {
 
  // specialize numeric_limits for Angle class
  template<>
  class numeric_limits<Angle> : public numeric_limits<double> {
  };
}
 
// More Arctangle stuff... NOTE: remember the Angle = operator
// DOESN'T set wrap_point, so only use these to e.g. pass directly
// to a function
//inline Angle Arctangle() {Angle a; a.makeArctangle(); return a;}
//inline Angle Arctangle(const Angle&a1) {Angle a(a1); a.makeArctangle(); return a;}
//inline Angle Arctangle(double rad) {Angle a(rad); a.makeArctangle(); return a;}
 
 
 
class AnglePair
{
 public:
  Angle angle1, angle2;
 
  // both HMS and DMS in one string
  int setHMSDMS (const char* coordstr);
 
  // HMS and DMS in separate strings
  int    setHMSDMS (const char *, const char *);
  int    setHMSDMS (const std::string&, const std::string&);
  void   getHMSDMS (char* s1, char* s2, int places1=3, int places2=2) const;
  std::string getHMSDMS (int places1 = 3, int places2 = 2, bool space = true) const;
 
  void   setDegrees (double, double);
  void   getDegrees (double *, double *) const;
  std::string getDegrees () const;
  void   setRadians (double, double);
  void   getRadians (double *, double *) const;
  std::string getRadians () const;
 
  Angle angularSeparation (const AnglePair& other) const;
 
  // set the angles in radians from 
  // az and ze in milliseconds of an hour
  void setRadMS(long int , long int );
  void getRadMS(long int*, long int*);
 
  AnglePair & operator=  (const AnglePair &);
  AnglePair & operator*= (const double);
 
  // default constructor
  AnglePair (const double = 0.0, const double = 0.0);
 
  // copy constructor
  AnglePair (const AnglePair&);
 
  AnglePair (const Angle &, const Angle &);
  AnglePair (const std::string& astr);
 
  friend AnglePair operator * (const AnglePair&, const double);
 
  friend int operator == (const AnglePair &, const AnglePair &);
  friend int operator != (const AnglePair &, const AnglePair &);
 
  friend std::ostream& operator<< (std::ostream&, const AnglePair&);
};
 
 
// redwards : trig functions
double sin(const Angle&);
double cos(const Angle&);
double tan(const Angle&);
// Unfortunately can't overload atan as only the return type differs
// from the standard library function
Angle arctan(double);
Angle arctan(double y, double x); // returns atan y/x
 
#endif //ANGLE_H