system.c

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <time.h>
#include <sys/time.h>
#include <sys/resource.h>
#include <unistd.h>
#include <signal.h>
#include <gsl/gsl_rng.h>
#include <mpi.h>

#include "allvars.h"
#include "proto.h"


double get_random_number(int id)
{
  return RndTable[(id % RNDTABLE)];
}


void set_random_numbers(void)
{
  int i;

  for(i = 0; i < RNDTABLE; i++)
    RndTable[i] = gsl_rng_uniform(random_generator);
}


double second(void)
{
#ifdef WALLCLOCK
  return MPI_Wtime();
#else
  return ((double) clock()) / CLOCKS_PER_SEC;
#endif
  
  /* note: on AIX and presumably many other 32bit systems, 
   * clock() has only a resolution of 10ms=0.01sec 
   */
}


double timediff(double t0, double t1)
{
  double dt;

  dt = t1 - t0;

  if(dt < 0)    /* overflow has occured (for systems with 32bit tick counter) */
    {
#ifdef WALLCLOCK
      dt = 0;
#else
      dt = t1 + pow(2, 32) / CLOCKS_PER_SEC - t0;
#endif
    }

  return dt;
}


double dmax(double x, double y)
{
  if(x > y)
    return x;
  else
    return y;
}

double dmin(double x, double y)
{
  if(x < y)
    return x;
  else
    return y;
}

int imax(int x, int y)
{
  if(x > y)
    return x;
  else
    return y;
}

int imin(int x, int y)
{
  if(x < y)
    return x;
  else
    return y;
}

---