Killing sphere for SPH simulation#

This simple example shows how to setup a killing sphere for sph simulations

 9 import shamrock
10
11 # If we use the shamrock executable to run this script instead of the python interpreter,
12 # we should not initialize the system as the shamrock executable needs to handle specific MPI logic
13 if not shamrock.sys.is_initialized():
14     shamrock.change_loglevel(1)
15     shamrock.sys.init("0:0")
16
17
18 def is_in_sphere(pt):
19     x, y, z = pt
20     return (x**2 + y**2 + z**2) < 1

Setup parameters

27 dr = 0.1
28 pmass = 1
29
30 C_cour = 0.3
31 C_force = 0.25
32
33 bsize = 4
34
35
36 render_gif = True
37
38 dump_folder = "_to_trash"
39 sim_name = "kill_particle_sphere"
40
41 import os
42
43 # Create the dump directory if it does not exist
44 if shamrock.sys.world_rank() == 0:
45     os.makedirs(dump_folder, exist_ok=True)

Setup

50 ctx = shamrock.Context()
51 ctx.pdata_layout_new()
52
53 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
54
55 cfg = model.gen_default_config()
56 # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2)
57 # cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV)
58 # cfg.set_eos_locally_isothermalLP07(cs0=cs0, q=q, r0=r0)
59 cfg.set_artif_viscosity_VaryingCD10(
60     alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
61 )
62 cfg.set_boundary_periodic()
63 cfg.set_eos_adiabatic(1.00001)

The important part to enable killing

67 cfg.add_kill_sphere(center=(0.0, 0.0, 0.0), radius=4.0)

Rest of the setup

 72 cfg.print_status()
 73 model.set_solver_config(cfg)
 74
 75 model.init_scheduler(int(1e7), 1)
 76
 77 bmin = (-bsize, -bsize, -bsize)
 78 bmax = (bsize, bsize, bsize)
 79 model.resize_simulation_box(bmin, bmax)
 80
 81 model.set_particle_mass(pmass)
 82
 83 setup = model.get_setup()
 84 lat = setup.make_generator_lattice_hcp(dr, (-bsize, -bsize, -bsize), (bsize, bsize, bsize))
 85
 86 thesphere = setup.make_modifier_filter(parent=lat, filter=is_in_sphere)
 87
 88 offset_sphere = setup.make_modifier_offset(
 89     parent=thesphere, offset_position=(0.0, 0.0, 0.0), offset_velocity=(-1.0, -1.0, -1.0)
 90 )
 91
 92 setup.apply_setup(offset_sphere)
 93
 94 model.set_value_in_a_box("uint", "f64", 1, bmin, bmax)
 95
 96 model.set_cfl_cour(C_cour)
 97 model.set_cfl_force(C_force)
 98
 99 model.change_htolerance(1.3)
