#ifndef EPHEMERIS_H
#define EPHEMERIS_H

#include <newmat.h>
#include <QtCore>
#include <stdio.h>
#include <string>
#include "bnctime.h"
#include "bncconst.h"
#include "t_prn.h"
#include "gnss.h"


class t_orbCorr;
class t_clkCorr;

class t_eph {
 public:
  enum e_type {unknown, GPS, QZSS, GLONASS, Galileo, SBAS, BDS, IRNSS};
  enum e_checkState {unchecked, ok, bad, outdated, unhealthy};

  t_eph();
  virtual ~t_eph();

  virtual e_type  type() const = 0;
  virtual QString toString(double version) const = 0;
  virtual unsigned int IOD() const = 0;
  virtual unsigned int isUnhealthy() const = 0;
  virtual int     slotNum() const {return 0;}
  bncTime TOC() const {return _TOC;}
  bool    isNewerThan(const t_eph* eph) const {return earlierTime(eph, this);}
  e_checkState checkState() const {return _checkState;}
  void    setCheckState(e_checkState checkState) {_checkState = checkState;}
  t_prn   prn() const {return _prn;}
  t_irc   getCrd(const bncTime& tt, ColumnVector& xc, ColumnVector& vv, bool useCorr) const;
  void    setOrbCorr(const t_orbCorr* orbCorr);
  void    setClkCorr(const t_clkCorr* clkCorr);
  const QDateTime& receptDateTime() const {return _receptDateTime;}
    const QString receptStaID() const {return _receptStaID;}
  static QString rinexDateStr(const bncTime& tt, const t_prn& prn, double version);
  static QString rinexDateStr(const bncTime& tt, const QString& prnStr, double version);
  static bool earlierTime(const t_eph* eph1, const t_eph* eph2) {return eph1->_TOC < eph2->_TOC;}
  static bool prnSort(const t_eph* eph1, const t_eph* eph2) {return eph1->prn() < eph2->prn();}

 protected:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const = 0;
  t_prn        _prn;
  bncTime      _TOC;
  QDateTime    _receptDateTime;
  QString      _receptStaID;
  e_checkState _checkState;
  t_orbCorr*   _orbCorr;
  t_clkCorr*   _clkCorr;
};


class t_ephGPS : public t_eph {
 friend class t_ephEncoder;
 friend class RTCM3Decoder;
 public:
  t_ephGPS() {
    _clock_bias       = 0.0;
    _clock_drift      = 0.0;
    _clock_driftrate  = 0.0;
    _IODE             = 0.0;
    _Crs              = 0.0;
    _Delta_n          = 0.0;
    _M0               = 0.0;
    _Cuc              = 0.0;
    _e                = 0.0;
    _Cus              = 0.0;
    _sqrt_A           = 0.0;
    _TOEsec           = 0.0;
    _Cic              = 0.0;
    _OMEGA0           = 0.0;
    _Cis              = 0.0;
    _i0               = 0.0;
    _Crc              = 0.0;
    _omega            = 0.0;
    _OMEGADOT         = 0.0;
    _IDOT             = 0.0;
    _L2Codes          = 0.0;
    _TOEweek          = 0.0;
    _L2PFlag          = 0.0;
    _ura              = 0.0;
    _health           = 0.0;
    _TGD              = 0.0;
    _IODC             = 0.0;
    _TOT              = 0.0;
    _fitInterval      = 0.0;
    _receptStaID      = "";
  }
  t_ephGPS(float rnxVersion, const QStringList& lines);
  virtual ~t_ephGPS() {}

    virtual e_type type() const {
    switch (_prn.system()) {
      case 'J':
        return t_eph::QZSS;
      case 'I':
        return t_eph::IRNSS;
    };
    return t_eph::GPS;
  }
  virtual QString toString(double version) const;
  virtual unsigned int  IOD() const { return static_cast<unsigned int>(_IODE); }
  virtual unsigned int  isUnhealthy() const { return static_cast<unsigned int>(_health); }
  double TGD() const {return _TGD;} // Timing Group Delay (P1-P2 DCB)

 private:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;

  double  _clock_bias;      // [s]
  double  _clock_drift;     // [s/s]
  double  _clock_driftrate; // [s/s^2]

