harmonic.h

Go to the documentation of this file.

 
#ifndef __HARMONIC__
#define __HARMONIC__
 
#include "transfer.h"
 
struct harmonic {
 
 
  int non_diag; 
 
 
  int md_size;           
  int index_md_scalars; 
  int * ic_size;         
  int * ic_ic_size;      
  short ** is_non_zero; 
 
 
  int has_tt; 
  int has_ee; 
  int has_te; 
  int has_bb; 
  int has_pp; 
  int has_tp; 
  int has_ep; 
  int has_dd; 
  int has_td; 
  int has_pd; 
  int has_ll; 
  int has_tl; 
  int has_dl; 
  int index_ct_tt; 
  int index_ct_ee; 
  int index_ct_te; 
  int index_ct_bb; 
  int index_ct_pp; 
  int index_ct_tp; 
  int index_ct_ep; 
  int index_ct_dd; 
  int index_ct_td; 
  int index_ct_pd; 
  int index_ct_ll; 
  int index_ct_tl; 
  int index_ct_dl; 
  int d_size;      
  int ct_size; 
 
 
  int * l_size;   
  int l_size_max; 
  double * l;    
  int ** l_max_ct;    
  int * l_max;    
  int l_max_tot; 
  double ** cl;   
  double ** ddcl; 
 
 
  struct fourier * pfo; 
  short harmonic_verbose; 
  ErrorMsg error_message; 
  short is_allocated; 
};
 
/*************************************************************************************************************/
/* @cond INCLUDE_WITH_DOXYGEN */
/*
 * Boilerplate for C++
 */
#ifdef __cplusplus
extern "C" {
#endif
 
  /* external functions (meant to be called from other modules) */
 
  int harmonic_cl_at_l(
                       struct harmonic * phr,
                       double l,
                       double * cl,
                       double ** cl_md,
                       double ** cl_md_ic
                       );
 
  /* internal functions */
 
  int harmonic_init(
                    struct precision * ppr,
                    struct background * pba,
                    struct perturbations * ppt,
                    struct primordial * ppm,
                    struct fourier *pfo,
                    struct transfer * ptr,
                    struct harmonic * phr
                    );
 
  int harmonic_free(
                    struct harmonic * phr
                    );
 
  int harmonic_indices(
                       struct background * pba,
                       struct perturbations * ppt,
                       struct transfer * ptr,
                       struct primordial * ppm,
                       struct harmonic * phr
                       );
 
  int harmonic_cls(
                   struct precision * ppr,
                   struct background * pba,
                   struct perturbations * ppt,
                   struct transfer * ptr,
                   struct primordial * ppm,
                   struct harmonic * phr
                   );
 
  int harmonic_compute_cl(
                          struct precision * ppr,
                          struct background * pba,
                          struct perturbations * ppt,
                          struct transfer * ptr,
                          struct primordial * ppm,
                          struct harmonic * phr,
                          int index_md,
                          int index_ic1,
                          int index_ic2,
                          int index_l,
                          int cl_integrand_num_columns,
                          double * cl_integrand,
                          double * cl_integrand_limber,
                          double * primordial_pk,
                          double * transfer_ic1,
                          double * transfer_ic2
                          );
 
  int harmonic_k_and_tau(
                         struct background * pba,
                         struct perturbations * ppt,
                         struct fourier *pfo,
                         struct harmonic * phr
                         );
 
  /* deprecated functions (since v2.8) */
 
  int harmonic_pk_at_z(
                       struct background * pba,
                       struct harmonic * phr,
                       enum linear_or_logarithmic mode,
                       double z,
                       double * output_tot,
                       double * output_ic,
                       double * output_cb_tot,
                       double * output_cb_ic
                       );
 
  int harmonic_pk_at_k_and_z(
                             struct background * pba,
                             struct primordial * ppm,
                             struct harmonic * phr,
                             double k,
                             double z,
                             double * pk,
                             double * pk_ic,
                             double * pk_cb,
                             double * pk_cb_ic
                             );
 
  int harmonic_pk_nl_at_z(
                          struct background * pba,
                          struct harmonic * phr,
                          enum linear_or_logarithmic mode,
                          double z,
                          double * output_tot,
                          double * output_cb_tot
                          );
 
  int harmonic_pk_nl_at_k_and_z(
                                struct background * pba,
                                struct primordial * ppm,
                                struct harmonic * phr,
                                double k,
                                double z,
                                double * pk_tot,
                                double * pk_cb_tot
                                );
 
  int harmonic_fast_pk_at_kvec_and_zvec(
                                        struct background * pba,
                                        struct harmonic * phr,
                                        double * kvec,
                                        int kvec_size,
                                        double * zvec,
                                        int zvec_size,
                                        double * pk_tot_out,
                                        double * pk_cb_tot_out,
                                        int nonlinear);
 
  int harmonic_sigma(
                     struct background * pba,
                     struct primordial * ppm,
                     struct harmonic * phr,
                     double R,
                     double z,
                     double *sigma
                     );
 
  int harmonic_sigma_cb(
                        struct background * pba,
                        struct primordial * ppm,
                        struct harmonic * phr,
                        double R,
                        double z,
                        double *sigma_cb
                        );
 
  /* deprecated functions (since v2.1) */
 
  int harmonic_tk_at_z(
                       struct background * pba,
                       struct harmonic * phr,
                       double z,
                       double * output
                       );
 
  int harmonic_tk_at_k_and_z(
                             struct background * pba,
                             struct harmonic * phr,
                             double k,
                             double z,
                             double * output
                             );
 
  /* end deprecated functions */
 
#ifdef __cplusplus
}
#endif
 
#endif
/* @endcond */