coord.h

Go to the documentation of this file.

#pragma once

#include <limits>
#include <map>
#include <string>
#include "Eigen/Dense"

namespace coordConv {
    
    const double MinParallax = 1e-7; // minimum parallax (arcsec)

    double distanceFromParallax(double parallax);

    double parallaxFromDistance(double dist);

    class Coord {
    public:
        explicit Coord(double equatAng, double polarAng, double parallax=0);

        explicit Coord(double equatAng, double polarAng, double parallax, double equatPM, double polarPM, double radVel);
        
        explicit Coord(Eigen::Vector3d const &pos);
        
        explicit Coord(Eigen::Vector3d const &pos, Eigen::Vector3d const &pm);
        
        explicit Coord();
        
        ~Coord() {};
    
        bool atInfinity() const { return _atInfinity; }

        bool atPole() const { return _atPole; }
        
        double getDistance() const { return _dist; }
        
        double getParallax() const;
    
        bool getSphPos(double &equatAng, double &polarAng) const;
        
        bool getPM(double &equatPM, double &polarPM) const;
        
        double getRadVel() const;

        Eigen::Vector3d const getVecPos() const { return _pos; }

        Eigen::Vector3d const getVecPM() const { return _pm; }
        
        bool isfinite() const;

        double angularSeparation(Coord const &coord) const;

        double orientationTo(Coord const &coord) const;
        
        Coord offset(double &toOrient, double fromOrient, double dist) const;

        bool operator==(Coord const &rhs) {
            return (this->getVecPos() == rhs.getVecPos()) && (this->getVecPM() == rhs.getVecPM());
        }

        bool operator!=(Coord const &rhs) {
            return !operator==(rhs);
        }
        
        std::string __repr__() const;

    private:
        Eigen::Vector3d _pos;   // vector position (AU)
        Eigen::Vector3d _pm;    // vector proper motion and radial velocity (AU/year)
        double _dist;           // distance (AU); a cache of _pos.norm()
        bool _atInfinity;       // true if distance far enough; a cached value
        bool _atPole;           // true if very near the pole; a cached value
        
        void _setPosFromSph(double equatAng, double polarAng, double parallax);
        
        void _setCache();
    };
    
    std::ostream &operator<<(std::ostream &out, Coord const &coord);

}