  double  _IODE;            // IODEC in case of IRNSS
  double  _Crs;             // [m]
  double  _Delta_n;         // [rad/s]
  double  _M0;              // [rad]

  double  _Cuc;             // [rad]
  double  _e;               //
  double  _Cus;             // [rad]
  double  _sqrt_A;          // [m^0.5]

  double  _TOEsec;          // [s]
  double  _Cic;             // [rad]
  double  _OMEGA0;          // [rad]
  double  _Cis;             // [rad]

  double  _i0;              // [rad]
  double  _Crc;             // [m]
  double  _omega;           // [rad]
  double  _OMEGADOT;        // [rad/s]

  double  _IDOT;            // [rad/s]
  double  _L2Codes;         // Codes on L2 channel  (not valid for IRNSS)
  double  _TOEweek;
  double  _L2PFlag;         // L2 P data flag (not valid for IRNSS and QZSS)

  mutable double  _ura;     // SV accuracy
  double  _health;          // SV health
  double  _TGD;             // [s]
  double  _IODC;            // (not valid for IRNSS)

  double  _TOT;             // Transmission time
  double  _fitInterval;     // Fit interval in hours (not valid for IRNSS)

};

class t_ephGlo : public t_eph {
 friend class t_ephEncoder;
 friend class RTCM3Decoder;
 public:
  t_ephGlo() {
    _xv.ReSize(6); _xv = 0.0;
    _gps_utc          = 0.0;
    _tau              = 0.0;
    _gamma            = 0.0;
    _tki              = 0.0;
    _x_pos            = 0.0;
    _x_velocity       = 0.0;
    _x_acceleration   = 0.0;
    _health           = 0.0;
    _y_pos            = 0.0;
    _y_velocity       = 0.0;
    _y_acceleration   = 0.0;
    _frequency_number = 0.0;
    _z_pos            = 0.0;
    _z_velocity       = 0.0;
    _z_acceleration   = 0.0;
    _E                = 0.0;
    _almanac_health   = 0.0;
    _almanac_health_availablility_indicator = 0.0;
    _additional_data_availability = 0.0;
    _tauC             = 0.0;
    _P1               = 0.0;
    _P2               = 0.0;
    _P3               = 0.0;
    _NA               = 0.0;
    _M_P              = 0.0;
    _M_l3             = 0.0;
    _M_delta_tau      = 0.0;
    _M_P4             = 0.0;
    _M_FT             = 0.0;
    _M_NT             = 0.0;
    _M_M              = 0.0;
    _M_N4             = 0.0;
    _M_tau_GPS        = 0.0;
    _M_l5             = 0.0;
    _receptStaID      = "";    
  }
  t_ephGlo(float rnxVersion, const QStringList& lines);
  virtual ~t_ephGlo() {}

  virtual e_type type() const {return t_eph::GLONASS;}
  virtual QString toString(double version) const;
  virtual unsigned int  IOD() const;
  virtual unsigned int isUnhealthy() const;
  virtual int slotNum() const {return int(_frequency_number);}

 private:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
  static ColumnVector glo_deriv(double /* tt */, const ColumnVector& xv, double* acc);

  mutable bncTime      _tt;  // time
  mutable ColumnVector _xv;  // status vector (position, velocity) at time _tt

  double  _gps_utc;
  double  _tau;              // [s]
  double  _gamma;            // [-]
  mutable double  _tki;      // message frame time

  double  _x_pos;            // [km]
  double  _x_velocity;       // [km/s]
  double  _x_acceleration;   // [km/s^2]
  double  _health;           // 0 = O.K. MSB of Bn word

  double  _y_pos;            // [km]
  double  _y_velocity;       // [km/s]
  double  _y_acceleration;   // [km/s^2]
  double  _frequency_number; // ICD-GLONASS data position

  double  _z_pos;            // [km]
  double  _z_velocity;       // [km/s]
  double  _z_acceleration;   // [km/s^2]
  double  _E;                // Age of Information [days]

  double _almanac_health;     // Cn word
  double _almanac_health_availablility_indicator;

  double _additional_data_availability;  //
  double _tauC;               // GLONASS time scale correction to UTC(SU) time [sec]
  double _P1;                 // flag of the immediate data updating [-]
  double _P2;                 // flag of oddness or evenness of the value of tb for intervals 30 or 60 minutes [-]
  double _P3;                 // flag indicating a number of satellites for which almanac is transmitted within given frame [-]
  double _NA;                 // calendar day number within the 4-year period [days]

