Homec4science
Diffusion Gear (master) Browse

Gear/src/system.c
b87f34a0ce7emaster

Last Change
Raw File

system.c
View Options

#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"
/*! \file system.c
* \brief contains miscellaneous routines, e.g. elapsed time measurements
*/
/*! This routine returns a random number taken from a table of random numbers,
* which is refilled every timestep. This method is used to allow random
* number application to particles independent of the number of processors
* used, and independent of the particular order the particles have. In order
* to work properly, the particle IDs should be set properly to unique
* integer values.
*/
double get_random_number(int id)
{
return RndTable[(id % RNDTABLE)];
}
#ifdef SFR
double get_StarFormation_random_number(int id)
{
return StarFormationRndTable[(id % RNDTABLE)];
}
#endif
#ifdef FEEDBACK_WIND
double get_FeedbackWind_random_number(int id)
{
return FeedbackWindRndTable[(id % RNDTABLE)];
}
#endif
#ifdef CHIMIE
double get_Chimie_random_number(int id)
{
return ChimieRndTable[(id % RNDTABLE)];
}
#endif
#ifdef CHIMIE_KINETIC_FEEDBACK
double get_ChimieKineticFeedback_random_number(int id)
{
return ChimieKineticFeedbackRndTable[(id % RNDTABLE)];
}
#endif
#ifdef GAS_ACCRETION
double get_gasAccretion_random_number(int id)
{
return gasAccretionRndTable[(id % RNDTABLE)];
}
#endif
#ifdef HOT_HALO
double get_Halo_random_number(int id)
{
return HaloRndTable[(id % RNDTABLE)];
}
#endif
/*! This routine fills the random number table.
*/
void set_random_numbers(void)
{
int i;
for(i = 0; i < RNDTABLE; i++)
RndTable[i] = gsl_rng_uniform(random_generator);
#ifdef SFR
for(i = 0; i < RNDTABLE; i++)
StarFormationRndTable[i] = gsl_rng_uniform(random_generator);
#endif
#ifdef FEEDBACK_WIND
for(i = 0; i < RNDTABLE; i++)
FeedbackWindRndTable[i] = gsl_rng_uniform(random_generator);
#endif
#ifdef CHIMIE
for(i = 0; i < RNDTABLE; i++)
ChimieRndTable[i] = gsl_rng_uniform(random_generator);
#endif
#ifdef CHIMIE_KINETIC_FEEDBACK
for(i = 0; i < RNDTABLE; i++)
ChimieKineticFeedbackRndTable[i] = gsl_rng_uniform(random_generator);
#endif
#ifdef GAS_ACCRETION
for(i = 0; i < RNDTABLE; i++)
gasAccretionRndTable[i] = gsl_rng_uniform(random_generator);
#endif
#ifdef HOT_HALO
for(i = 0; i < RNDTABLE; i++)
HaloRndTable[i] = gsl_rng_uniform(random_generator);
#endif
}
/*! returns the number of cpu-ticks in seconds that have elapsed, or the
* wall-clock time obtained with MPI_Wtime().
*/
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
*/
}
/*! returns the time difference between two measurements obtained with
* second(). The routine takes care of the possible overflow of the tick
* counter on 32bit systems, but depending on the system, this may not always
* work properly. Similarly, in some MPI implementations, the MPI_Wtime()
* function may also overflow, in which case a negative time difference would
* be returned. The routine returns instead a time difference equal to 0.
*/
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;
}
/*! returns the maximum of two double
*/
double dmax(double x, double y)
{
if(x > y)
return x;
else
return y;
}
/*! returns the minimum of two double
*/
double dmin(double x, double y)
{
if(x < y)
return x;
else
return y;
}
/*! returns the maximum of two integers
*/
int imax(int x, int y)
{
if(x > y)
return x;
else
return y;
}
/*! returns the minimum of two integers
*/
int imin(int x, int y)
{
if(x < y)
return x;
else
return y;
}
c4science · Help