wxmspherical.h

Go to the documentation of this file.

// WARPXM includes
#include "warpxm/warpxm_config.h" // for pi
 
// STL includes
#include <iostream>
#include <cmath>
#include <vector>
#include <sstream>
#include <limits>
 
// typedef double real;
 
// namespace wxm
// {
//     constexpr real pi = 3.14159265358979323846264338327950288419716939937510582;
// }
 
// forward declaration
class SphericalPositionVector;
 
class SphericalHemisphereConverter
{
public:
    SphericalHemisphereConverter();
 
    SphericalPositionVector convert(const SphericalPositionVector& vecObj) const;
 
    // protected:
    // const SphericalPositionVector& _vecObj;
};
 
class SphericalPositionVector
{
public:
    SphericalPositionVector(const std::vector<real> n);
 
    SphericalPositionVector(const real r, const real theta, const real phi);
 
    void CalculateFromRTHETAPHI(real r, real theta, real phi);
 
    // const std::vector<real>  get_vec() { return _vec; } const
 
    real get_r() const
    {
        return _r;
    }
 
    real get_theta() const
    {
        return _theta;
    }
 
    real get_phi() const
    {
        return _phi;
    }
 
    real get_x() const
    {
        return _vec[0];
    }
 
    real get_y() const
    {
        return _vec[1];
    }
 
    real get_z() const
    {
        return _vec[2];
    }
 
    // /**
    //  * setter of x
    //  */
    // const real set_x() const { return _vec[0]; }
 
    // /**
    //  * setter of y
    //  */
    // const real set_y() const { return _vec[1]; }
 
    // /**
    //  * setter of z
    //  */
    // const real set_z() const { return _vec[2]; }
 
    void print(std::stringstream& ss);
 
protected:
    void set_r(real r)
    {
        _r = r;
    }
 
    void set_theta(real theta)
    {
        _theta = theta;
    }
 
    void set_phi(real phi)
    {
        _phi = phi;
    }
 
    real calculate_r();
 
    real calculate_theta();
 
    real calculate_phi();
 
    real calculate_x();
 
    real calculate_y();
 
    real calculate_z();
 
    std::vector<real> calculate_vec();
 
    // attributes
    std::vector<real> _vec;
    real _r;
    real _theta;
    real _phi;
};