full run as add-on

Author: pancetta

docs/source/projects/doc_fput.rst

@@ -12,3 +12,5 @@ Results:
    :scale: 100 %
 .. image:: ../../../data/fput_energy.png
    :scale: 100 %
+.. image:: ../../../data/fput_energy_full.png
+   :scale: 100 %

pySDC/projects/Hamiltonian/README.rst

@@ -35,6 +35,7 @@ Fermi-Pasta-Ulam-Tsingou problem
 This is the famous FPUT problem, one of the first numerical simulation done in physics.
 The basis for this setup can be found `here <http://www.scholarpedia.org/article/Fermi-Pasta-Ulam_nonlinear_lattice_oscillations>`_ and we implement this in the ``FermiPastaUlamTsingou`` problem class.
 Due to time limitations in the CI environment, we only run the first few steps and not until Tend=10000.
+We refer to the left-most picture below for the full run.
 Note that we run MLSDC here, since PFASST does not seem to converge easily.
 
 .. include:: doc_fput.rst

pySDC/projects/Hamiltonian/data/fput_energy_full.png

@@ -106,10 +106,10 @@ def run_simulation():
     iter_counts = sort_stats(filtered_stats, sortby='time')
 
     # compute and print statistics
-    for item in iter_counts:
-        out = 'Number of iterations for time %4.2f: %2i' % item
-        f.write(out)
-        print(out)
+    # for item in iter_counts:
+    #     out = 'Number of iterations for time %4.2f: %2i' % item
+    #     f.write(out + '\n')
+    #     print(out)
 
     niters = np.array([item[1] for item in iter_counts])
     out = '   Mean number of iterations: %4.2f' % np.mean(niters)
@@ -140,7 +140,7 @@ def run_simulation():
     dill.dump(stats, f)
     f.close()
 
-    assert os.path.isfile(fname), 'Run for %s did not create stats file' % prob
+    assert os.path.isfile(fname), 'Run for %s did not create stats file'
 
 
 def show_results(cwd=''):