cosmology.h

Go to the documentation of this file.

/*******************************************************************************
 * This file is part of SWIFT.
 * Copyright (c) 2017 Matthieu Schaller (schaller@strw.leidenuniv.nl)
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published
 * by the Free Software Foundation, either version 3 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 Lesser General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 ******************************************************************************/
#ifndef SWIFT_COSMOLOGY_H
#define SWIFT_COSMOLOGY_H
 
/* Config parameters. */
#include <config.h>
 
/* Local includes. */
//#include "parser.h"
//#include "physical_constants.h"
//#include "timeline.h"
//#include "units.h"
 
struct cosmology {
 
  double a;
 
  double a_inv;
 
  double a2_inv;
 
  double a3_inv;
 
  double a_factor_internal_energy;
 
  double a_factor_pressure;
 
  double a_factor_sound_speed;
 
  double a_factor_mu;
 
  double a_factor_Balsara_eps;
 
  double a_factor_grav_accel;
 
  double a_factor_hydro_accel;
 
  double z;
 
  double H;
 
  double critical_density;
 
  double critical_density_0;
 
  double mean_density;
 
  double mean_density_Omega_m;
 
  double mean_density_Omega_b;
 
  double overdensity_BN98;
 
  double time_step_factor;
 
  double a_dot;
 
  double time;
 
  double lookback_time;
 
  double w;
 
  double a_old;
 
  double z_old;
 
  /*------------------------------------------------------------------ */
 
  double a_begin;
 
  double a_end;
 
  double time_begin;
 
  double time_end;
 
  double time_base;
 
  double time_base_inv;
 
  double h;
 
  double H0;
 
  double Hubble_time;
 
  double Omega_cdm;
 
  double Omega_b;
 
  double Omega_lambda;
 
  double Omega_r;
 
  double Omega_g;
 
  double Omega_nu;
 
  double Omega_nu_0;
 
  double Omega_ur;
 
  double Omega_k;
 
  double w_0;
 
  double w_a;
 
  double T_CMB_0;
 
  double T_CMB_0_K;
 
  double T_nu_0;
 
  /* Neutrino temperature at z = 0 (electron-volts) */
  double T_nu_0_eV;
 
  int N_nu;
 
  double N_ur;
 
  double N_eff;
 
  double *M_nu_eV;
 
  double *deg_nu;
 
  double deg_nu_tot;
 
  double log_a_long_begin;
 
  double log_a_long_mid;
 
  double log_a_long_end;
 
  double log_a_begin;
 
  double log_a_end;
 
  double const_speed_light_c;
 
  double *drift_fac_interp_table;
 
  double *grav_kick_fac_interp_table;
 
  double *hydro_kick_fac_interp_table;
 
  double *hydro_kick_corr_interp_table;
 
  double *time_interp_table;
 
  double *scale_factor_interp_table;
 
  double *comoving_distance_interp_table;
 
  double comoving_distance_interp_table_offset;
 
  double comoving_distance_start_to_end;
 
  double *comoving_distance_inverse_interp_table;
 
  double *neutrino_density_early_table;
 
  double *neutrino_density_late_table;
 
  double time_interp_table_offset;
 
  double time_interp_table_max;
 
  double universe_age_at_present_day;
};
 
//void cosmology_update(struct cosmology *c, const struct phys_const *phys_const,
//                      integertime_t ti_current);
//
//double cosmology_get_drift_factor(const struct cosmology *cosmo,
//                                  const integertime_t ti_start,
//                                  const integertime_t ti_end);
//double cosmology_get_grav_kick_factor(const struct cosmology *cosmo,
//                                      const integertime_t ti_start,
//                                      const integertime_t ti_end);
//double cosmology_get_hydro_kick_factor(const struct cosmology *cosmo,
//                                       const integertime_t ti_start,
//                                       const integertime_t ti_end);
//double cosmology_get_therm_kick_factor(const struct cosmology *cosmo,
//                                       const integertime_t ti_start,
//                                       const integertime_t ti_end);
//double cosmology_get_corr_kick_factor(const struct cosmology *cosmo,
//                                      const integertime_t ti_start,
//                                      const integertime_t ti_end);
//double cosmology_get_delta_time(const struct cosmology *c,
//                                const integertime_t ti_start,
//                                const integertime_t ti_end);
//double cosmology_get_neutrino_density(const struct cosmology *c, double a);
//
//double cosmology_get_delta_time_from_scale_factors(const struct cosmology *c,
//                                                   const double a_start,
//                                                   const double a_end);
//
//double cosmology_get_timebase(struct cosmology *c,
//                              const integertime_t ti_current);
//
//double cosmology_get_scale_factor(const struct cosmology *cosmo, double t);
//
//double cosmology_get_comoving_distance(const struct cosmology *c,
//                                       const double a);
//
//double cosmology_scale_factor_at_comoving_distance(const struct cosmology *c,
//                                                   double r);
//
//double cosmology_get_time_since_big_bang(const struct cosmology *c, double a);
//void cosmology_init(struct swift_params *params, const struct unit_system *us,
//                    const struct phys_const *phys_const, struct cosmology *c);
//
//void cosmology_init_no_cosmo(struct cosmology *c);
//
//void cosmology_print(const struct cosmology *c);
//void cosmology_clean(struct cosmology *c);
//
//#ifdef HAVE_HDF5
//void cosmology_write_model(hid_t h_grp, const struct cosmology *c);
//#endif
//
//void cosmology_struct_dump(const struct cosmology *cosmology, FILE *stream);
//void cosmology_struct_restore(int enabled, struct cosmology *cosmology,
//                              FILE *stream);
 
#endif /* SWIFT_COSMOLOGY_H */