X-Git-Url: https://git.mxchange.org/?a=blobdiff_plain;f=simgear%2Ftiming%2Fsg_time.cxx;h=d32302ee1eef9c5575a5595f82ea612ea1ce734d;hb=b47d1ad5fd8ed111cae99c1f65f5bb65a5371501;hp=2db16af87d7e7587f8b2589f117ab04500953cad;hpb=e2b4c0787a4398199ef2d35e05963acd58b7a04c;p=simgear.git diff --git a/simgear/timing/sg_time.cxx b/simgear/timing/sg_time.cxx index 2db16af8..d32302ee 100644 --- a/simgear/timing/sg_time.cxx +++ b/simgear/timing/sg_time.cxx @@ -2,7 +2,7 @@ // // Written by Curtis Olson, started August 1997. // -// Copyright (C) 1997 Curtis L. Olson - curt@infoplane.com +// Copyright (C) 1997 Curtis L. Olson - http://www.flightgear.org/~curt // // This library is free software; you can redistribute it and/or // modify it under the terms of the GNU Library General Public @@ -14,39 +14,39 @@ // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU // Library General Public License for more details. // -// You should have received a copy of the GNU Library General Public -// License along with this library; if not, write to the -// Free Software Foundation, Inc., 59 Temple Place - Suite 330, -// Boston, MA 02111-1307, USA. +// You should have received a copy of the GNU General Public License +// along with this program; if not, write to the Free Software +// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. // // $Id$ +#ifdef HAVE_CONFIG_H +# include +#endif + #include #include // for errno -#ifdef SG_HAVE_STD_INCLUDES -# include -# include -# include -# include -#else -# include -# include -# include -# include -#endif +#include +#include +#include +#include +#include + +#ifdef HAVE_SYS_TIME_H +# include // for get/setitimer, gettimeofday, struct timeval +#endif #ifdef HAVE_SYS_TIMEB_H # include // for ftime() and struct timeb #endif #ifdef HAVE_UNISTD_H # include // for gettimeofday() #endif -#ifdef HAVE_SYS_TIME_H -# include // for get/setitimer, gettimeofday, struct timeval -#endif + +#include // for NAN #include #include @@ -56,69 +56,74 @@ #include "timezone.h" #include "lowleveltime.h" - #define DEGHR(x) ((x)/15.) #define RADHR(x) DEGHR(x*SGD_RADIANS_TO_DEGREES) +using std::string; static const double MJD0 = 2415020.0; static const double J2000 = 2451545.0 - MJD0; static const double SIDRATE = 0.9972695677; -SGTime::SGTime( double lon, double lat, const string& root ) +void SGTime::init( double lon_rad, double lat_rad, + const string& root, time_t init_time ) { SG_LOG( SG_EVENT, SG_INFO, "Initializing Time" ); gst_diff = -9999.0; - cur_time = time(NULL); - cout << "Current greenwich mean time = " << asctime(gmtime(&cur_time)) - << endl; - cout << "Current local time = " - << asctime(localtime(&cur_time)) << endl; + if ( init_time ) { + cur_time = init_time; + } else { + cur_time = time(NULL); + } + + SG_LOG( SG_EVENT, SG_INFO, + "Current greenwich mean time = " << asctime(gmtime(&cur_time))); + SG_LOG( SG_EVENT, SG_INFO, + "Current local time = " << asctime(localtime(&cur_time))); - if ( root != (string)"" ) { + if ( !root.empty()) { SGPath zone( root ); zone.append( "zone.tab" ); SG_LOG( SG_EVENT, SG_INFO, "Reading timezone info from: " << zone.str() ); - tzContainer = new TimezoneContainer( zone.c_str() ); - - GeoCoord location( SGD_RADIANS_TO_DEGREES * lat, SGD_RADIANS_TO_DEGREES * lon ); - GeoCoord* nearestTz = tzContainer->getNearest(location); + tzContainer = new SGTimeZoneContainer( zone.c_str() ); + SGGeod location(SGGeod::fromRad(lon_rad, lat_rad)); + SGTimeZone* nearestTz = tzContainer->getNearest(location); SGPath name( root ); name.append( nearestTz->getDescription() ); - zonename = strdup( name.c_str() ); + zonename = name.str(); SG_LOG( SG_EVENT, SG_INFO, "Using zonename = " << zonename ); } else { SG_LOG( SG_EVENT, SG_INFO, "*** NO TIME ZONE NAME ***" ); tzContainer = NULL; - zonename = NULL; + zonename.erase(); } } +SGTime::SGTime( double lon_rad, double lat_rad, const string& root, + time_t init_time ) +{ + init( lon_rad, lat_rad, root, init_time ); +} + SGTime::SGTime( const string& root ) { - SGTime( 0.0, 0.0, root ); + init( 0.0, 0.0, root, 0 ); } SGTime::SGTime() { - SGTime( 0.0, 0.0, "" ); + init( 0.0, 0.0, "", 0 ); } SGTime::~SGTime() { - if ( tzContainer != NULL ) { - delete tzContainer; - } - - if ( zonename != NULL ) { - delete zonename; - } + delete tzContainer; } @@ -145,7 +150,7 @@ static double sidereal_precise( double mjd, double lng ) // return a courser but cheaper estimate of sidereal time -static double sidereal_course( time_t cur_time, struct tm *gmt, double lng ) +static double sidereal_course( time_t cur_time, const struct tm *gmt, double lng ) { time_t start_gmt, now; double diff, part, days, hours, lstTmp; @@ -182,31 +187,39 @@ static double sidereal_course( time_t cur_time, struct tm *gmt, double lng ) // Update the time related variables -void SGTime::update( double lon, double lat, long int warp ) { +void SGTime::update( double lon_rad, double lat_rad, + time_t ct, long int warp ) +{ double gst_precise, gst_course; + + tm * gmt = &m_gmt; + + SG_LOG( SG_EVENT, SG_DEBUG, "Updating time" ); // get current Unix calendar time (in seconds) // warp += warp_delta; - cur_time = time(NULL) + warp; + if ( ct ) { + cur_time = ct + warp; + } else { + cur_time = time(NULL) + warp; + } SG_LOG( SG_EVENT, SG_DEBUG, " Current Unix calendar time = " << cur_time << " warp = " << warp ); // get GMT break down for current time - gmt = gmtime(&cur_time); + + memcpy( gmt, gmtime(&cur_time), sizeof(tm) ); SG_LOG( SG_EVENT, SG_DEBUG, " Current GMT = " << gmt->tm_mon+1 << "/" - << gmt->tm_mday << "/" << gmt->tm_year << " " + << gmt->tm_mday << "/" << (1900 + gmt->tm_year) << " " << gmt->tm_hour << ":" << gmt->tm_min << ":" << gmt->tm_sec ); - // calculate modified Julian date - // t->mjd = cal_mjd ((int)(t->gmt->tm_mon+1), (double)t->gmt->tm_mday, - // (int)(t->gmt->tm_year + 1900)); - mjd = sgTimeCalcMJD( (int)(gmt->tm_mon+1), (double)gmt->tm_mday, - (int)(gmt->tm_year + 1900) ); + // calculate modified Julian date starting with current + mjd = sgTimeCurrentMJD( ct, warp ); // add in partial day mjd += (gmt->tm_hour / 24.0) + (gmt->tm_min / (24.0 * 60.0)) + @@ -228,40 +241,67 @@ void SGTime::update( double lon, double lat, long int warp ) { gst_diff = gst_precise - gst_course; - lst = sidereal_course( cur_time, gmt, -(lon * SGD_RADIANS_TO_DEGREES) ) + gst_diff; + lst = sidereal_course( cur_time, gmt, + -(lon_rad * SGD_RADIANS_TO_DEGREES) ) + gst_diff; } else { // course + difference should drift off very slowly gst = sidereal_course( cur_time, gmt, 0.00 ) + gst_diff; - lst = sidereal_course( cur_time, gmt, -(lon * SGD_RADIANS_TO_DEGREES) ) + gst_diff; + lst = sidereal_course( cur_time, gmt, + -(lon_rad * SGD_RADIANS_TO_DEGREES) ) + gst_diff; } SG_LOG( SG_EVENT, SG_DEBUG, " Current lon=0.00 Sidereal Time = " << gst ); SG_LOG( SG_EVENT, SG_DEBUG, " Current LOCAL Sidereal Time = " << lst << " (" - << sidereal_precise( mjd, -(lon * SGD_RADIANS_TO_DEGREES) ) + << sidereal_precise( mjd, -(lon_rad * SGD_RADIANS_TO_DEGREES) ) << ") (diff = " << gst_diff << ")" ); } // Given lon/lat, update timezone information and local_offset -void SGTime::updateLocal( double lon, double lat, const string& root ) -{ - time_t currGMT; - time_t aircraftLocalTime; - GeoCoord location( SGD_RADIANS_TO_DEGREES * lat, SGD_RADIANS_TO_DEGREES * lon ); - GeoCoord* nearestTz = tzContainer->getNearest(location); - SGPath zone( root ); - zone.append ( nearestTz->getDescription() ); - if ( zonename ) { - delete zonename; - } - zonename = strdup( zone.c_str() ); - currGMT = sgTimeGetGMT( gmtime(&cur_time) ); - aircraftLocalTime = sgTimeGetGMT( (fgLocaltime(&cur_time, zone.c_str())) ); - local_offset = aircraftLocalTime - currGMT; - // cout << "Using " << local_offset << " as local time offset Timezone is " - // << zonename << endl; +void SGTime::updateLocal( double lon_rad, double lat_rad, const string& root ) { + // sanity checking + if ( lon_rad < -SGD_PI || lon_rad> SGD_PI ) { + // not within -180 ... 180 + lon_rad = 0.0; + } + if ( lat_rad < -SGD_PI_2 || lat_rad > SGD_PI_2 ) { + // not within -90 ... 90 + lat_rad = 0.0; + } + if ( lon_rad != lon_rad ) { + // only true if lon_rad == nan + SG_LOG( SG_EVENT, SG_ALERT, + " Detected lon_rad == nan, resetting to 0.0" ); + lon_rad = 0.0; + } + if ( lat_rad != lat_rad ) { + // only true if lat_rad == nan + SG_LOG( SG_EVENT, SG_ALERT, + " Detected lat_rad == nan, resetting to 0.0" ); + lat_rad = 0.0; + } + time_t currGMT; + time_t aircraftLocalTime; + SGGeod location(SGGeod::fromRad(lon_rad, lat_rad)); + SGTimeZone* nearestTz = tzContainer->getNearest(location); + SGPath zone( root ); + zone.append ( nearestTz->getDescription() ); + zonename = zone.str(); + + //Avoid troubles when zone.tab hasn't got the right line endings + if (zonename[zonename.size()-1] == '\r') + { + zonename[zonename.size()-1]=0; + zone.set( zonename ); + } + + currGMT = sgTimeGetGMT( gmtime(&cur_time) ); + aircraftLocalTime = sgTimeGetGMT( (fgLocaltime(&cur_time, zone.c_str())) ); + local_offset = aircraftLocalTime - currGMT; + // cout << "Using " << local_offset << " as local time offset Timezone is " + // << zonename << endl; } @@ -316,6 +356,43 @@ double sgTimeCalcMJD(int mn, double dy, int yr) { } +// return the current modified Julian date (number of days elapsed +// since 1900 jan 0.5), mjd. +double sgTimeCurrentMJD( time_t ct, long int warp ) { + + struct tm m_gmt; // copy of system gmtime(&time_t) structure + struct tm *gmt = &m_gmt; + + // get current Unix calendar time (in seconds) + // warp += warp_delta; + time_t cur_time; + if ( ct ) { + cur_time = ct + warp; + } else { + cur_time = time(NULL) + warp; + } + SG_LOG( SG_EVENT, SG_DEBUG, + " Current Unix calendar time = " << cur_time + << " warp = " << warp ); + + // get GMT break down for current time + memcpy( gmt, gmtime(&cur_time), sizeof(tm) ); + SG_LOG( SG_EVENT, SG_DEBUG, + " Current GMT = " << gmt->tm_mon+1 << "/" + << gmt->tm_mday << "/" << (1900 + gmt->tm_year) << " " + << gmt->tm_hour << ":" << gmt->tm_min << ":" + << gmt->tm_sec ); + + // calculate modified Julian date + // t->mjd = cal_mjd ((int)(t->gmt->tm_mon+1), (double)t->gmt->tm_mday, + // (int)(t->gmt->tm_year + 1900)); + double mjd = sgTimeCalcMJD( (int)(gmt->tm_mon+1), (double)gmt->tm_mday, + (int)(gmt->tm_year + 1900) ); + + return mjd; +} + + // given an mjd, calculate greenwich mean sidereal time, gst double sgTimeCalcGST( double mjd ) { double gst; @@ -431,7 +508,7 @@ time_t sgTimeGetGMT(int year, int month, int day, int hour, int min, int sec) #else // ! defined ( MK_TIME_IS_GMT ) // timezone seems to work as a proper offset for Linux & Solaris -# if defined( __linux__ ) || defined( __sun__ ) ||defined(__CYGWIN__) +# if defined( __linux__ ) || defined(__sun) ||defined(__CYGWIN__) # define TIMEZONE_OFFSET_WORKS 1 # endif