.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples_mesh/example_meshconvergence_thacker2d.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_mesh_example_meshconvergence_thacker2d.py: ================================================================================ Thacker2d mesh convergence ================================================================================ In this example mesh convergence is carried out in the case of the thacker2d. H, QX and QY errors are computed by comparison to single layer analytical solution. .. GENERATED FROM PYTHON SOURCE LINES 10-19 .. code-block:: default import os import h5py import numpy as np import matplotlib.pyplot as plt import freshkiss3d as fk import freshkiss3d.extra.plots as fk_plt #sphinx_gallery_thumbnail_number = 7 .. GENERATED FROM PYTHON SOURCE LINES 20-22 Mesh parameters: -------------------- .. GENERATED FROM PYTHON SOURCE LINES 23-28 .. code-block:: default N_mesh = 5 #Number of mesh for mesh convergence First_mesh = 2 NL = [1 for r in range(N_mesh)] .. GENERATED FROM PYTHON SOURCE LINES 29-31 Case parameters: -------------------- .. GENERATED FROM PYTHON SOURCE LINES 32-41 .. code-block:: default PERIOD = 2.*np.pi/(np.sqrt(2.*9.81*0.1)) #une période: 4.485s FINAL_TIME = 1.*PERIOD MESH_SPLITTING = False PLOT_MESH = True PLOT_SOL = False PLOT_ERROR = True SAVE_ERROR = True .. GENERATED FROM PYTHON SOURCE LINES 42-44 Lists initialization: ----------------------- .. GENERATED FROM PYTHON SOURCE LINES 45-55 .. code-block:: default cases = [] simutime_list = [] triangular_mesh_list = [] NT = [] layer_list = [] primitives_list = [] problem_list = [] thacker2d_analytic_list = [] .. GENERATED FROM PYTHON SOURCE LINES 56-58 Meshes: -------------------- .. GENERATED FROM PYTHON SOURCE LINES 59-92 .. code-block:: default os.system('gmsh -2 -format msh2 inputs/thacker2d{}.geo -o inputs/thacker2d{}.msh'.format(First_mesh, First_mesh)) TG, vertex_labels, boundary_labels = fk.read_msh('inputs/thacker2d{}.msh'.format(First_mesh)) os.system('rm inputs/thacker2d{}.msh'.format(First_mesh)) x = np.asarray(TG.x) y = np.asarray(TG.y) trivtx = np.asarray(TG.trivtx) x *= 0.4 y *= 0.4 triangular_mesh_list.append( fk.TriangularMesh(TG, vertex_labels, boundary_labels) ) NT.append(triangular_mesh_list[0].NT) if PLOT_MESH: fk_plt.plot_mesh(triangular_mesh_list[0]) for M in range(1, N_mesh): print(M) if MESH_SPLITTING: triangular_mesh_list.append( triangular_mesh_list[M-1].refine_by_splitting() ) NT.append( triangular_mesh_list[M].NT) if PLOT_MESH: fk_plt.plot_mesh(triangular_mesh_list[M]) else: mesh_id = First_mesh+M os.system('gmsh -2 -format msh2 inputs/thacker2d{}.geo -o inputs/thacker2d{}.msh'.format(mesh_id, mesh_id)) TG, vertex_labels, boundary_labels = fk.read_msh('inputs/thacker2d{}.msh'.format(mesh_id)) os.system('rm inputs/thacker2d{}.msh'.format(mesh_id)) x = np.asarray(TG.x) y = np.asarray(TG.y) trivtx = np.asarray(TG.trivtx) x *= 0.4 y *= 0.4 triangular_mesh_list.append( fk.TriangularMesh(TG, vertex_labels, boundary_labels) ) NT.append(triangular_mesh_list[M].NT) if PLOT_MESH: fk_plt.plot_mesh(triangular_mesh_list[M]) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_001.png :alt: Mesh :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_001.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_002.png :alt: Mesh :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_002.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_003.png :alt: Mesh :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_003.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_004.png :alt: Mesh :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_004.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_005.png :alt: Mesh :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_005.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none 1 2 3 4 .. GENERATED FROM PYTHON SOURCE LINES 93-98 Cases: -------------------- Here we define cases based on a ``case`` class containing info such as space order or time order for each simulation we want to run. .. GENERATED FROM PYTHON SOURCE LINES 99-116 .. code-block:: default class case(): def __init__(self, NL=1, mesh_id=0, time_order=False, space_order=False): self.NL = NL self.mesh_id = mesh_id self.time_order = time_order self.space_order = space_order self.ipres = True if space_order is True: self.order = 2 else: self.order = 1 for M in range(N_mesh): cases.append(case(NL=NL[M], mesh_id=M, time_order=False, space_order=False)) cases.append(case(NL=NL[M], mesh_id=M, time_order=True, space_order=True)) .. GENERATED FROM PYTHON SOURCE LINES 117-119 Boundary conditions -------------------- .. GENERATED FROM PYTHON SOURCE LINES 120-123 .. code-block:: default fluvial_heights = [ fk.FluvialHeight(ref=r, height=0.0) for r in [1,2,3,4] ] .. GENERATED FROM PYTHON SOURCE LINES 124-126 Solving (loop over all cases): ------------------------------- .. GENERATED FROM PYTHON SOURCE LINES 127-187 .. code-block:: default for I, case in enumerate(cases): print(" ") print(" ") print(" SOLVING CASE {} ".format(I)) print(" ") print(" ") mesh_id = case.mesh_id time_second_order = case.time_order numerical_params = { 'space_second_order':case.space_order, 'ipres':case.ipres} simutime_list.append(fk.SimuTime( final_time=FINAL_TIME, time_iteration_max=10000, second_order=time_second_order)) thacker2d_analytic_list.append(fk.Thacker2D( triangular_mesh_list[mesh_id], a=1., h0=0.1, compute_error=True, error_type='L2', error_output='txt')) thacker2d_analytic_list[I](0.) layer_list.append(fk.Layer( case.NL, triangular_mesh_list[mesh_id], topography=thacker2d_analytic_list[I].topography)) primitives_list.append(fk.Primitive( triangular_mesh_list[mesh_id], layer_list[I], height=thacker2d_analytic_list[I].H, Uinit=thacker2d_analytic_list[I].U, Vinit=thacker2d_analytic_list[I].V)) if PLOT_SOL: create_figure = {'plot':fk_plt.plot_freesurface_3d_analytic_2} create_figure_scheduler = fk.schedules(times=[0.99*FINAL_TIME]) else: create_figure = None create_figure_scheduler = None problem_list.append(fk.Problem( simutime_list[I], triangular_mesh_list[mesh_id], layer_list[I], primitives_list[I], analytic_sol=thacker2d_analytic_list[I], numerical_parameters=numerical_params, fluvial_heights=fluvial_heights, custom_funct=create_figure, custom_funct_scheduler=create_figure_scheduler)) problem_list[I].solve() if PLOT_SOL: plt.show() .. rst-class:: sphx-glr-script-out .. code-block:: none SOLVING CASE 0 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=463, Layers=1, Triangles=848, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 4 ; Dt = 0.0248s ; Time = 0.10s ; ETA = 0.07s Iter = 8 ; Dt = 0.0248s ; Time = 0.20s ; ETA = 0.07s Iter = 12 ; Dt = 0.0250s ; Time = 0.30s ; ETA = 0.06s Iter = 15 ; Dt = 0.0250s ; Time = 0.37s ; ETA = 0.06s Iter = 19 ; Dt = 0.0252s ; Time = 0.47s ; ETA = 0.06s Iter = 22 ; Dt = 0.0254s ; Time = 0.55s ; ETA = 0.06s Iter = 26 ; Dt = 0.0255s ; Time = 0.65s ; ETA = 0.06s Iter = 30 ; Dt = 0.0257s ; Time = 0.75s ; ETA = 0.06s Iter = 33 ; Dt = 0.0259s ; Time = 0.83s ; ETA = 0.06s Iter = 37 ; Dt = 0.0262s ; Time = 0.94s ; ETA = 0.05s Iter = 40 ; Dt = 0.0263s ; Time = 1.01s ; ETA = 0.05s Iter = 44 ; Dt = 0.0265s ; Time = 1.12s ; ETA = 0.05s Iter = 47 ; Dt = 0.0267s ; Time = 1.20s ; ETA = 0.05s Iter = 51 ; Dt = 0.0269s ; Time = 1.31s ; ETA = 0.05s Iter = 54 ; Dt = 0.0271s ; Time = 1.39s ; ETA = 0.05s Iter = 57 ; Dt = 0.0273s ; Time = 1.47s ; ETA = 0.04s Iter = 61 ; Dt = 0.0275s ; Time = 1.58s ; ETA = 0.04s Iter = 64 ; Dt = 0.0277s ; Time = 1.66s ; ETA = 0.04s Iter = 67 ; Dt = 0.0279s ; Time = 1.75s ; ETA = 0.04s Iter = 71 ; Dt = 0.0281s ; Time = 1.86s ; ETA = 0.04s Iter = 74 ; Dt = 0.0282s ; Time = 1.94s ; ETA = 0.04s Iter = 77 ; Dt = 0.0283s ; Time = 2.03s ; ETA = 0.03s Iter = 80 ; Dt = 0.0285s ; Time = 2.11s ; ETA = 0.03s Iter = 83 ; Dt = 0.0286s ; Time = 2.20s ; ETA = 0.03s Iter = 87 ; Dt = 0.0287s ; Time = 2.31s ; ETA = 0.03s Iter = 90 ; Dt = 0.0287s ; Time = 2.40s ; ETA = 0.03s Iter = 93 ; Dt = 0.0288s ; Time = 2.49s ; ETA = 0.03s Iter = 96 ; Dt = 0.0288s ; Time = 2.57s ; ETA = 0.03s Iter = 99 ; Dt = 0.0288s ; Time = 2.66s ; ETA = 0.02s Iter = 103 ; Dt = 0.0289s ; Time = 2.77s ; ETA = 0.02s Iter = 106 ; Dt = 0.0289s ; Time = 2.86s ; ETA = 0.02s Iter = 109 ; Dt = 0.0290s ; Time = 2.95s ; ETA = 0.02s Iter = 112 ; Dt = 0.0290s ; Time = 3.03s ; ETA = 0.02s Iter = 115 ; Dt = 0.0290s ; Time = 3.12s ; ETA = 0.02s Iter = 118 ; Dt = 0.0291s ; Time = 3.21s ; ETA = 0.02s Iter = 121 ; Dt = 0.0292s ; Time = 3.30s ; ETA = 0.02s Iter = 125 ; Dt = 0.0294s ; Time = 3.41s ; ETA = 0.01s Iter = 128 ; Dt = 0.0295s ; Time = 3.50s ; ETA = 0.01s Iter = 131 ; Dt = 0.0296s ; Time = 3.59s ; ETA = 0.01s Iter = 134 ; Dt = 0.0298s ; Time = 3.68s ; ETA = 0.01s Iter = 137 ; Dt = 0.0299s ; Time = 3.77s ; ETA = 0.01s Iter = 140 ; Dt = 0.0301s ; Time = 3.86s ; ETA = 0.01s Iter = 143 ; Dt = 0.0302s ; Time = 3.95s ; ETA = 0.01s Iter = 146 ; Dt = 0.0304s ; Time = 4.04s ; ETA = 0.01s Iter = 149 ; Dt = 0.0306s ; Time = 4.13s ; ETA = 0.00s Iter = 152 ; Dt = 0.0308s ; Time = 4.22s ; ETA = 0.00s Iter = 155 ; Dt = 0.0309s ; Time = 4.32s ; ETA = 0.00s Iter = 158 ; Dt = 0.0311s ; Time = 4.41s ; ETA = 0.00s Iter = 161 ; Dt = 0.0130s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.06898927688598633s (wall time) SOLVING CASE 1 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=463, Layers=1, Triangles=848, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 4 ; Dt = 0.0248s ; Time = 0.10s ; ETA = 0.19s Iter = 8 ; Dt = 0.0248s ; Time = 0.20s ; ETA = 0.19s Iter = 12 ; Dt = 0.0248s ; Time = 0.30s ; ETA = 0.19s Iter = 15 ; Dt = 0.0249s ; Time = 0.37s ; ETA = 0.19s Iter = 19 ; Dt = 0.0248s ; Time = 0.47s ; ETA = 0.18s Iter = 23 ; Dt = 0.0249s ; Time = 0.57s ; ETA = 0.18s Iter = 26 ; Dt = 0.0250s ; Time = 0.65s ; ETA = 0.18s Iter = 30 ; Dt = 0.0250s ; Time = 0.75s ; ETA = 0.18s Iter = 34 ; Dt = 0.0250s ; Time = 0.85s ; ETA = 0.17s Iter = 37 ; Dt = 0.0250s ; Time = 0.92s ; ETA = 0.17s Iter = 41 ; Dt = 0.0251s ; Time = 1.02s ; ETA = 0.16s Iter = 45 ; Dt = 0.0251s ; Time = 1.12s ; ETA = 0.16s Iter = 48 ; Dt = 0.0233s ; Time = 1.19s ; ETA = 0.17s Iter = 53 ; Dt = 0.0202s ; Time = 1.30s ; ETA = 0.19s Iter = 57 ; Dt = 0.0204s ; Time = 1.38s ; ETA = 0.18s Iter = 62 ; Dt = 0.0212s ; Time = 1.49s ; ETA = 0.18s Iter = 66 ; Dt = 0.0197s ; Time = 1.57s ; ETA = 0.18s Iter = 71 ; Dt = 0.0181s ; Time = 1.66s ; ETA = 0.18s Iter = 76 ; Dt = 0.0177s ; Time = 1.75s ; ETA = 0.18s Iter = 82 ; Dt = 0.0151s ; Time = 1.85s ; ETA = 0.22s Iter = 88 ; Dt = 0.0155s ; Time = 1.94s ; ETA = 0.20s Iter = 93 ; Dt = 0.0161s ; Time = 2.02s ; ETA = 0.19s Iter = 99 ; Dt = 0.0173s ; Time = 2.12s ; ETA = 0.16s Iter = 104 ; Dt = 0.0183s ; Time = 2.21s ; ETA = 0.15s Iter = 109 ; Dt = 0.0183s ; Time = 2.30s ; ETA = 0.14s Iter = 114 ; Dt = 0.0199s ; Time = 2.40s ; ETA = 0.13s Iter = 118 ; Dt = 0.0181s ; Time = 2.47s ; ETA = 0.14s Iter = 124 ; Dt = 0.0157s ; Time = 2.57s ; ETA = 0.15s Iter = 130 ; Dt = 0.0161s ; Time = 2.67s ; ETA = 0.14s Iter = 135 ; Dt = 0.0178s ; Time = 2.75s ; ETA = 0.12s Iter = 140 ; Dt = 0.0197s ; Time = 2.84s ; ETA = 0.10s Iter = 144 ; Dt = 0.0236s ; Time = 2.93s ; ETA = 0.08s Iter = 148 ; Dt = 0.0223s ; Time = 3.03s ; ETA = 0.08s Iter = 153 ; Dt = 0.0206s ; Time = 3.13s ; ETA = 0.08s Iter = 157 ; Dt = 0.0201s ; Time = 3.21s ; ETA = 0.08s Iter = 162 ; Dt = 0.0198s ; Time = 3.31s ; ETA = 0.08s Iter = 166 ; Dt = 0.0191s ; Time = 3.39s ; ETA = 0.07s Iter = 171 ; Dt = 0.0185s ; Time = 3.48s ; ETA = 0.07s Iter = 176 ; Dt = 0.0175s ; Time = 3.57s ; ETA = 0.06s Iter = 182 ; Dt = 0.0188s ; Time = 3.68s ; ETA = 0.05s Iter = 186 ; Dt = 0.0200s ; Time = 3.76s ; ETA = 0.05s Iter = 190 ; Dt = 0.0228s ; Time = 3.85s ; ETA = 0.04s Iter = 195 ; Dt = 0.0204s ; Time = 3.95s ; ETA = 0.03s Iter = 199 ; Dt = 0.0198s ; Time = 4.03s ; ETA = 0.03s Iter = 204 ; Dt = 0.0189s ; Time = 4.13s ; ETA = 0.02s Iter = 209 ; Dt = 0.0180s ; Time = 4.22s ; ETA = 0.02s Iter = 214 ; Dt = 0.0168s ; Time = 4.31s ; ETA = 0.01s Iter = 220 ; Dt = 0.0158s ; Time = 4.41s ; ETA = 0.01s Iter = 226 ; Dt = 0.0054s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.278918981552124s (wall time) SOLVING CASE 2 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=1054, Layers=1, Triangles=1990, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 6 ; Dt = 0.0162s ; Time = 0.10s ; ETA = 0.15s Iter = 12 ; Dt = 0.0162s ; Time = 0.19s ; ETA = 0.16s Iter = 17 ; Dt = 0.0163s ; Time = 0.28s ; ETA = 0.14s Iter = 23 ; Dt = 0.0163s ; Time = 0.37s ; ETA = 0.14s Iter = 29 ; Dt = 0.0164s ; Time = 0.47s ; ETA = 0.13s Iter = 34 ; Dt = 0.0165s ; Time = 0.55s ; ETA = 0.13s Iter = 40 ; Dt = 0.0166s ; Time = 0.65s ; ETA = 0.13s Iter = 45 ; Dt = 0.0167s ; Time = 0.74s ; ETA = 0.13s Iter = 51 ; Dt = 0.0168s ; Time = 0.84s ; ETA = 0.12s Iter = 56 ; Dt = 0.0169s ; Time = 0.92s ; ETA = 0.12s Iter = 62 ; Dt = 0.0170s ; Time = 1.02s ; ETA = 0.12s Iter = 67 ; Dt = 0.0171s ; Time = 1.11s ; ETA = 0.11s Iter = 72 ; Dt = 0.0172s ; Time = 1.19s ; ETA = 0.11s Iter = 78 ; Dt = 0.0173s ; Time = 1.30s ; ETA = 0.10s Iter = 83 ; Dt = 0.0174s ; Time = 1.38s ; ETA = 0.10s Iter = 88 ; Dt = 0.0175s ; Time = 1.47s ; ETA = 0.10s Iter = 93 ; Dt = 0.0176s ; Time = 1.56s ; ETA = 0.09s Iter = 99 ; Dt = 0.0177s ; Time = 1.66s ; ETA = 0.10s Iter = 104 ; Dt = 0.0177s ; Time = 1.75s ; ETA = 0.09s Iter = 109 ; Dt = 0.0178s ; Time = 1.84s ; ETA = 0.09s Iter = 114 ; Dt = 0.0179s ; Time = 1.93s ; ETA = 0.08s Iter = 119 ; Dt = 0.0179s ; Time = 2.02s ; ETA = 0.08s Iter = 124 ; Dt = 0.0180s ; Time = 2.11s ; ETA = 0.08s Iter = 129 ; Dt = 0.0180s ; Time = 2.20s ; ETA = 0.08s Iter = 134 ; Dt = 0.0180s ; Time = 2.29s ; ETA = 0.08s Iter = 139 ; Dt = 0.0180s ; Time = 2.38s ; ETA = 0.07s Iter = 145 ; Dt = 0.0180s ; Time = 2.49s ; ETA = 0.06s Iter = 150 ; Dt = 0.0180s ; Time = 2.58s ; ETA = 0.06s Iter = 155 ; Dt = 0.0181s ; Time = 2.67s ; ETA = 0.06s Iter = 160 ; Dt = 0.0181s ; Time = 2.76s ; ETA = 0.06s Iter = 165 ; Dt = 0.0181s ; Time = 2.85s ; ETA = 0.05s Iter = 170 ; Dt = 0.0181s ; Time = 2.94s ; ETA = 0.05s Iter = 175 ; Dt = 0.0181s ; Time = 3.03s ; ETA = 0.05s Iter = 180 ; Dt = 0.0182s ; Time = 3.12s ; ETA = 0.05s Iter = 185 ; Dt = 0.0182s ; Time = 3.21s ; ETA = 0.04s Iter = 190 ; Dt = 0.0183s ; Time = 3.30s ; ETA = 0.04s Iter = 195 ; Dt = 0.0183s ; Time = 3.40s ; ETA = 0.03s Iter = 200 ; Dt = 0.0184s ; Time = 3.49s ; ETA = 0.03s Iter = 205 ; Dt = 0.0185s ; Time = 3.58s ; ETA = 0.03s Iter = 210 ; Dt = 0.0185s ; Time = 3.67s ; ETA = 0.03s Iter = 215 ; Dt = 0.0186s ; Time = 3.77s ; ETA = 0.03s Iter = 220 ; Dt = 0.0187s ; Time = 3.86s ; ETA = 0.02s Iter = 225 ; Dt = 0.0188s ; Time = 3.95s ; ETA = 0.02s Iter = 229 ; Dt = 0.0189s ; Time = 4.03s ; ETA = 0.02s Iter = 234 ; Dt = 0.0190s ; Time = 4.12s ; ETA = 0.01s Iter = 239 ; Dt = 0.0191s ; Time = 4.22s ; ETA = 0.01s Iter = 244 ; Dt = 0.0192s ; Time = 4.31s ; ETA = 0.01s Iter = 249 ; Dt = 0.0193s ; Time = 4.41s ; ETA = 0.00s Iter = 253 ; Dt = 0.0173s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.15381264686584473s (wall time) SOLVING CASE 3 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=1054, Layers=1, Triangles=1990, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 6 ; Dt = 0.0162s ; Time = 0.10s ; ETA = 0.54s Iter = 12 ; Dt = 0.0162s ; Time = 0.19s ; ETA = 0.54s Iter = 17 ; Dt = 0.0162s ; Time = 0.28s ; ETA = 0.53s Iter = 23 ; Dt = 0.0162s ; Time = 0.37s ; ETA = 0.52s Iter = 29 ; Dt = 0.0162s ; Time = 0.47s ; ETA = 0.52s Iter = 34 ; Dt = 0.0162s ; Time = 0.55s ; ETA = 0.51s Iter = 40 ; Dt = 0.0162s ; Time = 0.65s ; ETA = 0.50s Iter = 46 ; Dt = 0.0162s ; Time = 0.75s ; ETA = 0.49s Iter = 51 ; Dt = 0.0162s ; Time = 0.83s ; ETA = 0.48s Iter = 57 ; Dt = 0.0162s ; Time = 0.93s ; ETA = 0.47s Iter = 63 ; Dt = 0.0156s ; Time = 1.02s ; ETA = 0.48s Iter = 69 ; Dt = 0.0156s ; Time = 1.11s ; ETA = 0.46s Iter = 75 ; Dt = 0.0141s ; Time = 1.20s ; ETA = 0.50s Iter = 81 ; Dt = 0.0130s ; Time = 1.28s ; ETA = 0.52s Iter = 89 ; Dt = 0.0119s ; Time = 1.38s ; ETA = 0.55s Iter = 96 ; Dt = 0.0118s ; Time = 1.46s ; ETA = 0.54s Iter = 105 ; Dt = 0.0107s ; Time = 1.57s ; ETA = 0.58s Iter = 113 ; Dt = 0.0107s ; Time = 1.65s ; ETA = 0.56s Iter = 122 ; Dt = 0.0108s ; Time = 1.75s ; ETA = 0.54s Iter = 130 ; Dt = 0.0099s ; Time = 1.83s ; ETA = 0.58s Iter = 139 ; Dt = 0.0111s ; Time = 1.93s ; ETA = 0.50s Iter = 147 ; Dt = 0.0121s ; Time = 2.02s ; ETA = 0.45s Iter = 154 ; Dt = 0.0124s ; Time = 2.11s ; ETA = 0.43s Iter = 162 ; Dt = 0.0123s ; Time = 2.20s ; ETA = 0.40s Iter = 169 ; Dt = 0.0132s ; Time = 2.30s ; ETA = 0.39s Iter = 176 ; Dt = 0.0134s ; Time = 2.39s ; ETA = 0.35s Iter = 183 ; Dt = 0.0120s ; Time = 2.48s ; ETA = 0.38s Iter = 191 ; Dt = 0.0112s ; Time = 2.57s ; ETA = 0.39s Iter = 199 ; Dt = 0.0121s ; Time = 2.66s ; ETA = 0.34s Iter = 207 ; Dt = 0.0119s ; Time = 2.76s ; ETA = 0.34s Iter = 214 ; Dt = 0.0146s ; Time = 2.85s ; ETA = 0.25s Iter = 220 ; Dt = 0.0143s ; Time = 2.94s ; ETA = 0.25s Iter = 226 ; Dt = 0.0140s ; Time = 3.02s ; ETA = 0.24s Iter = 233 ; Dt = 0.0143s ; Time = 3.12s ; ETA = 0.23s Iter = 239 ; Dt = 0.0152s ; Time = 3.21s ; ETA = 0.20s Iter = 245 ; Dt = 0.0134s ; Time = 3.30s ; ETA = 0.20s Iter = 253 ; Dt = 0.0122s ; Time = 3.40s ; ETA = 0.21s Iter = 260 ; Dt = 0.0132s ; Time = 3.49s ; ETA = 0.17s Iter = 267 ; Dt = 0.0130s ; Time = 3.58s ; ETA = 0.16s Iter = 274 ; Dt = 0.0125s ; Time = 3.67s ; ETA = 0.15s Iter = 282 ; Dt = 0.0114s ; Time = 3.76s ; ETA = 0.15s Iter = 289 ; Dt = 0.0145s ; Time = 3.85s ; ETA = 0.10s Iter = 296 ; Dt = 0.0128s ; Time = 3.94s ; ETA = 0.10s Iter = 303 ; Dt = 0.0134s ; Time = 4.04s ; ETA = 0.08s Iter = 310 ; Dt = 0.0133s ; Time = 4.13s ; ETA = 0.06s Iter = 317 ; Dt = 0.0123s ; Time = 4.22s ; ETA = 0.05s Iter = 324 ; Dt = 0.0135s ; Time = 4.31s ; ETA = 0.03s Iter = 331 ; Dt = 0.0139s ; Time = 4.40s ; ETA = 0.01s Iter = 338 ; Dt = 0.0023s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.7573263645172119s (wall time) SOLVING CASE 4 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=1848, Layers=1, Triangles=3538, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 8 ; Dt = 0.0121s ; Time = 0.10s ; ETA = 0.27s Iter = 16 ; Dt = 0.0121s ; Time = 0.19s ; ETA = 0.27s Iter = 23 ; Dt = 0.0121s ; Time = 0.28s ; ETA = 0.26s Iter = 31 ; Dt = 0.0121s ; Time = 0.37s ; ETA = 0.27s Iter = 38 ; Dt = 0.0122s ; Time = 0.46s ; ETA = 0.25s Iter = 46 ; Dt = 0.0122s ; Time = 0.56s ; ETA = 0.25s Iter = 53 ; Dt = 0.0123s ; Time = 0.64s ; ETA = 0.25s Iter = 61 ; Dt = 0.0123s ; Time = 0.74s ; ETA = 0.24s Iter = 68 ; Dt = 0.0124s ; Time = 0.83s ; ETA = 0.23s Iter = 76 ; Dt = 0.0124s ; Time = 0.93s ; ETA = 0.24s Iter = 83 ; Dt = 0.0125s ; Time = 1.01s ; ETA = 0.22s Iter = 90 ; Dt = 0.0126s ; Time = 1.10s ; ETA = 0.21s Iter = 97 ; Dt = 0.0126s ; Time = 1.19s ; ETA = 0.21s Iter = 105 ; Dt = 0.0127s ; Time = 1.29s ; ETA = 0.21s Iter = 112 ; Dt = 0.0127s ; Time = 1.38s ; ETA = 0.20s Iter = 119 ; Dt = 0.0128s ; Time = 1.47s ; ETA = 0.19s Iter = 126 ; Dt = 0.0129s ; Time = 1.56s ; ETA = 0.19s Iter = 133 ; Dt = 0.0129s ; Time = 1.65s ; ETA = 0.18s Iter = 140 ; Dt = 0.0130s ; Time = 1.74s ; ETA = 0.18s Iter = 147 ; Dt = 0.0130s ; Time = 1.83s ; ETA = 0.17s Iter = 154 ; Dt = 0.0130s ; Time = 1.92s ; ETA = 0.16s Iter = 161 ; Dt = 0.0131s ; Time = 2.01s ; ETA = 0.16s Iter = 168 ; Dt = 0.0131s ; Time = 2.11s ; ETA = 0.15s Iter = 175 ; Dt = 0.0131s ; Time = 2.20s ; ETA = 0.15s Iter = 182 ; Dt = 0.0131s ; Time = 2.29s ; ETA = 0.14s Iter = 189 ; Dt = 0.0131s ; Time = 2.38s ; ETA = 0.14s Iter = 196 ; Dt = 0.0131s ; Time = 2.47s ; ETA = 0.13s Iter = 203 ; Dt = 0.0131s ; Time = 2.56s ; ETA = 0.13s Iter = 210 ; Dt = 0.0131s ; Time = 2.66s ; ETA = 0.12s Iter = 217 ; Dt = 0.0131s ; Time = 2.75s ; ETA = 0.12s Iter = 224 ; Dt = 0.0131s ; Time = 2.84s ; ETA = 0.11s Iter = 231 ; Dt = 0.0131s ; Time = 2.93s ; ETA = 0.11s Iter = 238 ; Dt = 0.0131s ; Time = 3.02s ; ETA = 0.10s Iter = 245 ; Dt = 0.0132s ; Time = 3.11s ; ETA = 0.09s Iter = 252 ; Dt = 0.0132s ; Time = 3.21s ; ETA = 0.08s Iter = 259 ; Dt = 0.0132s ; Time = 3.30s ; ETA = 0.08s Iter = 266 ; Dt = 0.0133s ; Time = 3.39s ; ETA = 0.07s Iter = 273 ; Dt = 0.0133s ; Time = 3.48s ; ETA = 0.06s Iter = 280 ; Dt = 0.0133s ; Time = 3.58s ; ETA = 0.06s Iter = 287 ; Dt = 0.0134s ; Time = 3.67s ; ETA = 0.05s Iter = 294 ; Dt = 0.0135s ; Time = 3.77s ; ETA = 0.05s Iter = 300 ; Dt = 0.0135s ; Time = 3.85s ; ETA = 0.04s Iter = 307 ; Dt = 0.0136s ; Time = 3.94s ; ETA = 0.04s Iter = 314 ; Dt = 0.0136s ; Time = 4.04s ; ETA = 0.03s Iter = 321 ; Dt = 0.0137s ; Time = 4.13s ; ETA = 0.02s Iter = 327 ; Dt = 0.0137s ; Time = 4.21s ; ETA = 0.02s Iter = 334 ; Dt = 0.0138s ; Time = 4.31s ; ETA = 0.01s Iter = 341 ; Dt = 0.0138s ; Time = 4.41s ; ETA = 0.01s Iter = 347 ; Dt = 0.0096s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.2960951328277588s (wall time) SOLVING CASE 5 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=1848, Layers=1, Triangles=3538, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 8 ; Dt = 0.0121s ; Time = 0.10s ; ETA = 1.23s Iter = 16 ; Dt = 0.0121s ; Time = 0.19s ; ETA = 1.19s Iter = 23 ; Dt = 0.0121s ; Time = 0.28s ; ETA = 1.17s Iter = 31 ; Dt = 0.0121s ; Time = 0.37s ; ETA = 1.17s Iter = 38 ; Dt = 0.0121s ; Time = 0.46s ; ETA = 1.16s Iter = 46 ; Dt = 0.0121s ; Time = 0.56s ; ETA = 1.13s Iter = 54 ; Dt = 0.0121s ; Time = 0.65s ; ETA = 1.10s Iter = 61 ; Dt = 0.0121s ; Time = 0.74s ; ETA = 1.08s Iter = 69 ; Dt = 0.0121s ; Time = 0.83s ; ETA = 1.06s Iter = 76 ; Dt = 0.0119s ; Time = 0.92s ; ETA = 1.04s Iter = 84 ; Dt = 0.0115s ; Time = 1.01s ; ETA = 1.05s Iter = 93 ; Dt = 0.0100s ; Time = 1.11s ; ETA = 1.22s Iter = 102 ; Dt = 0.0093s ; Time = 1.19s ; ETA = 1.24s Iter = 112 ; Dt = 0.0084s ; Time = 1.28s ; ETA = 1.32s Iter = 124 ; Dt = 0.0082s ; Time = 1.38s ; ETA = 1.31s Iter = 135 ; Dt = 0.0078s ; Time = 1.47s ; ETA = 1.37s Iter = 147 ; Dt = 0.0075s ; Time = 1.56s ; ETA = 1.37s Iter = 159 ; Dt = 0.0068s ; Time = 1.65s ; ETA = 1.47s Iter = 172 ; Dt = 0.0082s ; Time = 1.74s ; ETA = 1.22s Iter = 184 ; Dt = 0.0075s ; Time = 1.83s ; ETA = 1.25s Iter = 194 ; Dt = 0.0092s ; Time = 1.92s ; ETA = 0.99s Iter = 204 ; Dt = 0.0093s ; Time = 2.02s ; ETA = 0.94s Iter = 214 ; Dt = 0.0090s ; Time = 2.11s ; ETA = 0.96s Iter = 224 ; Dt = 0.0086s ; Time = 2.20s ; ETA = 0.98s Iter = 235 ; Dt = 0.0078s ; Time = 2.29s ; ETA = 1.06s Iter = 246 ; Dt = 0.0094s ; Time = 2.38s ; ETA = 0.83s Iter = 256 ; Dt = 0.0098s ; Time = 2.48s ; ETA = 0.75s Iter = 265 ; Dt = 0.0095s ; Time = 2.57s ; ETA = 0.75s Iter = 275 ; Dt = 0.0095s ; Time = 2.66s ; ETA = 0.72s Iter = 284 ; Dt = 0.0100s ; Time = 2.75s ; ETA = 0.64s Iter = 293 ; Dt = 0.0109s ; Time = 2.85s ; ETA = 0.57s Iter = 302 ; Dt = 0.0101s ; Time = 2.94s ; ETA = 0.57s Iter = 310 ; Dt = 0.0103s ; Time = 3.02s ; ETA = 0.54s Iter = 320 ; Dt = 0.0093s ; Time = 3.12s ; ETA = 0.56s Iter = 330 ; Dt = 0.0099s ; Time = 3.21s ; ETA = 0.50s Iter = 339 ; Dt = 0.0087s ; Time = 3.30s ; ETA = 0.51s Iter = 349 ; Dt = 0.0101s ; Time = 3.39s ; ETA = 0.41s Iter = 358 ; Dt = 0.0104s ; Time = 3.48s ; ETA = 0.36s Iter = 367 ; Dt = 0.0095s ; Time = 3.57s ; ETA = 0.37s Iter = 377 ; Dt = 0.0090s ; Time = 3.66s ; ETA = 0.34s Iter = 388 ; Dt = 0.0088s ; Time = 3.76s ; ETA = 0.31s Iter = 399 ; Dt = 0.0076s ; Time = 3.85s ; ETA = 0.32s Iter = 411 ; Dt = 0.0093s ; Time = 3.94s ; ETA = 0.22s Iter = 420 ; Dt = 0.0098s ; Time = 4.03s ; ETA = 0.18s Iter = 430 ; Dt = 0.0099s ; Time = 4.13s ; ETA = 0.14s Iter = 439 ; Dt = 0.0102s ; Time = 4.22s ; ETA = 0.10s Iter = 448 ; Dt = 0.0104s ; Time = 4.31s ; ETA = 0.06s Iter = 456 ; Dt = 0.0107s ; Time = 4.40s ; ETA = 0.03s Iter = 466 ; Dt = 0.0009s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 1.703859806060791s (wall time) SOLVING CASE 6 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=2914, Layers=1, Triangles=5630, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 10 ; Dt = 0.0096s ; Time = 0.10s ; ETA = 0.46s Iter = 20 ; Dt = 0.0096s ; Time = 0.19s ; ETA = 0.48s Iter = 29 ; Dt = 0.0096s ; Time = 0.28s ; ETA = 0.45s Iter = 39 ; Dt = 0.0097s ; Time = 0.37s ; ETA = 0.45s Iter = 48 ; Dt = 0.0097s ; Time = 0.46s ; ETA = 0.45s Iter = 58 ; Dt = 0.0097s ; Time = 0.56s ; ETA = 0.44s Iter = 67 ; Dt = 0.0097s ; Time = 0.65s ; ETA = 0.45s Iter = 76 ; Dt = 0.0098s ; Time = 0.73s ; ETA = 0.42s Iter = 86 ; Dt = 0.0098s ; Time = 0.83s ; ETA = 0.41s Iter = 95 ; Dt = 0.0098s ; Time = 0.92s ; ETA = 0.41s Iter = 104 ; Dt = 0.0099s ; Time = 1.01s ; ETA = 0.40s Iter = 114 ; Dt = 0.0099s ; Time = 1.11s ; ETA = 0.39s Iter = 123 ; Dt = 0.0100s ; Time = 1.20s ; ETA = 0.38s Iter = 132 ; Dt = 0.0100s ; Time = 1.29s ; ETA = 0.36s Iter = 141 ; Dt = 0.0100s ; Time = 1.38s ; ETA = 0.36s Iter = 150 ; Dt = 0.0101s ; Time = 1.47s ; ETA = 0.35s Iter = 159 ; Dt = 0.0101s ; Time = 1.56s ; ETA = 0.33s Iter = 168 ; Dt = 0.0102s ; Time = 1.65s ; ETA = 0.34s Iter = 177 ; Dt = 0.0102s ; Time = 1.74s ; ETA = 0.32s Iter = 186 ; Dt = 0.0102s ; Time = 1.83s ; ETA = 0.31s Iter = 195 ; Dt = 0.0103s ; Time = 1.93s ; ETA = 0.30s Iter = 204 ; Dt = 0.0103s ; Time = 2.02s ; ETA = 0.30s Iter = 213 ; Dt = 0.0103s ; Time = 2.11s ; ETA = 0.27s Iter = 222 ; Dt = 0.0103s ; Time = 2.20s ; ETA = 0.27s Iter = 231 ; Dt = 0.0103s ; Time = 2.30s ; ETA = 0.26s Iter = 240 ; Dt = 0.0103s ; Time = 2.39s ; ETA = 0.25s Iter = 249 ; Dt = 0.0103s ; Time = 2.48s ; ETA = 0.24s Iter = 258 ; Dt = 0.0103s ; Time = 2.57s ; ETA = 0.23s Iter = 266 ; Dt = 0.0103s ; Time = 2.66s ; ETA = 0.23s Iter = 275 ; Dt = 0.0103s ; Time = 2.75s ; ETA = 0.21s Iter = 284 ; Dt = 0.0103s ; Time = 2.84s ; ETA = 0.20s Iter = 293 ; Dt = 0.0103s ; Time = 2.93s ; ETA = 0.19s Iter = 302 ; Dt = 0.0103s ; Time = 3.03s ; ETA = 0.17s Iter = 311 ; Dt = 0.0103s ; Time = 3.12s ; ETA = 0.16s Iter = 320 ; Dt = 0.0103s ; Time = 3.21s ; ETA = 0.16s Iter = 329 ; Dt = 0.0103s ; Time = 3.30s ; ETA = 0.14s Iter = 338 ; Dt = 0.0104s ; Time = 3.40s ; ETA = 0.13s Iter = 346 ; Dt = 0.0104s ; Time = 3.48s ; ETA = 0.12s Iter = 355 ; Dt = 0.0104s ; Time = 3.57s ; ETA = 0.11s Iter = 364 ; Dt = 0.0105s ; Time = 3.67s ; ETA = 0.10s Iter = 373 ; Dt = 0.0105s ; Time = 3.76s ; ETA = 0.09s Iter = 381 ; Dt = 0.0105s ; Time = 3.85s ; ETA = 0.08s Iter = 390 ; Dt = 0.0106s ; Time = 3.94s ; ETA = 0.07s Iter = 399 ; Dt = 0.0106s ; Time = 4.04s ; ETA = 0.05s Iter = 407 ; Dt = 0.0106s ; Time = 4.12s ; ETA = 0.04s Iter = 416 ; Dt = 0.0107s ; Time = 4.22s ; ETA = 0.03s Iter = 425 ; Dt = 0.0107s ; Time = 4.31s ; ETA = 0.02s Iter = 433 ; Dt = 0.0107s ; Time = 4.40s ; ETA = 0.01s Iter = 442 ; Dt = 0.0010s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.5304944515228271s (wall time) SOLVING CASE 7 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=2914, Layers=1, Triangles=5630, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 10 ; Dt = 0.0096s ; Time = 0.10s ; ETA = 2.32s Iter = 20 ; Dt = 0.0096s ; Time = 0.19s ; ETA = 2.31s Iter = 29 ; Dt = 0.0096s ; Time = 0.28s ; ETA = 2.31s Iter = 39 ; Dt = 0.0096s ; Time = 0.37s ; ETA = 2.31s Iter = 48 ; Dt = 0.0096s ; Time = 0.46s ; ETA = 2.24s Iter = 58 ; Dt = 0.0096s ; Time = 0.56s ; ETA = 2.19s Iter = 67 ; Dt = 0.0096s ; Time = 0.64s ; ETA = 2.19s Iter = 77 ; Dt = 0.0096s ; Time = 0.74s ; ETA = 2.10s Iter = 86 ; Dt = 0.0096s ; Time = 0.83s ; ETA = 2.05s Iter = 96 ; Dt = 0.0093s ; Time = 0.92s ; ETA = 2.09s Iter = 106 ; Dt = 0.0084s ; Time = 1.01s ; ETA = 2.25s Iter = 118 ; Dt = 0.0076s ; Time = 1.11s ; ETA = 2.37s Iter = 130 ; Dt = 0.0069s ; Time = 1.19s ; ETA = 2.53s Iter = 144 ; Dt = 0.0064s ; Time = 1.29s ; ETA = 2.69s Iter = 159 ; Dt = 0.0060s ; Time = 1.38s ; ETA = 2.75s Iter = 174 ; Dt = 0.0058s ; Time = 1.47s ; ETA = 2.84s Iter = 190 ; Dt = 0.0056s ; Time = 1.56s ; ETA = 2.88s Iter = 206 ; Dt = 0.0055s ; Time = 1.65s ; ETA = 2.83s Iter = 221 ; Dt = 0.0069s ; Time = 1.74s ; ETA = 2.16s Iter = 234 ; Dt = 0.0071s ; Time = 1.83s ; ETA = 2.09s Iter = 246 ; Dt = 0.0077s ; Time = 1.93s ; ETA = 1.89s Iter = 258 ; Dt = 0.0074s ; Time = 2.02s ; ETA = 1.86s Iter = 271 ; Dt = 0.0072s ; Time = 2.11s ; ETA = 1.83s Iter = 284 ; Dt = 0.0069s ; Time = 2.20s ; ETA = 1.84s Iter = 297 ; Dt = 0.0072s ; Time = 2.29s ; ETA = 1.72s Iter = 310 ; Dt = 0.0065s ; Time = 2.38s ; ETA = 1.81s Iter = 323 ; Dt = 0.0072s ; Time = 2.47s ; ETA = 1.57s Iter = 337 ; Dt = 0.0072s ; Time = 2.57s ; ETA = 1.55s Iter = 349 ; Dt = 0.0081s ; Time = 2.66s ; ETA = 1.31s Iter = 360 ; Dt = 0.0075s ; Time = 2.75s ; ETA = 1.52s Iter = 372 ; Dt = 0.0075s ; Time = 2.84s ; ETA = 1.41s Iter = 385 ; Dt = 0.0077s ; Time = 2.93s ; ETA = 1.21s Iter = 397 ; Dt = 0.0077s ; Time = 3.02s ; ETA = 1.17s Iter = 409 ; Dt = 0.0074s ; Time = 3.11s ; ETA = 1.18s Iter = 422 ; Dt = 0.0073s ; Time = 3.21s ; ETA = 1.12s Iter = 434 ; Dt = 0.0075s ; Time = 3.30s ; ETA = 1.04s Iter = 446 ; Dt = 0.0084s ; Time = 3.39s ; ETA = 0.85s Iter = 456 ; Dt = 0.0085s ; Time = 3.48s ; ETA = 0.71s Iter = 468 ; Dt = 0.0083s ; Time = 3.58s ; ETA = 0.64s Iter = 479 ; Dt = 0.0077s ; Time = 3.67s ; ETA = 0.61s Iter = 491 ; Dt = 0.0076s ; Time = 3.76s ; ETA = 0.55s Iter = 503 ; Dt = 0.0071s ; Time = 3.85s ; ETA = 0.53s Iter = 516 ; Dt = 0.0071s ; Time = 3.94s ; ETA = 0.44s Iter = 529 ; Dt = 0.0075s ; Time = 4.03s ; ETA = 0.39s Iter = 541 ; Dt = 0.0083s ; Time = 4.13s ; ETA = 0.29s Iter = 552 ; Dt = 0.0080s ; Time = 4.22s ; ETA = 0.22s Iter = 564 ; Dt = 0.0077s ; Time = 4.31s ; ETA = 0.16s Iter = 576 ; Dt = 0.0070s ; Time = 4.40s ; ETA = 0.08s Iter = 589 ; Dt = 0.0019s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 3.4363200664520264s (wall time) SOLVING CASE 8 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=4155, Layers=1, Triangles=8072, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 12 ; Dt = 0.0080s ; Time = 0.10s ; ETA = 0.88s Iter = 23 ; Dt = 0.0080s ; Time = 0.18s ; ETA = 0.88s Iter = 35 ; Dt = 0.0080s ; Time = 0.28s ; ETA = 0.84s Iter = 46 ; Dt = 0.0080s ; Time = 0.37s ; ETA = 0.72s Iter = 58 ; Dt = 0.0080s ; Time = 0.46s ; ETA = 1.06s Iter = 69 ; Dt = 0.0080s ; Time = 0.55s ; ETA = 0.75s Iter = 81 ; Dt = 0.0081s ; Time = 0.65s ; ETA = 0.72s Iter = 92 ; Dt = 0.0081s ; Time = 0.74s ; ETA = 0.66s Iter = 103 ; Dt = 0.0081s ; Time = 0.83s ; ETA = 0.66s Iter = 114 ; Dt = 0.0081s ; Time = 0.92s ; ETA = 0.65s Iter = 126 ; Dt = 0.0082s ; Time = 1.01s ; ETA = 0.63s Iter = 137 ; Dt = 0.0082s ; Time = 1.10s ; ETA = 0.61s Iter = 148 ; Dt = 0.0082s ; Time = 1.19s ; ETA = 0.59s Iter = 159 ; Dt = 0.0083s ; Time = 1.28s ; ETA = 0.58s Iter = 170 ; Dt = 0.0083s ; Time = 1.38s ; ETA = 0.59s Iter = 181 ; Dt = 0.0083s ; Time = 1.47s ; ETA = 0.56s Iter = 192 ; Dt = 0.0084s ; Time = 1.56s ; ETA = 0.54s Iter = 203 ; Dt = 0.0084s ; Time = 1.65s ; ETA = 0.53s Iter = 214 ; Dt = 0.0084s ; Time = 1.74s ; ETA = 0.52s Iter = 225 ; Dt = 0.0084s ; Time = 1.84s ; ETA = 0.51s Iter = 236 ; Dt = 0.0084s ; Time = 1.93s ; ETA = 0.50s Iter = 247 ; Dt = 0.0085s ; Time = 2.02s ; ETA = 0.48s Iter = 257 ; Dt = 0.0085s ; Time = 2.11s ; ETA = 0.45s Iter = 268 ; Dt = 0.0084s ; Time = 2.20s ; ETA = 0.44s Iter = 279 ; Dt = 0.0084s ; Time = 2.29s ; ETA = 0.43s Iter = 290 ; Dt = 0.0084s ; Time = 2.39s ; ETA = 0.41s Iter = 301 ; Dt = 0.0084s ; Time = 2.48s ; ETA = 0.39s Iter = 312 ; Dt = 0.0084s ; Time = 2.57s ; ETA = 0.38s Iter = 323 ; Dt = 0.0084s ; Time = 2.66s ; ETA = 0.37s Iter = 333 ; Dt = 0.0084s ; Time = 2.75s ; ETA = 0.35s Iter = 344 ; Dt = 0.0084s ; Time = 2.84s ; ETA = 0.35s Iter = 355 ; Dt = 0.0084s ; Time = 2.93s ; ETA = 0.32s Iter = 366 ; Dt = 0.0084s ; Time = 3.03s ; ETA = 0.31s Iter = 377 ; Dt = 0.0084s ; Time = 3.12s ; ETA = 0.29s Iter = 388 ; Dt = 0.0085s ; Time = 3.21s ; ETA = 0.27s Iter = 398 ; Dt = 0.0085s ; Time = 3.30s ; ETA = 0.25s Iter = 409 ; Dt = 0.0085s ; Time = 3.39s ; ETA = 0.22s Iter = 420 ; Dt = 0.0085s ; Time = 3.48s ; ETA = 0.20s Iter = 431 ; Dt = 0.0085s ; Time = 3.58s ; ETA = 0.19s Iter = 441 ; Dt = 0.0086s ; Time = 3.66s ; ETA = 0.17s Iter = 452 ; Dt = 0.0086s ; Time = 3.76s ; ETA = 0.15s Iter = 463 ; Dt = 0.0086s ; Time = 3.85s ; ETA = 0.14s Iter = 473 ; Dt = 0.0086s ; Time = 3.94s ; ETA = 0.11s Iter = 484 ; Dt = 0.0087s ; Time = 4.03s ; ETA = 0.10s Iter = 494 ; Dt = 0.0087s ; Time = 4.12s ; ETA = 0.07s Iter = 505 ; Dt = 0.0087s ; Time = 4.22s ; ETA = 0.06s Iter = 515 ; Dt = 0.0088s ; Time = 4.30s ; ETA = 0.04s Iter = 526 ; Dt = 0.0088s ; Time = 4.40s ; ETA = 0.02s Iter = 536 ; Dt = 0.0072s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 0.8902091979980469s (wall time) SOLVING CASE 9 =================================================================== | INITIALIZATION | =================================================================== Problem size: Nodes=4155, Layers=1, Triangles=8072, Iter = 0 ; Dt = 0.0000s ; Time = 0.00s ; ETA = 0.00s =================================================================== | TIME LOOP | =================================================================== Iter = 12 ; Dt = 0.0080s ; Time = 0.10s ; ETA = 4.19s Iter = 23 ; Dt = 0.0080s ; Time = 0.18s ; ETA = 4.16s Iter = 35 ; Dt = 0.0080s ; Time = 0.28s ; ETA = 4.04s Iter = 46 ; Dt = 0.0080s ; Time = 0.37s ; ETA = 3.99s Iter = 58 ; Dt = 0.0080s ; Time = 0.46s ; ETA = 3.89s Iter = 69 ; Dt = 0.0080s ; Time = 0.55s ; ETA = 3.82s Iter = 81 ; Dt = 0.0080s ; Time = 0.65s ; ETA = 3.75s Iter = 92 ; Dt = 0.0080s ; Time = 0.73s ; ETA = 3.68s Iter = 104 ; Dt = 0.0080s ; Time = 0.83s ; ETA = 3.55s Iter = 116 ; Dt = 0.0073s ; Time = 0.92s ; ETA = 3.73s Iter = 129 ; Dt = 0.0067s ; Time = 1.01s ; ETA = 3.98s Iter = 143 ; Dt = 0.0061s ; Time = 1.10s ; ETA = 4.26s Iter = 158 ; Dt = 0.0056s ; Time = 1.19s ; ETA = 4.55s Iter = 176 ; Dt = 0.0049s ; Time = 1.28s ; ETA = 5.15s Iter = 195 ; Dt = 0.0047s ; Time = 1.37s ; ETA = 5.26s Iter = 214 ; Dt = 0.0046s ; Time = 1.46s ; ETA = 5.19s Iter = 234 ; Dt = 0.0050s ; Time = 1.56s ; ETA = 4.68s Iter = 251 ; Dt = 0.0052s ; Time = 1.65s ; ETA = 4.36s Iter = 268 ; Dt = 0.0061s ; Time = 1.74s ; ETA = 3.64s Iter = 283 ; Dt = 0.0065s ; Time = 1.84s ; ETA = 3.30s Iter = 296 ; Dt = 0.0067s ; Time = 1.92s ; ETA = 3.05s Iter = 310 ; Dt = 0.0065s ; Time = 2.01s ; ETA = 3.12s Iter = 326 ; Dt = 0.0060s ; Time = 2.11s ; ETA = 3.35s Iter = 341 ; Dt = 0.0061s ; Time = 2.20s ; ETA = 3.10s Iter = 356 ; Dt = 0.0056s ; Time = 2.29s ; ETA = 3.28s Iter = 372 ; Dt = 0.0057s ; Time = 2.38s ; ETA = 3.05s Iter = 388 ; Dt = 0.0060s ; Time = 2.47s ; ETA = 2.82s Iter = 403 ; Dt = 0.0064s ; Time = 2.56s ; ETA = 2.49s Iter = 417 ; Dt = 0.0065s ; Time = 2.66s ; ETA = 2.40s Iter = 431 ; Dt = 0.0063s ; Time = 2.75s ; ETA = 2.35s Iter = 447 ; Dt = 0.0056s ; Time = 2.84s ; ETA = 2.47s Iter = 462 ; Dt = 0.0063s ; Time = 2.93s ; ETA = 2.05s Iter = 477 ; Dt = 0.0062s ; Time = 3.03s ; ETA = 1.97s Iter = 491 ; Dt = 0.0063s ; Time = 3.12s ; ETA = 1.87s Iter = 505 ; Dt = 0.0064s ; Time = 3.21s ; ETA = 1.68s Iter = 519 ; Dt = 0.0075s ; Time = 3.30s ; ETA = 1.38s Iter = 530 ; Dt = 0.0077s ; Time = 3.39s ; ETA = 1.22s Iter = 543 ; Dt = 0.0076s ; Time = 3.49s ; ETA = 1.15s Iter = 555 ; Dt = 0.0072s ; Time = 3.57s ; ETA = 1.09s Iter = 569 ; Dt = 0.0065s ; Time = 3.67s ; ETA = 1.07s Iter = 583 ; Dt = 0.0060s ; Time = 3.75s ; ETA = 1.05s Iter = 599 ; Dt = 0.0060s ; Time = 3.85s ; ETA = 0.91s Iter = 614 ; Dt = 0.0061s ; Time = 3.94s ; ETA = 0.78s Iter = 629 ; Dt = 0.0063s ; Time = 4.03s ; ETA = 0.60s Iter = 642 ; Dt = 0.0073s ; Time = 4.12s ; ETA = 0.43s Iter = 655 ; Dt = 0.0074s ; Time = 4.22s ; ETA = 0.31s Iter = 668 ; Dt = 0.0070s ; Time = 4.31s ; ETA = 0.22s Iter = 681 ; Dt = 0.0063s ; Time = 4.40s ; ETA = 0.12s Iter = 696 ; Dt = 0.0037s ; Time = 4.49s ; ETA = 0.00s =================================================================== | END | =================================================================== Problem.solve() completed in 5.732445955276489s (wall time) .. GENERATED FROM PYTHON SOURCE LINES 188-192 Errors: -------------------- Errors computed during solving loop are stored in ``error_*`` matrix. .. GENERATED FROM PYTHON SOURCE LINES 193-210 .. code-block:: default error_H = np.zeros((N_mesh, 2)) error_QX = np.zeros((N_mesh, 2)) error_QY = np.zeros((N_mesh, 2)) error_H_cumul = np.zeros((N_mesh, 2)) error_QX_cumul = np.zeros((N_mesh, 2)) error_QY_cumul = np.zeros((N_mesh, 2)) for I, case in enumerate(cases): sol = problem_list[I].analytic_sol error_H[case.mesh_id, case.order-1] = sol.error_H error_QX[case.mesh_id, case.order-1] = sol.error_QX error_QY[case.mesh_id, case.order-1] = sol.error_QY error_H_cumul[case.mesh_id, case.order-1] = sol.error_H_cumul/FINAL_TIME error_QX_cumul[case.mesh_id, case.order-1] = sol.error_QX_cumul/FINAL_TIME error_QY_cumul[case.mesh_id, case.order-1] = sol.error_QY_cumul/FINAL_TIME .. GENERATED FROM PYTHON SOURCE LINES 211-212 Errors are save in h5 format in \data repertory. .. GENERATED FROM PYTHON SOURCE LINES 213-227 .. code-block:: default if SAVE_ERROR: os.system('rm data/thacker2d_errors.h5') REF_FILE = "data/thacker2d_errors.h5" with h5py.File(REF_FILE, 'w') as f: f.create_dataset('NT', data=NT) f.create_dataset('NL', data=NL) f.create_dataset('error_H', data=error_H) f.create_dataset('error_QX', data=error_QX) f.create_dataset('error_QY', data=error_QY) f.create_dataset('error_H_cumul', data=error_H_cumul) f.create_dataset('error_QX_cumul', data=error_QX_cumul) f.create_dataset('error_QY_cumul', data=error_QY_cumul) .. GENERATED FROM PYTHON SOURCE LINES 228-231 Errors plots: -------------------- .. GENERATED FROM PYTHON SOURCE LINES 232-288 .. code-block:: default def plot_error(ax, error): err = np.zeros((N_mesh)) ab = np.zeros((N_mesh)) for order in range(2): for M in range(N_mesh): err[M] = np.log(error[M, order]) ab[M] = 0.5*np.log(NT[M]/NT[0]) ax.plot(ab, err, color=colors[order], marker=markers[order], label=labels[order]) def plot_reference(ax, error): xmax = 0.5*np.log(NT[N_mesh-1]/NT[0]) ax.plot([0.0, xmax], [np.log(error[0, 0]), -1.0*xmax+np.log(error[0, 0])],\ color='k', ls=':', lw=0.5) ax.plot([0.0, xmax], [np.log(error[0, 1]), -2.0*xmax+np.log(error[0, 1])],\ color='k', ls='-.', lw=0.5) if PLOT_ERROR == True: plt.rcParams["figure.figsize"] = [10, 5] colors = ['blue','red','green','purple'] markers = ['s','*','o','v'] labels = ['$1^{st}$ order','$2^{nd}$ order'] fig = plt.figure() ax0 = fig.add_subplot(111) plot_error(ax0, error_H_cumul) plot_reference(ax0, error_H_cumul) ax0.set_xlabel('$log(l_0/l_i)$') ax0.set_ylabel('Error of $h$ in $L_2$ norm') ax0.set_title("Convergence of $h$") plt.legend(loc=1) plt.show() fig = plt.figure() ax1 = fig.add_subplot(121) plot_error(ax1, error_QX_cumul) plot_reference(ax1, error_QX_cumul) ax1.set_xlabel('$log(l_0/l_i)$') ax1.set_ylabel('Error of $hu$ in $L_2$ norm') ax1.set_title("Convergence of $hu$") plt.legend(loc=3) ax2 = fig.add_subplot(122) plot_error(ax2, error_QY_cumul) plot_reference(ax2, error_QY_cumul) ax2.set_xlabel('$log(l_0/l_i)$') ax2.set_ylabel('Error of $hv$ in $L_2$ norm') ax2.set_title("Convergence of $hv$") plt.legend(loc=3) plt.show() os.system('rm -r outputs') .. rst-class:: sphx-glr-horizontal * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_006.png :alt: Convergence of $h$ :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_006.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_007.png :alt: Convergence of $hu$, Convergence of $hv$ :srcset: /auto_examples_mesh/images/sphx_glr_example_meshconvergence_thacker2d_007.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none 0 .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 17.535 seconds) .. _sphx_glr_download_auto_examples_mesh_example_meshconvergence_thacker2d.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: example_meshconvergence_thacker2d.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: example_meshconvergence_thacker2d.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_