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
Setup parameters
27 dr = 0.1
28 pmass = 1
29
30 C_cour = 0.3
31 C_force = 0.25
32
33 bsize = 4
34
35
36 render_gif = True
37
38 dump_folder = "_to_trash"
39 sim_name = "kill_particle_sphere"
40
41 import os
42
43 # Create the dump directory if it does not exist
44 if shamrock.sys.world_rank() == 0:
45 os.makedirs(dump_folder, exist_ok=True)
Setup
50 ctx = shamrock.Context()
51 ctx.pdata_layout_new()
52
53 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
54
55 cfg = model.gen_default_config()
56 # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2)
57 # cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV)
58 # cfg.set_eos_locally_isothermalLP07(cs0=cs0, q=q, r0=r0)
59 cfg.set_artif_viscosity_VaryingCD10(
60 alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
61 )
62 cfg.set_boundary_periodic()
63 cfg.set_eos_adiabatic(1.00001)
The important part to enable killing
67 cfg.add_kill_sphere(center=(0.0, 0.0, 0.0), radius=4.0)
Rest of the setup
72 cfg.print_status()
73 model.set_solver_config(cfg)
74
75 model.init_scheduler(int(1e7), 1)
76
77 bmin = (-bsize, -bsize, -bsize)
78 bmax = (bsize, bsize, bsize)
79 model.resize_simulation_box(bmin, bmax)
80
81 model.set_particle_mass(pmass)
82
83 setup = model.get_setup()
84 lat = setup.make_generator_lattice_hcp(dr, (-bsize, -bsize, -bsize), (bsize, bsize, bsize))
85
86 thesphere = setup.make_modifier_filter(parent=lat, filter=is_in_sphere)
87
88 offset_sphere = setup.make_modifier_offset(
89 parent=thesphere, offset_position=(0.0, 0.0, 0.0), offset_velocity=(-1.0, -1.0, -1.0)
90 )
91
92 setup.apply_setup(offset_sphere)
93
94 model.set_value_in_a_box("uint", "f64", 1, bmin, bmax)
95
96 model.set_cfl_cour(C_cour)
97 model.set_cfl_force(C_force)
98
99 model.change_htolerance(1.3)
100 model.timestep()
101 model.change_htolerance(1.1)
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"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,
"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,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"time_state": {
"cfl_multiplier": 0.01,
"dt_sph": 0.0,
"time": 0.0
},
"tree_reduction_level": 3,
"type_id": "sycl::vec<f64,3>",
"unit_sys": null,
"use_two_stage_search": true
}
]
------------------------------------
Info: pushing data in scheduler, N = 735 [DataInserterUtility][rank=0]
Info: reattributing data ... [DataInserterUtility][rank=0]
Info: reattributing data done in 1746.19 us [DataInserterUtility][rank=0]
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.04 us (50.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1363.00 ns (0.3%)
patch tree reduce : 732.00 ns (0.2%)
gen split merge : 641.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 892.00 ns (0.2%)
LB compute : 414.96 us (97.5%)
LB move op cnt : 0
LB apply : 3.33 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.55 us (64.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1623.00 ns (0.4%)
patch tree reduce : 331.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.1%)
LB compute : 375.60 us (98.0%)
LB move op cnt : 0
LB apply : 1834.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.50 us (66.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1242.00 ns (0.3%)
patch tree reduce : 331.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 352.98 us (98.0%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Final load balancing step 0 of 3 [SPH setup][rank=0]
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 : 2.50 us (60.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1372.00 ns (0.4%)
patch tree reduce : 331.00 ns (0.1%)
gen split merge : 360.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.1%)
LB compute : 356.18 us (98.0%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.40 us (62.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1313.00 ns (0.4%)
patch tree reduce : 311.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 335.23 us (97.9%)
LB move op cnt : 0
LB apply : 1873.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.42 us (65.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1222.00 ns (0.3%)
patch tree reduce : 331.00 ns (0.1%)
gen split merge : 361.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 385.52 us (98.2%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Final load balancing step 1 of 3 [SPH setup][rank=0]
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 : 2.27 us (64.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1212.00 ns (0.3%)
patch tree reduce : 311.00 ns (0.1%)
gen split merge : 361.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.1%)
LB compute : 353.93 us (98.0%)
LB move op cnt : 0
LB apply : 1753.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.25 us (66.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1282.00 ns (0.4%)
patch tree reduce : 320.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 310.00 ns (0.1%)
LB compute : 319.84 us (97.8%)
LB move op cnt : 0
LB apply : 1773.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.46 us (66.3%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1202.00 ns (0.3%)
patch tree reduce : 321.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.1%)
LB compute : 364.21 us (98.1%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Final load balancing step 2 of 3 [SPH setup][rank=0]
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 : 2.28 us (63.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1543.00 ns (0.4%)
patch tree reduce : 331.00 ns (0.1%)
gen split merge : 360.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.1%)
LB compute : 412.48 us (98.2%)
LB move op cnt : 0
LB apply : 1894.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.24 us (63.8%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1263.00 ns (0.3%)
patch tree reduce : 311.00 ns (0.1%)
gen split merge : 351.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.1%)
LB compute : 399.26 us (98.3%)
LB move op cnt : 0
LB apply : 1814.00 ns (0.4%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.50 us (66.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 1784.00 ns (0.5%)
patch tree reduce : 360.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 335.43 us (97.7%)
LB move op cnt : 0
LB apply : 1893.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: the setup took : 0.08388714200000001 s [SPH setup][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
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.92 us (0.8%)
patch tree reduce : 381.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 371.00 ns (0.1%)
LB compute : 346.01 us (97.4%)
LB move op cnt : 0
LB apply : 1813.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.32 us (70.1%)
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
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735,-735,-735)
sum a = (-5.590417451878382e-18,1.904130065427667e-18,2.2090619264392153e-18)
sum e = 1837.5
sum de = -9.75781955236954e-19
Info: CFL hydro = 0.1305431643141645 sink sink = inf [SPH][rank=0]
Info: cfl dt = 0.1305431643141645 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.0642e+04 | 735 | 6.907e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
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
Draw utilities#
107 import matplotlib.pyplot as plt
108 import numpy as np
109
110
111 def plot_state(iplot):
112
113 pos = ctx.collect_data()["xyz"]
114
115 if shamrock.sys.world_rank() == 0:
116 X = pos[:, 0]
117 Y = pos[:, 1]
118 Z = pos[:, 2]
119
120 plt.cla()
121
122 ax.set_xlim3d(bmin[0], bmax[0])
123 ax.set_ylim3d(bmin[1], bmax[1])
124 ax.set_zlim3d(bmin[2], bmax[2])
125
126 ax.scatter(X, Y, Z, s=1)
127
128 ax.set_title(f"t = {model.get_time():.2f} ")
129
130 plt.savefig(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.png"))
Run the simulation#
136 nstop = 28 # To be increased when epmty simulations will be fixed
137 dt_stop = 0.1
138
139 t_stop = [i * dt_stop for i in range(nstop + 1)]
140
141 # Init MPL
142 fig = plt.figure(dpi=120)
143 ax = fig.add_subplot(111, projection="3d")
144
145 iplot = 0
146 istop = 0
147 for ttarg in t_stop:
148
149 model.evolve_until(ttarg)
150
151 # if do_plots:
152 plot_state(iplot)
153
154 iplot += 1
155 istop += 1
156
157 plt.close(fig)
Info: iteration since start : 1 [SPH][rank=0]
Info: time since start : 200.923104614 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0, dt = 0.1 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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 : 7.55 us (1.6%)
patch tree reduce : 822.00 ns (0.2%)
gen split merge : 651.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1032.00 ns (0.2%)
LB compute : 440.01 us (95.7%)
LB move op cnt : 0
LB apply : 3.19 us (0.7%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.62 us (68.7%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999987,-735.0000000000001,-734.9999999999998)
sum a = (3.211948935988307e-18,-2.6495190590114515e-18,-8.131516293641283e-19)
sum e = 1837.5000000099578
sum de = -6.478107980600889e-18
Info: CFL hydro = 4.4384722909924195 sink sink = inf [SPH][rank=0]
Info: cfl dt = 4.4384722909924195 cfl multiplier : 0.34 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 3.9339e+04 | 735 | 1.868e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 19268.223899544897 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 2 [SPH][rank=0]
Info: time since start : 201.134162053 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.1, dt = 0.1 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.40 us (0.9%)
patch tree reduce : 451.00 ns (0.1%)
gen split merge : 390.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 451.00 ns (0.1%)
LB compute : 358.57 us (97.1%)
LB move op cnt : 0
LB apply : 2.22 us (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (68.3%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999992,-734.9999999999995,-735.0000000000003)
sum a = (3.2255014631443757e-18,-1.5585406229479126e-19,1.6398557858843255e-18)
sum e = 1837.500000009958
sum de = -1.3552527156068805e-19
Info: CFL hydro = 7.310448194870652 sink sink = inf [SPH][rank=0]
Info: cfl dt = 7.310448194870652 cfl multiplier : 0.56 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 3.9148e+04 | 735 | 1.877e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 19174.515790213754 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 3 [SPH][rank=0]
Info: time since start : 201.32926797300001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.2, dt = 0.10000000000000003 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.17 us (0.9%)
patch tree reduce : 350.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 355.18 us (97.2%)
LB move op cnt : 0
LB apply : 1934.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.46 us (69.7%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.999999999999,-735.0000000000013,-735.0000000000005)
sum a = (-8.321251673826247e-18,-3.889575293791747e-18,-1.802486111757151e-18)
sum e = 1837.5000000099571
sum de = 2.846030702774449e-18
Info: CFL hydro = 9.225138276950844 sink sink = inf [SPH][rank=0]
Info: cfl dt = 9.225138276950844 cfl multiplier : 0.7066666666666667 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2443e+04 | 735 | 1.732e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20788.527319936133 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 4 [SPH][rank=0]
Info: time since start : 201.522054335 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.30000000000000004, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.20 us (0.9%)
patch tree reduce : 391.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 451.00 ns (0.1%)
LB compute : 362.60 us (97.3%)
LB move op cnt : 0
LB apply : 1933.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.63 us (70.1%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000007,-735,-734.9999999999999)
sum a = (3.2729353081906165e-18,-5.082197683525802e-19,4.607859233063394e-19)
sum e = 1837.5000000099567
sum de = 6.2341624917916505e-19
Info: CFL hydro = 10.501650417394147 sink sink = inf [SPH][rank=0]
Info: cfl dt = 10.501650417394147 cfl multiplier : 0.8044444444444444 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2171e+04 | 735 | 1.743e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20654.966699029934 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 5 [SPH][rank=0]
Info: time since start : 201.71465175400002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.4, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.18 us (0.9%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 461.00 ns (0.1%)
LB compute : 353.22 us (97.3%)
LB move op cnt : 0
LB apply : 1844.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.69 us (69.5%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735,-734.9999999999993,-735)
sum a = (5.1567365828841805e-18,9.0801931945661e-19,1.3145951341386741e-18)
sum e = 1837.5000000099558
sum de = -3.1712913545201005e-18
Info: CFL hydro = 11.352721173598153 sink sink = inf [SPH][rank=0]
Info: cfl dt = 11.352721173598153 cfl multiplier : 0.8696296296296296 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.3402e+04 | 735 | 1.693e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21258.079103438147 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 6 [SPH][rank=0]
Info: time since start : 201.906154726 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.5, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.21 us (0.8%)
patch tree reduce : 430.00 ns (0.1%)
gen split merge : 410.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 360.00 ns (0.1%)
LB compute : 382.27 us (97.5%)
LB move op cnt : 0
LB apply : 1893.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (69.6%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000003,-734.9999999999999,-735.0000000000003)
sum a = (-9.886568560352194e-18,1.6601845766184287e-18,-8.809142651444724e-19)
sum e = 1837.5000000099567
sum de = 7.589415207398531e-19
Info: CFL hydro = 11.92017386487619 sink sink = inf [SPH][rank=0]
Info: cfl dt = 11.92017386487619 cfl multiplier : 0.9130864197530864 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2241e+04 | 735 | 1.740e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20689.277058039992 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 7 [SPH][rank=0]
Info: time since start : 202.098902877 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.6000000000000001, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.91 us (1.0%)
patch tree reduce : 411.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 366.38 us (97.2%)
LB move op cnt : 0
LB apply : 1863.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.36 us (70.9%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000002,-734.9999999999997,-735.0000000000002)
sum a = (7.020209066843641e-18,1.6195269951502222e-18,-2.846030702774449e-19)
sum e = 1837.5000000099558
sum de = -1.5178830414797062e-18
Info: CFL hydro = 12.29855686966423 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.29855686966423 cfl multiplier : 0.9420576131687243 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.1173e+04 | 735 | 1.785e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20166.33529456545 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 8 [SPH][rank=0]
Info: time since start : 202.292598459 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.7000000000000001, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.33 us (0.8%)
patch tree reduce : 371.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 451.00 ns (0.1%)
LB compute : 386.28 us (97.4%)
LB move op cnt : 0
LB apply : 1873.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.39 us (71.5%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999997,-734.9999999999994,-734.9999999999999)
sum a = (1.6059744679941534e-18,1.2265037076242269e-18,2.439454888092385e-19)
sum e = 1837.5000000099549
sum de = -5.908901840045999e-18
Info: CFL hydro = 12.550902239785634 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.550902239785634 cfl multiplier : 0.9613717421124829 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.1318e+04 | 735 | 1.779e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20237.412961557056 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 9 [SPH][rank=0]
Info: time since start : 202.48543671800002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.8, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.33 us (0.8%)
patch tree reduce : 391.00 ns (0.1%)
gen split merge : 400.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 360.00 ns (0.1%)
LB compute : 389.25 us (97.4%)
LB move op cnt : 0
LB apply : 1914.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.32 us (69.0%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000007,-735.0000000000002,-734.9999999999997)
sum a = (-8.239936510889834e-18,-1.3145951341386741e-18,-1.5585406229479126e-18)
sum e = 1837.5000000099556
sum de = 9.405453846311751e-18
Info: CFL hydro = 12.719231298543177 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.719231298543177 cfl multiplier : 0.9742478280749886 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2373e+04 | 735 | 1.735e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20754.12894747856 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 10 [SPH][rank=0]
Info: time since start : 202.67787244500002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0.9, dt = 0.09999999999999998 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.14 us (0.8%)
patch tree reduce : 361.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 362.28 us (97.4%)
LB move op cnt : 0
LB apply : 1793.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.38 us (70.8%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000008,-735.0000000000001,-735.0000000000001)
sum a = (-2.846030702774449e-18,-1.0503208545953324e-18,2.290377089375628e-18)
sum e = 1837.5000000099549
sum de = 7.020209066843641e-18
Info: CFL hydro = 12.831558291123326 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.831558291123326 cfl multiplier : 0.9828318853833258 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.3114e+04 | 735 | 1.705e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21117.095526525147 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 11 [SPH][rank=0]
Info: time since start : 202.870009423 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.44 us (0.9%)
patch tree reduce : 360.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 441.00 ns (0.1%)
LB compute : 357.03 us (97.2%)
LB move op cnt : 0
LB apply : 1783.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.60 us (67.9%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999998,-735.0000000000013,-734.9999999999998)
sum a = (4.174178364069192e-18,-4.811147140404426e-19,1.2468324983583301e-18)
sum e = 1837.5000000099544
sum de = -2.4936649967166602e-18
Info: CFL hydro = 12.906559440326541 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.906559440326541 cfl multiplier : 0.9885545902555505 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.1391e+04 | 735 | 1.776e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20273.025849346883 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 12 [SPH][rank=0]
Info: time since start : 203.06287592500001 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.1, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.43 us (0.9%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 461.00 ns (0.1%)
LB compute : 386.56 us (97.4%)
LB move op cnt : 0
LB apply : 2.10 us (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.46 us (67.7%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999995,-735.0000000000002,-735.0000000000001)
sum a = (-8.063753657860939e-19,3.2729353081906165e-18,-1.3552527156068805e-20)
sum e = 1837.5000000099535
sum de = -8.456776945386935e-18
Info: CFL hydro = 12.956685627293245 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.956685627293245 cfl multiplier : 0.9923697268370336 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2444e+04 | 735 | 1.732e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20789.116758782682 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 13 [SPH][rank=0]
Info: time since start : 203.255818529 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.2000000000000002, dt = 0.09999999999999987 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.32 us (0.9%)
patch tree reduce : 491.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 461.00 ns (0.1%)
LB compute : 358.63 us (97.3%)
LB move op cnt : 0
LB apply : 1853.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.32 us (71.2%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000005,-735.0000000000006,-735.0000000000005)
sum a = (-4.797594613248357e-18,-9.554531645028508e-19,-1.2739375526704677e-18)
sum e = 1837.5000000099526
sum de = 7.399679827213568e-18
Info: CFL hydro = 12.990237500990482 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.990237500990482 cfl multiplier : 0.9949131512246892 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.3177e+04 | 735 | 1.702e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21148.078670617648 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 14 [SPH][rank=0]
Info: time since start : 203.44780943400002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.3, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.30 us (0.8%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 632.00 ns (0.2%)
LB compute : 391.07 us (97.4%)
LB move op cnt : 0
LB apply : 1933.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.37 us (70.9%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-735.0000000000003,-735.0000000000008,-734.9999999999995)
sum a = (-2.9002408113987244e-18,-1.8566962203814263e-18,-2.168404344971009e-19)
sum e = 1837.500000009951
sum de = 9.595189226496714e-18
Info: CFL hydro = 13.012748880137543 sink sink = inf [SPH][rank=0]
Info: cfl dt = 13.012748880137543 cfl multiplier : 0.9966087674831261 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2223e+04 | 735 | 1.741e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20680.57117209954 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 15 [SPH][rank=0]
Info: time since start : 203.640810627 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.4000000000000001, dt = 0.09999999999999987 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.28 us (0.9%)
patch tree reduce : 360.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 356.96 us (97.3%)
LB move op cnt : 0
LB apply : 1854.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.49 us (68.2%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999993,-735.0000000000002,-735.0000000000001)
sum a = (2.1751806085490433e-18,4.607859233063394e-19,1.2739375526704677e-18)
sum e = 1837.5000000099508
sum de = -2.846030702774449e-18
Info: CFL hydro = 13.027909020178868 sink sink = inf [SPH][rank=0]
Info: cfl dt = 13.027909020178868 cfl multiplier : 0.997739178322084 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.2204e+04 | 735 | 1.742e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20671.292174532053 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 16 [SPH][rank=0]
Info: time since start : 203.833615865 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.5, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.22 us (0.8%)
patch tree reduce : 361.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 371.28 us (97.3%)
LB move op cnt : 0
LB apply : 1864.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (69.8%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999993,-735.0000000000014,-735.0000000000001)
sum a = (7.914675859144182e-18,-1.0909784360635388e-18,2.303929616531697e-19)
sum e = 1837.5000000099483
sum de = 6.613633252161577e-18
Info: CFL hydro = 13.038177458584492 sink sink = inf [SPH][rank=0]
Info: cfl dt = 13.038177458584492 cfl multiplier : 0.998492785548056 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.1404e+04 | 735 | 1.775e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20279.47824958361 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 17 [SPH][rank=0]
Info: time since start : 204.026720783 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.6, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.70 us (0.9%)
patch tree reduce : 452.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 551.00 ns (0.1%)
LB compute : 418.82 us (97.4%)
LB move op cnt : 0
LB apply : 1984.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.48 us (70.4%)
Info: conservation infos : [sph::Model][rank=0]
sum v = (-734.9999999999997,-735,-734.9999999999999)
sum a = (4.5672016515951874e-18,8.334804200982315e-19,-1.1655173354219173e-18)
sum e = 1837.5000000099494
sum de = 3.713392440762853e-18
Info: CFL hydro = 13.045193913020135 sink sink = inf [SPH][rank=0]
Info: cfl dt = 13.045193913020135 cfl multiplier : 0.9989951903653708 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.1130e+04 | 735 | 1.787e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20145.52117218731 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 18 [SPH][rank=0]
Info: time since start : 204.219200852 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.7000000000000002, dt = 0.09999999999999987 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 735 min = 735 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 735 min = 735 [LoadBalance][rank=0]
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.8%)
patch tree reduce : 641.00 ns (0.2%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 551.00 ns (0.1%)
LB compute : 386.30 us (97.3%)
LB move op cnt : 0
LB apply : 1903.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (70.7%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.2962445542109251 unconverged cnt = 4
Info: conservation infos : [sph::Model][rank=0]
sum v = (-723.9990621681544,-723.9989589588944,-723.999080836305)
sum a = (-1.8905775382715984e-18,-1.7279472123987727e-18,8.944667923005412e-19)
sum e = 1809.9971019165423
sum de = 4.851804721872632e-18
Info: CFL hydro = 12.89661811090234 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.89661811090234 cfl multiplier : 0.9993301269102473 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 2.1158e+04 | 724 | 3.422e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 10520.529805698365 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 19 [SPH][rank=0]
Info: time since start : 204.428614748 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.8, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 724 min = 724 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 724 min = 724 [LoadBalance][rank=0]
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 : 3.15 us (0.9%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 410.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 370.00 ns (0.1%)
LB compute : 357.31 us (97.3%)
LB move op cnt : 0
LB apply : 1934.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.30 us (69.9%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.30056167192620137 unconverged cnt = 19
Info: conservation infos : [sph::Model][rank=0]
sum v = (-681.9968714883288,-681.9966204670972,-681.9967286191298)
sum a = (-3.0696474008495844e-18,1.260385025514399e-18,2.710505431213761e-19)
sum e = 1704.9902204838572
sum de = -3.2526065174565133e-19
Info: CFL hydro = 12.789615550227142 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.789615550227142 cfl multiplier : 0.9995534179401648 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.9764e+04 | 682 | 3.451e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 10432.369253941582 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 20 [SPH][rank=0]
Info: time since start : 204.63764578700003 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 1.9000000000000001, dt = 0.09999999999999987 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 682 min = 682 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 682 min = 682 [LoadBalance][rank=0]
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 : 3.42 us (0.9%)
patch tree reduce : 390.00 ns (0.1%)
gen split merge : 400.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 501.00 ns (0.1%)
LB compute : 361.76 us (97.3%)
LB move op cnt : 0
LB apply : 1884.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.35 us (70.7%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.29625970614486274 unconverged cnt = 34
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.29755427911452514 unconverged cnt = 2
Info: conservation infos : [sph::Model][rank=0]
sum v = (-612.9949011129379,-612.9946819821328,-612.9947797392336)
sum a = (-5.210946691508456e-18,1.3078188705606397e-18,-1.4907779871675686e-18)
sum e = 1532.4843627578634
sum de = -4.174178364069192e-18
Info: CFL hydro = 12.861008218464518 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.861008218464518 cfl multiplier : 0.9997022786267765 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.3603e+04 | 613 | 4.506e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 7988.919014264137 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 21 [SPH][rank=0]
Info: time since start : 204.85511147000003 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 613 min = 613 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 613 min = 613 [LoadBalance][rank=0]
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 : 3.39 us (0.9%)
patch tree reduce : 350.00 ns (0.1%)
gen split merge : 400.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 681.00 ns (0.2%)
LB compute : 364.03 us (97.2%)
LB move op cnt : 0
LB apply : 1964.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.60 us (69.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.30671834896297384 unconverged cnt = 47
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.307306761640605 unconverged cnt = 2
Info: conservation infos : [sph::Model][rank=0]
sum v = (-532.9935533115207,-532.9935945920531,-532.9934668702186)
sum a = (-4.858580985450667e-18,7.928228386300251e-19,4.0657581468206416e-19)
sum e = 1332.480614721095
sum de = -4.20128341838133e-18
Info: CFL hydro = 12.685397983886743 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.685397983886743 cfl multiplier : 0.9998015190845176 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.0928e+04 | 533 | 4.877e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 7380.854857990513 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 22 [SPH][rank=0]
Info: time since start : 205.072708626 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.1, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 533 min = 533 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 533 min = 533 [LoadBalance][rank=0]
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 : 3.46 us (1.0%)
patch tree reduce : 361.00 ns (0.1%)
gen split merge : 390.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 471.00 ns (0.1%)
LB compute : 334.11 us (96.9%)
LB move op cnt : 0
LB apply : 1934.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.24 us (68.9%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.2962764478718123 unconverged cnt = 55
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.31069429926209974 unconverged cnt = 2
Info: conservation infos : [sph::Model][rank=0]
sum v = (-433.9927338736061,-433.99250861466425,-433.99249597858295)
sum a = (2.846030702774449e-19,2.0328790734103208e-19,-1.734723475976807e-18)
sum e = 1084.9777384174881
sum de = -8.971772977317549e-18
Info: CFL hydro = 12.731993623699823 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.731993623699823 cfl multiplier : 0.9998676793896785 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 9.9759e+03 | 434 | 4.350e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 8274.959432586038 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 23 [SPH][rank=0]
Info: time since start : 205.27938880000002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.2, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 434 min = 434 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 434 min = 434 [LoadBalance][rank=0]
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 : 3.21 us (0.9%)
patch tree reduce : 350.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 451.00 ns (0.1%)
LB compute : 364.19 us (97.3%)
LB move op cnt : 0
LB apply : 1954.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.25 us (69.2%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3109500729833232 unconverged cnt = 57
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3156077724953211 unconverged cnt = 4
Info: conservation infos : [sph::Model][rank=0]
sum v = (-341.99193142018476,-341.9921450989848,-341.9925044093816)
sum a = (6.098637220230962e-20,-8.605854744103691e-19,-4.336808689942018e-19)
sum e = 854.9765809341818
sum de = -3.1712913545201005e-18
Info: CFL hydro = 12.667051368296761 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.667051368296761 cfl multiplier : 0.9999117862597856 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 7.9167e+03 | 342 | 4.320e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 8333.38271634203 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 24 [SPH][rank=0]
Info: time since start : 205.481041572 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.3000000000000003, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 342 min = 342 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 342 min = 342 [LoadBalance][rank=0]
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 : 3.45 us (0.8%)
patch tree reduce : 532.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 379.23 us (92.9%)
LB move op cnt : 0
LB apply : 2.24 us (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.42 us (69.7%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3020184485298713 unconverged cnt = 56
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.328576480695004 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3614341287645044 unconverged cnt = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (-245.9921769555506,-245.99263707959494,-245.99229063406636)
sum a = (-4.87890977618477e-19,1.3552527156068805e-19,-3.9302328752599536e-19)
sum e = 614.9771046889904
sum de = -2.9544509200229996e-18
Info: CFL hydro = 12.502435326683461 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.502435326683461 cfl multiplier : 0.9999411908398571 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 4.6236e+03 | 246 | 5.321e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 6766.257583118734 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 25 [SPH][rank=0]
Info: time since start : 205.68646765300002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.4000000000000004, dt = 0.09999999999999964 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 246 min = 246 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 246 min = 246 [LoadBalance][rank=0]
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 : 3.04 us (0.8%)
patch tree reduce : 361.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 371.00 ns (0.1%)
LB compute : 380.43 us (97.4%)
LB move op cnt : 0
LB apply : 2.18 us (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.31 us (70.4%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.32774218784270714 unconverged cnt = 53
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.343573140831559 unconverged cnt = 10
Info: conservation infos : [sph::Model][rank=0]
sum v = (-166.99278662167782,-166.99320940422416,-166.99324858622313)
sum a = (1.6940658945086007e-19,-4.54009659728305e-19,5.149960319306146e-19)
sum e = 417.47924469434525
sum de = -2.1141942363467336e-18
Info: CFL hydro = 12.481349233596474 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.481349233596474 cfl multiplier : 0.999960793893238 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 3.9516e+03 | 167 | 4.226e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 8518.332990628147 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 26 [SPH][rank=0]
Info: time since start : 205.876591356 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.5, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 167 min = 167 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 167 min = 167 [LoadBalance][rank=0]
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 : 3.34 us (0.8%)
patch tree reduce : 432.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 441.00 ns (0.1%)
LB compute : 388.11 us (97.5%)
LB move op cnt : 0
LB apply : 1974.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.27 us (67.4%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3058811268159916 unconverged cnt = 40
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3364692394975908 unconverged cnt = 17
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.36630356642675027 unconverged cnt = 5
Info: conservation infos : [sph::Model][rank=0]
sum v = (-91.99460804607467,-91.99491947485237,-91.99457795712206)
sum a = (2.168404344971009e-19,1.7618285302889447e-19,2.0328790734103208e-19)
sum e = 229.98410556558173
sum de = 2.168404344971009e-19
Info: CFL hydro = 12.1850822359298 sink sink = inf [SPH][rank=0]
Info: cfl dt = 12.1850822359298 cfl multiplier : 0.999973862595492 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.7629e+03 | 92 | 5.219e-02 | 0 % | 2 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 6898.370673830553 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 27 [SPH][rank=0]
Info: time since start : 206.07181298700002 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.6, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 92 min = 92 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 92 min = 92 [LoadBalance][rank=0]
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 : 3.16 us (0.9%)
patch tree reduce : 371.00 ns (0.1%)
gen split merge : 410.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 352.44 us (97.3%)
LB move op cnt : 0
LB apply : 1863.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.35 us (71.1%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.33322739338798624 unconverged cnt = 28
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3665501327267849 unconverged cnt = 18
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.4032051459994634 unconverged cnt = 9
Info: conservation infos : [sph::Model][rank=0]
sum v = (-43.996389524961266,-43.996148049667234,-43.996195016207054)
sum a = (1.4907779871675686e-19,1.0842021724855044e-19,-1.6263032587282567e-19)
sum e = 109.9887327309434
sum de = 1.3552527156068805e-19
Info: CFL hydro = 11.902519640066306 sink sink = inf [SPH][rank=0]
Info: cfl dt = 11.902519640066306 cfl multiplier : 0.9999825750636614 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 8.4946e+02 | 44 | 5.180e-02 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 6950.092601296303 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 28 [SPH][rank=0]
Info: time since start : 206.261830897 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 2.7, dt = 0.10000000000000009 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 44 min = 44 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 44 min = 44 [LoadBalance][rank=0]
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 : 3.15 us (0.8%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 384.61 us (97.4%)
LB move op cnt : 0
LB apply : 1924.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.40 us (70.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3175808603008827 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.34933894633097096 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.3842728409640681 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.42270012506047494 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.46497013756652245 unconverged cnt = 8
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.5114671513231748 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.5610681639786606 unconverged cnt = 3
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.609277261524694 unconverged cnt = 3
Info: conservation infos : [sph::Model][rank=0]
sum v = (-10.998438229211487,-10.998627137787167,-10.998397802857001)
sum a = (-6.776263578034403e-21,0,2.710505431213761e-20)
sum e = 27.495463315833934
sum de = -1.3552527156068805e-20
Info: CFL hydro = 15.643421624235584 sink sink = inf [SPH][rank=0]
Info: cfl dt = 15.643421624235584 cfl multiplier : 0.9999883833757742 [sph::Model][rank=0]
Info: processing rate infos : [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank | rate (N.s^-1) | Nobj | t compute (s) | MPI | alloc | mem (max) |
---------------------------------------------------------------------------------------
| 0 | 1.0461e+02 | 11 | 1.052e-01 | 0 % | 1 % | 179.82 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 3423.5083531606797 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 29 [SPH][rank=0]
Info: time since start : 206.503478623 (s) [SPH][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
Convert PNG sequence to Image sequence in mpl#
162 import matplotlib.animation as animation
163
164
165 def show_image_sequence(glob_str):
166
167 if render_gif and shamrock.sys.world_rank() == 0:
168
169 import glob
170
171 files = sorted(glob.glob(glob_str))
172
173 from PIL import Image
174
175 image_array = []
176 for my_file in files:
177 image = Image.open(my_file)
178 image_array.append(image)
179
180 if not image_array:
181 raise RuntimeError(f"Warning: No images found for glob pattern: {glob_str}")
182
183 pixel_x, pixel_y = image_array[0].size
184
185 # Create the figure and axes objects
186 # Remove axes, ticks, and frame & set aspect ratio
187 dpi = 200
188 fig = plt.figure(dpi=dpi)
189 plt.gca().set_position((0, 0, 1, 1))
190 plt.gcf().set_size_inches(pixel_x / dpi, pixel_y / dpi)
191 plt.axis("off")
192
193 # Set the initial image with correct aspect ratio
194 im = plt.imshow(image_array[0], animated=True, aspect="auto")
195
196 def update(i):
197 im.set_array(image_array[i])
198 return (im,)
199
200 # Create the animation object
201 ani = animation.FuncAnimation(
202 fig,
203 update,
204 frames=len(image_array),
205 interval=50,
206 blit=True,
207 repeat_delay=10,
208 )
209
210 return ani
211
212
213 # If the animation is not returned only a static image will be shown in the doc
214 glob_str = os.path.join(dump_folder, f"{sim_name}_*.png")
215 ani = show_image_sequence(glob_str)
216
217 if render_gif and shamrock.sys.world_rank() == 0:
218 # To save the animation using Pillow as a gif
219 # writer = animation.PillowWriter(fps=15,
220 # metadata=dict(artist='Me'),
221 # bitrate=1800)
222 # ani.save('scatter.gif', writer=writer)
223
224 # Show the animation
225 plt.show()
Total running time of the script: (0 minutes 11.463 seconds)
Estimated memory usage: 122 MB