summaryrefslogtreecommitdiff
path: root/src/fracture_square.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/fracture_square.cpp')
-rw-r--r--src/fracture_square.cpp87
1 files changed, 0 insertions, 87 deletions
diff --git a/src/fracture_square.cpp b/src/fracture_square.cpp
deleted file mode 100644
index f73ed85..0000000
--- a/src/fracture_square.cpp
+++ /dev/null
@@ -1,87 +0,0 @@
-
-#include <random>
-#include <iostream>
-
-#include <cholmod.h>
-
-#include "randutils/randutils.hpp"
-
-#include <graph.hpp>
-#include <network.hpp>
-#include <hooks.hpp>
-#include "measurements.hpp"
-
-#include <csignal>
-#include <cstring>
-#include <atomic>
-
-std::atomic<bool> quit(false); // signal flag
-
-void got_signal(int) {
- quit.store(true);
-}
-
-int main(int argc, char* argv[]) {
- struct sigaction sa;
- memset( &sa, 0, sizeof(sa) );
- sa.sa_handler = got_signal;
- sigfillset(&sa.sa_mask);
- sigaction(SIGINT, &sa, NULL);
-
- int opt;
-
- unsigned N = 1;
- unsigned Lx = 16;
- unsigned Ly = 16;
- double beta = 0.5;
-
- while ((opt = getopt(argc, argv, "N:X:Y:b:")) != -1) {
- switch (opt) {
- case 'N':
- N = (unsigned)atof(optarg);
- break;
- case 'X':
- Lx = atoi(optarg);
- break;
- case 'Y':
- Ly = atoi(optarg);
- break;
- case 'b':
- beta = atof(optarg);
- break;
- default:
- exit(1);
- }
- }
-
- cholmod_common c;
- CHOL_F(start)(&c);
-
- ma meas(Lx, Ly, beta, 0.5);
- graph G(Lx, Ly);
- elastic_network perm_network(G, &c);
-
- randutils::auto_seed_128 seeds;
- std::mt19937 rng{seeds};
-
- for (unsigned trial = 0; trial < N; trial++) {
- while (true) {
- try {
- elastic_network tmp_network(perm_network);
- tmp_network.set_thresholds(beta, rng);
- tmp_network.fracture(meas);
- break;
- } catch (std::exception &e) {
- std::cout << e.what() << '\n';
- }
- }
-
- if (quit.load())
- break;
- }
-
- CHOL_F(finish)(&c);
-
- return 0;
-}
-