diff options
| -rw-r--r-- | src/fracture.cpp | 6 |
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}; |