  double _M_P;                // control segment parameter that indicates the satellite operation mode with respect of time parameters
  double _M_l3;               // health flag
  double _M_delta_tau;        // [sec]
  double _M_P4;               // flag to show that ephemeris parameters are present [-]
  double _M_FT;               // indicator for predicted satellite user range accuracy [-]
  double _M_NT;               // current date, calendar number of day within 4-year interval [days]
  double _M_M;                // type of satellite transmitting navigation signal: 0 = GLONASS, 1 = GLONASS-M satellite [-]
  double _M_N4;               // 4-year interval number starting from 1996
  double _M_tau_GPS;          // correction to GPS time relative to GLONASS time [days]
  double _M_l5;               // health flag
};

class t_ephGal : public t_eph {
 friend class t_ephEncoder;
 friend class RTCM3Decoder;
 public:
  t_ephGal() {
    _clock_bias      = 0.0;
    _clock_drift     = 0.0;
    _clock_driftrate = 0.0;
    _IODnav          = 0.0;
    _Crs             = 0.0;
    _Delta_n         = 0.0;
    _M0              = 0.0;
    _Cuc             = 0.0;
    _e               = 0.0;
    _Cus             = 0.0;
    _sqrt_A          = 0.0;
    _TOEsec          = 0.0;
    _Cic             = 0.0;
    _OMEGA0          = 0.0;
    _Cis             = 0.0;
    _i0              = 0.0;
    _Crc             = 0.0;
    _omega           = 0.0;
    _OMEGADOT        = 0.0;
    _IDOT            = 0.0;
    _TOEweek         = 0.0;
    _SISA            = 0.0;
    _E5aHS           = 0.0;
    _E5bHS           = 0.0;
    _E1_bHS          = 0.0;
    _BGD_1_5A        = 0.0;
    _BGD_1_5B        = 0.0;
    _TOT             = 0.0;
    _receptStaID     = "";    
  };
  t_ephGal(float rnxVersion, const QStringList& lines);
  virtual ~t_ephGal() {}

  virtual QString toString(double version) const;
  virtual e_type type() const {return t_eph::Galileo;}
  virtual unsigned int  IOD() const { return static_cast<unsigned long>(_IODnav); }
  virtual unsigned int  isUnhealthy() const;

 private:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;

  double  _clock_bias;       //  [s]
  double  _clock_drift;      //  [s/s]
  double  _clock_driftrate;  //  [s/s^2]

  double  _IODnav;
  double  _Crs;              //  [m]
  double  _Delta_n;          //  [rad/s]
  double  _M0;               //  [rad]

  double  _Cuc;              //  [rad]
  double  _e;                //
  double  _Cus;              //  [rad]
  double  _sqrt_A;           //  [m^0.5]

  double  _TOEsec;           //  [s]
  double  _Cic;              //  [rad]
  double  _OMEGA0;           //  [rad]
  double  _Cis;              //  [rad]

  double  _i0;               //  [rad]
  double  _Crc;              //  [m]
  double  _omega;            //  [rad]
  double  _OMEGADOT;         //  [rad/s]

  double  _IDOT;             //  [rad/s]
  double  _TOEweek;
  // spare

  mutable double  _SISA;     // Signal In Space Accuracy
  double  _E5aHS;            //  [0..3] E5a Health Status
  double  _E5bHS;            //  [0..3] E5b Health Status
  double  _E1_bHS;           //  [0..3] E1-b Health Status
  double  _BGD_1_5A;         //  group delay [s]
  double  _BGD_1_5B;         //  group delay [s]

  double  _TOT;              // [s]
  /** Data comes from I/NAV when <code>true</code> */
  bool    _inav;
  /** Data comes from F/NAV when <code>true</code> */
  bool    _fnav;
  /** EE Data is not valid */
  bool    _e1DataInValid;
  /** E5A Data is not valid */
  bool    _e5aDataInValid;
  /** E5B Data is not valid */
  bool    _e5bDataInValid;
};

