#ifndef MLMC_MC_HPP
#define MLMC_MC_HPP

#include "EmpiricMeasures.hpp"
#include "pdesolver/PDESolver.hpp"


class MonteCarlo {
protected:
    int verbose = 1;

    int plotting = 0;

    void method();

    void computeSampleSolution(int m, SampleID &id, SampleSolution &solution);

public:
    Level level;
    bool onlyFine;

    SampleCounter ctr;
    Sums sums;
    Averages avgs;
    Variances vars;
    Kurtosis kurtosis;

    PDESolver *pdeSolver;
    MatrixGraphs *solMGraphs;



    MonteCarlo(Level level, int dM, bool onlyFine,
               MatrixGraphs *solMGraphs, PDESolver *pdeSolver) :
        level(level), onlyFine(onlyFine),
        solMGraphs(solMGraphs), pdeSolver(pdeSolver) {
//        coarseSolution(SampleSolution(*solMGraphs, level.mGcoarse)),
//        fineSolution(SampleSolution(*solMGraphs, level.mGfine)) {

        config.get("MCPlotting", plotting);
        config.get("MCVerbose", verbose);

        ctr.dM = dM;
//        fineId.level = level;
//        fineId.coarse = false;
//        coarseId.level = level;
//        coarseId.coarse = true;
    }

    void Method();
};

#endif //MLMC_MC_HPP