100 model.timestep()
101 model.change_htolerance(1.1)
----- SPH Solver configuration -----
[
    {
        "artif_viscosity": {
            "alpha_max": 1.0,
            "alpha_min": 0.0,
            "alpha_u": 1.0,
            "av_type": "varying_cd10",
            "beta_AV": 2.0,
            "sigma_decay": 0.1
        },
        "boundary_config": {
            "bc_type": "periodic"
        },
        "cfl_config": {
            "cfl_cour": 0.0,
            "cfl_force": 0.0,
            "cfl_multiplier_stiffness": 2.0,
            "eta_sink": 0.05
        },
        "combined_dtdiv_divcurlv_compute": false,
        "debug_dump_filename": "",
        "do_debug_dump": false,
        "enable_particle_reordering": false,
        "eos_config": {
            "Tvec": "f64_3",
            "eos_type": "adiabatic",
            "gamma": 1.00001
        },
        "epsilon_h": 1e-06,
        "ext_force_config": {
            "force_list": []
        },
        "gpart_mass": 0.0,
        "h_iter_per_subcycles": 50,
        "h_max_subcycles_count": 100,
        "htol_up_coarse_cycle": 1.1,
        "htol_up_fine_cycle": 1.1,
        "kernel_id": "M4<f64>",
        "mhd_config": {
            "mhd_type": "none"
        },
        "particle_killing": [
            {
                "center": [
                    0.0,
                    0.0,
                    0.0
                ],
                "radius": 4.0,
                "type": "sphere"
            }
        ],
        "particle_reordering_step_freq": 1000,
        "self_grav_config": {
            "softening_length": 1e-09,
            "softening_mode": "plummer",
            "type": "none"
        },
        "show_neigh_stats": false,
        "smoothing_length_config": {
            "type": "density_based"
        },
        "time_state": {
            "cfl_multiplier": 0.01,
            "dt_sph": 0.0,
            "time": 0.0
        },
        "tree_reduction_level": 3,
        "type_id": "sycl::vec<f64,3>",
        "unit_sys": null,
        "use_two_stage_search": true
    }
]
------------------------------------
Info: pushing data in scheduler, N = 735                              [DataInserterUtility][rank=0]
Info: reattributing data ...                                          [DataInserterUtility][rank=0]
Info: reattributing data done in  2.11 ms                             [DataInserterUtility][rank=0]
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 4.09 us    (40.2%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1753.00 ns (0.4%)
   patch tree reduce : 661.00 ns  (0.1%)
   gen split merge   : 642.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 802.00 ns  (0.2%)
   LB compute        : 477.27 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 12.29 us   (2.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.73 us    (61.8%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1753.00 ns (0.5%)
   patch tree reduce : 361.00 ns  (0.1%)
   gen split merge   : 411.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 311.00 ns  (0.1%)
   LB compute        : 329.97 us  (97.6%)
   LB move op cnt    : 0
   LB apply          : 2.05 us    (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.75 us    (42.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1393.00 ns (0.4%)
   patch tree reduce : 340.00 ns  (0.1%)
   gen split merge   : 371.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 291.00 ns  (0.1%)
   LB compute        : 386.77 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 2.09 us    (0.5%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Final load balancing step 0 of 3                                          [SPH setup][rank=0]
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.85 us    (59.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1623.00 ns (0.4%)
   patch tree reduce : 381.00 ns  (0.1%)
   gen split merge   : 411.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 331.00 ns  (0.1%)
   LB compute        : 366.14 us  (97.8%)
   LB move op cnt    : 0
   LB apply          : 2.10 us    (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.75 us    (65.1%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1362.00 ns (0.4%)
   patch tree reduce : 331.00 ns  (0.1%)
   gen split merge   : 360.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 311.00 ns  (0.1%)
   LB compute        : 342.21 us  (97.9%)
   LB move op cnt    : 0
   LB apply          : 1984.00 ns (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.85 us    (66.4%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1673.00 ns (0.5%)
   patch tree reduce : 361.00 ns  (0.1%)
   gen split merge   : 460.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 300.00 ns  (0.1%)
   LB compute        : 357.42 us  (97.8%)
   LB move op cnt    : 0
   LB apply          : 2.12 us    (0.6%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Final load balancing step 1 of 3                                          [SPH setup][rank=0]
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.58 us    (66.1%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1203.00 ns (0.3%)
   patch tree reduce : 340.00 ns  (0.1%)
   gen split merge   : 360.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 251.00 ns  (0.1%)
   LB compute        : 368.37 us  (98.0%)
   LB move op cnt    : 0
   LB apply          : 2.11 us    (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (66.6%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1673.00 ns (0.5%)
   patch tree reduce : 381.00 ns  (0.1%)
   gen split merge   : 391.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 310.00 ns  (0.1%)
   LB compute        : 355.97 us  (97.8%)
   LB move op cnt    : 0
   LB apply          : 2.05 us    (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.65 us    (64.6%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1673.00 ns (0.5%)
   patch tree reduce : 340.00 ns  (0.1%)
   gen split merge   : 381.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 251.00 ns  (0.1%)
   LB compute        : 348.87 us  (97.8%)
   LB move op cnt    : 0
   LB apply          : 1944.00 ns (0.5%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Final load balancing step 2 of 3                                          [SPH setup][rank=0]
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.46 us    (67.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1724.00 ns (0.4%)
   patch tree reduce : 381.00 ns  (0.1%)
   gen split merge   : 411.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 251.00 ns  (0.1%)
   LB compute        : 422.65 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 2.03 us    (0.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.79 us    (68.2%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1243.00 ns (0.4%)
   patch tree reduce : 361.00 ns  (0.1%)
   gen split merge   : 370.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 320.00 ns  (0.1%)
   LB compute        : 330.02 us  (97.9%)
   LB move op cnt    : 0
   LB apply          : 2.00 us    (0.6%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.83 us    (67.2%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1232.00 ns (0.4%)
   patch tree reduce : 371.00 ns  (0.1%)
   gen split merge   : 371.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 240.00 ns  (0.1%)
   LB compute        : 328.43 us  (97.8%)
   LB move op cnt    : 0
   LB apply          : 2.06 us    (0.6%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: the setup took : 0.08152437 s                                             [SPH setup][rank=0]
Warning: .change_htolerance(val) is deprecated,                                       [SPH][rank=0]
    -> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
    see: https://shamrock-code.github.io/Shamrock/mkdocs/models/sph/smoothing_length_tolerance
---------------- t = 0, dt = 0 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.14 us    (0.9%)
   patch tree reduce : 421.00 ns  (0.1%)
   gen split merge   : 401.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 361.00 ns  (0.1%)
   LB compute        : 351.48 us  (97.3%)
   LB move op cnt    : 0
   LB apply          : 2.09 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.65 us    (70.6%)
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]
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 = (-5.590417451878382e-18,1.904130065427667e-18,2.2090619264392153e-18)
    sum e = 1837.5
    sum de = -9.75781955236954e-19
Info: CFL hydro = 0.1305431643141645 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 0.1305431643141645 cfl multiplier : 0.01                        [sph::Model][rank=0]
Info: 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.0142e+04 |  735 |      1 | 7.247e-02 | 0.1% |   1.6% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr)                                             [sph::Model][rank=0]
Warning: .change_htolerance(val) is deprecated,                                       [SPH][rank=0]
    -> calling this is replaced internally by .change_htolerances(coarse=val, fine=min(val, 1.1))
    see: https://shamrock-code.github.io/Shamrock/mkdocs/models/sph/smoothing_length_tolerance

Draw utilities#

107 import matplotlib.pyplot as plt
108 import numpy as np
109
110
111 def plot_state(iplot):
112
113     pos = ctx.collect_data()["xyz"]
114
115     if shamrock.sys.world_rank() == 0:
116         X = pos[:, 0]
117         Y = pos[:, 1]
118         Z = pos[:, 2]
119
120         plt.cla()
121
122         ax.set_xlim3d(bmin[0], bmax[0])
123         ax.set_ylim3d(bmin[1], bmax[1])
124         ax.set_zlim3d(bmin[2], bmax[2])
125
126         ax.scatter(X, Y, Z, s=1)
127
128         ax.set_title(f"t = {model.get_time():.2f} ")
129
130         plt.savefig(os.path.join(dump_folder, f"{sim_name}_{iplot:04}.png"))

Run the simulation#

136 nstop = 28  # To be increased when epmty simulations will be fixed
137 dt_stop = 0.1
138
139 t_stop = [i * dt_stop for i in range(nstop + 1)]
140
141 # Init MPL
142 fig = plt.figure(dpi=120)
143 ax = fig.add_subplot(111, projection="3d")
144
145 iplot = 0
146 istop = 0
147 for ttarg in t_stop:
148
149     model.evolve_until(ttarg)
150
151     # if do_plots:
152     plot_state(iplot)
153
154     iplot += 1
155     istop += 1
156
157 plt.close(fig)
Info: iteration since start : 1                                                       [SPH][rank=0]
Info: time since start : 1143.3721865730001 (s)                                       [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0, dt = 0.1 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 12.06 us   (2.3%)
   patch tree reduce : 1152.00 ns (0.2%)
   gen split merge   : 792.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1092.00 ns (0.2%)
   LB compute        : 491.69 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.86 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.57 us    (69.8%)
Warning: the unit system is not set                                           [sph::Config][rank=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.211948935988307e-18,-2.6495190590114515e-18,-8.131516293641283e-19)
    sum e = 1837.5000000099571
    sum de = -6.478107980600889e-18
Info: CFL hydro = 4.4384722909924195 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 4.4384722909924195 cfl multiplier : 0.34                        [sph::Model][rank=0]
Info: 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.8718e+04 |  735 |      1 | 1.898e-02 | 0.3% |   2.0% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 18963.85399876303 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 2                                                       [SPH][rank=0]
Info: time since start : 1143.5883551860002 (s)                                       [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.1, dt = 0.1 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.59 us    (1.6%)
   patch tree reduce : 802.00 ns  (0.2%)
   gen split merge   : 631.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1253.00 ns (0.3%)
   LB compute        : 447.21 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 3.17 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.03 us    (73.3%)
Warning: the unit system is not set                                           [sph::Config][rank=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 = (3.2255014631443757e-18,-1.5585406229479126e-19,1.6398557858843255e-18)
    sum e = 1837.5000000099576
    sum de = -1.3552527156068805e-19
Info: CFL hydro = 7.310448194870652 sink sink = inf                                   [SPH][rank=0]
Info: cfl dt = 7.310448194870652 cfl multiplier : 0.56                         [sph::Model][rank=0]
Info: 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    | 4.0507e+04 |  735 |      1 | 1.815e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19840.081229701464 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 3                                                       [SPH][rank=0]
Info: time since start : 1143.786211268 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.2, dt = 0.10000000000000003 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 5.48 us    (1.3%)
   patch tree reduce : 611.00 ns  (0.1%)
   gen split merge   : 612.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 561.00 ns  (0.1%)
   LB compute        : 416.12 us  (96.5%)
   LB move op cnt    : 0
   LB apply          : 2.90 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.62 us    (72.2%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-8.321251673826247e-18,-3.889575293791747e-18,-1.802486111757151e-18)
    sum e = 1837.500000009957
    sum de = 2.846030702774449e-18
Info: CFL hydro = 9.225138276950844 sink sink = inf                                   [SPH][rank=0]
Info: cfl dt = 9.225138276950844 cfl multiplier : 0.7066666666666667           [sph::Model][rank=0]
Info: 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.9164e+04 |  735 |      1 | 1.877e-02 | 0.1% |   1.2% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19182.365523491822 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 4                                                       [SPH][rank=0]
Info: time since start : 1143.983480436 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.30000000000000004, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.06 us   (1.9%)
   patch tree reduce : 1162.00 ns (0.2%)
   gen split merge   : 822.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.2%)
   LB compute        : 516.93 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.75 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.13 us    (74.6%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (3.2729353081906165e-18,-5.082197683525802e-19,4.607859233063394e-19)
    sum e = 1837.5000000099562
    sum de = 6.2341624917916505e-19
Info: CFL hydro = 10.501650417394147 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 10.501650417394147 cfl multiplier : 0.8044444444444444          [sph::Model][rank=0]
Info: 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    | 4.0772e+04 |  735 |      1 | 1.803e-02 | 0.1% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19969.905352633577 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 5                                                       [SPH][rank=0]
Info: time since start : 1144.180490807 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.4, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.03 us    (1.5%)
   patch tree reduce : 731.00 ns  (0.2%)
   gen split merge   : 551.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1132.00 ns (0.3%)
   LB compute        : 396.23 us  (96.2%)
   LB move op cnt    : 0
   LB apply          : 2.62 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.42 us    (69.5%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (5.1567365828841805e-18,9.0801931945661e-19,1.3145951341386741e-18)
    sum e = 1837.5000000099562
    sum de = -3.1712913545201005e-18
Info: CFL hydro = 11.352721173598153 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 11.352721173598153 cfl multiplier : 0.8696296296296296          [sph::Model][rank=0]
Info: 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    | 4.1069e+04 |  735 |      1 | 1.790e-02 | 0.1% |   1.2% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 20115.41667325539 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 6                                                       [SPH][rank=0]
Info: time since start : 1144.376138818 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.5, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 9.96 us    (1.9%)
   patch tree reduce : 1122.00 ns (0.2%)
   gen split merge   : 851.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1082.00 ns (0.2%)
   LB compute        : 500.23 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.62 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.41 us    (76.2%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-9.886568560352194e-18,1.6601845766184287e-18,-8.809142651444724e-19)
    sum e = 1837.5000000099565
    sum de = 7.589415207398531e-19
Info: CFL hydro = 11.92017386487619 sink sink = inf                                   [SPH][rank=0]
Info: cfl dt = 11.92017386487619 cfl multiplier : 0.9130864197530864           [sph::Model][rank=0]
Info: 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    | 4.0461e+04 |  735 |      1 | 1.817e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19817.81264636282 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 7                                                       [SPH][rank=0]
Info: time since start : 1144.572688233 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.6000000000000001, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.85 us    (1.5%)
   patch tree reduce : 942.00 ns  (0.2%)
   gen split merge   : 651.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 771.00 ns  (0.2%)
   LB compute        : 432.56 us  (96.1%)
   LB move op cnt    : 0
   LB apply          : 3.15 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.92 us    (69.3%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (7.020209066843641e-18,1.6195269951502222e-18,-2.846030702774449e-19)
    sum e = 1837.5000000099549
    sum de = -1.5178830414797062e-18
Info: CFL hydro = 12.29855686966423 sink sink = inf                                   [SPH][rank=0]
Info: cfl dt = 12.29855686966423 cfl multiplier : 0.9420576131687243           [sph::Model][rank=0]
Info: 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    | 4.1384e+04 |  735 |      1 | 1.776e-02 | 0.1% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 20269.843416585703 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 8                                                       [SPH][rank=0]
Info: time since start : 1144.768341513 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.7000000000000001, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 9.41 us    (2.0%)
   patch tree reduce : 1382.00 ns (0.3%)
   gen split merge   : 782.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1142.00 ns (0.2%)
   LB compute        : 459.80 us  (95.4%)
   LB move op cnt    : 0
   LB apply          : 3.83 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.04 us    (73.7%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (1.6059744679941534e-18,1.2265037076242269e-18,2.439454888092385e-19)
    sum e = 1837.5000000099542
    sum de = -5.908901840045999e-18
Info: CFL hydro = 12.550902239785634 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.550902239785634 cfl multiplier : 0.9613717421124829          [sph::Model][rank=0]
Info: 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    | 4.0154e+04 |  735 |      1 | 1.830e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19667.126240060104 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 9                                                       [SPH][rank=0]
Info: time since start : 1144.9649121280002 (s)                                       [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.8, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.04 us   (2.1%)
   patch tree reduce : 1173.00 ns (0.2%)
   gen split merge   : 791.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1113.00 ns (0.2%)
   LB compute        : 463.71 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.84 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 4.10 us    (76.9%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-8.239936510889834e-18,-1.3145951341386741e-18,-1.5585406229479126e-18)
    sum e = 1837.5000000099546
    sum de = 9.405453846311751e-18
Info: CFL hydro = 12.719231298543177 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.719231298543177 cfl multiplier : 0.9742478280749886          [sph::Model][rank=0]
Info: 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    | 4.0373e+04 |  735 |      1 | 1.821e-02 | 0.1% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19774.6676619916 (tsim/hr)                              [sph::Model][rank=0]
Info: iteration since start : 10                                                      [SPH][rank=0]
Info: time since start : 1145.162261166 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 0.9, dt = 0.09999999999999998 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 8.99 us    (1.9%)
   patch tree reduce : 1042.00 ns (0.2%)
   gen split merge   : 831.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1072.00 ns (0.2%)
   LB compute        : 461.17 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.38 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (70.9%)
Warning: the unit system is not set                                           [sph::Config][rank=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 = (-2.846030702774449e-18,-1.0503208545953324e-18,2.290377089375628e-18)
    sum e = 1837.5000000099558
    sum de = 7.020209066843641e-18
Info: CFL hydro = 12.831558291123326 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.831558291123326 cfl multiplier : 0.9828318853833258          [sph::Model][rank=0]
Info: 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    | 4.1002e+04 |  735 |      1 | 1.793e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 20082.53363475728 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 11                                                      [SPH][rank=0]
Info: time since start : 1145.359673072 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.65 us   (2.0%)
   patch tree reduce : 1473.00 ns (0.3%)
   gen split merge   : 792.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1172.00 ns (0.2%)
   LB compute        : 503.30 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 3.98 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.68 us    (71.3%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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.174178364069192e-18,-4.811147140404426e-19,1.2468324983583301e-18)
    sum e = 1837.5000000099549
    sum de = -2.4936649967166602e-18
Info: CFL hydro = 12.906559440326541 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.906559440326541 cfl multiplier : 0.9885545902555505          [sph::Model][rank=0]
Info: 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    | 4.0095e+04 |  735 |      1 | 1.833e-02 | 0.1% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19638.329611688052 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 12                                                      [SPH][rank=0]
Info: time since start : 1145.55767122 (s)                                            [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.1, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.58 us   (2.1%)
   patch tree reduce : 1132.00 ns (0.2%)
   gen split merge   : 1022.00 ns (0.2%)
   split / merge op  : 0/0
   apply split merge : 1111.00 ns (0.2%)
   LB compute        : 477.80 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.96 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.83 us    (70.3%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-8.063753657860939e-19,3.2729353081906165e-18,-1.3552527156068805e-20)
    sum e = 1837.5000000099526
    sum de = -8.456776945386935e-18
Info: CFL hydro = 12.956685627293245 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.956685627293245 cfl multiplier : 0.9923697268370336          [sph::Model][rank=0]
Info: 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    | 4.0381e+04 |  735 |      1 | 1.820e-02 | 0.1% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19778.291940115305 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 13                                                      [SPH][rank=0]
Info: time since start : 1145.7558136 (s)                                             [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.2000000000000002, dt = 0.09999999999999987 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.33 us   (2.0%)
   patch tree reduce : 1452.00 ns (0.3%)
   gen split merge   : 901.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.2%)
   LB compute        : 482.26 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.81 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.79 us    (72.8%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-4.797594613248357e-18,-9.554531645028508e-19,-1.2739375526704677e-18)
    sum e = 1837.500000009954
    sum de = 7.399679827213568e-18
Info: CFL hydro = 12.990237500990482 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.990237500990482 cfl multiplier : 0.9949131512246892          [sph::Model][rank=0]
Info: 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    | 4.0660e+04 |  735 |      1 | 1.808e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19915.107429004413 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 14                                                      [SPH][rank=0]
Info: time since start : 1145.954928267 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.3, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 9.94 us    (2.0%)
   patch tree reduce : 1533.00 ns (0.3%)
   gen split merge   : 1132.00 ns (0.2%)
   split / merge op  : 0/0
   apply split merge : 1212.00 ns (0.2%)
   LB compute        : 481.05 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.68 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.06 us    (73.5%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-2.9002408113987244e-18,-1.8566962203814263e-18,-2.168404344971009e-19)
    sum e = 1837.5000000099512
    sum de = 9.595189226496714e-18
Info: CFL hydro = 13.012748880137543 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 13.012748880137543 cfl multiplier : 0.9966087674831261          [sph::Model][rank=0]
Info: 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.9842e+04 |  735 |      1 | 1.845e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19514.494553450382 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 15                                                      [SPH][rank=0]
Info: time since start : 1146.152088409 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.4000000000000001, dt = 0.09999999999999987 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 8.90 us    (2.0%)
   patch tree reduce : 992.00 ns  (0.2%)
   gen split merge   : 772.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 922.00 ns  (0.2%)
   LB compute        : 431.54 us  (95.4%)
   LB move op cnt    : 0
   LB apply          : 3.57 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.83 us    (74.3%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
Warning: the unit system is not set                                           [sph::Config][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-734.9999999999993,-735.0000000000002,-735.0000000000001)
    sum a = (2.1751806085490433e-18,4.607859233063394e-19,1.2739375526704677e-18)
    sum e = 1837.5000000099517
    sum de = -2.846030702774449e-18
Info: CFL hydro = 13.027909020178868 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 13.027909020178868 cfl multiplier : 0.997739178322084           [sph::Model][rank=0]
Info: 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.8705e+04 |  735 |      1 | 1.899e-02 | 0.1% |   1.2% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 18957.63048576556 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 16                                                      [SPH][rank=0]
Info: time since start : 1146.349748803 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.5, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.88 us    (1.7%)
   patch tree reduce : 1062.00 ns (0.2%)
   gen split merge   : 801.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1013.00 ns (0.2%)
   LB compute        : 442.52 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.60 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.92 us    (67.5%)
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (7.914675859144182e-18,-1.0909784360635388e-18,2.303929616531697e-19)
    sum e = 1837.5000000099483
    sum de = 6.613633252161577e-18
Info: CFL hydro = 13.038177458584492 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 13.038177458584492 cfl multiplier : 0.998492785548056           [sph::Model][rank=0]
Info: 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    | 4.0011e+04 |  735 |      1 | 1.837e-02 | 0.1% |   1.9% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19597.21196997709 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 17                                                      [SPH][rank=0]
Info: time since start : 1146.546661501 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.6, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.42 us   (2.1%)
   patch tree reduce : 1252.00 ns (0.3%)
   gen split merge   : 791.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1122.00 ns (0.2%)
   LB compute        : 460.49 us  (94.9%)
   LB move op cnt    : 0
   LB apply          : 3.96 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.06 us    (73.5%)
Warning: the unit system is not set                                           [sph::Config][rank=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 = (4.5672016515951874e-18,8.334804200982315e-19,-1.1655173354219173e-18)
    sum e = 1837.5000000099478
    sum de = 3.713392440762853e-18
Info: CFL hydro = 13.045193913020135 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 13.045193913020135 cfl multiplier : 0.9989951903653708          [sph::Model][rank=0]
Info: 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.9615e+04 |  735 |      1 | 1.855e-02 | 0.1% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 19403.394138719734 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 18                                                      [SPH][rank=0]
Info: time since start : 1146.744763114 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.7000000000000002, dt = 0.09999999999999987 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 735 min = 735                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 735 min = 735                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 735
    max = 735
    avg = 735
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.93 us   (2.3%)
   patch tree reduce : 1192.00 ns (0.3%)
   gen split merge   : 882.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 447.40 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.97 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.37 us    (72.9%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.2962445542109251 unconverged cnt = 4
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-1.8905775382715984e-18,-1.7279472123987727e-18,8.944667923005412e-19)
    sum e = 1809.9971019165425
    sum de = 4.851804721872632e-18
Info: CFL hydro = 12.89661811090234 sink sink = inf                                   [SPH][rank=0]
Info: cfl dt = 12.89661811090234 cfl multiplier : 0.9993301269102473           [sph::Model][rank=0]
Info: 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.0767e+04 |  724 |      1 | 3.486e-02 | 0.1% |   1.5% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 10326.275601239495 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 19                                                      [SPH][rank=0]
Info: time since start : 1146.960209689 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.8, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 724 min = 724                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 724 min = 724                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 724
    max = 724
    avg = 724
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 9.97 us    (2.1%)
   patch tree reduce : 1102.00 ns (0.2%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 457.64 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.90 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.34 us    (71.8%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.30056167192620137 unconverged cnt = 19
Warning: the unit system is not set                                           [sph::Config][rank=0]
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.0696474008495844e-18,1.260385025514399e-18,2.710505431213761e-19)
    sum e = 1704.9902204838565
    sum de = -3.2526065174565133e-19
Info: CFL hydro = 12.789615550227142 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.789615550227142 cfl multiplier : 0.9995534179401648          [sph::Model][rank=0]
Info: 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.9104e+04 |  682 |      1 | 3.570e-02 | 0.1% |   1.5% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 10083.99237358862 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 20                                                      [SPH][rank=0]
Info: time since start : 1147.1748970150002 (s)                                       [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 1.9000000000000001, dt = 0.09999999999999987 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 682 min = 682                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 682 min = 682                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 682
    max = 682
    avg = 682
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.47 us    (1.6%)
   patch tree reduce : 751.00 ns  (0.2%)
   gen split merge   : 692.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 641.00 ns  (0.2%)
   LB compute        : 399.01 us  (96.1%)
   LB move op cnt    : 0
   LB apply          : 2.65 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.88 us    (71.2%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.29625970614486274 unconverged cnt = 34
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.29755427911452514 unconverged cnt = 2
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-5.210946691508456e-18,1.3078188705606397e-18,-1.4907779871675686e-18)
    sum e = 1532.484362757863
    sum de = -4.174178364069192e-18
Info: CFL hydro = 12.861008218464518 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.861008218464518 cfl multiplier : 0.9997022786267765          [sph::Model][rank=0]
Info: 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.3447e+04 |  613 |      1 | 4.558e-02 | 0.0% |   1.5% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 7897.352666988634 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 21                                                      [SPH][rank=0]
Info: time since start : 1147.396045507 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 613 min = 613                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 613 min = 613                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 613
    max = 613
    avg = 613
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 10.33 us   (2.1%)
   patch tree reduce : 1142.00 ns (0.2%)
   gen split merge   : 882.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1132.00 ns (0.2%)
   LB compute        : 458.14 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.53 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.88 us    (74.0%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.30671834896297384 unconverged cnt = 47
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.307306761640605 unconverged cnt = 2
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-4.858580985450667e-18,7.928228386300251e-19,4.0657581468206416e-19)
    sum e = 1332.4806147210948
    sum de = -4.20128341838133e-18
Info: CFL hydro = 12.685397983886743 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.685397983886743 cfl multiplier : 0.9998015190845176          [sph::Model][rank=0]
Info: 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.1175e+04 |  533 |      1 | 4.769e-02 | 0.0% |   1.5% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 7548.036648737151 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 22                                                      [SPH][rank=0]
Info: time since start : 1147.615834874 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.1, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 533 min = 533                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 533 min = 533                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 533
    max = 533
    avg = 533
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 9.31 us    (0.7%)
   patch tree reduce : 1242.00 ns (0.1%)
   gen split merge   : 852.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1021.00 ns (0.1%)
   LB compute        : 1233.85 us (98.1%)
   LB move op cnt    : 0
   LB apply          : 3.99 us    (0.3%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.88 us    (77.0%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.2962764478718123 unconverged cnt = 55
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.31069429926209974 unconverged cnt = 2
Warning: the unit system is not set                                           [sph::Config][rank=0]
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.846030702774449e-19,2.0328790734103208e-19,-1.734723475976807e-18)
    sum e = 1084.9777384174886
    sum de = -8.971772977317549e-18
Info: CFL hydro = 12.731993623699823 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.731993623699823 cfl multiplier : 0.9998676793896785          [sph::Model][rank=0]
Info: 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.6724e+03 |  434 |      1 | 4.487e-02 | 0.1% |   1.5% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 8023.147493687963 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 23                                                      [SPH][rank=0]
Info: time since start : 1147.828584718 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.2, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 434 min = 434                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 434 min = 434                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 434
    max = 434
    avg = 434
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 8.37 us    (1.9%)
   patch tree reduce : 812.00 ns  (0.2%)
   gen split merge   : 651.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 842.00 ns  (0.2%)
   LB compute        : 420.24 us  (95.8%)
   LB move op cnt    : 0
   LB apply          : 2.83 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.58 us    (70.3%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3109500729833232 unconverged cnt = 57
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3156077724953211 unconverged cnt = 4
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (6.098637220230962e-20,-8.605854744103691e-19,-4.336808689942018e-19)
    sum e = 854.9765809341825
    sum de = -3.1712913545201005e-18
Info: CFL hydro = 12.667051368296761 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.667051368296761 cfl multiplier : 0.9999117862597856          [sph::Model][rank=0]
Info: 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.9832e+03 |  342 |      1 | 4.897e-02 | 0.0% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 7350.70593729587 (tsim/hr)                              [sph::Model][rank=0]
Info: iteration since start : 24                                                      [SPH][rank=0]
Info: time since start : 1148.0391489840001 (s)                                       [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.3000000000000003, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 342 min = 342                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 342 min = 342                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 342
    max = 342
    avg = 342
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.95 us    (1.7%)
   patch tree reduce : 1102.00 ns (0.2%)
   gen split merge   : 581.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1302.00 ns (0.3%)
   LB compute        : 443.50 us  (95.5%)
   LB move op cnt    : 0
   LB apply          : 3.65 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.35 us    (70.8%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3020184485298713 unconverged cnt = 56
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.328576480695004 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3614341287645044 unconverged cnt = 1
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (-4.87890977618477e-19,1.3552527156068805e-19,-3.9302328752599536e-19)
    sum e = 614.9771046889904
    sum de = -2.9544509200229996e-18
Info: CFL hydro = 12.502435326683461 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.502435326683461 cfl multiplier : 0.9999411908398571          [sph::Model][rank=0]
Info: 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    | 4.6252e+03 |  246 |      1 | 5.319e-02 | 0.0% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6768.641897385702 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 25                                                      [SPH][rank=0]
Info: time since start : 1148.251672043 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.4000000000000004, dt = 0.09999999999999964 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 246 min = 246                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 246 min = 246                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 246
    max = 246
    avg = 246
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 4.46 us    (1.1%)
   patch tree reduce : 611.00 ns  (0.1%)
   gen split merge   : 561.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 551.00 ns  (0.1%)
   LB compute        : 404.27 us  (96.7%)
   LB move op cnt    : 0
   LB apply          : 3.17 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.77 us    (72.6%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.32774218784270714 unconverged cnt = 53
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.343573140831559 unconverged cnt = 10
Warning: the unit system is not set                                           [sph::Config][rank=0]
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.6940658945086007e-19,-4.54009659728305e-19,5.149960319306146e-19)
    sum e = 417.47924469434525
    sum de = -2.1141942363467336e-18
Info: CFL hydro = 12.481349233596474 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 12.481349233596474 cfl multiplier : 0.999960793893238           [sph::Model][rank=0]
Info: 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    | 4.0118e+03 |  167 |      1 | 4.163e-02 | 0.0% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 8648.119477229795 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 26                                                      [SPH][rank=0]
Info: time since start : 1148.44311106 (s)                                            [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.5, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 167 min = 167                              [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 167 min = 167                         [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 167
    max = 167
    avg = 167
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.11 us    (1.5%)
   patch tree reduce : 751.00 ns  (0.2%)
   gen split merge   : 731.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 951.00 ns  (0.2%)
   LB compute        : 402.64 us  (96.1%)
   LB move op cnt    : 0
   LB apply          : 2.73 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.96 us    (68.8%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3058811268159916 unconverged cnt = 40
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3364692394975908 unconverged cnt = 17
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.36630356642675027 unconverged cnt = 5
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (2.168404344971009e-19,1.7618285302889447e-19,2.0328790734103208e-19)
    sum e = 229.98410556558173
    sum de = 2.168404344971009e-19
Info: CFL hydro = 12.1850822359298 sink sink = inf                                    [SPH][rank=0]
Info: cfl dt = 12.1850822359298 cfl multiplier : 0.999973862595492             [sph::Model][rank=0]
Info: 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.7717e+03 |   92 |      1 | 5.193e-02 | 0.0% |   1.4% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6932.669604680697 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 27                                                      [SPH][rank=0]
Info: time since start : 1148.640262955 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.6, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 92 min = 92                                [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 92 min = 92                           [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 92
    max = 92
    avg = 92
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.37 us    (1.4%)
   patch tree reduce : 1011.00 ns (0.2%)
   gen split merge   : 711.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.2%)
   LB compute        : 432.64 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 3.31 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.96 us    (70.1%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.33322739338798624 unconverged cnt = 28
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3665501327267849 unconverged cnt = 18
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.4032051459994634 unconverged cnt = 9
Warning: the unit system is not set                                           [sph::Config][rank=0]
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 = (1.4907779871675686e-19,1.0842021724855044e-19,-1.6263032587282567e-19)
    sum e = 109.9887327309434
    sum de = 1.3552527156068805e-19
Info: CFL hydro = 11.902519640066306 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 11.902519640066306 cfl multiplier : 0.9999825750636614          [sph::Model][rank=0]
Info: 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.0638e+02 |   44 |      1 | 5.457e-02 | 0.0% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 6597.63149427777 (tsim/hr)                              [sph::Model][rank=0]
Info: iteration since start : 28                                                      [SPH][rank=0]
Info: time since start : 1148.835949436 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed
---------------- t = 2.7, dt = 0.10000000000000009 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 44 min = 44                                [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 44 min = 44                           [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 44
    max = 44
    avg = 44
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.61 us    (1.5%)
   patch tree reduce : 891.00 ns  (0.2%)
   gen split merge   : 692.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 651.00 ns  (0.1%)
   LB compute        : 428.25 us  (96.0%)
   LB move op cnt    : 0
   LB apply          : 3.43 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.62 us    (70.4%)
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3175808603008827 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.34933894633097096 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3842728409640681 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.42270012506047494 unconverged cnt = 11
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.46497013756652245 unconverged cnt = 8
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.5114671513231748 unconverged cnt = 7
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.5610681639786606 unconverged cnt = 3
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.609277261524694 unconverged cnt = 3
Warning: the unit system is not set                                           [sph::Config][rank=0]
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.776263578034403e-21,0,2.710505431213761e-20)
    sum e = 27.495463315833934
    sum de = -1.3552527156068805e-20
Info: CFL hydro = 15.643421624235584 sink sink = inf                                  [SPH][rank=0]
Info: cfl dt = 15.643421624235584 cfl multiplier : 0.9999883833757742          [sph::Model][rank=0]
Info: 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.0636e+02 |   11 |      1 | 1.034e-01 | 0.0% |   1.3% 0.0% |   935.14 MB |     5.04 MB |
+------+------------+------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 3480.7158335761887 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 29                                                      [SPH][rank=0]
Info: time since start : 1149.078960623 (s)                                           [SPH][rank=0]
Info: collected : 1 patches                                                [PatchScheduler][rank=0]
adding ->  xyz
adding ->  vxyz
adding ->  axyz
adding ->  axyz_ext
adding ->  hpart
adding ->  uint
adding ->  duint
adding ->  alpha_AV
adding ->  divv
adding ->  dtdivv
adding ->  curlv
adding ->  soundspeed

Convert PNG sequence to Image sequence in mpl#

162 import matplotlib.animation as animation
163
164
165 def show_image_sequence(glob_str):
166
167     if render_gif and shamrock.sys.world_rank() == 0:
168
169         import glob
170
171         files = sorted(glob.glob(glob_str))
172
173         from PIL import Image
174
175         image_array = []
176         for my_file in files:
177             image = Image.open(my_file)
178             image_array.append(image)
179
180         if not image_array:
181             raise RuntimeError(f"Warning: No images found for glob pattern: {glob_str}")
182
183         pixel_x, pixel_y = image_array[0].size
184
185         # Create the figure and axes objects
186         # Remove axes, ticks, and frame & set aspect ratio
187         dpi = 200
188         fig = plt.figure(dpi=dpi)
189         plt.gca().set_position((0, 0, 1, 1))
190         plt.gcf().set_size_inches(pixel_x / dpi, pixel_y / dpi)
191         plt.axis("off")
192
193         # Set the initial image with correct aspect ratio
194         im = plt.imshow(image_array[0], animated=True, aspect="auto")
195
196         def update(i):
197             im.set_array(image_array[i])
198             return (im,)
199
200         # Create the animation object
201         ani = animation.FuncAnimation(
202             fig,
203             update,
204             frames=len(image_array),
205             interval=50,
206             blit=True,
207             repeat_delay=10,
208         )
209
210         return ani
211
212
213 # If the animation is not returned only a static image will be shown in the doc
214 glob_str = os.path.join(dump_folder, f"{sim_name}_*.png")
215 ani = show_image_sequence(glob_str)
216
217 if render_gif and shamrock.sys.world_rank() == 0:
218     # To save the animation using Pillow as a gif
219     # writer = animation.PillowWriter(fps=15,
220     #                                 metadata=dict(artist='Me'),
221     #                                 bitrate=1800)
222     # ani.save('scatter.gif', writer=writer)
223
224     # Show the animation
225     plt.show()

Total running time of the script: (0 minutes 12.324 seconds)

Estimated memory usage: 113 MB

Gallery generated by Sphinx-Gallery