Note
Go to the end to download the full example code.
Testing 3D circularly polarised Alfven wave with SPMHD#
CI test for the 3D circularly polarised Alfven wave (Toth 2000, 3D setup of Gardiner & Stone 2008), following the description given in the Phantom code paper (Price et al. 2018, Section 5.6.2 / 5.7.1).
11 import os
12
13 import matplotlib.pyplot as plt
14 import numpy as np
15
16 import shamrock
17
18 # If we use the shamrock executable to run this script instead of the python interpreter,
19 # we should not initialize the system as the shamrock executable needs to handle specific MPI logic
20 if not shamrock.sys.is_initialized():
21 shamrock.change_loglevel(1)
22 shamrock.sys.init("0:0")
25 gamma = 5.0 / 3.0
26 rho0 = 1.0
27 P0 = 0.1
28
29 wavelength = 1.0
30 amplitude = 0.1
31
32 # Alfven speed set by B1 = 1, rho = 1 -> v_A = 1, period = wavelength / v_A = 1
33 v_A = 1.0
34 n_periods = 5
35 t_target = n_periods * wavelength / v_A
36
37 resol = 32
38
39 plot_extra_resols = [] # 64, 128, too much for the ci
40
41 do_plot = False
42 dump_folder = "_to_trash"
43 if do_plot and shamrock.sys.world_rank() == 0:
44 os.makedirs(dump_folder, exist_ok=True)
48 def best_transverse_counts(model, xcnt, target_ratio=0.5, search_frac=0.18):
49
50 lo = max(2, int(xcnt * (target_ratio - search_frac)))
51 hi = int(xcnt * (target_ratio + search_frac)) + 1
52 lo_y = lo + lo % 2
53 lo_z = max(3, lo + (-lo) % 3)
54 best = None
55 for ycnt in range(lo_y, hi, 2):
56 for zcnt in range(lo_z, hi, 3):
57 xs, ys, zs = model.get_box_dim_true_fcc_3d(1, xcnt, ycnt, zcnt)
58 score = abs(ys / xs - target_ratio) + abs(zs / xs - target_ratio)
59 if best is None or score < best[0]:
60 best = (score, ycnt, zcnt)
61 return best[1], best[2]
62
63
64 def wave_basis(xs, ys, zs):
65 # Athena-style: r_hat proportional to (1/Lx, 1/Ly, 1/Lz) and lambda = 1/|(1/Lx, 1/Ly, 1/Lz)|,
66 # so a periodic translation along any axis shifts x1 by exactly one
67 # wavelength. A lattice-derived box is only approximately 2:1:1, and
68 # keeping the paper's fixed angles would leave a phase seam at the y/z
69 # boundaries. For an exact 3 x 1.5 x 1.5 box this reduces to
70 # sin(a) = 2/3, sin(b) = 2/sqrt(5), lambda = 1.
71 k = np.array([1.0 / xs, 1.0 / ys, 1.0 / zs])
72 lam = 1.0 / np.linalg.norm(k)
73 r_hat = k * lam
74
75 sin_a = r_hat[2]
76 cos_a = np.hypot(r_hat[0], r_hat[1])
77 sin_b = r_hat[1] / cos_a
78 cos_b = r_hat[0] / cos_a
79
80 e2_hat = np.array([-sin_b, cos_b, 0.0])
81 e3_hat = np.array([-sin_a * cos_b, -sin_a * sin_b, cos_a])
82 return lam, r_hat, e2_hat, e3_hat
83
84
85 def run_alfven_wave(resol):
86 ctx = shamrock.Context()
87 ctx.pdata_layout_new()
88
89 codeu = shamrock.UnitSystem(unit_time=1.0, unit_length=1.0, unit_mass=1.2566370621219e-06)
90
91 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
92
93 cfg = model.gen_default_config()
94 cfg.set_units(codeu)
95 cfg.set_artif_viscosity_None()
96 cfg.set_IdealMHD(sigma_mhd=1, sigma_u=1)
97 cfg.set_boundary_periodic()
98 cfg.set_eos_adiabatic(gamma)
99 model.set_solver_config(cfg)
100
101 crit_split = int(1e7)
102 crit_merge = 1
103 model.init_scheduler(crit_split, crit_merge)
104
105 ycnt, zcnt = best_transverse_counts(model, resol)
106
107 # lambda scales linearly with dr: pick dr so that lambda == wavelength exactly
108 lam_unit = wave_basis(*model.get_box_dim_true_fcc_3d(1, resol, ycnt, zcnt))[0]
109 dr = wavelength / lam_unit
110 (xs, ys, zs) = model.get_box_dim_true_fcc_3d(dr, resol, ycnt, zcnt)
111 lam, r_hat, e2_hat, e3_hat = wave_basis(xs, ys, zs)
112 rotmat = np.column_stack([r_hat, e2_hat, e3_hat])
113 print(f"[resol={resol}] Box dims: xs={xs} ys={ys} zs={zs}")
114 print(f"[resol={resol}] r_hat={r_hat} (paper: [1/3, 2/3, 2/3]) lambda={lam}")
115
116 def vel_func(r):
117 x1 = np.dot(r, r_hat)
118 v_wave = np.array(
119 [
120 0.0,
121 amplitude * np.sin(2.0 * np.pi * x1 / lam),
122 amplitude * np.cos(2.0 * np.pi * x1 / lam),
123 ]
124 )
125 return tuple(rotmat @ v_wave)
126
127 def mag_func(r):
128 x1 = np.dot(r, r_hat)
129 B_wave = np.array(
130 [
131 1.0,
132 amplitude * np.sin(2.0 * np.pi * x1 / lam),
133 amplitude * np.cos(2.0 * np.pi * x1 / lam),
134 ]
135 )
136 # field is stored as B/rho in SPMHD
137 return tuple(rotmat @ B_wave / rho0)
138
139 box_min = (0.0, 0.0, 0.0)
140 box_max = (xs, ys, zs)
141
142 model.resize_simulation_box(box_min, box_max)
143
144 setup = model.get_setup()
145 gen = setup.make_generator_lattice_fcc(dr, box_min, box_max)
146 setup.apply_setup(gen)
147
148 gam1 = gamma - 1.0
149 uuzero = P0 / (gam1 * rho0)
150 model.set_value_in_a_box("uint", "f64", uuzero, box_min, box_max)
151
152 model.set_field_value_lambda_f64_3("vxyz", vel_func)
153 model.set_field_value_lambda_f64_3("B/rho", mag_func)
154
155 vol_b = xs * ys * zs
156 totmass = rho0 * vol_b
157 pmass = model.total_mass_to_part_mass(totmass)
158 model.set_particle_mass(pmass)
159
160 print(f"[resol={resol}] Total mass :", totmass)
161 print(f"[resol={resol}] Current part mass :", pmass)
162
163 model.set_cfl_cour(0.3)
164 model.set_cfl_force(0.25)
165
166 model.timestep()
167
168 model.evolve_until(t_target)
169
170 # %%
171 # Compare the transverse field component B2 against the exact solution.
172 #
173 # t_target is an integer number of wave periods, the exact
174 # (undamped) solution coincides with the initial condition:
175 # B2_exact(x1) = amplitude * sin(2*pi*x1/wavelength)
176
177 data = ctx.collect_data()
178
179 xyz = data["xyz"]
180 B_on_rho = data["B/rho"]
181
182 x1 = xyz @ r_hat
183 # SPMHD stores B/rho: recover the physical field before projecting it on e2
184 rho = pmass * (model.get_hfact() / data["hpart"]) ** 3
185 B2 = (B_on_rho * rho[:, None]) @ e2_hat
186
187 return x1, B2
188
189
190 x1, B2 = run_alfven_wave(resol)
191
192 B2_exact = amplitude * np.sin(2.0 * np.pi * x1 / wavelength)
193
194 l2_err_B2 = np.sqrt(np.mean((B2 - B2_exact) ** 2))
195
196 print(f"L2 error on B2 : {l2_err_B2}")
[resol=32] Box dims: xs=2.949065100472731 ys=1.4366054781132147 zs=1.5801854062767053
[resol=32] r_hat=[0.33909051 0.69608533 0.63283713] (paper: [1/3, 2/3, 2/3]) lambda=1.0
SPH setup: generating particles ...
SPH setup: Nstep = 12096 ( 1.2e+04 ) Ntotal = 12096 ( 1.2e+04 rank min = 4.4e+06 max = 1.2e+04) rate = 1.209600e+04 N.s^-1
SPH setup: the generation step took : 0.008611486 s
SPH setup: final particle count = 12096 beginning injection ...
SPH setup: injected 12096 / 12096 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005715532000000001 s
SPH setup: the setup took : 0.018698562000000002 s
[resol=32] Total mass : 6.694681564488508
[resol=32] Current part mass : 0.0005534624309266292
---------------- t = 0, dt = 0 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.09031406717174761 unconverged cnt = 12096
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0, dt = 0.00027288947639721585 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.00027288947639721585, dt = 0.009278239649430232 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.009551129125827448, dt = 0.01528146050713228 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.02483258963295973, dt = 0.019282955190414144 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.044115544823373876, dt = 0.021949977978129098 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.06606552280150298, dt = 0.02372746075860039 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01110471585823701
---------------- t = 0.08979298356010337, dt = 0.012457067709741993 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.10225005126984536, dt = 0.017399207981155653 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.11964925925100102, dt = 0.02069329082211636 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.14034255007311738, dt = 0.022888800739862047 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010274610793354726
---------------- t = 0.16323135081297943, dt = 0.012176957594496085 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.1754083084074755, dt = 0.017211595348119738 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.19261990375559523, dt = 0.020567503047719204 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.21318740680331444, dt = 0.0228044187791328 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010143191023244746
---------------- t = 0.23599182558244725, dt = 0.01214857279245679 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.24814039837490404, dt = 0.017192225836487243 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.26533262421139125, dt = 0.020554349059937513 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.2858869732713288, dt = 0.022795515208774937 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010095245723794371
---------------- t = 0.3086824884801037, dt = 0.012145575330693919 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.32082806381079765, dt = 0.017190155565015224 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.3380182193758129, dt = 0.020552948137204517 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.35857116751301743, dt = 0.02279452064244343 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010100271168618394
---------------- t = 0.38136568815546085, dt = 0.012145101569391985 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.3935107897248528, dt = 0.01718936369196458 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.41070015341681737, dt = 0.020551564984792436 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.4312517184016098, dt = 0.02279251356645679 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010196964357094555
---------------- t = 0.4540442319680666, dt = 0.012143873902564231 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.4661881058706308, dt = 0.017187564442663356 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.48337567031329415, dt = 0.020549391470152424 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.5039250617834465, dt = 0.022789916628284307 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010242752694283722
---------------- t = 0.5267149784117309, dt = 0.012142304783340262 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.5388572831950711, dt = 0.017185335818788564 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.5560426190138597, dt = 0.020546706316079286 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.5765893253299389, dt = 0.022786925055819442 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010239534108047151
---------------- t = 0.5993762503857584, dt = 0.012140796349555595 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.6115170467353139, dt = 0.017183253197540334 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.6287002999328543, dt = 0.020544180856042135 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.6492444807888964, dt = 0.022784258254016344 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01019122893249363
---------------- t = 0.6720287390429127, dt = 0.012139435467851232 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.684168174510764, dt = 0.017181461586580033 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.701349636097344, dt = 0.020542311662598928 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.7218919477599429, dt = 0.022782640259317907 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010257515344900296
---------------- t = 0.7446745880192608, dt = 0.012139046059818682 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.7568136340790794, dt = 0.017181115686591714 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.7739947497656712, dt = 0.020542438304868946 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.7945371880705401, dt = 0.022783646180816696 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010359526656103334
---------------- t = 0.8173208342513568, dt = 0.012139955296903834 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.8294607895482606, dt = 0.01718275916156706 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.8466435487098277, dt = 0.020544823495888872 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.8671883722057165, dt = 0.02278673311711123 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01041937531040456
---------------- t = 0.8899751053228278, dt = 0.012141836260984556 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.9021169415838123, dt = 0.017185566499418294 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.9193025080832306, dt = 0.020548409337643816 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.9398509174208743, dt = 0.022790035667174847 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010509561097303359
---------------- t = 0.9626409530880492, dt = 0.012142821842361283 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.9747837749304105, dt = 0.01718619307503258 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 0.991969968005443, dt = 0.02054792082997508 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.012517888835418, dt = 0.02278855063240848 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010482040358084144
---------------- t = 1.0353064394678266, dt = 0.012141849883276335 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.047448289351103, dt = 0.017184711581472115 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.064633000932575, dt = 0.020546004551614153 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.0851790054841892, dt = 0.02278638708655861 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010432900718965164
---------------- t = 1.1079653925707478, dt = 0.012140602863987786 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.1201059954347357, dt = 0.017183054615445278 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.1372890500501809, dt = 0.020544193447606593 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.1578332434977874, dt = 0.022784597706691674 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010367376225040003
---------------- t = 1.180617841204479, dt = 0.012139805129656836 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.1927576463341358, dt = 0.017182165490016575 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.2099398118241524, dt = 0.020543471336045973 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.2304832831601984, dt = 0.022784302361676615 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010367534699371475
---------------- t = 1.2532675855218751, dt = 0.01213996938564839 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.2654075549075234, dt = 0.017182632273443715 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.282590187180967, dt = 0.020544346831625534 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.3031345340125926, dt = 0.022785738897343543 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010332260379935888
---------------- t = 1.3259202729099362, dt = 0.012141002189155896 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.3380612750990921, dt = 0.017184159069600486 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.3552454341686926, dt = 0.02054629138520502 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.3757917255538976, dt = 0.022788004097926303 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010357399126342457
---------------- t = 1.3985797296518239, dt = 0.012141985956846131 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.41072171560867, dt = 0.01718554790248004 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.42790726351115, dt = 0.020547776154954827 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.448455039666105, dt = 0.02278906424850468 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010385881611186868
---------------- t = 1.4712441039146096, dt = 0.01214162709882678 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.4833857310134364, dt = 0.017184404945444494 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.500570135958881, dt = 0.02054550391714589 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.5211156398760268, dt = 0.02278548375816867 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010434180899419582
---------------- t = 1.5439011236341955, dt = 0.012139665162450127 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.5560407887966456, dt = 0.01718159961604677 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.5732223884126924, dt = 0.020542199361526042 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.5937645877742184, dt = 0.022781965435594603 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010307150651296764
---------------- t = 1.6165465532098129, dt = 0.012137896442055816 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.6286844496518686, dt = 0.01717918979386756 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.6458636394457362, dt = 0.020539406054073298 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.6664030454998096, dt = 0.022779140546878653 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01031625035768625
---------------- t = 1.6891821860466882, dt = 0.012136573761204061 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.7013187598078923, dt = 0.01717754917079953 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.7184963089786918, dt = 0.02053786528928687 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.7390341742679787, dt = 0.02277787005066303 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010430285040673476
---------------- t = 1.7618120443186418, dt = 0.012136341354781073 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.7739483856734228, dt = 0.017177428720364207 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.7911258143937872, dt = 0.02053815028648102 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.8116639646802681, dt = 0.022778844760790853 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010551735996272037
---------------- t = 1.834442809441059, dt = 0.012137131383070012 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.846579940824129, dt = 0.01717886544371424 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.8637588062678432, dt = 0.020540184833545624 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.8842989911013888, dt = 0.02278167071455595 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01062063729071743
---------------- t = 1.9070806618159448, dt = 0.012138786161492996 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.9192194479774378, dt = 0.017181222348229927 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.9364006703256678, dt = 0.020543085808831348 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.9569437561344991, dt = 0.022784689391636862 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010664643309562367
---------------- t = 1.979728445526136, dt = 0.012140206332478932 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 1.9918686518586148, dt = 0.01718300300629441 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.0090516548649093, dt = 0.020544741067584877 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.029596395932494, dt = 0.022785927370986406 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01074141642686474
---------------- t = 2.0523823233034806, dt = 0.012140521033186344 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.064522844336667, dt = 0.01718331234262497 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.081706156679292, dt = 0.02054496440228715 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.102251121081579, dt = 0.022785889421026457 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01074213979176944
---------------- t = 2.1250370105026053, dt = 0.012140434208877416 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.1371774447114826, dt = 0.017183111486337624 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.1543605561978203, dt = 0.020544763823388067 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.1749053200212085, dt = 0.022785671284892537 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01067199003182933
---------------- t = 2.197690991306101, dt = 0.012140458363574826 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.2098314496696756, dt = 0.0171833352674863 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.227014784937162, dt = 0.020545116204467962 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.24755990114163, dt = 0.022786357772218678 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010633246115204144
---------------- t = 2.2703462589138486, dt = 0.012140773310025242 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.282487032223874, dt = 0.017183711464364246 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.299670743688238, dt = 0.020545554919175472 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.3202162986074133, dt = 0.022786985780561368 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010643969004529732
---------------- t = 2.3430032843879744, dt = 0.0121412641175124 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.355144548505487, dt = 0.01718462638458557 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.3723291748900723, dt = 0.020546844137396863 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.392876019027469, dt = 0.02278854280811747 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010627225170461543
---------------- t = 2.4156645618355865, dt = 0.012141430947824822 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.4278059927834112, dt = 0.017184257294622585 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.4449902500780336, dt = 0.020545566977628615 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.4655358170556623, dt = 0.022785680269965802 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01061608563217777
---------------- t = 2.4883214973256282, dt = 0.012139733553090145 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.5004612308787184, dt = 0.017181775877899135 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.5176430067566176, dt = 0.020542571643712877 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.5381855784003307, dt = 0.022782682828237262 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010616311108057821
---------------- t = 2.560968261228568, dt = 0.012138226518319235 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.5731064877468874, dt = 0.017179658144721232 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.590286145891609, dt = 0.020540043973902567 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.6108261898655116, dt = 0.02277979658543629 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010551559902820746
---------------- t = 2.633605986450948, dt = 0.012136721930447042 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.6457427083813947, dt = 0.01717761010773223 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.662920318489127, dt = 0.020537846299863497 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.6834581647889904, dt = 0.022777554309911314 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010474774268052992
---------------- t = 2.7062357190989017, dt = 0.012135602324201352 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.718371321423103, dt = 0.01717619695926395 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.735547518382367, dt = 0.02053620085719553 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.756083719239563, dt = 0.022776126863893884 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010575542680858491
---------------- t = 2.778859846103457, dt = 0.012135116090206535 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.7909949621936634, dt = 0.01717572678249497 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.8081706889761584, dt = 0.020536112242200793 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.828706801218359, dt = 0.022776504304731604 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010709071336104064
---------------- t = 2.851483305523091, dt = 0.01213578699655678 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.8636190925196474, dt = 0.01717685191160913 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.8807959444312567, dt = 0.020537816330570613 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.9013337607618275, dt = 0.0227790521683539 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010815897153191636
---------------- t = 2.9241128129301814, dt = 0.012137284027188963 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.9362500969573704, dt = 0.01717917728754395 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.953429274244914, dt = 0.02054069360416613 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 2.9739699678490803, dt = 0.022782165005090647 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010905977549613149
---------------- t = 2.996752132854171, dt = 0.012138803788495663 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.0088909366426666, dt = 0.017181184260630745 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.0260721209032972, dt = 0.020542851727405744 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.046614972630703, dt = 0.022783928661843672 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010894140421946516
---------------- t = 3.069398901292547, dt = 0.012139439224908358 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.0815383405174552, dt = 0.01718196013248529 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.0987203006499406, dt = 0.02054350269712267 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.119263803347063, dt = 0.02278454263058838 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010894612271109553
---------------- t = 3.142048345977652, dt = 0.012139699269845686 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.1541880452474973, dt = 0.017182252085169315 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.1713702973326665, dt = 0.02054345360499463 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.191913750937661, dt = 0.0227837991455493 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010851434799443312
---------------- t = 3.2146975500832102, dt = 0.012139144458574122 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.2268366945417846, dt = 0.017181213302734023 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.2440179078445186, dt = 0.02054237533589853 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.2645602831804172, dt = 0.022783055056684674 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010804523620575246
---------------- t = 3.287343338237102, dt = 0.012138877865921365 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.2994822161030233, dt = 0.017181042949129602 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.3166632590521528, dt = 0.02054238811281646 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.337205647164969, dt = 0.02278348682838712 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010763898371740768
---------------- t = 3.359989133993356, dt = 0.01213932350066549 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.372128457494022, dt = 0.01718184228971813 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.38931029978374, dt = 0.0205437738853608 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.4098540736691008, dt = 0.022785217078554198 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010778924312433895
---------------- t = 3.432639290747655, dt = 0.01214037516704 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.444779665914695, dt = 0.01718335547514034 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.461963021389835, dt = 0.020545493179455147 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.4825085145692904, dt = 0.022786021650719614 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010778612546969542
---------------- t = 3.50529453622001, dt = 0.012139708360007433 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.5174342445800173, dt = 0.0171816400966244 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.534615884676642, dt = 0.020542138923803143 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.5551580236004447, dt = 0.022781957364870847 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010815995129877956
---------------- t = 3.5779399809653154, dt = 0.012137658732914327 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.59007763969823, dt = 0.017178844053182617 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.6072564837514127, dt = 0.0205389200572602 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.627795403808673, dt = 0.022778252218076295 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01068759737603632
---------------- t = 3.6505736560267494, dt = 0.012135747938992211 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.6627094039657417, dt = 0.017176266927135933 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.6798856708928778, dt = 0.020536070820224 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.7004217417131016, dt = 0.02277562383228766 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010656390148623076
---------------- t = 3.7231973655453894, dt = 0.012134469485510545 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.7353318350309, dt = 0.01717447479801987 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.7525063098289198, dt = 0.020534230293473933 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.773040540122394, dt = 0.02277372509231848 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010817005664556227
---------------- t = 3.7958142652147124, dt = 0.012133802162122509 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.807948067376835, dt = 0.017173820825718927 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.825121888202554, dt = 0.020533671075944927 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.845655559278499, dt = 0.02277366171586928 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010971954730966073
---------------- t = 3.868429220994368, dt = 0.012134076717892794 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.8805632977122606, dt = 0.01717445700898573 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.897737754721246, dt = 0.02053498959583514 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.9182727443170813, dt = 0.022775577014467335 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01107038366704029
---------------- t = 3.941048321331549, dt = 0.012135474218885094 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.953183795550434, dt = 0.017176555187733956 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.970360350738168, dt = 0.02053771148756275 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 3.990898062225731, dt = 0.022778904872491873 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01113376068021861
---------------- t = 4.013676967098223, dt = 0.012137141094787687 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.0258141081930106, dt = 0.017178885720894616 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.042992993913905, dt = 0.02054031053543786 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.063533304449343, dt = 0.02278141760571148 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011205615849429553
---------------- t = 4.086314722055054, dt = 0.012138154617897374 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.098452876672951, dt = 0.017180113021213285 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.115632989694165, dt = 0.020541385860262027 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.136174375554426, dt = 0.022782229980149137 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01114378451315007
---------------- t = 4.158956605534575, dt = 0.012138500646551042 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.171095106181126, dt = 0.017180517414856867 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.188275623595983, dt = 0.020541888405594226 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.208817512001577, dt = 0.0227827037510785 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01112276939849549
---------------- t = 4.231600215752655, dt = 0.012138809729315674 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.243739025481971, dt = 0.017181102139898548 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.260920127621869, dt = 0.02054251043518322 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.281462638057052, dt = 0.022783229334345484 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01111006610596835
---------------- t = 4.304245867391398, dt = 0.01213893890929165 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.3163848063006895, dt = 0.01718107778277307 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.3335658840834625, dt = 0.02054247778755508 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.354108361871018, dt = 0.022783443030451444 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011118516982419387
---------------- t = 4.37689180490147, dt = 0.012139222425031343 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.389031027326501, dt = 0.017181723968169804 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.406212751294671, dt = 0.02054340458501601 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.426756155879687, dt = 0.022784900246852525 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011098778130617958
---------------- t = 4.44954105612654, dt = 0.01214009849389664 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.461681154620437, dt = 0.017182981861247743 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.478864136481684, dt = 0.020545017581531197 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.4994091540632155, dt = 0.02278654467488556 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011053723657372283
---------------- t = 4.522195698738101, dt = 0.012140384901996636 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.534336083640098, dt = 0.01718267159500131 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.551518755235099, dt = 0.02054343545585442 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.5720621906909535, dt = 0.022783200578992538 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.01100410711629281
---------------- t = 4.5948453912699465, dt = 0.01213823211011815 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.6069836233800645, dt = 0.01717975265001283 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.6241633760300775, dt = 0.02054031270963158 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.6447036887397095, dt = 0.022780020207983146 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010924100291445854
---------------- t = 4.667483708947692, dt = 0.01213649474324536 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.679620203690938, dt = 0.01717719648764403 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.696797400178582, dt = 0.020537181196582963 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.717334581375165, dt = 0.022776898841349945 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.010913544291358015
---------------- t = 4.740111480216515, dt = 0.012135113864984276 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.752246594081499, dt = 0.01717557593784818 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.7694221700193475, dt = 0.020535510524454532 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.789957680543802, dt = 0.022775051974173798 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011037137171844645
---------------- t = 4.812732732517976, dt = 0.01213435897586751 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.824867091493843, dt = 0.01717458076041436 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.842041672254258, dt = 0.02053472182363671 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.8625763940778945, dt = 0.022774827081882148 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011198513433843614
---------------- t = 4.885351221159777, dt = 0.01213442235312067 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.897485643512898, dt = 0.017174879385049878 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.914660522897948, dt = 0.020535379133033833 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.935195902030982, dt = 0.02277601921858994 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
Warning: the corrector tolerance are broken the step will be re rerunned [BasicGasSPH][rank=0]
eps_v = 0.011327481325217678
---------------- t = 4.957971921249571, dt = 0.012135431951602988 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.970107353201175, dt = 0.017176503953030928 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
---------------- t = 4.987283857154205, dt = 0.01271614284579492 ----------------
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
L2 error on B2 : 0.077212856332666
201 # The box spans 3 wavelengths along x1, and the solution only depends on
202 # x1 modulo the wavelength, so fold x1 into a single period to overlay them.
203 def fold_x1(x1):
204 return np.mod(x1 + 0.5 * wavelength, wavelength) - 0.5 * wavelength
205
206
207 if do_plot:
208 all_resols = sorted(set([resol] + list(plot_extra_resols)))
209 results = {res: (x1, B2) if res == resol else run_alfven_wave(res) for res in all_resols}
210
211 # Fig. 26: B2 vs x1 (folded into one wavelength) for each resolution
212 fig, ax = plt.subplots(figsize=(8, 5))
213
214 colors = plt.cm.viridis(np.linspace(0.0, 0.8, len(all_resols)))
215
216 for color, res in zip(colors, all_resols):
217 x1_r, B2_r = results[res]
218 ax.plot(fold_x1(x1_r), B2_r, ".", color=color, markersize=1.5, label=f"resol={res}")
219
220 x1_exact = np.linspace(-0.5 * wavelength, 0.5 * wavelength, 500)
221 B2_exact_line = amplitude * np.sin(2.0 * np.pi * x1_exact / wavelength)
222 ax.plot(x1_exact, B2_exact_line, "-", color="red", linewidth=1.5, label="exact solution")
223
224 ax.set_xlabel("$x_1$")
225 ax.set_ylabel("$B_2$")
226 ax.set_title(f"3D circularly polarised Alfven wave, $t={t_target:.0f}$ periods")
227 ax.legend()
228 fig.tight_layout()
229 fig.savefig(os.path.join(dump_folder, "mhd_alfven_wave_b2_vs_x1.png"), dpi=150)
230 plt.close(fig)
231
232 nxs = np.array(all_resols, dtype=float)
233 l1_errs = np.array(
234 [
235 np.mean(np.abs(B2_r - amplitude * np.sin(2.0 * np.pi * x1_r / wavelength)))
236 for x1_r, B2_r in (results[res] for res in all_resols)
237 ]
238 )
239 for res, err in zip(all_resols, l1_errs):
240 print(f"L1 error on B2 (resol={res}) : {err}")
241
242 fig, ax = plt.subplots(figsize=(6, 5))
243 ax.loglog(nxs, l1_errs, "o-", color="black", label="Shamrock")
244 if len(nxs) > 1:
245 # Local order between consecutive resolutions: a single fit over all points
246 # would be biased by coarse runs that are not yet in the asymptotic regime.
247 orders = -np.diff(np.log(l1_errs)) / np.diff(np.log(nxs))
248 for lo, hi, order in zip(all_resols[:-1], all_resols[1:], orders):
249 print(f"L1 convergence order between resol={lo} and {hi} : {order}")
250 ax.loglog(
251 nxs,
252 l1_errs[-1] * (nxs / nxs[-1]) ** -2,
253 "--",
254 color="gray",
255 label="second order",
256 )
257 ax.set_title(f"Convergence, order between the two finest = {orders[-1]:.2f}")
258 ax.set_xlabel("number of particles in x")
259 ax.set_ylabel("L1 error on $B_2$")
260 ax.legend()
261 fig.tight_layout()
262 fig.savefig(os.path.join(dump_folder, "mhd_alfven_wave_convergence.png"), dpi=150)
263 plt.close(fig)
264
265 test_pass = True
266 err_log = ""
267
268 expect_l2_err_B2 = 0.077212856332666
269 tol = 1e-5 # too generous for now
270
271
272 def float_equal(val1, val2, prec):
273 return abs(val1 - val2) < prec
274
275
276 if not float_equal(l2_err_B2, expect_l2_err_B2, tol * expect_l2_err_B2):
277 err_log += "error on the L2 norm of B2 is outside of tolerances:\n"
278 err_log += f" expected error = {expect_l2_err_B2} +- {tol * expect_l2_err_B2}\n"
279 err_log += f" obtained error = {l2_err_B2} (relative error = {(l2_err_B2 - expect_l2_err_B2) / expect_l2_err_B2})\n"
280 test_pass = False
281
282 if test_pass == False:
283 exit("Test did not pass L2 margins : \n" + err_log)
Total running time of the script: (0 minutes 33.987 seconds)
Estimated memory usage: 162 MB