#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"
extern "C" {
#include "rtcm3torinex.h"
}

namespace BNC_PPP {
  class t_orbCorr;
  class t_clkCorr;
}

class t_eph {
 public:

  enum e_type {unknown, GPS, GLONASS, Galileo};

  t_eph();
  virtual ~t_eph() {};

  static bool earlierTime(const t_eph* eph1, const t_eph* eph2) {
    return eph1->_TOC < eph2->_TOC;
  }

  virtual e_type type() const = 0;
  virtual QString toString(double version) const = 0;
  virtual void position(int GPSweek, double GPSweeks, 
                        double* xc, double* vv) const = 0;
  virtual int  IOD() const = 0;
  virtual int  RTCM3(unsigned char *) = 0;

  bool ok() const {return _ok;}
  bncTime TOC() const {return _TOC;}
  bool isNewerThan(const t_eph* eph) const {
    return earlierTime(eph, this);
  }
  t_prn prn() const {return _prn;}
  const QDateTime& receptDateTime() const {return _receptDateTime;}

  void position(int GPSweek, double GPSweeks, 
                double& xx, double& yy, double& zz, double& cc) const {
    double tmp_xx[4];
    double tmp_vv[4];
    position(GPSweek, GPSweeks, tmp_xx, tmp_vv);

    xx = tmp_xx[0];
    yy = tmp_xx[1];
    zz = tmp_xx[2];
    cc = tmp_xx[3];
  }

  t_irc getCrd(const bncTime& tt, ColumnVector& xc,
               ColumnVector& vv, bool useCorr) const;
  void setOrbCorr(const BNC_PPP::t_orbCorr* orbCorr);
  void setClkCorr(const BNC_PPP::t_clkCorr* clkCorr);
  virtual int slotNum() const {return 0;}

  static QString rinexDateStr(const bncTime& tt, const t_prn& prn,
                              double version);

  static QString rinexDateStr(const bncTime& tt, const QString& prnStr,
                              double version);

 protected:  
  t_prn           _prn;
  bncTime         _TOC;
  QDateTime       _receptDateTime;
  bool            _ok;
  BNC_PPP::t_orbCorr* _orbCorr;
  BNC_PPP::t_clkCorr* _clkCorr;
};


class t_ephGPS : public t_eph {
 public:
  t_ephGPS() { }
  t_ephGPS(float rnxVersion, const QStringList& lines);
  virtual ~t_ephGPS() {}

  virtual e_type type() const {return t_eph::GPS;}

  virtual QString toString(double version) const;

  void set(const gpsephemeris* ee);

  virtual void position(int GPSweek, double GPSweeks, 
                        double* xc,
                        double* vv) const;

  virtual int  IOD() const { return static_cast<int>(_IODC); }

  virtual int  RTCM3(unsigned char *);

  double TGD() const {return _TGD;} // Timing Group Delay (P1-P2 DCB)

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

  double  _IODE;            
  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 
  double  _TOEweek;
  double  _L2PFlag;         // L2 P data flag

  double  _ura;             // SV accuracy
  double  _health;          // SV health
  double  _TGD;             // [s]    
  double  _IODC;            

  double  _TOT;             // Transmisstion time 
  double  _fitInterval;     // Fit interval
};

class t_ephGlo : public t_eph {
 public:
  t_ephGlo() { _xv.ReSize(6); }
  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 void position(int GPSweek, double GPSweeks, 
                        double* xc,
                        double* vv) const;

  virtual int  IOD() const;

  virtual int  RTCM3(unsigned char *);

  void set(const glonassephemeris* ee);

  virtual int slotNum() const {return int(_frequency_number);}

 private:
  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;            //          
  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. 

  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]   
};

class t_ephGal : public t_eph {
 public:
  t_ephGal() { }
  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;}

  void set(const galileoephemeris* ee);

  virtual void position(int GPSweek, double GPSweeks, 
                        double* xc,
                        double* vv) const;

  virtual int  IOD() const { return static_cast<int>(_IODnav); }

  virtual int  RTCM3(unsigned char *);

 private:
  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

  double  _SISA;             //  Signal In Space Accuracy
  double  _E5aHS;            //  E5a Health Status
  double  _E5bHS;            //  E5a Health Status
  double  _BGD_1_5A;         //  group delay [s] 
  double  _BGD_1_5B;         //  group delay [s] 

  double _TOT;               // [s]

  int    _flags;
};

#endif