summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--src/fracture.cpp6
1 files changed, 5 insertions, 1 deletions
diff --git a/src/fracture.cpp b/src/fracture.cpp
index ec00b89..595cf1d 100644
--- a/src/fracture.cpp
+++ b/src/fracture.cpp
@@ -57,7 +57,11 @@ int main(int argc, char* argv[]) {
cholmod_common c;
CHOL_F(start)(&c);
- ma meas(Lx, Ly, 4*ceil(Lx), 4*ceil(Ly), beta, 4);
+ // fourier transforms of powers of two are faster
+ unsigned Mx = pow(2, ceil(log2(4*Lx)));
+ unsigned My = pow(2, ceil(log2(4*Ly)));
+
+ ma meas(Lx, Ly, Mx, My, beta, 4);
randutils::auto_seed_128 seeds;
std::mt19937 rng{seeds};