From dc87673506bcd0abcfb731d4b90c5040c0aa730c Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@nians.uni-stuttgart.de>
Date: Fri, 6 Sep 2019 14:50:10 +0200
Subject: [PATCH] continue setting up the mesh study example
---
.../TP-TP-2-patch-pure-dd-convergence-study.py | 15 ++++++++-------
1 file changed, 8 insertions(+), 7 deletions(-)
diff --git a/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py b/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
index 441c93b..a4d2130 100755
--- a/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
+++ b/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
@@ -562,11 +562,12 @@ for mesh_resolution in resolutions:
simulation.initialise()
# simulation.write_exact_solution_to_xdmf()
- errornorms = simulation.run(analyse_condition=analyse_condition)
- for subdomain_index in errornorms.keys():
- for phase, different_errornorm in errornorms[subdomain_index].items():
- for errortype, errornorm in errornorms[subdomain_index][phase].items():
- eoc_filename = "{}_error".format(output_string)
- eocfile = open("eoc", "a")
- eocfile.write( str(dx) + " " + str(err) + "\n" )
+ output = simulation.run(analyse_condition=analyse_condition)
+ for subdomain_index in output.keys():
+ mesh_h = output['mesh_size']
+ for phase, different_errornorms in output[subdomain_index]['errornorm'].items():
+ for errortype, errornorm in different_errornorms.items():
+ eoc_filename = "{}_{}-errornorm-for-{}".format(output_string, errortype, phase)
+ eocfile = open("eoc_filename", "a")
+ eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
eocfile.close()
--
GitLab