file 1.17.1/1.17.1/polychord.hpp

[No description available] More…

Namespaces

Name
Gambit
TODO: see if we can use this one:
Gambit::PolyChord_1_17_1

Classes

Name
structSettings
classGambit::PolyChord_1_17_1::LogLikeWrapper
Class to connect PolyChord log-likelihood function and ScannerBit likelihood function.

Functions

Name
voidrun_polychord(double(*)(double *, int, double , int) loglikelihood, void()(int, int, int, double *, double *, double *, double, double) dumper, Settings )

Detailed Description

Author:

Date:

  • 2015 January
  • 2018 May
  • May 2020
  • June 2022

ScannerBit interface to PolyChord 1.17.1

Header file

Authors (add name and date if you modify):

Functions Documentation

function run_polychord 

void run_polychord(
    double(*)(double *, int, double *, int) loglikelihood,
    void(*)(int, int, int, double *, double *, double *, double, double) dumper,
    Settings 
)

Source code 

//  GAMBIT: Global and Modular BSM Inference Tool
//  *********************************************
///  \file
///
///  ScannerBit interface to PolyChord 1.17.1
///
///  Header file
///
///  *********************************************
///
///  Authors (add name and date if you modify):
///
///  \author Ben Farmer
///          (ben.farmer@gmail.com)
///  \date 2015 January
///
///  \author Will Handley
///          (wh260@cam.ac.uk)
///  \date 2018 May
///
///  \author Patrick Stoecker
///          (stoecker@physik.rwth-aachen.de)
///  \date May 2020
///
///  \author Anders Kvellestad
///          (anders.kvellestad@fys.uio.no)
///  \date June 2022
///
///  *********************************************

#ifndef __polychord_hpp__
#define __polychord_hpp__

#include "gambit/ScannerBit/scanner_plugin.hpp"

// Auxilliary classes and functions needed by polychord

struct Settings
{
   int nDims;
   int nDerived;
   int nlive;
   int num_repeats;
   int nprior;
   int nfail;
   bool do_clustering;
   int feedback;
   double precision_criterion;
   double logzero;
   int max_ndead;
   double boost_posterior;
   bool posteriors;
   bool equals;
   bool cluster_posteriors;
   bool write_resume;
   bool write_paramnames;
   bool read_resume;
   bool write_stats;
   bool write_live;
   bool write_dead;
   bool write_prior;
   bool maximise;
   double compression_factor;
   std::string base_dir;
   std::string file_root;
   std::vector<double> grade_frac;
   std::vector<int> grade_dims;
   std::vector<double> loglikes;
   std::vector<int> nlives;
   int seed;

   Settings(int _nDims=0,int _nDerived=0);
};

void run_polychord( 
         double (*loglikelihood)(double*,int,double*,int),
         void (*dumper)(int,int,int,double*,double*,double*,double,double), 
         Settings);


namespace Gambit
{

   namespace PolyChord_1_17_1
   {

      /// Typedef for the ScannerBit pointer to the external loglikelihood function
      typedef Gambit::Scanner::like_ptr scanPtr;

      /// Bring printer_interface and printer names into this scope
      using Gambit::Scanner::printer_interface;
      using Gambit::Scanner::printer;

      /// Class to connect PolyChord log-likelihood function and ScannerBit likelihood function
      class LogLikeWrapper
      {
         private:
            /// Scanner pointer (points to the ScannerBit provided log-likelihood function)
            scanPtr boundLogLike;

            /// Reference to a printer_interface object
            printer_interface& boundPrinter;

         public:
            /// Constructor
            LogLikeWrapper(scanPtr, printer_interface&);
   
            /// Main interface function from PolyChord to ScannerBit-supplied loglikelihood function 
            double LogLike(double*, int, double*, int);

            /// Main interface to PolyChord dumper routine   
            void dumper(int, int, int, double*, double*, double*, double, double);

            // PolyChord Hypercube index location for each parameter
            // gambit_hypercube[paramname] = polychord_hypercube[index_map[paramname]]
            // This is required because PolyChord needs parameters to be
            // ordered slow-to-fast to exploit a speed hierarchy
            std::unordered_map<std::string, int> index_map;

            /// copy of the settings in use.
            Settings boundSettings;

            /// Disable printing for speeds greater and equal than
            /// Speeds start at 0
            /// A value of -1 means that all evaluations are printed
            int printer_speed_threshold = -1;
      };

      ///@{ C-functions to pass to PolyChord for the callbacks
      double callback_loglike(double*, int, double*, int);
      void callback_dumper(int, int, int, double*, double*, double*, double, double);
      ///@}      

   } // End PolyChord namespace

} // End Gambit namespace

#endif

Updated on 2023-06-26 at 21:36:53 +0000