From c2ff5c7a8f818d4cf7c8af9df90bc2ab432687ce Mon Sep 17 00:00:00 2001 From: Lukas Riedel Date: Thu, 20 Apr 2017 17:12:04 +0200 Subject: [PATCH] ODE evaulation script now writes raw residual data into file --- python/testtools/dorie/testtools/dorie_run/ode.py | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/python/testtools/dorie/testtools/dorie_run/ode.py b/python/testtools/dorie/testtools/dorie_run/ode.py index bacd2e69..c5bc7db4 100644 --- a/python/testtools/dorie/testtools/dorie_run/ode.py +++ b/python/testtools/dorie/testtools/dorie_run/ode.py @@ -108,6 +108,11 @@ def evaluate(iniinfo,runtime): plt.savefig("{}/flux_residual.png".format(iniinfo["output.outputPath"])) plt.close() + # print data into file + filename = "{}/residuals.txt".format(iniinfo["output.outputPath"]) + f1,f2,f3 = res_flux.T + np.savetxt(filename,np.vstack((y,res_head,f1,f2,f3)).T,header="y,head,flux1,flux2,flux3") + head_tol = 1E-5 if not "_ode.head_abstol" in iniinfo else float(iniinfo["_ode.head_abstol"]) flux_tol = abs(influx) * 1e-5 if not "_ode.flux_abstol" in iniinfo else float(iniinfo["_ode.flux_abstol"]) print("tolerances: head L2 error: {:.2e}, flux L2 error: {:.2e}".format(head_tol,flux_tol)) -- GitLab