diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-10-02 10:46:17 -0400 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-10-02 10:46:17 -0400 |
commit | 47a7c96d1848e3a3edbcfb17e645790d31c01a67 (patch) | |
tree | e45bf1cc6c57338cdd064acf0aa8e7a210d306fb /src/fracture_elastic.cpp | |
parent | 172d281943371a3aaa04a1a4d564731178e782e7 (diff) | |
download | fuse_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.cpp | 87 |
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; -} - |