blob: a4ed5d740a13e85102f4af943e67e139e83f1dd9 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
#include <cstring>
#include <fstream>
#include <string>
#include <cinttypes>
#include <sstream>
#include <functional>
#include <iostream>
#include <array>
#include <hooks.hpp>
#include "analysis_tools.hpp"
class ma : public hooks {
// need:
// - interface for turning on and off specific measurements
//
private:
unsigned Mx;
unsigned My;
unsigned num;
uint64_t Nc;
uint64_t NC;
uint64_t Nl;
uint64_t Nm;
uint64_t Ns;
uint64_t Nn;
uint64_t Nb;
uint64_t Na;
uint64_t Ne;
uint64_t NA;
uint64_t NE;
uint64_t Nq;
// measurement storage
std::vector<uint64_t> sc;
std::vector<uint64_t> sC;
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<uint64_t> sb; // final avalanche size distribution
std::vector<uint64_t> sd; // final avalanche size distribution
std::vector<uint64_t> sD; // final avalanche size distribution
std::vector<uint64_t> sf; // final avalanche size distribution
std::vector<uint64_t> sa; // non-final avalanche size distribution
std::vector<uint64_t> sA; // non-final avalanche size distribution
std::vector<uint64_t> se; // non-final avalanche size distribution
std::vector<uint64_t> sE; // non-final avalanche size distribution
std::vector<uint64_t> si;
std::vector<uint64_t> sI;
std::vector<uint64_t> sy;
std::vector<uint64_t> sY;
std::vector<uint64_t> sz;
std::vector<uint64_t> sZ;
std::vector<uint64_t> cc;
std::vector<uint64_t> cC;
std::vector<uint64_t> cn;
std::vector<uint64_t> cs;
std::vector<uint64_t> cm;
std::vector<uint64_t> cl;
std::vector<uint64_t> cb;
std::vector<uint64_t> ca;
std::vector<uint64_t> ce;
std::vector<uint64_t> cA;
std::vector<uint64_t> cE;
std::vector<uint64_t> cq;
ClusterTree last_clusters;
ClusterTree last_clusters_v;
public:
long double lc;
long double lv;
long double lm;
std::list<double> moduliA;
std::list<double> moduliE;
std::list<std::list<unsigned>> avalanches;
std::list<std::list<unsigned>> evalanches;
std::list<unsigned> last_avalanche;
std::list<unsigned> last_evalanche;
std::string model_string;
ma(unsigned Lx, unsigned Ly, double beta, double w);
ma(unsigned n, double a, double beta, double w);
~ma();
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;
};
|