source: Sophya/trunk/SophyaLib/BaseTools/timestamp.cc@ 2826

#include "sopnamsp.h"
#include "machdefs.h"
#include <stdio.h>
#include <time.h>
#include <ctype.h>
#include <math.h>
#include "timestamp.h"
#include "pexceptions.h"
#include <iostream>

/*!
   \class SOPHYA::TimeStamp
   \ingroup BaseTools
   A simple class for representing date and time. Simple operations 
   on date/time are also provided. 
   
   \code
   // Create a object with the current date and time and prints it to cout
   TimeStamp ts; 
   cout << ts << endl;
   // Create an object with a specified date and time 
   TimeStamp ts2("01/01/1905","00:00:00");
   // Get the number of days since 0 Jan 1901
   cout << ts2.ToDays() << endl;
   \endcode
*/

#ifdef Linux 
// La fonction trunc non declaree ds math.h sous Linux
extern "C" { double trunc   (double x); }
#endif
//-------------------------------------------------------------------------
//--------------------------  Classe TimeStamp   --------------------------
//-------------------------------------------------------------------------

enum Jour {jour_Lundi=0, jour_Mardi, jour_Mercredi, jour_Jeudi, jour_Vendredi, jour_Samedi, jour_Dimanche};
enum Mois {mois_Janvier=1, mois_Fevrier, mois_Mars, mois_Avril, mois_Mai, mois_Juin, mois_Juillet,
                   mois_Aout, mois_Septembre, mois_Octobre, mois_Novembre, mois_Decembre};

  /*
  //! Day of the week enum
  enum WeekDay {day_Monday=0, day_Tuesday, day_Wednesday, day_Thursday, day_Friday, day_Saturday, day_Sunday};
  */
 
TimeStamp::TimeStamp()
{
  mSeconds = 0.;
  mDays = 0;
  SetNow();
}

TimeStamp::TimeStamp(TimeStamp const & ts)
{
  Set(ts);
  mSeconds = ts.mSeconds;
  mDays = ts.mDays;
}

TimeStamp::TimeStamp(int year, int month, int day, int hour, int min, double sec)
{
  SetDate(year, month, day);
  SetHour(hour, min, sec);
}


TimeStamp::TimeStamp(double days)
{
  Set(days);
}

TimeStamp::TimeStamp(int_8 days, r_8 seconds)
{
  Set(days, seconds);
}

TimeStamp::TimeStamp(string& date, string& hour)
{
  SetDate(date);
  SetHour(hour);
}

TimeStamp::TimeStamp(const char* date, const char* hour)
{
  SetDate(date);
  SetHour(hour);
}


void TimeStamp::Set(TimeStamp const & ts)
{
  mSeconds = ts.mSeconds;
  mDays = ts.mDays;
}

void TimeStamp::Set(double days)
{
  if (days >= 0.) {
    mDays = trunc(days);
    mSeconds = (days-trunc(days))*86400.;
  }
  else {
    if ( (trunc(days)-days) > 0.) {
      mDays = trunc(days)-1;
      mSeconds = (days-mDays)*86400.;
    }
    else {
      mDays = trunc(days);
      mSeconds = 0.;
    }
  }
}

void TimeStamp::Set(int_8 days, r_8 seconds)
{
  if ( (seconds < 0.) || (seconds > 86400.) ) 
    throw ParmError("TimeStamp::Set(int_8, r_8) seconds<0 or seconds>86400.");
  mDays = days;
}

void TimeStamp::SetNow()
{
  time_t t = time(NULL);
  struct tm* TM = gmtime(&t);

  int JJ,MM,AA;
  int hh,mm;
  double ss; 
  
  AA = TM->tm_year + 1900;
  MM = TM->tm_mon+1;
  JJ = TM->tm_mday;
  hh = TM->tm_hour;
  mm = TM->tm_min;
  ss = TM->tm_sec;
  SetDate(AA,MM,JJ);
  SetHour(hh,mm,ss);
}

void TimeStamp::SetDate(int year, int month, int day)
{
  mDays = ConvertToDays(year, month, day);
}

