test_CENamespace.cpp

Go to the documentation of this file.

/***************************************************************************
 *  test_CENamespace.cpp: CppEphem                                         *
 * ----------------------------------------------------------------------- *
 *  Copyright © 2019 JCardenzana                                           *
 * ----------------------------------------------------------------------- *
 *                                                                         *
 *  This program is free software: you can redistribute it and/or modify   *
 *  it under the terms of the GNU General Public License as published by   *
 *  the Free Software Foundation, either version 3 of the License, or      *
 *  (at your option) any later version.                                    *
 *                                                                         *
 *  This program 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 General Public License for more details.                           *
 *                                                                         *
 *  You should have received a copy of the GNU General Public License      *
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.  *
 *                                                                         *
 ***************************************************************************/
 
#include <iostream>
 
#include "test_CENamespace.h"
#include "CENamespace.h"
 
 
/**********************************************************************/
test_CENamespace::test_CENamespace() :
    CETestSuite()
{}
 
 
/**********************************************************************/
test_CENamespace::~test_CENamespace()
{}
 
 
/**********************************************************************/
bool test_CENamespace::runtests()
{
    std::cout << "\nTesting CENamespace:\n";
 
    // Run each of the tests
    test_SeaLevelVals();
    test_Conversions();
    test_Corrections();
    test_StrOpt();
 
    return pass();
}
 
 
/**********************************************************************/
bool test_CENamespace::test_SeaLevelVals(void)
{
    // Check the temperature at sea level
    double sea_level_temp_K = 288.2;
    test_double(CppEphem::SeaLevelTemp_K(), sea_level_temp_K, __func__, __LINE__);
 
    // Convert 
    double sea_level_temp_C = sea_level_temp_K - 273.15;
    test_double(CppEphem::SeaLevelTemp_C(), sea_level_temp_C, __func__, __LINE__);
 
    double sea_level_temp_F = (9.0/5.0) * sea_level_temp_C + 32.0;
    test_double(CppEphem::SeaLevelTemp_F(), sea_level_temp_F, __func__, __LINE__);
 
    return pass();
}
 
 
/**********************************************************************/
bool test_CENamespace::test_Conversions()
{
    // Test julian date 2000
    test_double(CppEphem::julian_date_J2000(), DJ00, __func__, __LINE__);
 
    // Speed of light (m/sec)
    test_double(CppEphem::c(), DC * DAU / DAYSEC, __func__, __LINE__);
 
    // Speed of light (AU/day)
    test_double(CppEphem::c_au_per_day(), DC, __func__, __LINE__);
 
    // Meters per AU
    test_double(CppEphem::m_per_au(), DAU, __func__, __LINE__);
 
    // Seconds per day
    test_double(CppEphem::sec_per_day(), DAYSEC, __func__, __LINE__);
 
    // Test temperature conversions
    double temp_K(0.0);
    double temp_C(-273.15);
    double temp_F(-459.67);
    // Reduce accuracy for the tests
    double old_tol = DblTol();
    SetDblTol(1.0e-11);
    test_double(CppEphem::Temp_C2F(temp_C), temp_F, __func__, __LINE__);
    test_double(CppEphem::Temp_C2K(temp_C), temp_K, __func__, __LINE__);
    test_double(CppEphem::Temp_F2C(temp_F), temp_C, __func__, __LINE__);
    test_double(CppEphem::Temp_F2K(temp_F), temp_K, __func__, __LINE__);
    test_double(CppEphem::Temp_K2C(temp_K), temp_C, __func__, __LINE__);
    test_double(CppEphem::Temp_K2F(temp_K), temp_F, __func__, __LINE__);
    SetDblTol(old_tol);
 
    return pass();
}
 
 
/**********************************************************************/
bool test_CENamespace::test_Corrections()
{
    double mjd(51544.5);
 
    // Get dut1 for J2000 (and test interpolation)
    CppEphem::CorrectionsInterp(false);
    test_double(CppEphem::dut1(mjd), 0.355499, __func__, __LINE__);
    CppEphem::CorrectionsInterp(true);
    test_double(CppEphem::dut1(mjd), 0.355066, __func__, __LINE__);
 
    // Turn off interpolation for future tests
    CppEphem::CorrectionsInterp(false);
 
    test_double(CppEphem::dut1Error(mjd), 0.0, __func__, __LINE__);
    test_double(CppEphem::dut1Calc(mjd), 0.0, __func__, __LINE__);
    
    // Test x,y polar motion
    test_double(CppEphem::xp(mjd), 0.043190 * DAS2R, __func__, __LINE__);
    test_double(CppEphem::yp(mjd), 0.377700 * DAS2R, __func__, __LINE__);
 
    // Test dpsi,deps corrections
    test_double(CppEphem::dpsi(mjd), -0.252*DMAS2R, __func__, __LINE__);
    test_double(CppEphem::deps(mjd), -0.119*DMAS2R, __func__, __LINE__);
 
    // Estimate altitude/pressure based on pressure
    test_double(CppEphem::EstimateAltitude_m(1013.25), 0.0, __func__, __LINE__);
    test_double(CppEphem::EstimatePressure_hPa(0.0), 1013.25, __func__, __LINE__);
 
    // Test the name of the corrections files
    std::string file_path(CECORRFILEPATH);
    std::string nut_file = file_path + "/nutation.txt";
    std::string ttut1_hist = file_path + "/ttut1_historic.txt";
    std::string ttut1_pred = file_path + "/ttut1_predicted.txt";
    
    test_string(CppEphem::NutationFile(), nut_file, __func__, __LINE__);
    test_string(CppEphem::TtUt1HistFile(), ttut1_hist, __func__, __LINE__);
    test_string(CppEphem::TtUt1PredFile(), ttut1_pred, __func__, __LINE__);
 
    // Test setting the corrections filename
    std::string newfilename("testfilename.txt");
    CppEphem::SetNutationFile(newfilename);
    CppEphem::SetTtUt1HistFile(newfilename);
    CppEphem::SetTtUt1PredFile(newfilename);
    test_string(CppEphem::NutationFile(), newfilename, __func__, __LINE__);
    test_string(CppEphem::TtUt1HistFile(), newfilename, __func__, __LINE__);
    test_string(CppEphem::TtUt1PredFile(), newfilename, __func__, __LINE__);
 
    // Reset the filenames
    CppEphem::SetNutationFile(nut_file);
    CppEphem::SetTtUt1HistFile(ttut1_hist);
    CppEphem::SetTtUt1PredFile(ttut1_pred);
 
    return pass();
}
 
 
/**********************************************************************/
bool test_CENamespace::test_StrOpt(void)
{
    // Create a string for parsing
    std::string test1 = "hello world!";
    test_vect(CppEphem::StrOpt::split(test1, ' '), {"hello", "world!"}, __func__, __LINE__);
 
    // See if we can pass the vector and get the same results back
    std::vector<std::string> test2;
    CppEphem::StrOpt::split(test1, ' ', &test2);
    test_vect(test2, {"hello", "world!"}, __func__, __LINE__);
 
    // Check if we can put it back together
    std::string test3 = CppEphem::StrOpt::join<std::string>(test2, ',');
    test_string(test3, "hello,world!", __func__, __LINE__);
 
    // See if we can join a vector of angle components into a legit string
    std::vector<double> angles = {123, 34, 52, 0.542};
    std::string angle_str = CppEphem::StrOpt::join_angle(angles, ':');
    test_string(angle_str, "123:34:52.54200000", __func__, __LINE__);
 
    return pass();
}
 
 
/**********************************************************************/
int main(int argc, char** argv) 
{    
    test_CENamespace tester;
    return (!tester.runtests());
}