nonlinear.h

Go to the documentation of this file.

 
#include <vector>
 
#include "primordial.h"
#include "trigonometric_integrals.h"
 
#ifndef __NONLINEAR__
#define __NONLINEAR__
 
#define _M_EV_TOO_BIG_FOR_HALOFIT_ \
  10. 
#define _M_SUN_ 1.98847e30 
#define _MAX_NUM_EXTRAPOLATION_ 100000
 
enum non_linear_method { nl_none, nl_halofit, nl_HMcode };
enum pk_outputs { pk_linear, pk_nonlinear };
 
enum source_extrapolation {
  extrap_zero,
  extrap_only_max,
  extrap_only_max_units,
  extrap_max_scaled,
  extrap_hmcode,
  extrap_user_defined
};
 
enum halofit_integral_type { halofit_integral_one, halofit_integral_two, halofit_integral_three };
 
enum hmcode_baryonic_feedback_model {
  nl_emu_dmonly,
  nl_owls_dmonly,
  nl_owls_ref,
  nl_owls_agn,
  nl_owls_dblim,
  nl_user_defined
};
enum out_sigmas { out_sigma, out_sigma_prime, out_sigma_disp };
 
struct nonlinear {
 
  enum non_linear_method method =
      nl_none; 
  enum source_extrapolation extrapolation_method =
      extrap_max_scaled; 
  enum hmcode_baryonic_feedback_model feedback =
      nl_emu_dmonly; 
  double c_min; 
  double eta_0; 
  double z_infinity =
      10.; 
 
  bool has_pk_eq = false; 
 
  short nonlinear_verbose = 0; 
};
 
struct nonlinear_workspace {
 
  std::vector<double> rtab;   
  std::vector<double> stab;   
  std::vector<double> ddstab; 
  std::vector<double> growtable;
  std::vector<double> ztable;
  std::vector<double> tautable;
 
  std::vector<double*> sigma_8;        /* row pointers into sigma_8_storage */
  std::vector<double*> sigma_disp;     /* row pointers into sigma_disp_storage */
  std::vector<double*> sigma_disp_100; /* row pointers into sigma_disp_100_storage */
  std::vector<double*> sigma_prime;    /* row pointers into sigma_prime_storage */
 
  std::vector<std::vector<double>> sigma_8_storage;
  std::vector<std::vector<double>> sigma_disp_storage;
  std::vector<std::vector<double>> sigma_disp_100_storage;
  std::vector<std::vector<double>> sigma_prime_storage;
 
  double dark_energy_correction; 
};
 
#endif