summaryrefslogtreecommitdiff
path: root/src/fracture_elastic.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2019-10-02 10:46:17 -0400
committerJaron Kent-Dobias <jaron@kent-dobias.com>2019-10-02 10:46:17 -0400
commit47a7c96d1848e3a3edbcfb17e645790d31c01a67 (patch)
treee45bf1cc6c57338cdd064acf0aa8e7a210d306fb /src/fracture_elastic.cpp
parent172d281943371a3aaa04a1a4d564731178e782e7 (diff)
downloadfuse_networks-47a7c96d1848e3a3edbcfb17e645790d31c01a67.tar.gz
fuse_networks-47a7c96d1848e3a3edbcfb17e645790d31c01a67.tar.bz2
fuse_networks-47a7c96d1848e3a3edbcfb17e645790d31c01a67.zip
cleaned up libraries and reduced number of binaries
Diffstat (limited to 'src/fracture_elastic.cpp')
-rw-r--r--src/fracture_elastic.cpp87
1 files changed, 0 insertions, 87 deletions
diff --git a/src/fracture_elastic.cpp b/src/fracture_elastic.cpp
deleted file mode 100644
index dec4a2c..0000000
--- a/src/fracture_elastic.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, 2*Lx, 2*Ly, beta);
- graph G(Lx, Ly);
- network perm_network(G, &c);
-
- randutils::auto_seed_128 seeds;
- std::mt19937 rng{seeds};
-
- for (unsigned trial = 0; trial < N; trial++) {
- while (true) {
- try {
- 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;
-}
-