class t_ephSBAS : public t_eph {
 friend class t_ephEncoder;
 friend class RTCM3Decoder;
 public:
  t_ephSBAS() {
    _IODN           = 0;
    _TOT            = 0.0;
    _agf0           = 0.0;
    _agf1           = 0.0;
    _x_pos          = 0.0;
    _x_velocity     = 0.0;
    _x_acceleration = 0.0;
    _y_pos          = 0.0;
    _y_velocity     = 0.0;
    _y_acceleration = 0.0;
    _z_pos          = 0.0;
    _z_velocity     = 0.0;
    _z_acceleration = 0.0;
    _ura            = 0.0;
    _health         = 0.0;
    _receptStaID    = "";    
  }
  t_ephSBAS(float rnxVersion, const QStringList& lines);
  virtual ~t_ephSBAS() {}

  virtual e_type  type() const {return t_eph::SBAS;}
  virtual unsigned int IOD() const;
  virtual unsigned int  isUnhealthy() const { return static_cast<unsigned int>(_health); }
  virtual QString toString(double version) const;

 private:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;

  int    _IODN;
  double _TOT;            // not used (set to  0.9999e9)
  double _agf0;           // [s]    clock correction
  double _agf1;           // [s/s]  clock correction drift

  double _x_pos;          // [m]
  double _x_velocity;     // [m/s]
  double _x_acceleration; // [m/s^2]

  double _y_pos;          // [m]
  double _y_velocity;     // [m/s]
  double _y_acceleration; // [m/s^2]

  double _z_pos;          // [m]
  double _z_velocity;     // [m/s]
  double _z_acceleration; // [m/s^2]

  mutable double _ura;
  double _health;
};

class t_ephBDS : public t_eph {
 friend class t_ephEncoder;
 friend class RTCM3Decoder;
 public:
 t_ephBDS() {
   _TOT             = 0.0;
   _AODE            = 0;
   _AODC            = 0;
   _URAI            = 0;
   _URA             = 0.0;
   _clock_bias      = 0.0;
   _clock_drift     = 0.0;
   _clock_driftrate = 0.0;
   _Crs             = 0.0;
   _Delta_n         = 0.0;
   _M0              = 0.0;
   _Cuc             = 0.0;
   _e               = 0.0;
   _Cus             = 0.0;
   _sqrt_A          = 0.0;
   _Cic             = 0.0;
   _OMEGA0          = 0.0;
   _Cis             = 0.0;
   _i0              = 0.0;
   _Crc             = 0.0;
   _omega           = 0.0;
   _OMEGADOT        = 0.0;
   _IDOT            = 0.0;
   _TGD1            = 0.0;
   _TGD2            = 0.0;
   _SatH1           = 0.0;
   _TOW             = 0.0;
   _TOEsec          = 0.0;
   _TOEweek         =-1.0;
    _receptStaID    = "";   
 }
 t_ephBDS(float rnxVersion, const QStringList& lines);
  virtual ~t_ephBDS() {}

  virtual e_type  type() const {return t_eph::BDS;}
  virtual unsigned int IOD() const;
  virtual unsigned int  isUnhealthy() const { return static_cast<unsigned int>(_SatH1); }
  virtual QString toString(double version) const;

 private:
  virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;

  double  _TOT;
  bncTime _TOE;
  int     _AODE;
  int     _AODC;
  int     _URAI;             //  [0..15] index from RTCM stream
  mutable double  _URA;      //  user range accuracy
  double  _clock_bias;       //  [s]
  double  _clock_drift;      //  [s/s]
  double  _clock_driftrate;  //  [s/s^2]
  double  _Crs;              //  [m]
  double  _Delta_n;          //  [rad/s]
  double  _M0;               //  [rad]
  double  _Cuc;              //  [rad]
  double  _e;                //
  double  _Cus;              //  [rad]
  double  _sqrt_A;           //  [m^0.5]
  double  _Cic;              //  [rad]
  double  _OMEGA0;           //  [rad]
  double  _Cis;              //  [rad]
  double  _i0;               //  [rad]
  double  _Crc;              //  [m]
  double  _omega;            //  [rad]
  double  _OMEGADOT;         //  [rad/s]
  double  _IDOT;             //  [rad/s]
  double  _TGD1;             //  [s]
  double  _TGD2;             //  [s]
  int     _SatH1;            //
  double  _TOW;              //  [s] of BDT week
  double  _TOEsec;           //  [s] of BDT week
  double  _TOEweek;          //  BDT week will be set only in case of RINEX file input
};

#endif