transfer.h

Go to the documentation of this file.

 
#ifndef __TRANSFER__
#define __TRANSFER__
 
#include <string>
#include <vector>
 
#include "hyperspherical.h"
#include "nonlinear.h"
 
/* macro: test if index_tt is in the range between index and index+num, while the flag is true */
#define _index_tt_in_range_(index, num, flag) \
  (flag) && (index_tt >= index) && (index_tt < index + num)
/* macro: test if index_tt corresponds to an integrated nCl/sCl contribution */
#define _integrated_ncl_                                                                         \
  (_index_tt_in_range_(index_tt_lensing_, ppt->selection_num, ppt->has_cl_lensing_potential)) || \
      (_index_tt_in_range_(index_tt_nc_lens_, ppt->selection_num, ppt->has_nc_lens)) ||          \
      (_index_tt_in_range_(index_tt_nc_g4_, ppt->selection_num, ppt->has_nc_gr)) ||              \
      (_index_tt_in_range_(index_tt_nc_g5_, ppt->selection_num, ppt->has_nc_gr))
/* macro: test if index_tt corresponds to an non-integrated nCl/sCl contribution */
#define _nonintegrated_ncl_                                                            \
  (_index_tt_in_range_(index_tt_density_, ppt->selection_num, ppt->has_nc_density)) || \
      (_index_tt_in_range_(index_tt_rsd_, ppt->selection_num, ppt->has_nc_rsd)) ||     \
      (_index_tt_in_range_(index_tt_d0_, ppt->selection_num, ppt->has_nc_rsd)) ||      \
      (_index_tt_in_range_(index_tt_d1_, ppt->selection_num, ppt->has_nc_rsd)) ||      \
      (_index_tt_in_range_(index_tt_nc_g1_, ppt->selection_num, ppt->has_nc_gr)) ||    \
      (_index_tt_in_range_(index_tt_nc_g2_, ppt->selection_num, ppt->has_nc_gr)) ||    \
      (_index_tt_in_range_(index_tt_nc_g3_, ppt->selection_num, ppt->has_nc_gr))
/* macro: bin number associated to particular redshift bin and selection function for non-integrated contributions*/
#define _get_bin_nonintegrated_ncl_(index_tt)                                          \
  if (_index_tt_in_range_(index_tt_density_, ppt->selection_num, ppt->has_nc_density)) \
    bin = index_tt - index_tt_density_;                                                \
  if (_index_tt_in_range_(index_tt_rsd_, ppt->selection_num, ppt->has_nc_rsd))         \
    bin = index_tt - index_tt_rsd_;                                                    \
  if (_index_tt_in_range_(index_tt_d0_, ppt->selection_num, ppt->has_nc_rsd))          \
    bin = index_tt - index_tt_d0_;                                                     \
  if (_index_tt_in_range_(index_tt_d1_, ppt->selection_num, ppt->has_nc_rsd))          \
    bin = index_tt - index_tt_d1_;                                                     \
  if (_index_tt_in_range_(index_tt_nc_g1_, ppt->selection_num, ppt->has_nc_gr))        \
    bin = index_tt - index_tt_nc_g1_;                                                  \
  if (_index_tt_in_range_(index_tt_nc_g2_, ppt->selection_num, ppt->has_nc_gr))        \
    bin = index_tt - index_tt_nc_g2_;                                                  \
  if (_index_tt_in_range_(index_tt_nc_g3_, ppt->selection_num, ppt->has_nc_gr))        \
    bin = index_tt - index_tt_nc_g3_;
/* macro: bin number associated to particular redshift bin and selection function for integrated contributions*/
#define _get_bin_integrated_ncl_(index_tt)                                                       \
  if (_index_tt_in_range_(index_tt_lensing_, ppt->selection_num, ppt->has_cl_lensing_potential)) \
    bin = index_tt - index_tt_lensing_;                                                          \
  if (_index_tt_in_range_(index_tt_nc_lens_, ppt->selection_num, ppt->has_nc_lens))              \
    bin = index_tt - index_tt_nc_lens_;                                                          \
  if (_index_tt_in_range_(index_tt_nc_g4_, ppt->selection_num, ppt->has_nc_gr))                  \
    bin = index_tt - index_tt_nc_g4_;                                                            \
  if (_index_tt_in_range_(index_tt_nc_g5_, ppt->selection_num, ppt->has_nc_gr))                  \
    bin = index_tt - index_tt_nc_g5_;
struct transfers {
 
  double lcmb_rescale = 1.;  
  double lcmb_tilt    = 0.;  
  double lcmb_pivot   = 0.1; 
  double selection_bias[_SELECTION_NUM_MAX_] = {
      1.}; 
  double selection_magnification_bias[_SELECTION_NUM_MAX_] = {
      0.}; 
  bool has_nz_file = false; 
  bool has_nz_analytic =
      false;                
  std::string nz_file_name; 
  bool has_nz_evo_file = false; 
  bool has_nz_evo_analytic =
      false; 
  std::string nz_evo_file_name; 
 
 
  bool initialise_HIS_cache =
      false; 
  short transfer_verbose =
      0; 
};
 
struct transfer_workspace {
 
  HyperInterpStruct
      HIS; 
  HyperInterpStruct*
      pBIS; 
  int l_size; 
 
 
  int tau_size;     
  int tau_size_max; 
  std::vector<double> interpolated_sources; 
  std::vector<double> sources;              
  std::vector<double> tau0_minus_tau; 
  std::vector<double>
      w_trapz; 
  std::vector<double> chi;     
  std::vector<double> cscKgen; 
  std::vector<double> cotKgen; 
  std::vector<double> Phi;        
  std::vector<double> dPhi;       
  std::vector<double> d2Phi;      
  std::vector<double> chireverse; 
  std::vector<double> rescale_function; 
  std::vector<double>
      radial_function; 
  std::vector<double> chi_full_reverse; 
 
 
  double K; 
  int sgnK; 
 
  double
      tau0_minus_tau_cut; 
  short
      neglect_late_source; 
};
 
typedef enum {
  SCALAR_TEMPERATURE_0,
  SCALAR_TEMPERATURE_1,
  SCALAR_TEMPERATURE_2,
  SCALAR_POLARISATION_E,
  VECTOR_TEMPERATURE_1,
  VECTOR_TEMPERATURE_2,
  VECTOR_POLARISATION_E,
  VECTOR_POLARISATION_B,
  TENSOR_TEMPERATURE_2,
  TENSOR_POLARISATION_E,
  TENSOR_POLARISATION_B,
  NC_RSD
} radial_function_type;
 
enum Hermite_Interpolation_Order { HERMITE4, HERMITE6 };
 
#endif