lensing.h

Go to the documentation of this file.

 
#ifndef __LENSING__
#define __LENSING__
 
#include "harmonic.h"
 
struct lensing {
 
 
  short has_lensed_cls; 
 
 
  int has_tt; 
  int has_ee; 
  int has_te; 
  int has_bb; 
  int has_pp; 
  int has_tp; 
  int has_dd; 
  int has_td; 
  int has_ll; 
  int has_tl; 
  int index_lt_tt; 
  int index_lt_ee; 
  int index_lt_te; 
  int index_lt_bb; 
  int index_lt_pp; 
  int index_lt_tp; 
  int index_lt_dd; 
  int index_lt_td; 
  int index_lt_ll; 
  int index_lt_tl; 
  int lt_size; 
 
 
  int l_unlensed_max;    
  int l_lensed_max;    
  /* interpolable version: */
 
  int l_size;       
  int * l_max_lt;    
  double * l;       
  double * cl_lens; 
  double * ddcl_lens; 
 
 
  short lensing_verbose; 
  ErrorMsg error_message; 
  short is_allocated; 
};
 
/*************************************************************************************************************/
/* @cond INCLUDE_WITH_DOXYGEN */
/*
 * Boilerplate for C++
 */
#ifdef __cplusplus
extern "C" {
#endif
 
  int lensing_cl_at_l(
                      struct lensing * ple,
                      int l,
                      double * cl_lensed
                      );
 
  int lensing_init(
                   struct precision * ppr,
                   struct perturbations * ppt,
                   struct harmonic * phr,
                   struct fourier * pfo,
                   struct lensing * ple
                   );
 
  int lensing_free(
                   struct lensing * ple
                   );
 
  int lensing_indices(
                      struct precision * ppr,
                      struct harmonic * phr,
                      struct lensing * ple
                      );
 
  int lensing_lensed_cl_tt(
                           double *ksi,
                           double **d00,
                           double *w8,
                           int nmu,
                           struct lensing * ple
                           );
 
  int lensing_lensed_cl_te(
                           double *ksiX,
                           double **d20,
                           double *w8,
                           int nmu,
                           struct lensing * ple
                           );
 
  int lensing_lensed_cl_ee_bb(
                              double *ksip,
                              double *ksim,
                              double **d22,
                              double **d2m2,
                              double *w8,
                              int nmu,
                              struct lensing * ple
                              );
  int lensing_addback_cl_tt(
                            struct lensing *ple,
                            double *cl_tt
                            );
 
  int lensing_addback_cl_te(
                            struct lensing *ple,
                            double *cl_te
                            );
 
  int lensing_addback_cl_ee_bb(
                               struct lensing *ple,
                               double *cl_ee,
                               double *cl_bb
                               );
 
 
  int lensing_X000(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X000
                   );
 
  int lensing_Xp000(
                    double * mu,
                    int num_mu,
                    int lmax,
                    double * sigma2,
                    double ** Xp000
                    );
 
  int lensing_X220(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X220
                   );
 
  int lensing_X022(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X022
                   );
 
  int lensing_Xp022(
                    double * mu,
                    int num_mu,
                    int lmax,
                    double * sigma2,
                    double ** Xp022
                    );
 
  int lensing_X121(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X121
                   );
 
  int lensing_X132(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X132
                   );
 
  int lensing_X242(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double * sigma2,
                   double ** X242
                   );
 
  int lensing_d00(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d00
                  );
 
  int lensing_d11(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d11
                  );
 
  int lensing_d1m1(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d1m1
                   );
 
  int lensing_d2m2(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d2m2
                   );
 
  int lensing_d22(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d22
                  );
 
  int lensing_d20(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d20
                  );
 
  int lensing_d31(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d3m1
                  );
 
  int lensing_d3m1(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d3m1
                   );
 
  int lensing_d3m3(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d3m3
                   );
 
  int lensing_d40(
                  double * mu,
                  int num_mu,
                  int lmax,
                  double ** d40
                  );
 
  int lensing_d4m2(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d4m2
                   );
 
  int lensing_d4m4(
                   double * mu,
                   int num_mu,
                   int lmax,
                   double ** d4m4
                   );
 
#ifdef __cplusplus
}
#endif
 
#endif
/* @endcond */