CAMB.h

Go to the documentation of this file.

/*******************************************************************
 *  Copyright (C) 2022 by Federico Marulli and Giorgio Lesci       *
 *  federico.marulli3@unibo.it                                     *
 *                                                                 *
 *  This program is free software; you can redistribute it and/or  *
 *  modify it under the terms of the GNU General Public License as *
 *  published by the Free Software Foundation; either version 2 of *
 *  the License, or (at your option) any later version.            *
 *                                                                 *
 *  This program is distributed in the hope that it will be useful,*
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the  *
 *  GNU General Public License for more details.                   *
 *                                                                 *
 *  You should have received a copy of the GNU General Public      *
 *  License along with this program; if not, write to the Free     *
 *  Software Foundation, Inc.,                                     *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.      *
 *******************************************************************/
 
#ifndef __CAMB__
#define __CAMB__
 
#include "Kernel.h"
 
namespace cbl {
 
  namespace wrapper {
 
    namespace camb {
 
      extern"C"
      {
 
    void* GetCAMBparams ();
    
    void* GetCAMBdata ();
    
    void SetCAMBparams (void* params, const double ombh2, const double omch2, const double omnuh2, const double massless_nu, const int massive_nu, const int neutrino_hierarchy, const double omk, const double H0, const int dark_energy_model, const double w, const double wa, const double tau, const double cs2_lam, const double T_cmb, const double helium_fraction);
 
    void SetCAMBPk (void* params, const double redshift, const double ns, const double As, const double pivot_scalar, const bool accurate_massive_nu, const double kmax, const bool nonlinear);
 
    void GetCAMBresults (void* params, void* data);
 
    void GetCAMBPk (void* data, double *Pk, const double mink, const double dlnkh, const int npoints);
 
    void* ReleaseCAMBparams (void* params);
 
    void* ReleaseCAMBdata (void* data);
    
      }
      
      std::vector<double> Pk_CAMB (const bool nonlinear, const double redshift, const double kmin, const double kmax, const int npoints, const double ombh2, const double omch2, const double omnuh2, const double massless_nu, const int massive_nu, const double omk, const double H0, const double ns, const double As, const double pivot_scalar=0.05, const double w=-1., const double wa=0., const double tau=2.1e-9, const bool accurate_massive_nu=false, const int neutrino_hierarchy=3, const int dark_energy_model=1, const double cs2_lam=1., const double T_cmb=2.7255, const double helium_fraction=0.24);
      
    }
  }
}
 
#endif