ParametricGeoMagneticField.h

Go to the documentation of this file.

#ifndef _fwk_ParametricGeoMagneticField_h_
#define _fwk_ParametricGeoMagneticField_h_

#include <utl/SVector.h>
#include <utl/Point.h>
#include <utl/Vector.h>

#include <valarray>


namespace fwk {

  class ParametricGeoMagneticField {
  public:
    ParametricGeoMagneticField(const double baseyear, const size_t order, const double g[], const double h[]);

    utl::Vector Get(const utl::Point& where) const;
    double GetDeclination(const utl::Point& where) const;

    double GetBaseYear() const { return fBaseYear; }

    ParametricGeoMagneticField Interpolate(const double year, const ParametricGeoMagneticField& later) const;

    static ParametricGeoMagneticField IGRF2010();

  protected:
    double fBaseYear = 0;
    size_t fOrder = 0;
    std::valarray<double> fG;
    std::valarray<double> fH;

    utl::SVector<3, double> GetXYZ(const double latitude, const double longitude, const double elevation) const;

    void Grow(const size_t order);

    friend class ParametricGeoMagChange;
  };

}


#endif