Note
Go to the end to download the full example code.
Killing sphere for SPH simulation#
This simple example shows how to setup a killing sphere for sph simulations
9 import shamrock
10
11 # If we use the shamrock executable to run this script instead of the python interpreter,
12 # we should not initialize the system as the shamrock executable needs to handle specific MPI logic
13 if not shamrock.sys.is_initialized():
14 shamrock.change_loglevel(1)
15 shamrock.sys.init("0:0")
16
17
18 def is_in_sphere(pt):
19 x, y, z = pt
20 return (x**2 + y**2 + z**2) < 1
Use shamrock documentation style for matplotlib
26 shamrock.matplotlib.set_shamrock_mpl_style()
Setup parameters
32 dr = 0.1
33 pmass = 1
34
35 C_cour = 0.3
36 C_force = 0.25
37
38 bsize = 4
39
40
41 render_gif = True
42
43 dump_folder = "_to_trash"
44 sim_name = "kill_particle_sphere"
45
46 import os
47
48 # Create the dump directory if it does not exist
49 if shamrock.sys.world_rank() == 0:
50 os.makedirs(dump_folder, exist_ok=True)
Setup
55 ctx = shamrock.Context()
56 ctx.pdata_layout_new()
57
58 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
59
60 cfg = model.gen_default_config()
61 # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2)
62 # cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV)
63 # cfg.set_eos_locally_isothermalLP07(cs0=cs0, q=q, r0=r0)
64 cfg.set_artif_viscosity_VaryingCD10(
65 alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
66 )
67 cfg.set_boundary_periodic()
68 cfg.set_eos_adiabatic(1.00001)
The important part to enable killing
72 cfg.add_kill_sphere(center=(0.0, 0.0, 0.0), radius=4.0)
Rest of the setup
77 cfg.print_status()
78 model.set_solver_config(cfg)
79
80 model.init_scheduler(int(1e7), 1)
81
82 bmin = (-bsize, -bsize, -bsize)
83 bmax = (bsize, bsize, bsize)
84 model.resize_simulation_box(bmin, bmax)
85
86 model.set_particle_mass(pmass)
87
88 setup = model.get_setup()
89 lat = setup.make_generator_lattice_hcp(dr, (-bsize, -bsize, -bsize), (bsize, bsize, bsize))
90
91 thesphere = setup.make_modifier_filter(parent=lat, filter=is_in_sphere)
92
93 offset_sphere = setup.make_modifier_offset(
94 parent=thesphere, offset_position=(0.0, 0.0, 0.0), offset_velocity=(-1.0, -1.0, -1.0)
95 )
96
97 setup.apply_setup(offset_sphere)
98
99 model.set_value_in_a_box("uint", "f64", 1, bmin, bmax)
100
101 model.set_cfl_cour(C_cour)
102 model.set_cfl_force(C_force)
103
104 model.change_htolerance(1.3)
105 model.timestep()
106 model.change_htolerance(1.1)
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.00001
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [
{
"center": [
0.0,
0.0,
0.0
],
"radius": 4.0,
"type": "sphere"
}
],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"self_grav_config": {
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "none"
},
"show_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"tree_reduction_level": 3,
"type_id": "sycl::vec<f64,3>",
"unit_sys": null,
"use_two_stage_search": true
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 735 ( 7.4e+02 ) Ntotal = 735 ( 7.4e+02 rank min = 9.5e+03 max = 7.4e+02) rate = 7.350000e+02 N.s^-1
SPH setup: the generation step took : 0.07837843500000001 s
SPH setup: final particle count = 735 beginning injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.69 us (54.9%)
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.17 us (0.4%)
patch tree reduce : 1.36 us (0.2%)
gen split merge : 761.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 892.00 ns (0.2%)
LB compute : 543.39 us (97.7%)
LB move op cnt : 0
LB apply : 3.76 us (0.7%)
Info: patch count stable after 1 runs npatch = 1 [DataInserterUtility][rank=0]
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 735 / 735 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.0069654700000000005 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 2.8% 0.0% | 337.11 MB | 460.80 kB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.08678398500000001 s
Warning: .change_htolerance(val) is deprecated, [SPH][rank=0]
-> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
see: https://shamrock-code.github.io/Shamrock/mkdocs/models/sph/smoothing_length_tolerance
---------------- t = 0, dt = 0 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.23 us (0.7%)
patch tree reduce : 521.00 ns (0.1%)
gen split merge : 531.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 391.00 ns (0.1%)
LB compute : 427.04 us (97.7%)
LB move op cnt : 0
LB apply : 1.94 us (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.85 us (67.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.13 unconverged cnt = 735
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.169 unconverged cnt = 735
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.21970000000000003 unconverged cnt = 334
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.28561000000000003 unconverged cnt = 2
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735,-735,-735)
sum a = (-4.2690460541616737e-19,1.0503208545953324e-18,-1.1519648082658485e-18)
sum e = 1837.5
sum de = -8.131516293641283e-20
Info: cfl dt = 0.13054316431416452 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 8.6447e+03 | 735 | 1 | 8.502e-02 | 0.1% | 1.9% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Warning: .change_htolerance(val) is deprecated, [SPH][rank=0]
-> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
see: https://shamrock-code.github.io/Shamrock/mkdocs/models/sph/smoothing_length_tolerance
Show the current solvergraph
111 print(model.get_solver_dot_graph())
subgraph cluster_1999 {
n_1998 [label="MapEdge"];
e_6604 -> n_1998 [style="dashed", color=green];
e_6604 [label="dt",shape=rect, style=filled];
n_1998 -> e_6615 [style="dashed", color=red];
e_6615 [label="dt_half",shape=rect, style=filled];
n_1976 [label="SetEdge"];
n_1976 -> e_6616 [style="dashed", color=red];
e_6616 [label="m",shape=rect, style=filled];
subgraph cluster_1985 {
n_1977 [label="SetEdge"];
n_1977 -> e_6613 [style="dashed", color=red];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1978 [label="GetObjCntFromLayer"];
e_6613 -> n_1978 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1978 -> e_6614 [style="dashed", color=red];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1979 [label="GetFieldRefFromLayer"];
e_6613 -> n_1979 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1979 -> e_6617 [style="dashed", color=red];
e_6617 [label="xyz",shape=rect, style=filled];
n_1980 [label="GetFieldRefFromLayer"];
e_6613 -> n_1980 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1980 -> e_6618 [style="dashed", color=red];
e_6618 [label="vxyz",shape=rect, style=filled];
n_1981 [label="GetFieldRefFromLayer"];
e_6613 -> n_1981 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1981 -> e_6619 [style="dashed", color=red];
e_6619 [label="axyz",shape=rect, style=filled];
n_1982 [label="GetFieldRefFromLayer"];
e_6613 -> n_1982 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1982 -> e_6620 [style="dashed", color=red];
e_6620 [label="uint",shape=rect, style=filled];
n_1983 [label="GetFieldRefFromLayer"];
e_6613 -> n_1983 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1983 -> e_6621 [style="dashed", color=red];
e_6621 [label="duint",shape=rect, style=filled];
n_1984 [label="GetFieldRefFromLayer"];
e_6613 -> n_1984 [style="dashed", color=green];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1984 -> e_6622 [style="dashed", color=red];
e_6622 [label="hpart",shape=rect, style=filled];
n_1977 -> n_1978 [weight=3];
n_1978 -> n_1979 [weight=3];
n_1979 -> n_1980 [weight=3];
n_1980 -> n_1981 [weight=3];
n_1981 -> n_1982 [weight=3];
n_1982 -> n_1983 [weight=3];
n_1983 -> n_1984 [weight=3];
label = "attach fields";
}
subgraph cluster_1993 {
subgraph cluster_1988 {
n_1986 [label="ForwardEuler"];
e_6615 -> n_1986 [style="dashed", color=green];
e_6615 [label="dt_half",shape=rect, style=filled];
e_6619 -> n_1986 [style="dashed", color=green];
e_6619 [label="axyz",shape=rect, style=filled];
e_6614 -> n_1986 [style="dashed", color=green];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1986 -> e_6618 [style="dashed", color=red];
e_6618 [label="vxyz",shape=rect, style=filled];
n_1987 [label="ForwardEuler"];
e_6615 -> n_1987 [style="dashed", color=green];
e_6615 [label="dt_half",shape=rect, style=filled];
e_6621 -> n_1987 [style="dashed", color=green];
e_6621 [label="duint",shape=rect, style=filled];
e_6614 -> n_1987 [style="dashed", color=green];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1987 -> e_6620 [style="dashed", color=red];
e_6620 [label="uint",shape=rect, style=filled];
n_1986 -> n_1987 [weight=3];
label = "half step";
}
n_1992 [label="ForwardEuler"];
e_6604 -> n_1992 [style="dashed", color=green];
e_6604 [label="dt",shape=rect, style=filled];
e_6618 -> n_1992 [style="dashed", color=green];
e_6618 [label="vxyz",shape=rect, style=filled];
e_6614 -> n_1992 [style="dashed", color=green];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1992 -> e_6617 [style="dashed", color=red];
e_6617 [label="xyz",shape=rect, style=filled];
subgraph cluster_1991 {
n_1989 [label="ForwardEuler"];
e_6615 -> n_1989 [style="dashed", color=green];
e_6615 [label="dt_half",shape=rect, style=filled];
e_6619 -> n_1989 [style="dashed", color=green];
e_6619 [label="axyz",shape=rect, style=filled];
e_6614 -> n_1989 [style="dashed", color=green];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1989 -> e_6618 [style="dashed", color=red];
e_6618 [label="vxyz",shape=rect, style=filled];
n_1990 [label="ForwardEuler"];
e_6615 -> n_1990 [style="dashed", color=green];
e_6615 [label="dt_half",shape=rect, style=filled];
e_6621 -> n_1990 [style="dashed", color=green];
e_6621 [label="duint",shape=rect, style=filled];
e_6614 -> n_1990 [style="dashed", color=green];
e_6614 [label="Npart_patch",shape=rect, style=filled];
n_1990 -> e_6620 [style="dashed", color=red];
e_6620 [label="uint",shape=rect, style=filled];
n_1989 -> n_1990 [weight=3];
label = "half step";
}
n_1987 -> n_1992 [weight=3];
n_1992 -> n_1989 [weight=3];
label = "leapfrog predictor";
}
subgraph cluster_1997 {
n_1994 [label="FreeAlloc"];
n_1994 -> e_6623 [style="dashed", color=red];
e_6623 [label="part_to_remove",shape=rect, style=filled];
n_1995 [label="GetParticlesOutsideSphere"];
e_6617 -> n_1995 [style="dashed", color=green];
e_6617 [label="xyz",shape=rect, style=filled];
n_1995 -> e_6623 [style="dashed", color=red];
e_6623 [label="part_to_remove",shape=rect, style=filled];
n_1996 [label="KillParticles"];
e_6623 -> n_1996 [style="dashed", color=green];
e_6623 [label="part_to_remove",shape=rect, style=filled];
n_1996 -> e_6613 [style="dashed", color=red];
e_6613 [label="patchdatas",shape=rect, style=filled];
n_1994 -> n_1995 [weight=3];
n_1995 -> n_1996 [weight=3];
label = "part killing step";
}
n_1998 -> n_1976 [weight=3];
n_1976 -> n_1977 [weight=3];
n_1984 -> n_1986 [weight=3];
n_1990 -> n_1994 [weight=3];
label = "time step";
}
Draw utilities#
117 import matplotlib.pyplot as plt
118 import numpy as np
119
120
121 def plot_state(iplot):
122 pos = ctx.collect_data()["xyz"]
123
124 if shamrock.sys.world_rank() == 0:
125 X = pos[:, 0]
126 Y = pos[:, 1]
127 Z = pos[:, 2]
128
129 plt.cla()
130
131 ax.set_xlim3d(bmin[0], bmax[0])
132 ax.set_ylim3d(bmin[1], bmax[1])
133 ax.set_zlim3d(bmin[2], bmax[2])
134
135 ax.scatter(X, Y, Z, s=1)
136
137 ax.minorticks_off()
138
139 ax.set_title(f"t = {model.get_time():.2f} ")
140
141 plt.savefig(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.png"))
Run the simulation#
147 nstop = 28 # To be increased when empty simulations will be fixed
148 dt_stop = 0.1
149
150 t_stop = [i * dt_stop for i in range(nstop + 1)]
151
152 # Init MPL
153 fig = plt.figure(dpi=120)
154 ax = fig.add_subplot(111, projection="3d")
155
156 iplot = 0
157 istop = 0
158 for ttarg in t_stop:
159 model.evolve_until(ttarg)
160
161 # if do_plots:
162 plot_state(iplot)
163
164 iplot += 1
165 istop += 1
166
167 plt.close(fig)
Info: evolve_until (target_time = 0.00s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
Info: iteration since start : 1 [SPH][rank=0]
Info: time since start : 536.127937064 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.10s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0, dt = 0.1 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.53 us (2.4%)
patch tree reduce : 1.65 us (0.4%)
gen split merge : 1.07 us (0.2%)
split / merge op : 0/0
apply split merge : 1.13 us (0.3%)
LB compute : 409.09 us (93.8%)
LB move op cnt : 0
LB apply : 3.93 us (0.9%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.80 us (69.7%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999987,-735.0000000000001,-734.9999999999998)
sum a = (3.0560948736935156e-18,-6.057979638762756e-18,9.486769009248164e-20)
sum e = 1837.5000000099571
sum de = -5.692061405548898e-19
Info: cfl dt = 4.4384722909924195 cfl multiplier : 0.34 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.6761e+04 | 735 | 1 | 2.747e-02 | 0.2% | 1.6% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13107.412992082322 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 2 [SPH][rank=0]
Info: time since start : 536.341176259 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.20s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.1, dt = 0.1 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.63 us (2.9%)
patch tree reduce : 1.75 us (0.5%)
gen split merge : 952.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.25 us (0.3%)
LB compute : 347.15 us (93.3%)
LB move op cnt : 0
LB apply : 4.12 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.30 us (69.0%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999992,-734.9999999999995,-735.0000000000003)
sum a = (-2.4191260973582818e-18,2.7037291676357267e-18,1.3688052427629493e-18)
sum e = 1837.5000000099576
sum de = -4.851804721872632e-18
Info: cfl dt = 7.310448194870652 cfl multiplier : 0.56 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.2799e+04 | 735 | 1 | 2.241e-02 | 0.1% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16064.998985896938 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 3 [SPH][rank=0]
Info: time since start : 536.534436258 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.30s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.2, dt = 0.10000000000000003 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.42 us (2.2%)
patch tree reduce : 2.31 us (0.5%)
gen split merge : 912.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.14 us (0.2%)
LB compute : 444.18 us (94.4%)
LB move op cnt : 0
LB apply : 4.52 us (1.0%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.89 us (67.8%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.999999999999,-735.0000000000013,-735.0000000000005)
sum a = (-1.1221492485224971e-17,1.7618285302889447e-19,-9.486769009248164e-20)
sum e = 1837.500000009957
sum de = 2.0057740190981832e-18
Info: cfl dt = 9.225138276950846 cfl multiplier : 0.7066666666666667 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.6718e+04 | 735 | 1 | 2.751e-02 | 0.1% | 1.2% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13086.60478762531 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 4 [SPH][rank=0]
Info: time since start : 536.738030268 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.40s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.30000000000000004, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.12 us (3.2%)
patch tree reduce : 2.18 us (0.6%)
gen split merge : 932.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.21 us (0.3%)
LB compute : 320.43 us (92.5%)
LB move op cnt : 0
LB apply : 4.47 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.09 us (67.9%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000007,-735,-734.9999999999999)
sum a = (6.776263578034403e-21,-5.082197683525802e-19,-3.1170812458958252e-19)
sum e = 1837.5000000099562
sum de = -7.15573433840433e-18
Info: cfl dt = 10.501650417394146 cfl multiplier : 0.8044444444444444 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.9150e+04 | 735 | 1 | 2.521e-02 | 0.1% | 1.8% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14277.518310322333 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 5 [SPH][rank=0]
Info: time since start : 536.93831514 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.50s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.4, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.69 us (2.7%)
patch tree reduce : 2.52 us (0.6%)
gen split merge : 942.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.25 us (0.3%)
LB compute : 367.87 us (93.5%)
LB move op cnt : 0
LB apply : 3.96 us (1.0%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.67 us (67.5%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735,-734.9999999999993,-735)
sum a = (3.8624702394796095e-18,4.472333961502706e-19,1.260385025514399e-18)
sum e = 1837.5000000099562
sum de = 6.80336863234654e-18
Info: cfl dt = 11.352721173598157 cfl multiplier : 0.8696296296296296 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.9010e+04 | 735 | 1 | 2.534e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14209.136624835093 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 6 [SPH][rank=0]
Info: time since start : 537.1425231210001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.60s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.5, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.79 us (2.6%)
patch tree reduce : 1.68 us (0.4%)
gen split merge : 1.04 us (0.2%)
split / merge op : 0/0
apply split merge : 1.09 us (0.2%)
LB compute : 427.97 us (94.1%)
LB move op cnt : 0
LB apply : 4.17 us (0.9%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.63 us (64.9%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000003,-734.9999999999999,-735.0000000000003)
sum a = (-1.0611628763201875e-17,3.9979955110402976e-18,-1.4501204056993622e-18)
sum e = 1837.5000000099565
sum de = 9.703609443745265e-18
Info: cfl dt = 11.920173864876189 cfl multiplier : 0.9130864197530864 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.1960e+04 | 735 | 1 | 2.300e-02 | 0.1% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 15653.813479855358 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 7 [SPH][rank=0]
Info: time since start : 537.3423696780001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.70s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.6000000000000001, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.84 us (3.4%)
patch tree reduce : 2.20 us (0.6%)
gen split merge : 993.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.45 us (0.4%)
LB compute : 322.90 us (92.0%)
LB move op cnt : 0
LB apply : 4.70 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.16 us (66.5%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000002,-734.9999999999997,-735.0000000000002)
sum a = (5.265156800132731e-18,3.1374100366299285e-18,-7.860465750519907e-19)
sum e = 1837.5000000099549
sum de = -1.235990476633475e-17
Info: cfl dt = 12.298556869664226 cfl multiplier : 0.9420576131687243 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.8110e+04 | 735 | 1 | 2.615e-02 | 0.1% | 1.3% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13768.103526347864 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 8 [SPH][rank=0]
Info: time since start : 537.546209354 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.80s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.7000000000000001, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.44 us (3.1%)
patch tree reduce : 1.78 us (0.5%)
gen split merge : 922.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.38 us (0.4%)
LB compute : 314.20 us (92.6%)
LB move op cnt : 0
LB apply : 3.89 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (64.3%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999997,-734.9999999999994,-734.9999999999999)
sum a = (3.3542504711270293e-18,-1.6805133673525319e-18,-1.260385025514399e-18)
sum e = 1837.500000009954
sum de = -5.827586677109586e-18
Info: cfl dt = 12.55090223978563 cfl multiplier : 0.9613717421124829 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.3016e+04 | 735 | 1 | 2.226e-02 | 0.1% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16171.234688592102 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 9 [SPH][rank=0]
Info: time since start : 537.7454711830001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 0.90s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.8, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.72 us (3.0%)
patch tree reduce : 1.96 us (0.6%)
gen split merge : 882.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.19 us (0.3%)
LB compute : 329.45 us (92.8%)
LB move op cnt : 0
LB apply : 4.35 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.13 us (65.9%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000007,-735.0000000000002,-734.9999999999997)
sum a = (-2.7037291676357267e-18,1.7211709488207383e-18,6.776263578034403e-19)
sum e = 1837.5000000099546
sum de = -9.839134715305953e-18
Info: cfl dt = 12.719231298543182 cfl multiplier : 0.9742478280749886 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.9951e+04 | 735 | 1 | 2.454e-02 | 0.1% | 1.2% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14669.788560670024 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 10 [SPH][rank=0]
Info: time since start : 537.943456156 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.00s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 0.9, dt = 0.09999999999999998 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.26 us (3.2%)
patch tree reduce : 2.32 us (0.7%)
gen split merge : 992.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.10 us (0.3%)
LB compute : 323.31 us (92.6%)
LB move op cnt : 0
LB apply : 4.02 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.31 us (67.1%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000008,-735.0000000000001,-735.0000000000001)
sum a = (1.2807138162485021e-18,1.917682592583736e-18,-9.893344823930228e-19)
sum e = 1837.5000000099558
sum de = 6.179952383167375e-18
Info: cfl dt = 12.831558291123326 cfl multiplier : 0.9828318853833258 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.3151e+04 | 735 | 1 | 2.217e-02 | 0.1% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16237.014954696693 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 11 [SPH][rank=0]
Info: time since start : 538.140247167 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.10s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.62 us (3.1%)
patch tree reduce : 2.15 us (0.6%)
gen split merge : 1.19 us (0.3%)
split / merge op : 0/0
apply split merge : 1.21 us (0.4%)
LB compute : 315.61 us (92.5%)
LB move op cnt : 0
LB apply : 4.29 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.12 us (64.1%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999998,-735.0000000000013,-734.9999999999998)
sum a = (4.851804721872632e-18,4.004771774618332e-18,3.279711571768651e-18)
sum e = 1837.5000000099549
sum de = 3.333921680392926e-18
Info: cfl dt = 12.906559440326545 cfl multiplier : 0.9885545902555505 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.3185e+04 | 735 | 1 | 2.215e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16253.796898242803 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 12 [SPH][rank=0]
Info: time since start : 538.336128435 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.20s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.1, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.85 us (2.9%)
patch tree reduce : 2.14 us (0.6%)
gen split merge : 25.36 us (6.8%)
split / merge op : 0/0
apply split merge : 1.14 us (0.3%)
LB compute : 323.21 us (86.3%)
LB move op cnt : 0
LB apply : 4.85 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.50 us (68.8%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999995,-735.0000000000002,-735.0000000000001)
sum a = (-2.4123498337802474e-18,-7.386127300057499e-19,-4.2012834183813297e-19)
sum e = 1837.5000000099526
sum de = -4.689174395999807e-18
Info: cfl dt = 12.956685627293247 cfl multiplier : 0.9923697268370336 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.7932e+04 | 735 | 1 | 2.631e-02 | 0.1% | 1.1% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13681.096676709616 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 13 [SPH][rank=0]
Info: time since start : 538.5359173950001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.30s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.2000000000000002, dt = 0.09999999999999987 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.43 us (3.1%)
patch tree reduce : 1.75 us (0.5%)
gen split merge : 1.26 us (0.4%)
split / merge op : 0/0
apply split merge : 1.10 us (0.3%)
LB compute : 314.04 us (92.7%)
LB move op cnt : 0
LB apply : 3.97 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.53 us (65.7%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000005,-735.0000000000006,-735.0000000000005)
sum a = (-3.9302328752599536e-18,6.369687763352339e-19,-1.6805133673525319e-18)
sum e = 1837.500000009954
sum de = -1.9786689647860456e-18
Info: cfl dt = 12.990237500990501 cfl multiplier : 0.9949131512246892 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.1738e+04 | 735 | 1 | 2.316e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 15545.068391823792 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 14 [SPH][rank=0]
Info: time since start : 538.731936169 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.40s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.3, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.89 us (3.2%)
patch tree reduce : 1.83 us (0.5%)
gen split merge : 972.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.16 us (0.3%)
LB compute : 314.98 us (92.7%)
LB move op cnt : 0
LB apply : 3.99 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.06 us (65.4%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000003,-735.0000000000008,-734.9999999999995)
sum a = (4.235164736271502e-18,9.554531645028508e-19,-5.014435047745458e-19)
sum e = 1837.5000000099512
sum de = 8.673617379884035e-19
Info: cfl dt = 13.012748880137558 cfl multiplier : 0.9966087674831261 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.2195e+04 | 735 | 1 | 2.283e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 15768.947770529234 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 15 [SPH][rank=0]
Info: time since start : 538.9286068260001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.50s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.4000000000000001, dt = 0.09999999999999987 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.30 us (3.0%)
patch tree reduce : 1.93 us (0.6%)
gen split merge : 1.00 us (0.3%)
split / merge op : 0/0
apply split merge : 1.09 us (0.3%)
LB compute : 323.42 us (92.8%)
LB move op cnt : 0
LB apply : 3.92 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.27 us (69.0%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999993,-735.0000000000002,-735.0000000000001)
sum a = (-4.655293078109635e-18,-2.2768245622195593e-18,-8.402566836762659e-19)
sum e = 1837.5000000099517
sum de = 2.168404344971009e-19
Info: cfl dt = 13.027909020178885 cfl multiplier : 0.997739178322084 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.2833e+04 | 735 | 1 | 2.239e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16081.240281737946 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 16 [SPH][rank=0]
Info: time since start : 539.126146561 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.60s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.5, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.10 us (2.9%)
patch tree reduce : 1.60 us (0.5%)
gen split merge : 1.03 us (0.3%)
split / merge op : 0/0
apply split merge : 1.12 us (0.3%)
LB compute : 325.06 us (92.9%)
LB move op cnt : 0
LB apply : 4.39 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.14 us (65.8%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999993,-735.0000000000014,-735.0000000000001)
sum a = (1.2739375526704677e-18,-6.030874584450618e-19,2.303929616531697e-19)
sum e = 1837.5000000099483
sum de = 1.5178830414797062e-18
Info: cfl dt = 13.0381774585845 cfl multiplier : 0.998492785548056 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.6999e+04 | 735 | 1 | 2.722e-02 | 0.1% | 1.2% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13224.174142978018 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 17 [SPH][rank=0]
Info: time since start : 539.326905555 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.70s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.6, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 33.86 us (9.2%)
patch tree reduce : 1.99 us (0.5%)
gen split merge : 932.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.53 us (0.4%)
LB compute : 319.16 us (86.5%)
LB move op cnt : 0
LB apply : 4.35 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.15 us (66.0%)
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999997,-735,-734.9999999999999)
sum a = (2.2632720350634905e-18,1.6872896309305663e-18,-7.724940478959219e-19)
sum e = 1837.5000000099478
sum de = 1.3552527156068805e-18
Info: cfl dt = 13.04519391302014 cfl multiplier : 0.9989951903653708 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 2.7401e+04 | 735 | 1 | 2.682e-02 | 0.1% | 1.1% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13420.742129759949 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 18 [SPH][rank=0]
Info: time since start : 539.529463767 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.80s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.7000000000000002, dt = 0.09999999999999987 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 735.0 min = 735.0 factor = 1
- strategy "round robin" : max = 698.2 min = 698.2 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 735
max = 735
avg = 735
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.57 us (3.1%)
patch tree reduce : 2.00 us (0.6%)
gen split merge : 982.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.22 us (0.4%)
LB compute : 313.33 us (92.5%)
LB move op cnt : 0
LB apply : 4.14 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.60 us (65.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.29624455421092516 unconverged cnt = 4
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-723.9990621681544,-723.9989589588944,-723.999080836305)
sum a = (-4.221612209115433e-18,-2.7376104855258987e-18,1.938011383317839e-18)
sum e = 1809.9971019165425
sum de = 9.107298248878237e-18
Info: cfl dt = 12.89661811090233 cfl multiplier : 0.9993301269102473 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 1.3934e+04 | 724 | 1 | 5.196e-02 | 0.0% | 1.3% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6928.733701164399 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 19 [SPH][rank=0]
Info: time since start : 539.756090913 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 1.90s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.8, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 724.0 min = 724.0 factor = 1
- strategy "round robin" : max = 687.8 min = 687.8 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 724
max = 724
avg = 724
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.11 us (2.9%)
patch tree reduce : 1.87 us (0.5%)
gen split merge : 1.24 us (0.3%)
split / merge op : 0/0
apply split merge : 1.14 us (0.3%)
LB compute : 353.09 us (93.2%)
LB move op cnt : 0
LB apply : 4.40 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.60 us (66.4%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3005616719262014 unconverged cnt = 19
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-681.9968714883288,-681.9966204670972,-681.9967286191298)
sum a = (-3.23227772672241e-18,1.260385025514399e-18,1.1248597539537109e-18)
sum e = 1704.9902204838565
sum de = -9.812029660993815e-18
Info: cfl dt = 12.78961555022715 cfl multiplier : 0.9995534179401648 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 1.3152e+04 | 682 | 1 | 5.185e-02 | 0.1% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6942.486111604751 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 20 [SPH][rank=0]
Info: time since start : 539.982904145 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.00s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 1.9000000000000001, dt = 0.09999999999999987 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 682.0 min = 682.0 factor = 1
- strategy "round robin" : max = 647.9 min = 647.9 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 682
max = 682
avg = 682
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.04 us (3.2%)
patch tree reduce : 1.95 us (0.6%)
gen split merge : 922.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.40 us (0.4%)
LB compute : 317.37 us (92.5%)
LB move op cnt : 0
LB apply : 4.48 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (67.1%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.2962597061448628 unconverged cnt = 34
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.29755427911452514 unconverged cnt = 2
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-612.9949011129379,-612.9946819821328,-612.9947797392336)
sum a = (-1.043544591017298e-18,-3.3406979439709605e-18,1.8973538018496328e-19)
sum e = 1532.484362757863
sum de = -1.0842021724855044e-19
Info: cfl dt = 12.861008218464523 cfl multiplier : 0.9997022786267765 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 9.1730e+03 | 613 | 1 | 6.683e-02 | 0.0% | 1.6% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5387.094864091893 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 21 [SPH][rank=0]
Info: time since start : 540.220826588 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.10s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 613.0 min = 613.0 factor = 1
- strategy "round robin" : max = 582.4 min = 582.4 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 613
max = 613
avg = 613
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.10 us (3.0%)
patch tree reduce : 2.09 us (0.6%)
gen split merge : 942.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.07 us (0.3%)
LB compute : 309.73 us (92.5%)
LB move op cnt : 0
LB apply : 4.33 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.42 us (68.3%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3067183489629737 unconverged cnt = 47
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.307306761640605 unconverged cnt = 2
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-532.9935533115207,-532.9935945920531,-532.9934668702186)
sum a = (-3.3881317890172014e-19,-1.8973538018496328e-19,-3.7947076036992655e-19)
sum e = 1332.4806147210948
sum de = -2.6291902682773483e-18
Info: cfl dt = 12.68539798388675 cfl multiplier : 0.9998015190845176 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 8.0873e+03 | 533 | 1 | 6.591e-02 | 0.0% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5462.319410978682 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 22 [SPH][rank=0]
Info: time since start : 540.452853201 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.20s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.1, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 533.0 min = 533.0 factor = 1
- strategy "round robin" : max = 506.3 min = 506.3 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 533
max = 533
avg = 533
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.26 us (2.8%)
patch tree reduce : 1.90 us (0.5%)
gen split merge : 741.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.20 us (0.3%)
LB compute : 343.62 us (93.3%)
LB move op cnt : 0
LB apply : 4.20 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.25 us (68.6%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.2962764478718124 unconverged cnt = 55
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.31069429926209985 unconverged cnt = 2
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-433.9927338736061,-433.99250861466425,-433.99249597858295)
sum a = (2.4868887331386258e-18,-2.7308342219478643e-18,-7.453889935837843e-19)
sum e = 1084.9777384174886
sum de = 3.06287113727155e-18
Info: cfl dt = 12.731993623699791 cfl multiplier : 0.9998676793896785 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 7.2844e+03 | 434 | 1 | 5.958e-02 | 0.0% | 1.3% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6042.397894721158 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 23 [SPH][rank=0]
Info: time since start : 540.670514655 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.30s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.2, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 434.0 min = 434.0 factor = 1
- strategy "round robin" : max = 412.3 min = 412.3 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 434
max = 434
avg = 434
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 11.76 us (3.5%)
patch tree reduce : 1.70 us (0.5%)
gen split merge : 1.13 us (0.3%)
split / merge op : 0/0
apply split merge : 1.36 us (0.4%)
LB compute : 309.91 us (91.9%)
LB move op cnt : 0
LB apply : 4.52 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.03 us (64.0%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.310950072983323 unconverged cnt = 57
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.31560777249532096 unconverged cnt = 4
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-341.99193142018476,-341.9921450989848,-341.9925044093816)
sum a = (2.3784685158900754e-18,-1.0164395367051604e-19,-5.963111948670274e-19)
sum e = 854.9765809341825
sum de = 5.5294310796760726e-18
Info: cfl dt = 12.667051368296768 cfl multiplier : 0.9999117862597856 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 5.5526e+03 | 342 | 1 | 6.159e-02 | 0.0% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5844.805414575788 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 24 [SPH][rank=0]
Info: time since start : 540.882422434 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.40s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.3000000000000003, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 342.0 min = 342.0 factor = 1
- strategy "round robin" : max = 324.9 min = 324.9 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 342
max = 342
avg = 342
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.21 us (2.6%)
patch tree reduce : 1.69 us (0.4%)
gen split merge : 851.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.23 us (0.3%)
LB compute : 374.25 us (93.9%)
LB move op cnt : 0
LB apply : 3.98 us (1.0%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.13 us (68.3%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.30201844852987103 unconverged cnt = 56
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.32857648069500406 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3614341287645045 unconverged cnt = 1
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-245.9921769555506,-245.99263707959494,-245.99229063406636)
sum a = (-1.0638733817514012e-18,-1.6398557858843255e-18,-2.710505431213761e-20)
sum e = 614.9771046889905
sum de = 1.4365678785432934e-18
Info: cfl dt = 12.50243532668346 cfl multiplier : 0.9999411908398571 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.3027e+03 | 246 | 1 | 7.448e-02 | 0.0% | 1.6% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 4833.230002903567 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 25 [SPH][rank=0]
Info: time since start : 541.100864191 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.50s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.4000000000000004, dt = 0.09999999999999964 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 246.0 min = 246.0 factor = 1
- strategy "round robin" : max = 233.7 min = 233.7 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 246
max = 246
avg = 246
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.93 us (3.4%)
patch tree reduce : 1.67 us (0.5%)
gen split merge : 841.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 1.04 us (0.3%)
LB compute : 300.84 us (92.3%)
LB move op cnt : 0
LB apply : 3.92 us (1.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.47 us (66.4%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3277421878427072 unconverged cnt = 53
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.34357314083155893 unconverged cnt = 10
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-166.99278662167782,-166.99320940422416,-166.99324858622313)
sum a = (-1.4840017235895342e-18,-7.182839392716467e-19,5.55653613398821e-19)
sum e = 417.47924469434525
sum de = 2.3852447794681098e-18
Info: cfl dt = 12.481349233596452 cfl multiplier : 0.999960793893238 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 3.3413e+03 | 167 | 1 | 4.998e-02 | 0.0% | 1.5% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 7202.721476282572 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 26 [SPH][rank=0]
Info: time since start : 541.2880832120001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.60s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.5, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 167.0 min = 167.0 factor = 1
- strategy "round robin" : max = 158.7 min = 158.7 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 167
max = 167
avg = 167
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.43 us (2.8%)
patch tree reduce : 2.03 us (0.5%)
gen split merge : 861.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1.16 us (0.3%)
LB compute : 348.97 us (93.2%)
LB move op cnt : 0
LB apply : 4.23 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.67 us (71.9%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.30588112681599167 unconverged cnt = 40
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.33646923949759083 unconverged cnt = 17
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.36630356642675005 unconverged cnt = 5
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-91.99460804607467,-91.99491947485237,-91.99457795712206)
sum a = (-4.472333961502706e-19,-2.642742795433417e-19,-8.131516293641283e-20)
sum e = 229.98410556558173
sum de = 2.439454888092385e-19
Info: cfl dt = 12.185082235929789 cfl multiplier : 0.999973862595492 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 1.3370e+03 | 92 | 1 | 6.881e-02 | 0.0% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5231.58657003481 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 27 [SPH][rank=0]
Info: time since start : 541.492475938 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.70s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.6, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 92.0 min = 92.0 factor = 1
- strategy "round robin" : max = 87.4 min = 87.4 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 92
max = 92
avg = 92
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.26 us (3.2%)
patch tree reduce : 2.09 us (0.6%)
gen split merge : 1.22 us (0.4%)
split / merge op : 0/0
apply split merge : 1.14 us (0.4%)
LB compute : 298.47 us (92.3%)
LB move op cnt : 0
LB apply : 4.18 us (1.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.10 us (65.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3332273933879859 unconverged cnt = 28
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.36655013272678455 unconverged cnt = 18
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.40320514599946305 unconverged cnt = 9
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-43.996389524961266,-43.996148049667234,-43.996195016207054)
sum a = (-9.486769009248164e-20,1.3552527156068805e-19,1.3552527156068805e-20)
sum e = 109.9887327309434
sum de = 2.981555974335137e-19
Info: cfl dt = 11.902519640066302 cfl multiplier : 0.9999825750636614 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 7.0891e+02 | 44 | 1 | 6.207e-02 | 0.0% | 1.6% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5800.1868497970445 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 28 [SPH][rank=0]
Info: time since start : 541.682884351 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Info: evolve_until (target_time = 2.80s, niter_max = -1, max_walltime = -1.00s) [SPH][rank=0]
---------------- t = 2.7, dt = 0.10000000000000009 ----------------
Info: Summary (strategy = round robin): [LoadBalance][rank=0]
- strategy "psweep" : max = 44.0 min = 44.0 factor = 1
- strategy "round robin" : max = 41.8 min = 41.8 factor = 0.95
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 44
max = 44
avg = 44
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.22 us (2.9%)
patch tree reduce : 1.68 us (0.5%)
gen split merge : 1.05 us (0.3%)
split / merge op : 0/0
apply split merge : 1.15 us (0.3%)
LB compute : 325.79 us (93.2%)
LB move op cnt : 0
LB apply : 3.94 us (1.1%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1.99 us (67.9%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.31758086030088317 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3493389463309715 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3842728409640687 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.4227001250604756 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.4649701375665232 unconverged cnt = 8
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.5114671513231756 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.5610681639786609 unconverged cnt = 3
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.609277261524694 unconverged cnt = 3
Warning: the unit system is not set [sph::Config][rank=0]
Info: conservation infos : [sph::Model][rank=0]
sum v = (-10.998438229211487,-10.998627137787167,-10.998397802857001)
sum a = (6.098637220230962e-20,-3.3881317890172014e-21,0)
sum e = 27.495463315833934
sum de = -7.453889935837843e-20
Info: cfl dt = 15.643421624235623 cfl multiplier : 0.9999883833757742 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+======+========+===========+======+=============+=============+=============+
| 0 | 8.6830e+01 | 11 | 1 | 1.267e-01 | 0.0% | 1.4% 0.0% | 337.11 MB | 460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 2841.7108785570726 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 29 [SPH][rank=0]
Info: time since start : 541.935583909 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Convert PNG sequence to Image sequence in mpl#
173 from matplotlib import animation
174 from shamrock.utils.plot import show_image_sequence
175
176 # If the animation is not returned only a static image will be shown in the doc
177 glob_str = os.path.join(dump_folder, f"{sim_name}_*.png")
178 ani = show_image_sequence(glob_str, render_gif=render_gif)
179
180 if render_gif and shamrock.sys.world_rank() == 0:
181 # To save the animation using Pillow as a gif
182 # writer = animation.PillowWriter(fps=15,
183 # metadata=dict(artist='Me'),
184 # bitrate=1800)
185 # ani.save('scatter.gif', writer=writer)
186
187 # Show the animation
188 plt.show()
Total running time of the script: (0 minutes 12.912 seconds)
Estimated memory usage: 160 MB