diff options
Diffstat (limited to 'langevin.cpp')
| -rw-r--r-- | langevin.cpp | 28 |
1 files changed, 12 insertions, 16 deletions
diff --git a/langevin.cpp b/langevin.cpp index e1464ec..e466860 100644 --- a/langevin.cpp +++ b/langevin.cpp @@ -182,30 +182,26 @@ int main(int argc, char* argv[]) { Rope<Complex> stokes(10, z1, z2); - std::cout << stokes.cost(Jp) << std::endl; + std::cout << stokes.cost(Jp) << " " << stokes.length() << " " << stokes.error(Jp) << std::endl; - stokes.relax(Jp, 10000, 1e-4); + stokes.relax(Jp, 1e5, 1e-1); - std::cout << stokes.cost(Jp) << std::endl; + std::cout << stokes.cost(Jp) << " " << stokes.length() << " " << stokes.error(Jp) << std::endl; - Rope<Complex> stokes2 = stokes.interpolate(); + stokes = stokes.interpolate(); + stokes.relax(Jp, 1e5, 1e-1); - stokes2.relax(Jp, 10000, 1e-4); + std::cout << stokes.cost(Jp) << " " << stokes.length() << " " << stokes.error(Jp) << std::endl; - std::cout << stokes2.cost(Jp) << std::endl; + stokes = stokes.interpolate(); + stokes.relax(Jp, 1e5, 1e-1); - stokes2 = stokes2.interpolate(); + std::cout << stokes.cost(Jp) << " " << stokes.length() << " " << stokes.error(Jp) << std::endl; - stokes2.relax(Jp, 10000, 1e-4); - - std::cout << stokes2.cost(Jp) << std::endl; - - stokes2 = stokes2.interpolate(); - - stokes2.relax(Jp, 10000, 1e-4); - - std::cout << stokes2.cost(Jp) << std::endl; + stokes = stokes.interpolate(); + stokes.relax(Jp, 1e5, 1e-1); + std::cout << stokes.cost(Jp) << " " << stokes.length() << " " << stokes.error(Jp) << std::endl; return 0; } |