GalINavISC.cpp

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.
//
//==============================================================================
#include <cmath>
#include "GalINavISC.hpp"
#include "FreqConv.hpp"
 
using namespace std;
 
namespace gnsstk
{
   GalINavISC ::
   GalINavISC()
         : bgdE1E5a(std::numeric_limits<double>::quiet_NaN()),
           bgdE1E5b(std::numeric_limits<double>::quiet_NaN())
   {
      weekFmt = "%4L(%4l)";
      msgLenSec = 2.0;
   }
 
 
   void GalINavISC ::
   dumpCorrections(std::ostream& s) const
   {
      const ios::fmtflags oldFlags = s.flags();
      s << "           CORRECTION"
        << endl << endl
        << scientific << setprecision(8) << setfill(' ')
        << setw(20) << left << "BGD(E1,E5a):" << setw(15) << bgdE1E5a << endl
        << setw(20) << left << "BGD(E1,E5b):" << setw(15) << bgdE1E5b << endl;
      s.flags(oldFlags);
   }
 
 
   bool GalINavISC ::
   validate() const
   {
      return !(std::isnan(bgdE1E5a) || std::isnan(bgdE1E5b));
   }
 
 
   bool GalINavISC ::
   getISC(const ObsID& oid, double& corrOut)
      const
   {
      if ((oid.band == CarrierBand::L1) && // same as E1
          ((oid.code == TrackingCode::E1B) ||
           (oid.code == TrackingCode::E1C)))
      {
         if (std::isnan(bgdE1E5b))
            return false;
         corrOut = -bgdE1E5b;
         return true;
      }
      else if ((oid.band == CarrierBand::L5) && // same as E5a
               ((oid.code == TrackingCode::E5aI) ||
                (oid.code == TrackingCode::E5aQ)))
      {
         if (std::isnan(bgdE1E5a))
            return false;
         corrOut = -(getGamma(CarrierBand::L1,oid.band) * bgdE1E5a);
         return true;
      }
      else if ((oid.band == CarrierBand::E5b) &&
               ((oid.code == TrackingCode::E5bI) ||
                (oid.code == TrackingCode::E5bQ)))
      {
         if (std::isnan(bgdE1E5b))
            return false;
         corrOut = -(getGamma(CarrierBand::L1,oid.band) * bgdE1E5b);
         return true;
      }
      return false;
   }
}