test/sct_test_accuracy.c

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2007 by Aleksander Morgado Juez                         *
 *   scitime@aleksander_morgado.mm.st                                      *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Library General Public License as       *
 *   published by the Free Software Foundation; either version 2 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 Library General Public     *
 *   License along with this program; if not, write to the                 *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

#include <stdio.h>
#include "scitime.h"

#define MAX_DATES   6
#define MAX_SECS    8

int main(int argc, char** argv)
{
    SCT_ACCURACY accuracy;
    double dates [MAX_DATES];
    double secs [MAX_SECS];
    sct_time    epoch1;
    sct_time    epoch2;
    sct_time    delta; 
    int ind;

    printf("\n\n---------------------------Test Getting accuracy from JD...\n");
    dates[0] = 1.0;
    dates[1] = 10.0;
    dates[2] = 100.0;
    dates[3] = 1.0E20;
    dates[4] = -1.0E20;
    dates[5] = 1.0E-20;
    for(ind = 0; ind < MAX_DATES; ind++)
    {
        accuracy = sct_get_accuracy_from_jd(dates[ind]);
        printf("For date '%.20lE', accuracy is '%d'\n",dates[ind],(int)accuracy);
    }

    printf("\n\n-------------------------Test Getting accuracy from secs...\n");
    secs[0] = 1.0;
    secs[1] = 1234.2;
    secs[2] = 1.1234567890;
    secs[3] = 1.123456789012345;
    secs[4] = 1.12345678901234567890;
    secs[5] = 1234567890.1234567890;
    secs[6] = 1E20;
    secs[7] = 1E-20;
    for(ind = 0; ind < MAX_SECS; ind++)
    {
        accuracy = sct_get_accuracy_from_secs(secs[ind]);
        printf("For secs '%.20lE', accuracy is '%d'\n",secs[ind],(int)accuracy);
    }

    printf("\n\n-------------------------Test scitime elements diff...\n");
    printf("Diff between a MJD2000 stored with automatic accuracy and a MJD2000"
        " with explicit accuracy set to SECS...\n");
    sct_set_mjd2000(&epoch1, 0.0);
    sct_set_mjd2000(&epoch2, 0.5);
    sct_set_accuracy(&epoch2, SCT_ACCURACY_SECS);
    sct_diff(&delta, &epoch1, &epoch2);

    printf("\n\n-------------------------Test scitime elements add...\n");
    printf("MJD2000 stored with automatic accuracy being added 15 femtosecs\n");
    sct_set_mjd2000(&epoch1, 0.0);
    sct_add_deltaDSF(&epoch1, &epoch1, 0, 0, 15);

    printf("\n\n");
    return 0;
}

---