BrcClockCorrection_T.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 "BrcClockCorrection.hpp"
#include "TestUtil.hpp"
#include "GPSWeekZcount.hpp"
 
using namespace std;
using namespace gnsstk;
 
#ifdef _MSC_VER
#define LDEXP(x,y) ldexp(x,y)
#else
#define LDEXP(x,y) std::ldexp(x,y)
#endif
 
class BrcClockCorrection_T
{
public:
   void fill(BrcClockCorrection& orbit);
   void fill2(BrcClockCorrection& orbit);
   unsigned initializationTest();
   unsigned equalityTest();
   unsigned svClockBiasTest();
};
 
 
void BrcClockCorrection_T ::
fill(BrcClockCorrection& orbit)
{
   ObsID oi(ObservationType::NavMsg, CarrierBand::L5, TrackingCode::Y);
   orbit.loadData("GPS", oi, 31, GPSWeekZcount(1886, 398400),
                  GPSWeekZcount(1887, 0), 5, 7, 9, true,
                  LDEXP(double( int32_t(0xfffff91d)), -31), // af0
                  LDEXP(double( int16_t(0xffed)),     -43), // af1
                  1.0e-15 // af2
                  );
}
 
 
void BrcClockCorrection_T ::
fill2(BrcClockCorrection& orbit)
{
      // By rules of Clock Correction, this must be week of Toc
   short weeknum   = 1638;
   short uraoc     = 5;
   short uraoc1    = 7;
   short uraoc2    = 7;
   double toc      = 388800.0;
   double top      = 302400.0;
   double af0      = 7.23189674E-04;
   double af1      = 5.11590770E-12;
   double af2      = 0.00000000E+00;
   CommonTime tocCT = GPSWeekSecond(weeknum, toc, TimeSystem::GPS);
   CommonTime topCT = GPSWeekSecond(weeknum, top, TimeSystem::GPS);
   ObsID oi(ObservationType::NavMsg, CarrierBand::L5, TrackingCode::Y);
   orbit.loadData("GPS", oi, 31, tocCT, topCT, uraoc, uraoc1, uraoc2, true,
                  af0, af1, af2);
}
 
 
unsigned BrcClockCorrection_T ::
initializationTest()
{
   BrcClockCorrection empty;
   CommonTime emptyTime;
   ObsID emptyObsID;
   TUDEF("BrcClockCorrection", "Default Constructor");
   TUASSERTE(bool, false, empty.dataLoaded);
   TUASSERTE(std::string, "", empty.satSys);
   TUASSERTE(ObsID, emptyObsID, empty.obsID);
   TUASSERTE(short, 0, empty.PRNID);
   TUASSERTE(CommonTime, emptyTime, empty.Toc);
   TUASSERTE(CommonTime, emptyTime, empty.Top);
   TUASSERTE(short, -16, empty.URAoc);
   TUASSERTE(short, 0, empty.URAoc1);
   TUASSERTE(short, 0, empty.URAoc2);
   TUASSERTE(bool, false, empty.healthy);
   TUASSERTE(double, 0, empty.af0);
   TUASSERTE(double, 0, empty.af1);
   TUASSERTE(double, 0, empty.af2);
   TURETURN();
}
 
 
unsigned BrcClockCorrection_T ::
equalityTest()
{
   TUDEF("BrcClockCorrection", "operator== / !=");
   BrcClockCorrection clock;
   fill(clock);
   BrcClockCorrection clockCopy(clock);
      // make sure our copy reports as being the same
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
      // Tweak each data member one by one and compare.
      // Yes the comments are fairly redundant but it helps to
      // visually separate each batch of statements per data
      // member.
   clockCopy.dataLoaded = false;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // satSys
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.satSys = "twaffle";
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // obsID
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.obsID = ObsID(ObservationType::NavMsg, CarrierBand::L1, TrackingCode::P);
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // PRNID
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.PRNID = 93;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // Toc
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.Toc = GPSWeekZcount(1234,56789);
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // Top
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.Top = GPSWeekZcount(1234,56789);
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // URAoc
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.URAoc = 943;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // URAoc1
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.URAoc1 = 943;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // URAoc2
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.URAoc2 = 943;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // healthy
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.healthy = false;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // af0
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.af0 = 1.5e-12;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // af1
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.af1 = 1.5e-12;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
      // af2
   TUCATCH(clockCopy = clock);
   TUASSERTE(BrcClockCorrection, clock, clockCopy);
   clockCopy.af2 = 1.5e-12;
   TUASSERT(clockCopy != clock);
   TUASSERT(!(clockCopy == clock));
   TURETURN();
}
 
 
unsigned BrcClockCorrection_T ::
svClockBiasTest()
{
   TUDEF("BrcClockCorrection", "svClockBias");
   BrcClockCorrection co1;
   CommonTime t = CivilTime(2011, 6, 2, 12, 14, 44.0, TimeSystem::GPS);
   fill2(co1);
   cout << setprecision(20)
        << co1.svClockDrift(t)
        << endl;
   TUASSERTFE(7.2319419646240680997e-4, co1.svClockBias(t));
   TUCSM("svClockBiasM");
   TUASSERTFE(216808.16576879983768, co1.svClockBiasM(t));
   TUCSM("svClockDrift");
   TUASSERTFE(5.1159076999999996399e-12, co1.svClockDrift(t));
   TURETURN();
}
 
 
int main() //Main function to initialize and run all tests above
{
   using namespace std;
   BrcClockCorrection_T testClass;
   unsigned errorTotal = 0;
 
   errorTotal += testClass.initializationTest();
   errorTotal += testClass.equalityTest();
   errorTotal += testClass.svClockBiasTest();
 
   cout << "Total Failures for " << __FILE__ << ": " << errorTotal << endl;
 
   return errorTotal; // Return the total number of errors
}