diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-02-24 17:23:29 -0500 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-02-24 17:23:29 -0500 |
commit | 67a11a04b1a16767e6ed30ed022ed9a0d3b07f76 (patch) | |
tree | ee4da8dc3a90f99c8857c8507299da3ad800e8d5 /src | |
parent | f861195c1416220c4039bda4d1cbf8c3aab07528 (diff) | |
download | fuse_networks-67a11a04b1a16767e6ed30ed022ed9a0d3b07f76.tar.gz fuse_networks-67a11a04b1a16767e6ed30ed022ed9a0d3b07f76.tar.bz2 fuse_networks-67a11a04b1a16767e6ed30ed022ed9a0d3b07f76.zip |
fixed problem computing Lx and Ly from a and n
Diffstat (limited to 'src')
-rw-r--r-- | src/fracture.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/fracture.cpp b/src/fracture.cpp index b1bde45..2651668 100644 --- a/src/fracture.cpp +++ b/src/fracture.cpp @@ -73,8 +73,8 @@ int main(int argc, char* argv[]) { std::mt19937 rng{seeds}; if (use_aN) { - unsigned Mx = pow(2, ceil(log2(4*2*n*a))); - unsigned My = pow(2, ceil(log2(4*2*n/a))); + unsigned Mx = pow(2, ceil(log2(4*sqrt(2*n*a)))); + unsigned My = pow(2, ceil(log2(4*sqrt(2*n/a)))); ma meas(n, a, Mx, My, beta); |