From d595c2fc02fc0bdcf9b5b0b3dfe1d2b14d602212 Mon Sep 17 00:00:00 2001 From: David Seus <david.seus@ians.uni-stuttgart.de> Date: Fri, 16 Aug 2019 16:42:21 +0200 Subject: [PATCH] update TP-one-patch-purely-postive --- .../TP-one-patch-purely-postive-pc.py | 27 ++++++++++--------- 1 file changed, 14 insertions(+), 13 deletions(-) diff --git a/TP-one-patch/TP-one-patch-purely-postive-pc.py b/TP-one-patch/TP-one-patch-purely-postive-pc.py index 916f716..5638ef0 100755 --- a/TP-one-patch/TP-one-patch-purely-postive-pc.py +++ b/TP-one-patch/TP-one-patch-purely-postive-pc.py @@ -14,16 +14,16 @@ import helpers as hlp sym.init_printing() -solver_tol = 5E-6 +solver_tol = 1E-6 ############ GRID ####################### -use_case="TP-one-patch-purely-positive-pc" -mesh_resolution = 30 -timestep_size = 0.0005 -number_of_timesteps = 2500 + +mesh_resolution = 40 +timestep_size = 0.001 +number_of_timesteps = 10 # decide how many timesteps you want analysed. Analysed means, that we write out # subsequent errors of the L-iteration within the timestep. -number_of_timesteps_to_analyse = 5 +number_of_timesteps_to_analyse = 1 starttime = 0 Lw = 0.25 #/timestep_size @@ -32,18 +32,19 @@ Lnw=Lw l_param_w = 40 l_param_nw = 40 -include_gravity = True -debugflag = False +include_gravity = False +debugflag = True analyse_condition = True -output_string = "./output/nondirichlet_pure_positive_pc_number_of_timesteps{}_".format(number_of_timesteps) +use_case="TP-one-patch-purely-positive-pc" +output_string = "./output/pure_positive_pc_number_of_timesteps{}_".format(number_of_timesteps) ##### Domain and Interface #### # global simulation domain domain -sub_domain0_vertices = [df.Point(-1.0,-1.0), # - df.Point(1.0,-1.0),# - df.Point(1.0,1.0),# - df.Point(-1.0,1.0)] +sub_domain0_vertices = [df.Point(-1.0, -1.0), # + df.Point(1.0, -1.0), # + df.Point(1.0, 1.0), # + df.Point(-1.0, 1.0)] subdomain0_outer_boundary_verts = { 0: [sub_domain0_vertices[0], -- GitLab