#include <cstring>
#include <fstream>
#include <string>
#include <cinttypes>
#include <sstream>
#include <functional>
#include <iostream>
#include <array>

#include <fftw3.h>

#include <hooks.hpp>
#include <boost/graph/graph_utility.hpp>
#include <boost/graph/incremental_components.hpp>
#include <boost/pending/disjoint_sets.hpp>
#include <boost/foreach.hpp>

#include "analysis_tools.hpp"
typedef boost::graph_traits<Graph>::vertex_descriptor Vertex;
typedef boost::graph_traits<Graph>::vertices_size_type VertexIndex;


class pm : public hooks {
  private:
    unsigned p;
    Graph G;
    std::vector<VertexIndex> rank;
    std::vector<Vertex> parent;
    boost::disjoint_sets<VertexIndex*, Vertex*> ds;

    // measurement storage
    std::vector<uint64_t> sc; // non-spanning cluster size distribution
    std::vector<std::vector<uint64_t>> sn; // non-spanning cluster size distribution
    std::vector<uint64_t> ss; // minimal spanning cluster size distribution
    std::vector<uint64_t> sm; // spanning cluster size distribution
    std::vector<uint64_t> sl; // final avalanche size distribution
    std::vector<std::vector<uint64_t>> sb; // final avalanche size distribution
    std::vector<uint64_t> sp;

  public:
    std::string model_string;

    pm(unsigned Lx, unsigned Ly);
    pm(unsigned n, double a);
    ~pm();

    void pre_fracture(const network &) override;
    void bond_broken(const network& net, const current_info& cur, unsigned i) override;
    void post_fracture(network &n) override;
};