primordial_module.h

#ifndef PRIMORDIAL_MODULE_H
#define PRIMORDIAL_MODULE_H
 
#include "base_module.h"
#include "input_module.h"
#include "perturbations_module.h"
 
class PrimordialModule : public BaseModule {
 public:
  PrimordialModule(InputModulePtr input_module, PerturbationsModulePtr perturbation_module);
  ~PrimordialModule();
 
  void primordial_spectrum_at_k(int index_md,
                                enum linear_or_logarithmic mode,
                                double k,
                                double* pk) const;
  void primordial_output_titles(std::string& titles) const;
  void primordial_output_data(int number_of_titles, double* data) const;
 
  std::vector<int>
      ic_size_; 
  std::vector<int>
      ic_ic_size_; 
  std::vector<std::vector<short>>
      is_non_zero_; 
  int lnk_size_;    
  double phi_pivot_; 
  double phi_min_;   
  double phi_max_;   
  double phi_stop_;  
  double A_s_;
  double n_s_;
  double alpha_s_;
  double beta_s_;
  double r_;
  double n_t_;
  double alpha_t_;
 
 private:
  void primordial_init();
  void primordial_indices();
  void primordial_get_lnk_list(double kmin, double kmax, double k_per_decade);
  void primordial_analytic_spectrum_init();
  void primordial_analytic_spectrum(int index_md, int index_ic1_ic2, double k, double* pk) const;
  void primordial_inflation_potential(double phi, double* V, double* dV, double* ddV) const;
  void primordial_inflation_hubble(
      double phi, double* H, double* dH, double* ddH, double* dddH) const;
  void primordial_inflation_indices();
  void primordial_inflation_solve_inflation();
  void primordial_inflation_analytic_spectra(double* y_ini);
  void primordial_inflation_spectra(double* y_ini);
  void primordial_inflation_one_wavenumber(double* y_ini, int index_k);
  void primordial_inflation_one_k(
      double k, double* y, double* dy, double* curvature, double* tensor);
  void primordial_inflation_find_attractor(
      double phi_0, double precision, double* y, double* dy, double* H_0, double* dphidt_0);
  void primordial_inflation_evolve_background(double* y,
                                              double* dy,
                                              enum target_quantity target,
                                              double stop,
                                              bool check_epsilon,
                                              enum integration_direction direction,
                                              enum time_definition time);
  void primordial_inflation_check_potential(double phi, double* V, double* dV, double* ddV);
  void primordial_inflation_check_hubble(
      double phi, double* H, double* dH, double* ddH, double* dddH);
  void primordial_inflation_get_epsilon(double phi, double* epsilon);
  void primordial_inflation_find_phi_pivot(double* y, double* dy);
  void primordial_inflation_derivs_member(double tau,
                                          double* y,
                                          double* dy,
                                          void* parameters_and_workspace) const;
  static void primordial_inflation_derivs(double tau,
                                          double* y,
                                          double* dy,
                                          void* parameters_and_workspace);
  void primordial_external_spectrum_init();
 
  PerturbationsModulePtr perturbations_module_;
 
  int md_size_;             
  std::vector<double> lnk_; 
  std::vector<std::vector<double>>
      lnpk_; 
  std::vector<std::vector<double>>
      ddlnpk_; 
 
 
  std::vector<std::vector<double>>
      amplitude_; 
  std::vector<std::vector<double>>
      tilt_; 
  std::vector<std::vector<double>>
      running_; 
 
 
  int index_in_a_;       
  int index_in_phi_;     
  int index_in_dphi_;    
  int index_in_ksi_re_;  
  int index_in_ksi_im_;  
  int index_in_dksi_re_; 
  int index_in_dksi_im_; 
  int index_in_ah_re_;   
  int index_in_ah_im_;   
  int index_in_dah_re_;  
  int index_in_dah_im_;  
  int in_bg_size_;       
  int in_size_;          
};
 
#endif  //PRIMORDIAL_MODULE_H