diff options
Diffstat (limited to 'collectStokesData.hpp')
| -rw-r--r-- | collectStokesData.hpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/collectStokesData.hpp b/collectStokesData.hpp index 42500c0..b48c6bc 100644 --- a/collectStokesData.hpp +++ b/collectStokesData.hpp @@ -5,8 +5,8 @@ using Complex = std::complex<Real>; template<int ...ps, class Generator, typename... T> -void collectStokesData(std::string tag, unsigned N, Generator& r, double ε, double δz, bool minimum, T... μs) { - unsigned nIts = 7; +void collectStokesData(std::string tag, unsigned N, Generator& r, Real ε, Real δz, bool minimum, T... μs) { + unsigned nIts = 5; unsigned nGs = 4; unsigned coDim = 16; Real newSaddleThres = 1e-4; @@ -83,11 +83,11 @@ void collectStokesData(std::string tag, unsigned N, Generator& r, double ε, dou std::vector<Vector<Complex>> oldGs; for (int j = 0; j < nIts; j++) { - Cord c(M, zMin, zSaddle, nGs + j * 2); + Cord c(M, zMin, zSaddle, nGs + pow(2, j)); for (unsigned i = 0; i < oldGs.size(); i++) { c.gs[i] = oldGs[i]; } - c.relaxNewton(c.gs.size() * 4, 1, 1e-10, 1e4); + c.relaxNewton(c.gs.size() * 2, 1, 1e-10, 1e4); oldGs = c.gs; Real reConstraintError = 0; |