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"
            },
            "evol_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"
        },
        "neigh_cache_strategy": "two_stage",
        "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
    }
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 735 ( 7.4e+02 ) Ntotal = 735 ( 7.4e+02 rank min = 9.4e+03 max = 7.4e+02) rate = 7.350000e+02 N.s^-1
SPH setup: the generation step took : 0.07885734700000001 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.50 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     : 1.70 us    (0.3%)
   patch tree reduce : 1.34 us    (0.2%)
   gen split merge   : 741.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 811.00 ns  (0.1%)
   LB compute        : 563.57 us  (96.4%)
   LB move op cnt    : 0
   LB apply          : 13.00 us   (2.2%)
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.007876017 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.1% 0.0% |   337.11 MB |   460.80 kB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.088267209 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/sphinx/user_guide/sph/smoothing_length_tolerance.html
---------------- 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.24 us    (0.8%)
   patch tree reduce : 441.00 ns  (0.1%)
   gen split merge   : 381.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 371.00 ns  (0.1%)
   LB compute        : 391.65 us  (97.6%)
   LB move op cnt    : 0
   LB apply          : 1.95 us    (0.5%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.22 us    (66.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
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    | 7.9317e+03 |  735 |      1 | 9.267e-02 | 0.0% |   1.3% 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/sphinx/user_guide/sph/smoothing_length_tolerance.html

Show the current solvergraph

G cluster_2214 time step cluster_2192 sink accretion cluster_2191 if_has_sinks cluster_2213 point mass accretion cluster_2212 if_has_ext_force_accretion cluster_2211 ext force accretion cluster_2203 sink predictor cluster_2202 sink predictor body cluster_2198 sink ext force cluster_2197 sink ext force body cluster_2167 attach fields cluster_2175 leapfrog predictor cluster_2170 half step cluster_2173 half step cluster_2179 part killing step n_2184 SetEdge e_7162 has_sinks n_2184->e_7162 n_2158 SetEdge n_2184->n_2158 n_2203 sink predictor e_7162->n_2203 n_2198 sink ext force e_7162->n_2198 e_7162->n_2198 e_7144 m n_2158->e_7144 n_2158->e_7144 n_2159 SetEdge n_2158->n_2159 n_2158->n_2159 n_2189 SinkParticlesAccreteQuantities e_7144->n_2189 e_7141 patchdatas n_2159->e_7141 n_2159->e_7141 n_2159->e_7141 n_2159->e_7141 n_2160 GetObjCntFromLayer n_2159->n_2160 n_2159->n_2160 n_2159->n_2160 n_2159->n_2160 e_7141->n_2160 e_7141->n_2160 e_7141->n_2160 e_7141->n_2160 n_2161 GetFieldRefFromLayer e_7141->n_2161 e_7141->n_2161 e_7141->n_2161 e_7141->n_2161 n_2162 GetFieldRefFromLayer e_7141->n_2162 e_7141->n_2162 e_7141->n_2162 n_2163 GetFieldRefFromLayer e_7141->n_2163 e_7141->n_2163 e_7141->n_2163 n_2164 GetFieldRefFromLayer e_7141->n_2164 e_7141->n_2164 n_2165 GetFieldRefFromLayer e_7141->n_2165 e_7141->n_2165 n_2166 GetFieldRefFromLayer e_7141->n_2166 e_7141->n_2166 e_7142 Npart_patch n_2160->e_7142 n_2160->e_7142 n_2160->e_7142 n_2160->e_7142 n_2160->n_2161 n_2160->n_2161 n_2160->n_2161 n_2160->n_2161 n_2188 SinkParticlesFlagAccreteHard e_7142->n_2188 e_7142->n_2189 n_2190 SinkParticlesEvictAccretedParticles e_7142->n_2190 n_2209 SinkParticlesFlagAccreteHard e_7142->n_2209 n_2210 SinkParticlesEvictAccretedParticles e_7142->n_2210 n_2168 ForwardEuler e_7142->n_2168 n_2169 ForwardEuler e_7142->n_2169 n_2174 ForwardEuler e_7142->n_2174 n_2171 ForwardEuler e_7142->n_2171 n_2172 ForwardEuler e_7142->n_2172 e_7145 xyz n_2161->e_7145 n_2161->e_7145 n_2161->e_7145 n_2161->e_7145 n_2161->n_2162 n_2161->n_2162 n_2161->n_2162 n_2161->n_2209 e_7145->n_2188 e_7145->n_2189 e_7145->n_2209 n_2177 GetParticlesOutsideSphere e_7145->n_2177 e_7146 vxyz n_2162->e_7146 n_2162->e_7146 n_2162->e_7146 n_2162->n_2163 n_2162->n_2163 n_2162->n_2163 e_7146->n_2189 e_7146->n_2174 e_7147 axyz n_2163->e_7147 n_2163->e_7147 n_2163->e_7147 n_2163->n_2188 n_2163->n_2164 n_2163->n_2164 e_7147->n_2189 e_7147->n_2168 e_7147->n_2171 e_7163 sink_accretion_table n_2188->e_7163 n_2188->n_2189 e_7134 sink_pos e_7134->n_2188 n_2196 SinkSelfGravityHost e_7134->n_2196 e_7134->n_2196 e_7140 sink_accretion_radius e_7140->n_2188 e_7163->n_2190 n_2189->e_7134 n_2189->e_7163 e_7135 sink_vel n_2189->e_7135 e_7136 sink_acc_sph n_2189->e_7136 e_7139 sink_angular_momentum n_2189->e_7139 e_7138 sink_mass n_2189->e_7138 n_2189->n_2190 e_7132 dt e_7132->n_2189 n_2199 MapEdge e_7132->n_2199 n_2201 ForwardEulerHost e_7132->n_2201 n_2180 MapEdge e_7132->n_2180 e_7132->n_2174 e_7135->n_2201 n_2200 ForwardEulerHost2Deriv e_7136->n_2200 e_7138->n_2196 e_7138->n_2196 n_2190->e_7141 n_2185 FreeAlloc n_2190->n_2185 n_2185->e_7145 n_2185->e_7145 n_2186 FreeAlloc n_2185->n_2186 n_2212_end n_2185->n_2212_end n_2186->e_7146 n_2187 FreeAlloc n_2186->n_2187 n_2187->e_7147 n_2206 SetEdge n_2187->n_2206 e_7167 ext_force_accretion_pos n_2206->e_7167 n_2207 SetEdge n_2206->n_2207 n_2208 MapEdge e_7167->n_2208 e_7167->n_2209 e_7168 ext_force_accretion_racc n_2207->e_7168 n_2207->n_2208 e_7168->n_2209 e_7170 has_ext_force_accretion n_2208->e_7170 n_2212 if_has_ext_force_accretion n_2208->n_2212 e_7170->n_2212 n_2212->n_2159 true n_2212->n_2212_end false e_7169 ext_force_accretion_table n_2209->e_7169 n_2209->n_2210 e_7169->n_2210 n_2210->e_7141 n_2210->n_2185 n_2212_end->n_2203 n_2203->n_2198 true n_2203_end n_2203->n_2203_end false n_2193 SetEdge n_2198->n_2193 true n_2198->n_2193 true n_2198_end n_2198->n_2198_end false n_2198->n_2198_end false e_7164 G n_2193->e_7164 n_2193->e_7164 n_2194 SetEdge n_2193->n_2194 n_2193->n_2194 e_7164->n_2196 e_7164->n_2196 e_7165 epsilon_grav_sink n_2194->e_7165 n_2194->e_7165 n_2195 SetEdge n_2194->n_2195 n_2194->n_2195 e_7165->n_2196 e_7165->n_2196 e_7137 sink_acc_ext n_2195->e_7137 n_2195->e_7137 n_2195->n_2196 n_2195->n_2196 e_7137->n_2200 n_2196->e_7137 n_2196->e_7137 n_2196->n_2198_end n_2196->n_2198_end n_2198_end->n_2199 e_7166 dt_half n_2199->e_7166 n_2199->n_2200 e_7166->n_2200 n_2200->e_7135 n_2200->n_2201 n_2201->e_7134 n_2201->n_2203_end n_2203_end->n_2180 n_2180->n_2158 e_7143 dt_half n_2180->e_7143 e_7143->n_2168 e_7143->n_2169 e_7143->n_2171 e_7143->n_2172 e_7148 uint n_2164->e_7148 n_2164->e_7148 n_2164->n_2165 n_2164->n_2165 e_7149 duint n_2165->e_7149 n_2165->e_7149 n_2165->n_2166 n_2165->n_2166 e_7149->n_2169 e_7149->n_2172 n_2166->n_2198 e_7150 hpart n_2166->e_7150 n_2166->e_7150 n_2166->n_2168 n_2168->e_7146 n_2168->n_2169 n_2169->e_7148 n_2169->n_2174 n_2174->e_7145 n_2174->n_2171 n_2171->e_7146 n_2171->n_2172 n_2172->e_7148 n_2176 FreeAlloc n_2172->n_2176 e_7151 part_to_remove n_2176->e_7151 n_2176->n_2177 n_2178 KillParticles e_7151->n_2178 n_2177->e_7151 n_2177->n_2178 n_2178->n_2159 n_2178->e_7141


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 : 640.528694573 (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     : 21.98 us   (3.8%)
   patch tree reduce : 1.44 us    (0.3%)
   gen split merge   : 802.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1.10 us    (0.2%)
   LB compute        : 536.07 us  (93.1%)
   LB move op cnt    : 0
   LB apply          : 4.02 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.76 us    (71.0%)
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.8868e+04 |  735 |      1 | 2.546e-02 | 0.2% |   1.5% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14139.234620450117 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 2                                                       [SPH][rank=0]
Info: time since start : 640.7364729750001 (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     : 9.51 us    (2.4%)
   patch tree reduce : 1.64 us    (0.4%)
   gen split merge   : 821.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.00 us    (0.3%)
   LB compute        : 371.95 us  (94.1%)
   LB move op cnt    : 0
   LB apply          : 4.02 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.59 us    (66.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    | 2.8337e+04 |  735 |      1 | 2.594e-02 | 0.1% |   1.3% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13879.147248229667 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 3                                                       [SPH][rank=0]
Info: time since start : 640.9309582740001 (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.40 us   (2.9%)
   patch tree reduce : 1.38 us    (0.4%)
   gen split merge   : 911.00 ns  (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.13 us    (0.3%)
   LB compute        : 332.18 us  (93.3%)
   LB move op cnt    : 0
   LB apply          : 3.68 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.71 us    (65.2%)
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.8048e+04 |  735 |      1 | 2.621e-02 | 0.1% |   1.2% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13737.761324179764 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 4                                                       [SPH][rank=0]
Info: time since start : 641.1284269470001 (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     : 10.49 us   (2.7%)
   patch tree reduce : 1.36 us    (0.3%)
   gen split merge   : 722.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.10 us    (0.3%)
   LB compute        : 366.22 us  (93.7%)
   LB move op cnt    : 0
   LB apply          : 4.19 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.33 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.8542e+04 |  735 |      1 | 2.575e-02 | 0.1% |   1.2% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13979.69596724215 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 5                                                       [SPH][rank=0]
Info: time since start : 641.3244125990001 (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.37 us   (2.6%)
   patch tree reduce : 1.69 us    (0.4%)
   gen split merge   : 1.05 us    (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.01 us    (0.3%)
   LB compute        : 377.78 us  (93.8%)
   LB move op cnt    : 0
   LB apply          : 4.01 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.76 us    (71.4%)
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.8283e+04 |  735 |      1 | 2.599e-02 | 0.1% |   1.2% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13852.752018499883 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 6                                                       [SPH][rank=0]
Info: time since start : 641.5216450280001 (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     : 9.88 us    (2.9%)
   patch tree reduce : 1.41 us    (0.4%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 982.00 ns  (0.3%)
   LB compute        : 323.53 us  (93.4%)
   LB move op cnt    : 0
   LB apply          : 3.38 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.14 us    (67.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    | 2.8284e+04 |  735 |      1 | 2.599e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13853.307481009244 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 7                                                       [SPH][rank=0]
Info: time since start : 641.71816946 (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     : 10.80 us   (2.9%)
   patch tree reduce : 1.66 us    (0.4%)
   gen split merge   : 912.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 912.00 ns  (0.2%)
   LB compute        : 346.17 us  (93.2%)
   LB move op cnt    : 0
   LB apply          : 4.17 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.14 us    (64.9%)
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.8761e+04 |  735 |      1 | 2.556e-02 | 0.1% |   1.2% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14087.117237372304 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 8                                                       [SPH][rank=0]
Info: time since start : 641.9144020450001 (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.32 us   (2.9%)
   patch tree reduce : 1.38 us    (0.4%)
   gen split merge   : 902.00 ns  (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.25 us    (0.4%)
   LB compute        : 332.15 us  (93.1%)
   LB move op cnt    : 0
   LB apply          : 3.57 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.01 us    (70.6%)
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    | 2.8724e+04 |  735 |      1 | 2.559e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14068.970627350247 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 9                                                       [SPH][rank=0]
Info: time since start : 642.111679179 (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     : 11.56 us   (3.2%)
   patch tree reduce : 1.43 us    (0.4%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 992.00 ns  (0.3%)
   LB compute        : 337.39 us  (93.1%)
   LB move op cnt    : 0
   LB apply          : 3.48 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.12 us    (63.5%)
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    | 3.2806e+04 |  735 |      1 | 2.240e-02 | 0.1% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16068.080815483661 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 10                                                      [SPH][rank=0]
Info: time since start : 642.3034854450001 (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     : 9.86 us    (2.8%)
   patch tree reduce : 1.38 us    (0.4%)
   gen split merge   : 992.00 ns  (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.02 us    (0.3%)
   LB compute        : 334.30 us  (93.4%)
   LB move op cnt    : 0
   LB apply          : 3.93 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.25 us    (68.0%)
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    | 2.7998e+04 |  735 |      1 | 2.625e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 13713.297834430286 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 11                                                      [SPH][rank=0]
Info: time since start : 642.50027832 (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     : 11.33 us   (2.9%)
   patch tree reduce : 1.44 us    (0.4%)
   gen split merge   : 1.27 us    (0.3%)
   split / merge op  : 0/0
   apply split merge : 912.00 ns  (0.2%)
   LB compute        : 363.10 us  (93.3%)
   LB move op cnt    : 0
   LB apply          : 4.67 us    (1.2%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.96 us    (69.7%)
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    | 2.8741e+04 |  735 |      1 | 2.557e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14077.398396670362 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 12                                                      [SPH][rank=0]
Info: time since start : 642.696138106 (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.58 us   (2.6%)
   patch tree reduce : 1.36 us    (0.3%)
   gen split merge   : 811.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.15 us    (0.3%)
   LB compute        : 379.40 us  (93.8%)
   LB move op cnt    : 0
   LB apply          : 4.09 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.65 us    (71.7%)
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.9246e+04 |  735 |      1 | 2.513e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14324.366770477222 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 13                                                      [SPH][rank=0]
Info: time since start : 642.8907948770001 (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.02 us   (2.6%)
   patch tree reduce : 1.79 us    (0.5%)
   gen split merge   : 802.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.12 us    (0.3%)
   LB compute        : 333.98 us  (87.6%)
   LB move op cnt    : 0
   LB apply          : 26.69 us   (7.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.71 us    (71.6%)
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.2930e+04 |  735 |      1 | 2.232e-02 | 0.1% |   1.5% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16128.772839182358 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 14                                                      [SPH][rank=0]
Info: time since start : 643.082700649 (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     : 8.21 us    (2.3%)
   patch tree reduce : 1.33 us    (0.4%)
   gen split merge   : 811.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.13 us    (0.3%)
   LB compute        : 329.41 us  (94.0%)
   LB move op cnt    : 0
   LB apply          : 3.62 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.00 us    (65.1%)
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.1331e+04 |  735 |      1 | 2.346e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 15345.605412906561 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 15                                                      [SPH][rank=0]
Info: time since start : 643.2755345710001 (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.97 us   (2.8%)
   patch tree reduce : 1.44 us    (0.4%)
   gen split merge   : 821.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.01 us    (0.3%)
   LB compute        : 363.98 us  (93.6%)
   LB move op cnt    : 0
   LB apply          : 3.91 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.52 us    (69.2%)
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.3516e+04 |  735 |      1 | 2.193e-02 | 0.1% |   1.2% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 16415.817771290065 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 16                                                      [SPH][rank=0]
Info: time since start : 643.468337225 (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.59 us   (2.8%)
   patch tree reduce : 1.42 us    (0.4%)
   gen split merge   : 842.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 882.00 ns  (0.2%)
   LB compute        : 349.21 us  (93.6%)
   LB move op cnt    : 0
   LB apply          : 3.82 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.05 us    (64.3%)
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.8613e+04 |  735 |      1 | 2.569e-02 | 0.1% |   1.1% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14014.706020324527 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 17                                                      [SPH][rank=0]
Info: time since start : 643.663784879 (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     : 9.20 us    (2.5%)
   patch tree reduce : 1.32 us    (0.4%)
   gen split merge   : 992.00 ns  (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.46 us    (0.4%)
   LB compute        : 339.70 us  (93.7%)
   LB move op cnt    : 0
   LB apply          : 3.63 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.44 us    (69.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.8937e+04 |  735 |      1 | 2.540e-02 | 0.1% |   1.0% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 14173.06429602745 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 18                                                      [SPH][rank=0]
Info: time since start : 643.859286543 (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.04 us   (2.8%)
   patch tree reduce : 1.74 us    (0.5%)
   gen split merge   : 801.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.08 us    (0.3%)
   LB compute        : 328.86 us  (93.2%)
   LB move op cnt    : 0
   LB apply          : 3.54 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.60 us    (71.0%)
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.6427e+04 |  724 |      1 | 4.407e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 8167.923800894322 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 19                                                      [SPH][rank=0]
Info: time since start : 644.0732288830001 (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     : 18.55 us   (4.9%)
   patch tree reduce : 1.67 us    (0.4%)
   gen split merge   : 802.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.46 us    (0.4%)
   LB compute        : 340.35 us  (90.7%)
   LB move op cnt    : 0
   LB apply          : 3.88 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.20 us    (67.3%)
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.4454e+04 |  682 |      1 | 4.718e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 7629.629309373591 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 20                                                      [SPH][rank=0]
Info: time since start : 644.290283167 (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     : 10.34 us   (2.9%)
   patch tree reduce : 1.56 us    (0.4%)
   gen split merge   : 821.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.04 us    (0.3%)
   LB compute        : 327.74 us  (93.3%)
   LB move op cnt    : 0
   LB apply          : 3.46 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.42 us    (62.5%)
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.5968e+03 |  613 |      1 | 6.388e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5635.966092963406 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 21                                                      [SPH][rank=0]
Info: time since start : 644.5211095440001 (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.25 us   (2.8%)
   patch tree reduce : 1.50 us    (0.4%)
   gen split merge   : 1.16 us    (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.09 us    (0.3%)
   LB compute        : 341.76 us  (93.3%)
   LB move op cnt    : 0
   LB apply          : 3.56 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.88 us    (66.1%)
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.6288e+03 |  533 |      1 | 6.177e-02 | 0.0% |   1.5% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5828.092165125721 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 22                                                      [SPH][rank=0]
Info: time since start : 644.746039575 (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.08 us   (2.6%)
   patch tree reduce : 1.39 us    (0.4%)
   gen split merge   : 811.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.08 us    (0.3%)
   LB compute        : 364.68 us  (93.6%)
   LB move op cnt    : 0
   LB apply          : 4.10 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.01 us    (70.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.1608e+03 |  434 |      1 | 6.061e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5939.807377986544 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 23                                                      [SPH][rank=0]
Info: time since start : 644.965027585 (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     : 10.03 us   (2.9%)
   patch tree reduce : 1.37 us    (0.4%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 982.00 ns  (0.3%)
   LB compute        : 323.87 us  (93.3%)
   LB move op cnt    : 0
   LB apply          : 3.77 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.49 us    (65.7%)
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.7048e+03 |  342 |      1 | 5.995e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6005.009879659106 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 24                                                      [SPH][rank=0]
Info: time since start : 645.172008275 (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     : 9.90 us    (2.5%)
   patch tree reduce : 1.36 us    (0.4%)
   gen split merge   : 811.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.40 us    (0.4%)
   LB compute        : 364.33 us  (93.7%)
   LB move op cnt    : 0
   LB apply          : 4.45 us    (1.1%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.38 us    (67.9%)
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.3532e+03 |  246 |      1 | 7.336e-02 | 0.1% |   1.5% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 4907.140476887638 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 25                                                      [SPH][rank=0]
Info: time since start : 645.389924014 (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     : 8.51 us    (2.4%)
   patch tree reduce : 1.16 us    (0.3%)
   gen split merge   : 742.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.02 us    (0.3%)
   LB compute        : 330.01 us  (94.0%)
   LB move op cnt    : 0
   LB apply          : 3.57 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1.88 us    (62.5%)
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    | 2.9769e+03 |  167 |      1 | 5.610e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6417.174869779902 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 26                                                      [SPH][rank=0]
Info: time since start : 645.5849289470001 (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     : 8.06 us    (2.0%)
   patch tree reduce : 1.43 us    (0.4%)
   gen split merge   : 821.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1.11 us    (0.3%)
   LB compute        : 385.60 us  (94.9%)
   LB move op cnt    : 0
   LB apply          : 3.48 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.17 us    (67.7%)
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.3552e+03 |   92 |      1 | 6.789e-02 | 0.0% |   1.3% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5303.0539949150525 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 27                                                      [SPH][rank=0]
Info: time since start : 645.7854665340001 (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     : 8.76 us    (2.5%)
   patch tree reduce : 1.37 us    (0.4%)
   gen split merge   : 872.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 912.00 ns  (0.3%)
   LB compute        : 327.80 us  (93.8%)
   LB move op cnt    : 0
   LB apply          : 3.44 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.01 us    (68.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    | 6.3387e+02 |   44 |      1 | 6.941e-02 | 0.0% |   1.4% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 5186.223678456211 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 28                                                      [SPH][rank=0]
Info: time since start : 645.9817503820001 (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.05 us   (2.8%)
   patch tree reduce : 1.36 us    (0.4%)
   gen split merge   : 1.20 us    (0.3%)
   split / merge op  : 0/0
   apply split merge : 1.21 us    (0.3%)
   LB compute        : 328.82 us  (93.2%)
   LB move op cnt    : 0
   LB apply          : 3.44 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1.98 us    (64.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.1088e+01 |   11 |      1 | 1.357e-01 | 0.0% |   1.3% 0.0% |   337.11 MB |   460.80 kB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 2653.7804396080205 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 29                                                      [SPH][rank=0]
Info: time since start : 646.2435653690001 (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.042 seconds)

Estimated memory usage: 160 MB

Gallery generated by Sphinx-Gallery