void TimeStamp::SetDate(const char* date)
{
  int day,month,year;
  sscanf(date,"%d/%d/%d", &day, &month, &year);
  mDays = ConvertToDays(year, month, day);
}

void TimeStamp::SetHour(int hour, int min, double sec)
{
  mSeconds = hour*3600.+min*60+sec;
}

void TimeStamp::SetHour(const char* shour)
{
  int hour, min;
  double sec;
  sscanf(shour,"%d:%d:%lf",&hour, &min, &sec);
  mSeconds = hour*3600.+min*60+sec;

}

void TimeStamp::GetDate(int& year, int& month, int& day) const
{
  int_8 jours = mDays;
  // Recherche de l'annee
  if (jours < 0) {
    year = 1901;
    while(jours < 0)  { 
      year--;
      jours += YearDays(year); 
    }
  }
  else {
    year = 1901;
    while(jours > YearDays(year))  { 
      jours -= YearDays(year); 
      year++;
    }
  }
  // Recherche du mois
  month = 1;
  while(jours > MonthDays(year, month) )  { 
    jours -= MonthDays(year, month);
    month++;
  }
  day = jours;
}

void TimeStamp::GetHour(int& hour, int& min, double& sec) const
{
  double seconds = mSeconds;
  hour = trunc(seconds/3600.);
  seconds -= hour*3600;
  min = trunc(seconds/60.);
  while (min >= 60) { hour++; min -= 60; }
  sec = seconds-min*60;
  while (sec >= 60.) { min++; sec -= 60.; }
}

double TimeStamp::ToDays() const
{
  return((double)mDays + mSeconds/86400.);
}

/*!
\param fgday : if false, ignore the date (dd/mm/yy) part 
\param fghour : if false, ignore the hour (hh:mm:ss) part
*/
string TimeStamp::ToString(bool fgday, bool fghour) const
{
  char buff[128];
  int aa, mm, jj;
  int hh, min;
  double sec;
  GetDate(aa, mm, jj);
  GetHour(hh, min, sec);
  if (!fghour) 
    sprintf(buff,"%02d/%02d/%02d ", jj,mm,aa);
  else if (!fgday) 
    sprintf(buff,"%02d:%02d:%02.3f ", hh,min,sec);
  else 
    sprintf(buff,"%02d/%02d/%02d %02d:%02d:%02.1f GMT", jj,mm,aa,hh,min,sec);
  return buff;
}

/*!
\param fgday : if false, ignore the date (dd/mm/yy) part 
\param fghour : if false, ignore the hour (hh:mm:ss) part
*/
void TimeStamp::Print(ostream& os, bool fgday, bool fghour)  const
{
  os << " " << ToString(fgday, fghour) << " "; 
}

int TimeStamp::MonthDays(int annee, int mois) 
{
  if (mois<1 || mois>12) throw ParmError("TimeStamp::MonthDays month out of range");
  
  switch(mois) {
    case mois_Janvier: 
    case mois_Mars:
    case mois_Mai:
    case mois_Juillet:
    case mois_Aout:
    case mois_Octobre:
    case mois_Decembre:
      return 31;
    case mois_Avril:
    case mois_Juin:
    case mois_Septembre:
    case mois_Novembre:
      return 30;
    case mois_Fevrier:
     return (((annee%4 == 0) && (annee%100 != 0)) || (annee%400 == 0)) ? 29 : 28;
  }
  return -1;
}

int TimeStamp::YearDays(int annee)
//      Retourne le nombre de jours dans l'année
{
  return (((annee%4 == 0) && (annee%100 != 0)) || (annee%400 == 0)) ? 366 : 365;
}

int_8 TimeStamp::ConvertToDays(int AA, int MM, int JJ)
{
  int_8 t = 0;
  //   if (!UndetDate()) {
  int nban = AA-1901;
  if (nban >= 0)  
    t = nban*365 + (nban/4) - (nban/100) + ((nban+300)/400);
  else 
    t = nban*365 + (nban/4) - (nban/100) + ((nban-100)/400);
  for (int i=1; i<MM; i++)  t += TimeStamp::MonthDays(AA, i);
  t += JJ;
  // }
  return t;
}