Main Page | Alphabetical List | Data Structures | File List | Data Fields | Globals | Related Pages

driftfac.c File Reference

compute loop-up tables for prefactors in cosmological integration More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <mpi.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_integration.h>
#include "allvars.h"
#include "proto.h"

Go to the source code of this file.

Defines

#define WORKSIZE   100000

Functions

void init_drift_table (void)
double get_drift_factor (int time0, int time1)
double get_gravkick_factor (int time0, int time1)
double get_hydrokick_factor (int time0, int time1)
double drift_integ (double a, void *param)
double gravkick_integ (double a, void *param)
double hydrokick_integ (double a, void *param)
double growthfactor_integ (double a, void *param)


Detailed Description

compute loop-up tables for prefactors in cosmological integration

Definition in file driftfac.c.


Define Documentation

#define WORKSIZE   100000
 

Referenced by init_drift_table().


Function Documentation

double drift_integ double  a,
void *  param
 

Integration kernel for drift factor computation.

Definition at line 179 of file driftfac.c.

References All, global_data_all_processes::Hubble, global_data_all_processes::Omega0, and global_data_all_processes::OmegaLambda.

double get_drift_factor int  time0,
int  time1
 

This function integrates the cosmological prefactor for a drift step between time0 and time1. The value returned is *

\[ \int_{a_0}^{a_1} \frac{{\rm d}a}{H(a)} * \]

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

Definition at line 67 of file driftfac.c.

References All, DRIFT_TABLE_LENGTH, DriftTable, and global_data_all_processes::Timebase_interval.

Referenced by move_particles().

double get_gravkick_factor int  time0,
int  time1
 

This function integrates the cosmological prefactor for a kick step of the gravitational force.

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

Definition at line 105 of file driftfac.c.

References All, DRIFT_TABLE_LENGTH, GravKickTable, and global_data_all_processes::Timebase_interval.

Referenced by advance_and_find_timesteps(), compute_global_quantities_of_system(), fill_write_buffer(), and move_particles().

double get_hydrokick_factor int  time0,
int  time1
 

This function integrates the cosmological prefactor for a kick step of the hydrodynamical force.

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

Definition at line 142 of file driftfac.c.

References All, DRIFT_TABLE_LENGTH, HydroKickTable, and global_data_all_processes::Timebase_interval.

Referenced by advance_and_find_timesteps(), compute_global_quantities_of_system(), fill_write_buffer(), and move_particles().

double gravkick_integ double  a,
void *  param
 

Integration kernel for gravitational kick factor computation.

Definition at line 191 of file driftfac.c.

References All, global_data_all_processes::Hubble, global_data_all_processes::Omega0, and global_data_all_processes::OmegaLambda.

double growthfactor_integ double  a,
void *  param
 

Definition at line 214 of file driftfac.c.

References All, global_data_all_processes::Omega0, global_data_all_processes::OmegaLambda, and pow().

double hydrokick_integ double  a,
void *  param
 

Integration kernel for hydrodynamical kick factor computation.

< adiabatic index of simulated gas

Definition at line 204 of file driftfac.c.

References All, GAMMA_MINUS1, global_data_all_processes::Hubble, global_data_all_processes::Omega0, global_data_all_processes::OmegaLambda, and pow().

void init_drift_table void   ) 
 

This function computes look-up tables for factors needed in cosmological integrations. The (simple) integrations are carried out with the GSL library. Separate factors are computed for the "drift", and the gravitational and hydrodynamical "kicks". The lookup-table is used for reasons of speed.

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

< length of the lookup table used to hold the drift and kick factors

Definition at line 26 of file driftfac.c.

References All, DRIFT_TABLE_LENGTH, DriftTable, GravKickTable, global_data_all_processes::Hubble, HydroKickTable, global_data_all_processes::TimeBegin, global_data_all_processes::TimeMax, and WORKSIZE.

Referenced by begrun().


Generated on Sun May 22 17:33:30 2005 for GADGET-2 by  doxygen 1.3.9.1