vegas_read_pdf.cpp

#include "hep/mc.hpp"

#include <fstream>
#include <iostream>

int main(int argc, char* argv[])
{
    if (argc > 1)
    {
        // open file
        std::ifstream input(argv[1]);

        // we don't need to know the dimensions and bins beforehand!
        hep::vegas_pdf<double> grid(input);

        // print the (left bin) boundaries of the pdf in scientific format
        grid.serialize(std::cout);
        std::cout << '\n';
    }
    else
    {
        std::cout << "usage: " << argv[0] << " [p.d.f. file]\n\n";
        std::cout << "run `vegas_write_pdf > [p.d.f. file]` to generate a valid"
            " file\n";
    }

    return 0;
}