GPStoFloatTest.c

Go to the documentation of this file.

/*
*  Copyright (C) 2007 David Chin, Reinhard Prix
*
*  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 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 General Public License
*  along with with program; see the file COPYING. If not, write to the
*  Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
*  MA  02110-1301  USA
*/
 
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>
 
#include <lal/LALStdlib.h>
#include <lal/Date.h>
 
 
static int test_random_doubles(unsigned int seed)
{
        int i;
 
        srand(seed);
 
        for(i = 0; i < 100000000; i++) {
                double in;
                double out;
                LIGOTimeGPS gps;
 
                in = rand() * 2000000000.0 / RAND_MAX;
                in += (double) rand() / RAND_MAX;
 
                out = XLALGPSGetREAL8(XLALGPSSetREAL8(&gps, in));
 
                /* max allowed round-trip error is 1 ns */
                if(fabs(in - out) > 1e-9) {
                        fprintf(stderr, "XLALGPSSetREAL8() + XLALGPSGetREAL8() failed:  input = %.17g s, output = %.17g s, difference = %.16g ns ~ %.2g%% (seed was %u)\n", in, out, (in - out) * 1e9, fabs(in - out) / in * 100.0, seed);
                        return -1;
                }
        }
 
        return 0;
}
 
 
 
int main(void)
{
 
  /* 5 */
  if(test_random_doubles(time(NULL)))
      return 5;
 
  return 0;
}