LinearClockModel.hpp

Go to the documentation of this file.

//==============================================================================
//
//  This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
//
//  The GNSSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 3.0 of the License, or
//  any later version.
//
//  The GNSSTk is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GNSSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin.
//  Copyright 2004-2022, The Board of Regents of The University of Texas System
//
//==============================================================================
 
//==============================================================================
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin, under contract to an agency or agencies
//  within the U.S. Department of Defense. The U.S. Government retains all
//  rights to use, duplicate, distribute, disclose, or release this software.
//
//  Pursuant to DoD Directive 523024
//
//  DISTRIBUTION STATEMENT A: This software has been approved for public
//                            release, distribution is unlimited.
//
//==============================================================================
 
#ifndef LINEARCLOCKMODEL_HPP
#define LINEARCLOCKMODEL_HPP
 
#include <map>
 
#include "Exception.hpp"
 
#include "ObsClockModel.hpp"
#include "ORDEpoch.hpp"
 
 
 
namespace gnsstk
{
 
 
   class LinearClockModel : public ObsClockModel
   {
   public:
      LinearClockModel(double sigma = 2, double elmask = 0, SvMode mode = ALWAYS)
            :ObsClockModel(sigma, elmask, mode) {reset();};
 
      virtual double getOffset(const gnsstk::CommonTime& t) const
         noexcept
      {
         if (!isOffsetValid(t))
            return 0;
         else
            return clockModel.Slope()*(t-baseTime) + clockModel.Intercept();
      };
 
      virtual bool isOffsetValid(const gnsstk::CommonTime& t) const noexcept
      {return t >= startTime && t <= endTime && clockModel.N() > 1;};
 
      virtual void addEpoch(const ORDEpoch& oe);
 
      void reset() noexcept;
 
      void dump(std::ostream& s, short detail=1) const noexcept;
 
      friend std::ostream& operator<<(std::ostream& s, const LinearClockModel& r)
      { r.dump(s, 0); return s; };
 
   private:
         // x is time y is clock offset
      gnsstk::TwoSampleStats<double> clockModel;
 
      gnsstk::CommonTime startTime, endTime, baseTime;
 
      unsigned long tossCount;
 
         // This is were we store what SVs were used to compute the individual
         // clock observations
      std::map<gnsstk::CommonTime, SvStatusMap> prnStatus;
 
         // This is a store of the clock observations that were added into the
         // clockModel object
      std::multimap<double,double> clockObs;
   };
 
 
}
#endif