Production run: Circular disc & central potential

This example demonstrates how to run a smoothed particle hydrodynamics (SPH) simulation of a circular disc orbiting around a central point mass potential.

The simulation models:

• A central star with a given mass and accretion radius

• A gaseous disc with specified mass, inner/outer radii, and vertical structure

• Artificial viscosity for angular momentum transport

• Locally isothermal equation of state

Also this simulation feature rolling dumps (see purge_old_dumps function) to save disk space.

This example is the accumulation of 3 files in a single one to showcase the complete workflow.

• The actual run script (runscript.py)

• Plot generation (make_plots.py)

• Animation from the plots (plot_to_gif.py)

On a cluster or laptop, one can run the code as follows:

mpirun <your parameters> ./shamrock --sycl-cfg 0:0 --loglevel 1 --rscript runscript.py

then after the run is done (or while it is running), one can run the following to generate the plots:

python make_plots.py

Runscript (runscript.py)

The runscript is the actual simulation with on the fly analysis & rolling dumps

import glob
import json
import os  # for makedirs

import numpy as np

import shamrock

# If we use the shamrock executable to run this script instead of the python interpreter,
# we should not initialize the system as the shamrock executable needs to handle specific MPI logic
if not shamrock.sys.is_initialized():
    shamrock.change_loglevel(1)
    shamrock.sys.init("0:0")

Setup units

si = shamrock.UnitSystem()
sicte = shamrock.Constants(si)
codeu = shamrock.UnitSystem(
    unit_time=sicte.year(),
    unit_length=sicte.au(),
    unit_mass=sicte.sol_mass(),
)
ucte = shamrock.Constants(codeu)
G = ucte.G()

List parameters

# Resolution
Npart = 100000

# Domain decomposition parameters
scheduler_split_val = int(1.0e7)  # split patches with more than 1e7 particles
scheduler_merge_val = scheduler_split_val // 16

# Dump and plot frequency and duration of the simulation
dump_freq_stop = 2
plot_freq_stop = 1

dt_stop = 0.01
nstop = 30

# The list of times at which the simulation will pause for analysis / dumping
t_stop = [i * dt_stop for i in range(nstop + 1)]


# Sink parameters
center_mass = 1.0
center_racc = 0.1

# Disc parameter
disc_mass = 0.01  # sol mass
rout = 10.0  # au
rin = 1.0  # au
H_r_0 = 0.05
q = 0.5
p = 3.0 / 2.0
r0 = 1.0

# Viscosity parameter
alpha_AV = 1.0e-3 / 0.08
alpha_u = 1.0
beta_AV = 2.0

# Integrator parameters
C_cour = 0.3
C_force = 0.25

sim_folder = f"_to_trash/circular_disc_central_pot_{Npart}/"

dump_folder = sim_folder + "dump/"
analysis_folder = sim_folder + "analysis/"
plot_folder = analysis_folder + "plots/"

dump_prefix = dump_folder + "dump_"


# Disc profiles
def sigma_profile(r):
    sigma_0 = 1.0  # We do not care as it will be renormalized
    return sigma_0 * (r / r0) ** (-p)


def kep_profile(r):
    return (G * center_mass / r) ** 0.5


def omega_k(r):
    return kep_profile(r) / r


def cs_profile(r):
    cs_in = (H_r_0 * r0) * omega_k(r0)
    return ((r / r0) ** (-q)) * cs_in

Create the dump directory if it does not exist

if shamrock.sys.world_rank() == 0:
    os.makedirs(sim_folder, exist_ok=True)
    os.makedirs(dump_folder, exist_ok=True)
    os.makedirs(analysis_folder, exist_ok=True)
    os.makedirs(plot_folder, exist_ok=True)

Utility functions and quantities deduced from the base one

# Deduced quantities
pmass = disc_mass / Npart

bsize = rout * 2
bmin = (-bsize, -bsize, -bsize)
bmax = (bsize, bsize, bsize)

cs0 = cs_profile(r0)


def rot_profile(r):
    return ((kep_profile(r) ** 2) - (2 * p + q) * cs_profile(r) ** 2) ** 0.5


def H_profile(r):
    H = cs_profile(r) / omega_k(r)
    # fact = (2.**0.5) * 3. # factor taken from phantom, to fasten thermalizing
    fact = 1.0
    return fact * H

Start the context The context holds the data of the code We then init the layout of the field (e.g. the list of fields used by the solver)

ctx = shamrock.Context()
ctx.pdata_layout_new()

Attach a SPH model to the context

model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")

Dump handling

def get_dump_name(idump):
    return dump_prefix + f"{idump:07}" + ".sham"


def get_vtk_dump_name(idump):
    return dump_prefix + f"{idump:07}" + ".vtk"


def get_ph_dump_name(idump):
    return dump_prefix + f"{idump:07}" + ".phdump"


def get_last_dump():
    res = glob.glob(dump_prefix + "*.sham")

    num_max = -1

    for f in res:
        try:
            dump_num = int(f[len(dump_prefix) : -5])
            if dump_num > num_max:
                num_max = dump_num
        except ValueError:
            pass

    if num_max == -1:
        return None
    else:
        return num_max


def purge_old_dumps():
    if shamrock.sys.world_rank() == 0:
        res = glob.glob(dump_prefix + "*.sham")
        res.sort()

        # The list of dumps to remove (keep the first and last 3 dumps)
        to_remove = res[1:-3]

        for f in to_remove:
            os.remove(f)


idump_last_dump = get_last_dump()

if shamrock.sys.world_rank() == 0:
    print("Last dump:", idump_last_dump)

Last dump: None

Load the last dump if it exists, setup otherwise

if idump_last_dump is not None:
    model.load_from_dump(get_dump_name(idump_last_dump))
else:
    # Generate the default config
    cfg = model.gen_default_config()
    cfg.set_artif_viscosity_ConstantDisc(alpha_u=alpha_u, alpha_AV=alpha_AV, beta_AV=beta_AV)
    cfg.set_eos_locally_isothermalLP07(cs0=cs0, q=q, r0=r0)

    cfg.add_ext_force_point_mass(center_mass, center_racc)
    cfg.add_kill_sphere(center=(0, 0, 0), radius=bsize)  # kill particles outside the simulation box

    cfg.set_units(codeu)
    cfg.set_particle_mass(pmass)
    # Set the CFL
    cfg.set_cfl_cour(C_cour)
    cfg.set_cfl_force(C_force)

    # Enable this to debug the neighbor counts
    # cfg.set_show_neigh_stats(True)

    # Standard way to set the smoothing length (e.g. Price et al. 2018)
    cfg.set_smoothing_length_density_based()

    # Standard density based smoothing lenght but with a neighbor count limit
    # Use it if you have large slowdowns due to giant particles
    # I recommend to use it if you have a circumbinary discs as the issue is very likely to happen
    # cfg.set_smoothing_length_density_based_neigh_lim(500)

    # Set the solver config to be the one stored in cfg
    model.set_solver_config(cfg)

    # Print the solver config
    model.get_current_config().print_status()

    # Init the scheduler & fields
    model.init_scheduler(scheduler_split_val, scheduler_merge_val)

    # Set the simulation box size
    model.resize_simulation_box(bmin, bmax)

    # Create the setup

    setup = model.get_setup()
    gen_disc = setup.make_generator_disc_mc(
        part_mass=pmass,
        disc_mass=disc_mass,
        r_in=rin,
        r_out=rout,
        sigma_profile=sigma_profile,
        H_profile=H_profile,
        rot_profile=rot_profile,
        cs_profile=cs_profile,
        random_seed=666,
    )

    # Print the dot graph of the setup
    print(gen_disc.get_dot())

    # Apply the setup
    setup.apply_setup(gen_disc)

    # correct the momentum and barycenter of the disc to 0
    analysis_momentum = shamrock.model_sph.analysisTotalMomentum(model=model)
    total_momentum = analysis_momentum.get_total_momentum()

    if shamrock.sys.world_rank() == 0:
        print(f"disc momentum = {total_momentum}")

    model.apply_momentum_offset((-total_momentum[0], -total_momentum[1], -total_momentum[2]))

    # Correct the barycenter
    analysis_barycenter = shamrock.model_sph.analysisBarycenter(model=model)
    barycenter, disc_mass = analysis_barycenter.get_barycenter()

    if shamrock.sys.world_rank() == 0:
        print(f"disc barycenter = {barycenter}")

    model.apply_position_offset((-barycenter[0], -barycenter[1], -barycenter[2]))

    total_momentum = shamrock.model_sph.analysisTotalMomentum(model=model).get_total_momentum()

    if shamrock.sys.world_rank() == 0:
        print(f"disc momentum after correction = {total_momentum}")

    barycenter, disc_mass = shamrock.model_sph.analysisBarycenter(model=model).get_barycenter()

    if shamrock.sys.world_rank() == 0:
        print(f"disc barycenter after correction = {barycenter}")

    if not np.allclose(total_momentum, 0.0):
        raise RuntimeError("disc momentum is not 0")
    if not np.allclose(barycenter, 0.0):
        raise RuntimeError("disc barycenter is not 0")

    # Run a single step to init the integrator and smoothing length of the particles
    # Here the htolerance is the maximum factor of evolution of the smoothing length in each
    # Smoothing length iterations, increasing it affect the performance negatively but increse the
    # convergence rate of the smoothing length
    # this is why we increase it temporely to 1.3 before lowering it back to 1.1 (default value)
    # Note that both ``change_htolerances`` can be removed and it will work the same but would converge
    # more slowly at the first timestep

    model.change_htolerances(coarse=1.3, fine=1.1)
    model.timestep()
    model.change_htolerances(coarse=1.1, fine=1.1)

----- SPH Solver configuration -----
[
    {
        "artif_viscosity": {
            "alpha_AV": 0.0125,
            "alpha_u": 1.0,
            "av_type": "constant_disc",
            "beta_AV": 2.0
        },
        "boundary_config": {
            "bc_type": "free"
        },
        "cfl_config": {
            "cfl_cour": 0.3,
            "cfl_force": 0.25,
            "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",
            "cs0": 0.31415811727277826,
            "eos_type": "locally_isothermal_lp07",
            "q": 0.5,
            "r0": 1.0
        },
        "epsilon_h": 1e-06,
        "ext_force_config": {
            "force_list": [
                {
                    "Racc": 0.1,
                    "central_mass": 1.0,
                    "force_type": "point_mass"
                }
            ]
        },
        "gpart_mass": 1e-07,
        "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": 20.0,
                "type": "sphere"
            }
        ],
        "particle_reordering_step_freq": 1000,
        "show_neigh_stats": false,
        "smoothing_length_config": {
            "type": "density_based"
        },
        "time_state": {
            "cfl_multiplier": 0.01,
            "dt_sph": 0.0,
            "time": 0.0
        },
        "tree_reduction_level": 3,
        "type_id": "sycl::vec<f64,3>",
        "unit_sys": {
            "unit_current": 1.0,
            "unit_length": 149597870700.0,
            "unit_lumint": 1.0,
            "unit_mass": 1.98847e+30,
            "unit_qte": 1.0,
            "unit_temperature": 1.0,
            "unit_time": 31557600.0
        },
        "use_two_stage_search": true
    }
]
------------------------------------
digraph G {
rankdir=LR;
node_0 [label="GeneratorMCDisc"];
node_2 [label="Simulation"];
node_0 -> node_2;
}

Info: pushing data in scheduler, N = 100000                           [DataInserterUtility][rank=0]
Info: reattributing data ...                                          [DataInserterUtility][rank=0]
Info: reattributing data done in  11.14 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     : 13.38 us   (73.1%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1804.00 ns (0.3%)
   patch tree reduce : 1383.00 ns (0.2%)
   gen split merge   : 902.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 912.00 ns  (0.2%)
   LB compute        : 572.28 us  (97.9%)
   LB move op cnt    : 0
   LB apply          : 3.51 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.60 us    (64.2%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1623.00 ns (0.4%)
   patch tree reduce : 441.00 ns  (0.1%)
   gen split merge   : 400.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 280.00 ns  (0.1%)
   LB compute        : 422.47 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 2.12 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.67 us    (67.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1232.00 ns (0.3%)
   patch tree reduce : 401.00 ns  (0.1%)
   gen split merge   : 370.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 270.00 ns  (0.1%)
   LB compute        : 384.20 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 2.02 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    (62.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1633.00 ns (0.4%)
   patch tree reduce : 431.00 ns  (0.1%)
   gen split merge   : 390.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 270.00 ns  (0.1%)
   LB compute        : 396.48 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 1954.00 ns (0.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.77 us    (63.5%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1462.00 ns (0.3%)
   patch tree reduce : 441.00 ns  (0.1%)
   gen split merge   : 401.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 290.00 ns  (0.1%)
   LB compute        : 425.17 us  (98.2%)
   LB move op cnt    : 0
   LB apply          : 1964.00 ns (0.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.60 us    (64.4%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1413.00 ns (0.4%)
   patch tree reduce : 411.00 ns  (0.1%)
   gen split merge   : 391.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 251.00 ns  (0.1%)
   LB compute        : 384.10 us  (98.1%)
   LB move op cnt    : 0
   LB apply          : 1844.00 ns (0.5%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Final load balancing step 1 of 3                                          [SPH setup][rank=0]
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.67 us    (65.9%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1933.00 ns (0.4%)
   patch tree reduce : 441.00 ns  (0.1%)
   gen split merge   : 420.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 290.00 ns  (0.1%)
   LB compute        : 483.63 us  (98.3%)
   LB move op cnt    : 0
   LB apply          : 1954.00 ns (0.4%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (68.9%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1332.00 ns (0.3%)
   patch tree reduce : 420.00 ns  (0.1%)
   gen split merge   : 380.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 251.00 ns  (0.1%)
   LB compute        : 395.71 us  (98.2%)
   LB move op cnt    : 0
   LB apply          : 1924.00 ns (0.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (66.8%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1663.00 ns (0.4%)
   patch tree reduce : 461.00 ns  (0.1%)
   gen split merge   : 380.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 280.00 ns  (0.1%)
   LB compute        : 434.66 us  (98.2%)
   LB move op cnt    : 0
   LB apply          : 1924.00 ns (0.4%)
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.54 us    (66.7%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1623.00 ns (0.4%)
   patch tree reduce : 451.00 ns  (0.1%)
   gen split merge   : 371.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 260.00 ns  (0.1%)
   LB compute        : 392.07 us  (98.0%)
   LB move op cnt    : 0
   LB apply          : 1884.00 ns (0.5%)
Info: Compute load ...                                                [DataInserterUtility][rank=0]
Info: run scheduler step ...                                          [DataInserterUtility][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.56 us    (65.8%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1614.00 ns (0.4%)
   patch tree reduce : 431.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        : 385.96 us  (98.0%)
   LB move op cnt    : 0
   LB apply          : 2.00 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.54 us    (70.6%)
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 1643.00 ns (0.5%)
   patch tree reduce : 381.00 ns  (0.1%)
   gen split merge   : 411.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 281.00 ns  (0.1%)
   LB compute        : 318.53 us  (97.7%)
   LB move op cnt    : 0
   LB apply          : 1853.00 ns (0.6%)
Info: ---------------------------------------------                   [DataInserterUtility][rank=0]
Info: the setup took : 0.326703954 s                                            [SPH setup][rank=0]
disc momentum = (1.60288854742703e-05, 0.00013766309275351513, 0.0)
disc barycenter = (0.013578565247347756, -0.009662650014606458, 0.000369997836038512)
disc momentum after correction = (-1.8058742435461683e-18, 5.7293308552280875e-18, 0.0)
disc barycenter after correction = (-1.895998549134026e-15, -2.1277467635028025e-16, 1.4653669987499396e-17)
---------------- t = 0, dt = 0 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 4.82 us    (1.1%)
   patch tree reduce : 752.00 ns  (0.2%)
   gen split merge   : 512.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 571.00 ns  (0.1%)
   LB compute        : 419.31 us  (96.6%)
   LB move op cnt    : 0
   LB apply          : 2.49 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.65 us    (68.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.2909764157267044 unconverged cnt = 100000
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.3782693404447157 unconverged cnt = 100000
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.49175014257813043 unconverged cnt = 100000
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.6392751853515696 unconverged cnt = 99999
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.8310577409570404 unconverged cnt = 99994
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9111537095081728 unconverged cnt = 99974
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 99896
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 99423
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 92847
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 46390
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 1529
Warning: smoothing length is not converged, rerunning the iterator ...    [Smoothinglength][rank=0]
     largest h = 0.9808577320660896 unconverged cnt = 6
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-4.206412995699793e-18,1.1630163498921319e-17,0)
    sum a = (-0.0005273027013605336,2.0373113576734206e-05,5.92470562869727e-07)
    sum e = 0.05019233353443156
    sum de = -2.2021610725640874e-05
Info: CFL hydro = 7.455574059655493e-05 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 7.455574059655493e-05 cfl multiplier : 0.01                     [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    3.5663e+04    |      100000 |   2.804e+00   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 0 (tsim/hr)                                             [sph::Model][rank=0]

On the fly analysis

def save_rho_integ(ext, arr_rho, iplot):
    if shamrock.sys.world_rank() == 0:
        metadata = {"extent": [-ext, ext, -ext, ext], "time": model.get_time()}
        np.save(plot_folder + f"rho_integ_{iplot:07}.npy", arr_rho)

        with open(plot_folder + f"rho_integ_{iplot:07}.json", "w") as fp:
            json.dump(metadata, fp)


def save_analysis_data(filename, key, value, ianalysis):
    """Helper to save analysis data to a JSON file."""
    if shamrock.sys.world_rank() == 0:
        filepath = os.path.join(analysis_folder, filename)
        try:
            with open(filepath, "r") as fp:
                data = json.load(fp)
        except (FileNotFoundError, json.JSONDecodeError):
            data = {key: []}
        data[key] = data[key][:ianalysis]
        data[key].append({"t": model.get_time(), key: value})
        with open(filepath, "w") as fp:
            json.dump(data, fp, indent=4)


def analysis(ianalysis):

    ext = rout * 1.5
    nx = 1024
    ny = 1024

    arr_rho2 = model.render_cartesian_column_integ(
        "rho",
        "f64",
        center=(0.0, 0.0, 0.0),
        delta_x=(ext * 2, 0, 0.0),
        delta_y=(0.0, ext * 2, 0.0),
        nx=nx,
        ny=ny,
    )

    save_rho_integ(ext, arr_rho2, ianalysis)

    barycenter, disc_mass = shamrock.model_sph.analysisBarycenter(model=model).get_barycenter()

    total_momentum = shamrock.model_sph.analysisTotalMomentum(model=model).get_total_momentum()

    potential_energy = shamrock.model_sph.analysisEnergyPotential(
        model=model
    ).get_potential_energy()

    kinetic_energy = shamrock.model_sph.analysisEnergyKinetic(model=model).get_kinetic_energy()

    save_analysis_data("barycenter.json", "barycenter", barycenter, ianalysis)
    save_analysis_data("disc_mass.json", "disc_mass", disc_mass, ianalysis)
    save_analysis_data("total_momentum.json", "total_momentum", total_momentum, ianalysis)
    save_analysis_data("potential_energy.json", "potential_energy", potential_energy, ianalysis)
    save_analysis_data("kinetic_energy.json", "kinetic_energy", kinetic_energy, ianalysis)

    sim_time_delta = model.solver_logs_cumulated_step_time()
    scount = model.solver_logs_step_count()

    save_analysis_data("sim_time_delta.json", "sim_time_delta", sim_time_delta, ianalysis)
    save_analysis_data("sim_step_count_delta.json", "sim_step_count_delta", scount, ianalysis)

    model.solver_logs_reset_cumulated_step_time()
    model.solver_logs_reset_step_count()

    part_count = model.get_total_part_count()
    save_analysis_data("part_count.json", "part_count", part_count, ianalysis)

Evolve the simulation

model.solver_logs_reset_cumulated_step_time()
model.solver_logs_reset_step_count()

t_start = model.get_time()

idump = 0
iplot = 0
istop = 0
for ttarg in t_stop:

    if ttarg >= t_start:
        model.evolve_until(ttarg)

        if istop % dump_freq_stop == 0:
            model.do_vtk_dump(get_vtk_dump_name(idump), True)
            model.dump(get_dump_name(idump))

            # dump = model.make_phantom_dump()
            # dump.save_dump(get_ph_dump_name(idump))

            purge_old_dumps()

        if istop % plot_freq_stop == 0:
            analysis(iplot)

    if istop % dump_freq_stop == 0:
        idump += 1

    if istop % plot_freq_stop == 0:
        iplot += 1

    istop += 1

Info: iteration since start : 1                                                       [SPH][rank=0]
Info: time since start : 1159.012839101 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000000.vtk   [VTK Dump][rank=0]
              - took 29.00 ms, bandwidth = 184.18 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000000.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 22.16 us   (57.4%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000000.sham [Shamrock Dump][rank=0]
              - took 51.48 ms, bandwidth = 222.40 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1739.66 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0, dt = 7.455574059655493e-05 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.18 us    (1.2%)
   patch tree reduce : 1663.00 ns (0.3%)
   gen split merge   : 892.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1072.00 ns (0.2%)
   LB compute        : 556.55 us  (96.1%)
   LB move op cnt    : 0
   LB apply          : 3.87 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.83 us    (69.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.931344342444163e-08,1.5189325894056083e-09,4.4172081596407504e-11)
    sum a = (-0.0005272704310153461,2.0235915746764177e-05,5.923962390347488e-07)
    sum e = 0.050192333513713126
    sum de = -2.166125369795201e-05
Info: CFL hydro = 0.0025347941335259246 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0025347941335259246 cfl multiplier : 0.34                     [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8407e+05    |      100000 |   5.433e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 0.4940505611738664 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 7.455574059655493e-05, dt = 0.0025347941335259246 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.01 us    (1.7%)
   patch tree reduce : 1442.00 ns (0.4%)
   gen split merge   : 721.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.3%)
   LB compute        : 384.81 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.49 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.90 us    (70.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-1.3758342357707875e-06,5.2807698657585434e-08,1.5457718223901913e-09)
    sum a = (-0.0005261867586338642,1.5605172255863687e-05,5.894651218246707e-07)
    sum e = 0.05019415085247674
    sum de = -9.47171291465733e-06
Info: CFL hydro = 0.004171090667471787 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004171090667471787 cfl multiplier : 0.56                      [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8567e+05    |      100000 |   5.386e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 16.942775912794662 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.0026093498741224797, dt = 0.004171090667471787 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.16 us    (1.7%)
   patch tree reduce : 1763.00 ns (0.4%)
   gen split merge   : 601.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1403.00 ns (0.3%)
   LB compute        : 391.68 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.45 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.40 us    (69.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.5692334708551866e-06,1.1202929630595974e-07,4.000769401478875e-09)
    sum a = (-0.0005244648563417138,8.122689299741859e-06,5.829270401181027e-07)
    sum e = 0.05019732597389974
    sum de = 1.0421135604408686e-05
Info: CFL hydro = 0.005264873170000351 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005264873170000351 cfl multiplier : 0.7066666666666667        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8463e+05    |      100000 |   5.416e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 27.72450233586342 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.006780440541594266, dt = 0.003219559458405734 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.35 us    (1.9%)
   patch tree reduce : 1342.00 ns (0.4%)
   gen split merge   : 902.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.3%)
   LB compute        : 357.32 us  (94.5%)
   LB move op cnt    : 0
   LB apply          : 3.40 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.40 us    (72.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-5.254188154409121e-06,1.2257572004840824e-07,5.86390220125697e-09)
    sum a = (-0.0005231928181809365,2.458088718159633e-06,5.764137771167129e-07)
    sum e = 0.05019535034735672
    sum de = 2.6089100424208043e-05
Info: CFL hydro = 0.006002219842118081 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006002219842118081 cfl multiplier : 0.8044444444444444        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8189e+05    |      100000 |   5.498e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21.08125907598202 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 5                                                       [SPH][rank=0]
Info: time since start : 1163.06106137 (s)                                            [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1734.93 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.01, dt = 0.006002219842118081 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.71 us    (1.6%)
   patch tree reduce : 1994.00 ns (0.4%)
   gen split merge   : 802.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 472.50 us  (95.5%)
   LB move op cnt    : 0
   LB apply          : 3.95 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.81 us    (69.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.392458767696038e-06,1.2821094973718628e-07,9.313179492787109e-09)
    sum a = (-0.0005209716276955883,-7.873149596371461e-06,5.608535282440876e-07)
    sum e = 0.050203011263207774
    sum de = 5.425545961200534e-05
Info: CFL hydro = 0.006522612771878552 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006522612771878552 cfl multiplier : 0.8696296296296296        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8258e+05    |      100000 |   5.477e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 39.452786694913215 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.01600221984211808, dt = 0.003997780157881919 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.15 us    (1.7%)
   patch tree reduce : 1523.00 ns (0.4%)
   gen split merge   : 701.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 407.84 us  (95.0%)
   LB move op cnt    : 0
   LB apply          : 3.82 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.94 us    (67.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-1.0468522766913522e-05,6.573064669065048e-08,1.1508650582213521e-08)
    sum a = (-0.0005196092921708382,-1.4602182144318539e-05,5.480032879145675e-07)
    sum e = 0.05019747268794535
    sum de = 7.407505750808332e-05
Info: CFL hydro = 0.006884931426892988 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006884931426892988 cfl multiplier : 0.9130864197530864        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8427e+05    |      100000 |   5.427e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 26.51961195885436 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 7                                                       [SPH][rank=0]
Info: time since start : 1165.918101385 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000001.vtk   [VTK Dump][rank=0]
              - took 8.04 ms, bandwidth = 664.57 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000001.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.62 us    (56.0%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000001.sham [Shamrock Dump][rank=0]
              - took 12.07 ms, bandwidth = 948.63 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1739.89 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.02, dt = 0.006884931426892988 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.50 us    (1.4%)
   patch tree reduce : 1693.00 ns (0.3%)
   gen split merge   : 952.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 532.00 us  (96.0%)
   LB move op cnt    : 0
   LB apply          : 3.72 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.88 us    (70.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-1.404327395332422e-05,-4.825497245273027e-08,1.5255929423309498e-08)
    sum a = (-0.00051750661267413,-2.5958067567607326e-05,5.212101323248673e-07)
    sum e = 0.05020717640454584
    sum de = 0.0001061303438083744
Info: CFL hydro = 0.006885244077367825 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006885244077367825 cfl multiplier : 0.9420576131687243        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8219e+05    |      100000 |   5.489e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 45.15716383902913 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.026884931426892987, dt = 0.003115068573107012 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.13 us    (1.0%)
   patch tree reduce : 1433.00 ns (0.2%)
   gen split merge   : 741.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1113.00 ns (0.2%)
   LB compute        : 665.00 us  (96.8%)
   LB move op cnt    : 0
   LB apply          : 3.80 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.75 us    (70.7%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-1.5648104136769944e-05,-1.6820837915498428e-07,1.6787300207027865e-08)
    sum a = (-0.000516661258751475,-3.10100694430528e-05,5.071353251762623e-07)
    sum e = 0.05019663841010546
    sum de = 0.0001222986316945174
Info: CFL hydro = 0.006938554078195827 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006938554078195827 cfl multiplier : 0.9613717421124829        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8348e+05    |      100000 |   5.450e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20.575897487919395 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 9                                                       [SPH][rank=0]
Info: time since start : 1168.8055364440002 (s)                                       [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1738.98 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.03, dt = 0.006938554078195827 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.90 us    (1.5%)
   patch tree reduce : 1473.00 ns (0.3%)
   gen split merge   : 661.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1051.00 ns (0.2%)
   LB compute        : 507.60 us  (95.8%)
   LB move op cnt    : 0
   LB apply          : 3.72 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.62 us    (67.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-1.9231669553002078e-05,-3.9124208910521643e-07,2.0284164091014856e-08)
    sum a = (-0.0005150329837712661,-4.21258862967926e-05,4.7143043775888644e-07)
    sum e = 0.05020870827501933
    sum de = 0.0001545183553556456
Info: CFL hydro = 0.0068941487652737325 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0068941487652737325 cfl multiplier : 0.9742478280749886       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8232e+05    |      100000 |   5.485e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 45.541805778719 (tsim/hr)                               [sph::Model][rank=0]
---------------- t = 0.03693855407819582, dt = 0.0030614459218041776 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.13 us    (1.9%)
   patch tree reduce : 1433.00 ns (0.4%)
   gen split merge   : 771.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.3%)
   LB compute        : 346.60 us  (94.2%)
   LB move op cnt    : 0
   LB apply          : 3.67 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.28 us    (70.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-2.080276624376339e-05,-5.587720600814066e-07,2.1603552736006746e-08)
    sum a = (-0.0005144285222390843,-4.6979743504858516e-05,4.537569701530646e-07)
    sum e = 0.050198005375001215
    sum de = 0.00017058207853391368
Info: CFL hydro = 0.006869163173025313 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006869163173025313 cfl multiplier : 0.9828318853833258        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8379e+05    |      100000 |   5.441e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20.255964378983474 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 11                                                      [SPH][rank=0]
Info: time since start : 1171.668240126 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000002.vtk   [VTK Dump][rank=0]
              - took 7.87 ms, bandwidth = 678.38 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000002.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.45 us    (53.2%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000002.sham [Shamrock Dump][rank=0]
              - took 12.96 ms, bandwidth = 883.72 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1744.85 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.04, dt = 0.006869163173025313 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 8.07 us    (1.5%)
   patch tree reduce : 2.01 us    (0.4%)
   gen split merge   : 711.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.2%)
   LB compute        : 513.86 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.79 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.54 us    (69.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-2.4335534440732185e-05,-8.889134947234992e-07,2.4693430222221777e-08)
    sum a = (-0.0005133364968694082,-5.781335292487208e-05,4.098563796294286e-07)
    sum e = 0.050210214060396434
    sum de = 0.00020251218278697738
Info: CFL hydro = 0.006705856862918982 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006705856862918982 cfl multiplier : 0.9885545902555505        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8038e+05    |      100000 |   5.544e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 44.60504433198048 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.046869163173025315, dt = 0.003130836826974688 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.01 us    (1.4%)
   patch tree reduce : 1794.00 ns (0.3%)
   gen split merge   : 741.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 493.58 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 3.54 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (69.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-2.5938956599539688e-05,-1.1071265845874847e-06,2.582584350948746e-08)
    sum a = (-0.0005129604331325715,-6.27348160892453e-05,3.879004688671905e-07)
    sum e = 0.05020007596077756
    sum de = 0.00021897210799720127
Info: CFL hydro = 0.006665885871089845 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006665885871089845 cfl multiplier : 0.9923697268370336        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8186e+05    |      100000 |   5.499e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20.49705903194765 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 13                                                      [SPH][rank=0]
Info: time since start : 1174.572614963 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1743.29 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.05, dt = 0.006665885871089845 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.87 us    (1.7%)
   patch tree reduce : 1663.00 ns (0.4%)
   gen split merge   : 1012.00 ns (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.2%)
   LB compute        : 441.00 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.67 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.47 us    (70.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-2.935770360608013e-05,-1.5330138578364823e-06,2.8377173577304913e-08)
    sum a = (-0.000512416833263766,-7.32321086973426e-05,3.371409469070733e-07)
    sum e = 0.05021164794318477
    sum de = 0.0002500360143082106
Info: CFL hydro = 0.006971783192670175 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006971783192670175 cfl multiplier : 0.9949131512246892        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8062e+05    |      100000 |   5.536e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 43.344655668066586 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.056665885871089845, dt = 0.0033341141289101525 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.14 us    (1.7%)
   patch tree reduce : 1583.00 ns (0.4%)
   gen split merge   : 771.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 393.12 us  (94.5%)
   LB move op cnt    : 0
   LB apply          : 3.57 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.57 us    (70.8%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.1064348022418073e-05,-1.8121649433769335e-06,2.93320613816959e-08)
    sum a = (-0.0005122746667287697,-7.850359675322324e-05,3.097204216452498e-07)
    sum e = 0.050202887190107674
    sum de = 0.0002674255500393733
Info: CFL hydro = 0.006922497212059127 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006922497212059127 cfl multiplier : 0.9966087674831261        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7959e+05    |      100000 |   5.568e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21.555494269815426 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 15                                                      [SPH][rank=0]
Info: time since start : 1177.4585582240002 (s)                                       [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000003.vtk   [VTK Dump][rank=0]
              - took 8.43 ms, bandwidth = 633.76 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000003.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.76 us    (54.4%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000003.sham [Shamrock Dump][rank=0]
              - took 16.60 ms, bandwidth = 689.61 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1742.21 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.06, dt = 0.006922497212059127 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.62 us    (1.3%)
   patch tree reduce : 1533.00 ns (0.3%)
   gen split merge   : 882.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.2%)
   LB compute        : 493.96 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 3.68 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.67 us    (68.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.461033097493133e-05,-2.3643937444399484e-06,3.143038855670381e-08)
    sum a = (-0.0005122521258236162,-8.955119786091305e-05,2.485376414408002e-07)
    sum e = 0.050215489961617114
    sum de = 0.00029944870490404617
Info: CFL hydro = 0.006751420201152761 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006751420201152761 cfl multiplier : 0.997739178322084         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8044e+05    |      100000 |   5.542e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 44.96739813317225 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.06692249721205912, dt = 0.003077502787940889 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.94 us    (1.5%)
   patch tree reduce : 1483.00 ns (0.3%)
   gen split merge   : 822.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1072.00 ns (0.2%)
   LB compute        : 428.21 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.76 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.52 us    (72.8%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.618671030060484e-05,-2.6782262994559005e-06,3.198349502845023e-08)
    sum a = (-0.0005123568965733009,-9.451709098250894e-05,2.19527019097133e-07)
    sum e = 0.05020532804770916
    sum de = 0.0003158173058747517
Info: CFL hydro = 0.0066640622555288655 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0066640622555288655 cfl multiplier : 0.998492785548056        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8303e+05    |      100000 |   5.464e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 20.278031856803477 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 17                                                      [SPH][rank=0]
Info: time since start : 1180.361217914 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1740.52 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.07, dt = 0.0066640622555288655 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.20 us    (1.6%)
   patch tree reduce : 1452.00 ns (0.3%)
   gen split merge   : 722.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 424.00 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.67 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.02 us    (69.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-3.960124977255662e-05,-3.3157353529350386e-06,3.3401796614911745e-08)
    sum a = (-0.0005128164721434881,-0.0001054379121856383,1.529748715571714e-07)
    sum e = 0.05021764106753366
    sum de = 0.00034661353319604565
Info: CFL hydro = 0.006459611315434467 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006459611315434467 cfl multiplier : 0.9989951903653708        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7981e+05    |      100000 |   5.561e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 43.1383551383335 (tsim/hr)                              [sph::Model][rank=0]
---------------- t = 0.07666406225552887, dt = 0.0033359377444711347 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.11 us    (1.7%)
   patch tree reduce : 1533.00 ns (0.4%)
   gen split merge   : 752.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 403.81 us  (95.0%)
   LB move op cnt    : 0
   LB apply          : 3.52 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.42 us    (71.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-4.131350491807338e-05,-3.7038581800745844e-06,3.369035743567547e-08)
    sum a = (-0.0005131605317201785,-0.00011100137212919653,1.1778498062642918e-07)
    sum e = 0.05020920762201752
    sum de = 0.0003640672984259994
Info: CFL hydro = 0.006387073048379539 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006387073048379539 cfl multiplier : 0.9993301269102473        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8190e+05    |      100000 |   5.497e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 21.84524579190739 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 19                                                      [SPH][rank=0]
Info: time since start : 1183.240108535 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000004.vtk   [VTK Dump][rank=0]
              - took 8.56 ms, bandwidth = 623.78 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000004.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.50 us    (53.0%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000004.sham [Shamrock Dump][rank=0]
              - took 13.19 ms, bandwidth = 868.11 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1742.15 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.08, dt = 0.006387073048379539 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 8.10 us    (2.0%)
   patch tree reduce : 1583.00 ns (0.4%)
   gen split merge   : 741.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.3%)
   LB compute        : 385.03 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.36 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.96 us    (71.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-4.459167260037727e-05,-4.422111730355288e-06,3.4383963068245475e-08)
    sum a = (-0.0005140167040719078,-0.00012188033669855778,4.7012507099540834e-08)
    sum e = 0.05022021486103016
    sum de = 0.0003935871974945001
Info: CFL hydro = 0.006238835997503293 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006238835997503293 cfl multiplier : 0.9995534179401648        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7968e+05    |      100000 |   5.565e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 41.31436549071488 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.08638707304837955, dt = 0.0036129269516204515 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.73 us    (1.8%)
   patch tree reduce : 1493.00 ns (0.4%)
   gen split merge   : 691.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 363.49 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.15 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 29.55 us   (95.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-4.64515116217771e-05,-4.897198854375949e-06,3.4327801343093846e-08)
    sum a = (-0.0005146075529739466,-0.00012817820048114645,5.056437294315569e-09)
    sum e = 0.05021364845838637
    sum de = 0.0004121351170443447
Info: CFL hydro = 0.006176586488189481 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006176586488189481 cfl multiplier : 0.9997022786267765        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8195e+05    |      100000 |   5.496e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 23.665436440776862 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 21                                                      [SPH][rank=0]
Info: time since start : 1186.148079689 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1746.93 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.09, dt = 0.006176586488189481 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.62 us    (1.5%)
   patch tree reduce : 1824.00 ns (0.4%)
   gen split merge   : 991.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 489.45 us  (95.6%)
   LB move op cnt    : 0
   LB apply          : 3.77 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.87 us    (67.8%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-4.9631097027158165e-05,-5.700279456452369e-06,3.4283240757673e-08)
    sum a = (-0.0005157762417440167,-0.00013922828724897038,-6.97918806005775e-08)
    sum e = 0.05022355503536553
    sum de = 0.0004406017969016454
Info: CFL hydro = 0.006014234148924096 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.006014234148924096 cfl multiplier : 0.9998015190845176        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8016e+05    |      100000 |   5.551e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 40.06054091357933 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.09617658648818948, dt = 0.0038234135118105222 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.27 us    (1.8%)
   patch tree reduce : 1463.00 ns (0.4%)
   gen split merge   : 722.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.3%)
   LB compute        : 382.75 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.79 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.69 us    (70.4%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-5.160673213254498e-05,-6.266732679453991e-06,3.3785243983881525e-08)
    sum a = (-0.0005165885898612574,-0.00014626050960692965,-1.1803572942543844e-07)
    sum e = 0.050218467144607246
    sum de = 0.00045989728263834336
Info: CFL hydro = 0.005883355137233783 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005883355137233783 cfl multiplier : 0.9998676793896785        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8186e+05    |      100000 |   5.499e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 25.032147680838026 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 23                                                      [SPH][rank=0]
Info: time since start : 1189.032332487 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000005.vtk   [VTK Dump][rank=0]
              - took 9.12 ms, bandwidth = 585.89 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000005.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.75 us    (54.8%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000005.sham [Shamrock Dump][rank=0]
              - took 13.19 ms, bandwidth = 868.04 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1742.24 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.1, dt = 0.005883355137233783 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.98 us    (1.5%)
   patch tree reduce : 1603.00 ns (0.3%)
   gen split merge   : 651.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.2%)
   LB compute        : 449.17 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 3.76 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.62 us    (68.4%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-5.464755923792454e-05,-7.140678747010443e-06,3.299856977705906e-08)
    sum a = (-0.0005179482673430554,-0.0001574044137531624,-1.9503234697460673e-07)
    sum e = 0.05022705529668905
    sum de = 0.000486916158741424
Info: CFL hydro = 0.0056912242397347546 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0056912242397347546 cfl multiplier : 0.9999117862597856       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7993e+05    |      100000 |   5.558e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 38.10948871677426 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.10588335513723379, dt = 0.004116644862766208 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.07 us    (1.8%)
   patch tree reduce : 1784.00 ns (0.5%)
   gen split merge   : 1122.00 ns (0.3%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.3%)
   LB compute        : 363.51 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.31 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.40 us    (70.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-5.6783768044615633e-05,-7.821438591120946e-06,3.196919164511237e-08)
    sum a = (-0.0005189636147884855,-0.0001654484679289775,-2.508247767821456e-07)
    sum e = 0.05022398567597083
    sum de = 0.0005071865838681318
Info: CFL hydro = 0.005572057795169468 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005572057795169468 cfl multiplier : 0.9999411908398571        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8166e+05    |      100000 |   5.505e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 26.922424830917098 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 25                                                      [SPH][rank=0]
Info: time since start : 1191.94069383 (s)                                            [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1749.16 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.11, dt = 0.005572057795169468 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.79 us    (1.4%)
   patch tree reduce : 1984.00 ns (0.4%)
   gen split merge   : 781.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1312.00 ns (0.2%)
   LB compute        : 516.55 us  (95.7%)
   LB move op cnt    : 0
   LB apply          : 3.71 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.04 us    (71.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-5.967755321222929e-05,-8.759884273699336e-06,3.045674268264413e-08)
    sum a = (-0.0005203962055905046,-0.000176689414195831,-3.2876027707396e-07)
    sum e = 0.05023098329068514
    sum de = 0.0005326607127616463
Info: CFL hydro = 0.0054194647391847175 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0054194647391847175 cfl multiplier : 0.999960793893238        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7991e+05    |      100000 |   5.558e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 36.08844456203835 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.11557205779516946, dt = 0.004427942204830532 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.90 us    (1.7%)
   patch tree reduce : 1533.00 ns (0.4%)
   gen split merge   : 691.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.2%)
   LB compute        : 392.95 us  (94.9%)
   LB move op cnt    : 0
   LB apply          : 3.32 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.73 us    (70.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-6.198582877356893e-05,-9.573572389102761e-06,2.878388062056022e-08)
    sum a = (-0.0005215629192974633,-0.00018592640781379082,-3.9258869501171656e-07)
    sum e = 0.050230074825999865
    sum de = 0.0005538421357854508
Info: CFL hydro = 0.005310660532810927 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005310660532810927 cfl multiplier : 0.999973862595492         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8123e+05    |      100000 |   5.518e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 28.888430629095783 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 27                                                      [SPH][rank=0]
Info: time since start : 1194.829436575 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000006.vtk   [VTK Dump][rank=0]
              - took 8.46 ms, bandwidth = 631.03 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000006.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.17 us    (51.8%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000006.sham [Shamrock Dump][rank=0]
              - took 13.10 ms, bandwidth = 873.69 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1745.21 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.12, dt = 0.005310660532810927 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.91 us    (1.5%)
   patch tree reduce : 1543.00 ns (0.3%)
   gen split merge   : 741.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 435.37 us  (95.4%)
   LB move op cnt    : 0
   LB apply          : 3.59 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.87 us    (70.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-6.475825545489097e-05,-1.0581414862022281e-05,2.6557661059507012e-08)
    sum a = (-0.0005229675395632233,-0.0001973838394392645,-4.7125465869209574e-07)
    sum e = 0.050235587796690614
    sum de = 0.0005779338206053286
Info: CFL hydro = 0.005584624906323711 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005584624906323711 cfl multiplier : 0.9999825750636614        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8024e+05    |      100000 |   5.548e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 34.45990733135671 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.12531066053281092, dt = 0.00468933946718908 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.14 us    (0.5%)
   patch tree reduce : 1543.00 ns (0.1%)
   gen split merge   : 811.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.1%)
   LB compute        : 1466.82 us (98.5%)
   LB move op cnt    : 0
   LB apply          : 3.91 us    (0.3%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 3.02 us    (71.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-6.721435750892302e-05,-1.1537437955461454e-05,2.4138903875107607e-08)
    sum a = (-0.0005241883621506704,-0.00020785926441693745,-5.425330117746971e-07)
    sum e = 0.05023655135287883
    sum de = 0.0005996837844120799
Info: CFL hydro = 0.005380472025086704 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005380472025086704 cfl multiplier : 0.9999883833757742        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7837e+05    |      100000 |   5.606e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.112381455807245 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 29                                                      [SPH][rank=0]
Info: time since start : 1197.749657865 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1748.61 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.13, dt = 0.005380472025086704 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.16 us    (1.7%)
   patch tree reduce : 1573.00 ns (0.4%)
   gen split merge   : 1011.00 ns (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 395.82 us  (94.8%)
   LB move op cnt    : 0
   LB apply          : 3.66 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.54 us    (71.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-7.003760075312107e-05,-1.268038032470796e-05,2.1052695985433607e-08)
    sum a = (-0.0005255338692504452,-0.0002203134562204754,-6.26301056537293e-07)
    sum e = 0.0502418798148296
    sum de = 0.0006235979337532802
Info: CFL hydro = 0.00548132660073205 sink sink = inf                                 [SPH][rank=0]
Info: cfl dt = 0.00548132660073205 cfl multiplier : 0.9999922555838495         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7854e+05    |      100000 |   5.601e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 34.58298844193254 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.1353804720250867, dt = 0.004619527974913301 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.57 us    (1.6%)
   patch tree reduce : 1482.00 ns (0.4%)
   gen split merge   : 752.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 397.65 us  (94.9%)
   LB move op cnt    : 0
   LB apply          : 4.08 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.39 us    (71.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-7.246893889555169e-05,-1.3731629214268079e-05,1.7934124923318088e-08)
    sum a = (-0.0005266148831053623,-0.00023138908652180307,-6.998098053863577e-07)
    sum e = 0.05024260201566397
    sum de = 0.0006446058804304744
Info: CFL hydro = 0.0053629461079758136 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0053629461079758136 cfl multiplier : 0.9999948370558996       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8085e+05    |      100000 |   5.529e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.07665409362497 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 31                                                      [SPH][rank=0]
Info: time since start : 1200.6438020770001 (s)                                       [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000007.vtk   [VTK Dump][rank=0]
              - took 8.52 ms, bandwidth = 627.03 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000007.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.28 us    (53.1%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000007.sham [Shamrock Dump][rank=0]
              - took 13.18 ms, bandwidth = 868.98 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1745.64 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.14, dt = 0.0053629461079758136 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.00 us    (1.7%)
   patch tree reduce : 1612.00 ns (0.4%)
   gen split merge   : 761.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.2%)
   LB compute        : 397.80 us  (95.0%)
   LB move op cnt    : 0
   LB apply          : 3.67 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.43 us    (71.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-7.529564302017087e-05,-1.4998138507261225e-05,1.4011294790346048e-08)
    sum a = (-0.0005277479730181583,-0.0002447077048673572,-7.868784990340767e-07)
    sum e = 0.050248298076401646
    sum de = 0.0006677482085564423
Info: CFL hydro = 0.005234197396630076 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005234197396630076 cfl multiplier : 0.9999965580372665        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7950e+05    |      100000 |   5.571e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 34.65478465611507 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.14536294610797582, dt = 0.004637053892024179 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.25 us    (2.0%)
   patch tree reduce : 1362.00 ns (0.4%)
   gen split merge   : 812.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.3%)
   LB compute        : 348.27 us  (94.2%)
   LB move op cnt    : 0
   LB apply          : 3.46 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.37 us    (71.4%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-7.77458771625331e-05,-1.6168574838737725e-05,1.0129024427990345e-08)
    sum a = (-0.0005285925268273856,-0.00025663122837664117,-8.635611551153989e-07)
    sum e = 0.050249333483235015
    sum de = 0.0006882053030857375
Info: CFL hydro = 0.005130842104743839 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005130842104743839 cfl multiplier : 0.9999977053581777        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7900e+05    |      100000 |   5.587e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 29.881398477815665 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 33                                                      [SPH][rank=0]
Info: time since start : 1203.564929158 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1743.65 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.15, dt = 0.005130842104743839 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.64 us    (1.8%)
   patch tree reduce : 1583.00 ns (0.4%)
   gen split merge   : 771.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1013.00 ns (0.2%)
   LB compute        : 395.30 us  (94.8%)
   LB move op cnt    : 0
   LB apply          : 3.55 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.67 us    (70.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.045996007619687e-05,-1.7512954161226085e-05,5.5204376888784215e-09)
    sum a = (-0.0005293450690473402,-0.0002702764400576835,-9.498265930836721e-07)
    sum e = 0.050254341952120904
    sum de = 0.0007097328048546787
Info: CFL hydro = 0.005023588607798861 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005023588607798861 cfl multiplier : 0.9999984702387851        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8020e+05    |      100000 |   5.549e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 33.28563447610275 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.15513084210474384, dt = 0.00486915789525616 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.17 us    (1.9%)
   patch tree reduce : 1643.00 ns (0.4%)
   gen split merge   : 802.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1102.00 ns (0.3%)
   LB compute        : 355.79 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.38 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.31 us    (70.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.303935538611614e-05,-1.8863978536555705e-05,6.742748633911332e-10)
    sum a = (-0.0005298485886218948,-0.00028367102884085536,-1.0329747059671885e-06)
    sum e = 0.050257094261824904
    sum de = 0.0007301559541244692
Info: CFL hydro = 0.00489637326058909 sink sink = inf                                 [SPH][rank=0]
Info: cfl dt = 0.00489637326058909 cfl multiplier : 0.99999898015919           [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7984e+05    |      100000 |   5.560e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.524280388209657 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 35                                                      [SPH][rank=0]
Info: time since start : 1206.4521504000002 (s)                                       [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000008.vtk   [VTK Dump][rank=0]
              - took 8.50 ms, bandwidth = 628.03 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000008.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.62 us    (55.9%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000008.sham [Shamrock Dump][rank=0]
              - took 12.76 ms, bandwidth = 897.14 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1747.15 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.16, dt = 0.00489637326058909 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.53 us    (1.6%)
   patch tree reduce : 1463.00 ns (0.3%)
   gen split merge   : 831.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1071.00 ns (0.2%)
   LB compute        : 457.50 us  (95.4%)
   LB move op cnt    : 0
   LB apply          : 3.91 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.55 us    (70.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.563491770576288e-05,-2.0285547960831494e-05,-4.585985510931123e-09)
    sum a = (-0.0005301125706722303,-0.00029758330068667415,-1.1177611603867345e-06)
    sum e = 0.05026080857119104
    sum de = 0.00075009920493455
Info: CFL hydro = 0.004777000507066701 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004777000507066701 cfl multiplier : 0.9999993201061267        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7837e+05    |      100000 |   5.606e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.442045432897707 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.1648963732605891, dt = 0.004777000507066701 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.96 us    (1.8%)
   patch tree reduce : 1483.00 ns (0.4%)
   gen split merge   : 742.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.3%)
   LB compute        : 355.57 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.58 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.29 us    (69.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.816791200199055e-05,-2.1741163377033657e-05,-1.0133104205018821e-08)
    sum a = (-0.0005301017251845929,-0.0003115869634190173,-1.201535120954125e-06)
    sum e = 0.050264110969750565
    sum de = 0.00076919037432409
Info: CFL hydro = 0.004669682099885259 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004669682099885259 cfl multiplier : 0.9999995467374179        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7972e+05    |      100000 |   5.564e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.907071826059926 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.1696733737676558, dt = 0.0003266262323442237 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.25 us    (1.9%)
   patch tree reduce : 1563.00 ns (0.4%)
   gen split merge   : 702.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.3%)
   LB compute        : 366.22 us  (94.5%)
   LB move op cnt    : 0
   LB apply          : 3.71 us    (1.0%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.28 us    (69.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-8.834103122679975e-05,-2.1876383604935704e-05,-1.0725651220660872e-08)
    sum a = (-0.0005300905500711502,-0.000312559827779536,-1.2072986437373787e-06)
    sum e = 0.050257691021018513
    sum de = 0.0007715714622936871
Info: CFL hydro = 0.004664245044635256 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004664245044635256 cfl multiplier : 0.9999996978249452        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8443e+05    |      100000 |   5.422e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 2.168643643773828 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 38                                                      [SPH][rank=0]
Info: time since start : 1209.918574601 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1747.18 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.17, dt = 0.004664245044635256 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.68 us    (1.7%)
   patch tree reduce : 1903.00 ns (0.4%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 441.45 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 3.60 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.78 us    (70.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.08135016231346e-05,-2.3334398114318734e-05,-1.6357729195975144e-08)
    sum a = (-0.000529770046196032,-0.0003266723940748216,-1.2900902045387758e-06)
    sum e = 0.05026769825440093
    sum de = 0.0007886465478107171
Info: CFL hydro = 0.004683907411468476 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004683907411468476 cfl multiplier : 0.9999997985499635        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8008e+05    |      100000 |   5.553e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.237836623910628 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.17466424504463526, dt = 0.004683907411468476 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.47 us    (2.0%)
   patch tree reduce : 1443.00 ns (0.4%)
   gen split merge   : 681.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.3%)
   LB compute        : 352.17 us  (94.3%)
   LB move op cnt    : 0
   LB apply          : 3.32 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.30 us    (69.7%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.329414801457522e-05,-2.489741359575273e-05,-2.2593472330081627e-08)
    sum a = (-0.0005291257584395418,-0.0003412518142396393,-1.3740850337655386e-06)
    sum e = 0.050271499836533104
    sum de = 0.000806298798027388
Info: CFL hydro = 0.004616914709026796 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004616914709026796 cfl multiplier : 0.9999998656999756        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7841e+05    |      100000 |   5.605e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.083832045756562 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.17934815245610375, dt = 0.0006518475438962446 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.87 us    (1.7%)
   patch tree reduce : 1493.00 ns (0.4%)
   gen split merge   : 712.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1102.00 ns (0.3%)
   LB compute        : 376.50 us  (94.8%)
   LB move op cnt    : 0
   LB apply          : 3.40 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.37 us    (70.4%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.363754844853025e-05,-2.515400207980132e-05,-2.3685878286018067e-08)
    sum a = (-0.0005290088590044117,-0.0003433126042735348,-1.3858370344045017e-06)
    sum e = 0.05026569009315309
    sum de = 0.0008098245778131402
Info: CFL hydro = 0.004610325166581742 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004610325166581742 cfl multiplier : 0.9999999104666504        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.6071e+05    |      100000 |   6.223e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 3.771214675974877 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 41                                                      [SPH][rank=0]
Info: time since start : 1213.4377137440001 (s)                                       [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000009.vtk   [VTK Dump][rank=0]
              - took 8.92 ms, bandwidth = 598.71 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000009.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.61 us    (55.2%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000009.sham [Shamrock Dump][rank=0]
              - took 13.21 ms, bandwidth = 866.47 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1751.95 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.18, dt = 0.004610325166581742 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.21 us    (1.6%)
   patch tree reduce : 1202.00 ns (0.3%)
   gen split merge   : 842.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1142.00 ns (0.3%)
   LB compute        : 435.21 us  (95.4%)
   LB move op cnt    : 0
   LB apply          : 3.25 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.56 us    (71.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.607641320423566e-05,-2.673745647973603e-05,-3.007886789888913e-08)
    sum a = (-0.0005279754248150916,-0.0003581132490709239,-1.46938288062994e-06)
    sum e = 0.0502756058974416
    sum de = 0.0008255909627032738
Info: CFL hydro = 0.004784436480133553 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004784436480133553 cfl multiplier : 0.9999999403111003        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8049e+05    |      100000 |   5.540e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 29.956481370164706 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.18461032516658174, dt = 0.004784436480133553 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.93 us    (1.9%)
   patch tree reduce : 1493.00 ns (0.4%)
   gen split merge   : 691.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 344.24 us  (94.3%)
   LB move op cnt    : 0
   LB apply          : 3.31 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.77 us    (71.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.860009585350944e-05,-2.848494446521738e-05,-3.7301623714961996e-08)
    sum a = (-0.000526498218522436,-0.00037388102496961595,-1.5568070086083669e-06)
    sum e = 0.050280091921165886
    sum de = 0.0008423084775995933
Info: CFL hydro = 0.0050752443852714315 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0050752443852714315 cfl multiplier : 0.9999999602074002       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7996e+05    |      100000 |   5.557e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.9964831431036 (tsim/hr)                              [sph::Model][rank=0]
---------------- t = 0.1893947616467153, dt = 0.000605238353284715 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.95 us    (1.7%)
   patch tree reduce : 1472.00 ns (0.4%)
   gen split merge   : 832.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 379.29 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.78 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.50 us    (71.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-9.891521896845693e-05,-2.874895156219502e-05,-3.845300061881122e-08)
    sum a = (-0.0005262804608269346,-0.00037590451766720244,-1.5679133739083733e-06)
    sum e = 0.05027394853955663
    sum de = 0.0008456076116870714
Info: CFL hydro = 0.0053295410543162305 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0053295410543162305 cfl multiplier : 0.9999999734716001       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8399e+05    |      100000 |   5.435e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 4.00883866088997 (tsim/hr)                              [sph::Model][rank=0]
Info: iteration since start : 44                                                      [SPH][rank=0]
Info: time since start : 1216.904318108 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1761.93 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.19, dt = 0.0053295410543162305 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.32 us    (1.6%)
   patch tree reduce : 1302.00 ns (0.3%)
   gen split merge   : 851.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1083.00 ns (0.2%)
   LB compute        : 445.40 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 4.35 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.42 us    (70.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00010171998639286118,-3.075296246930687e-05,-4.68126203137864e-08)
    sum a = (-0.000524038734509224,-0.0003940048152523129,-1.666155952225764e-06)
    sum e = 0.050286798160610206
    sum de = 0.0008618809194869914
Info: CFL hydro = 0.005274488568402859 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005274488568402859 cfl multiplier : 0.9999999823144           [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7927e+05    |      100000 |   5.578e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 34.39540396588711 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.19532954105431624, dt = 0.00467045894568377 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.70 us    (1.8%)
   patch tree reduce : 1723.00 ns (0.5%)
   gen split merge   : 671.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.3%)
   LB compute        : 354.73 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.55 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.33 us    (68.8%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00010416151410211963,-3.2641378922887386e-05,-5.485612721298141e-08)
    sum a = (-0.0005215730334179641,-0.00041026750577451474,-1.752812610994611e-06)
    sum e = 0.050288924427520566
    sum de = 0.0008772756706436187
Info: CFL hydro = 0.005212350291676712 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005212350291676712 cfl multiplier : 0.9999999882095999        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7965e+05    |      100000 |   5.567e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.204898078405563 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 46                                                      [SPH][rank=0]
Info: time since start : 1219.813887196 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000010.vtk   [VTK Dump][rank=0]
              - took 8.50 ms, bandwidth = 628.39 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000010.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.77 us    (55.4%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000010.sham [Shamrock Dump][rank=0]
              - took 13.33 ms, bandwidth = 858.71 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1756.16 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.2, dt = 0.005212350291676712 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.20 us    (1.5%)
   patch tree reduce : 1703.00 ns (0.4%)
   gen split merge   : 872.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1182.00 ns (0.3%)
   LB compute        : 445.34 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 3.44 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.49 us    (72.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00010687437747712671,-3.481781399048031e-05,-6.419476372073181e-08)
    sum a = (-0.000518223325543522,-0.0004288440727620283,-1.8500445680602435e-06)
    sum e = 0.050295142869922864
    sum de = 0.0008927606220750189
Info: CFL hydro = 0.005074074704324341 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005074074704324341 cfl multiplier : 0.9999999921397332        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7803e+05    |      100000 |   5.617e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 33.40583836859793 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.20521235029167673, dt = 0.004787649708323266 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.95 us    (1.6%)
   patch tree reduce : 1523.00 ns (0.4%)
   gen split merge   : 752.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 412.39 us  (95.1%)
   LB move op cnt    : 0
   LB apply          : 3.69 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.54 us    (71.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00010934671930510366,-3.6919382977533404e-05,-7.330553256726832e-08)
    sum a = (-0.0005145528324359958,-0.00044628487037586966,-1.939732883117478e-06)
    sum e = 0.050298206362027455
    sum de = 0.0009068330404480632
Info: CFL hydro = 0.005231260896779528 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005231260896779528 cfl multiplier : 0.9999999947598223        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7929e+05    |      100000 |   5.578e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.90171825870456 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 48                                                      [SPH][rank=0]
Info: time since start : 1222.750188431 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1756.56 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.21, dt = 0.005231260896779528 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.59 us    (1.9%)
   patch tree reduce : 1432.00 ns (0.4%)
   gen split merge   : 1012.00 ns (0.2%)
   split / merge op  : 0/0
   apply split merge : 1031.00 ns (0.3%)
   LB compute        : 386.78 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.61 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.34 us    (69.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00011202969289912391,-3.9295765783559796e-05,-8.366747946663123e-08)
    sum a = (-0.0005098470822267781,-0.0004657329897699652,-2.0380351330775476e-06)
    sum e = 0.05030432767254323
    sum de = 0.0009206254119190595
Info: CFL hydro = 0.005119973485917868 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005119973485917868 cfl multiplier : 0.9999999965065482        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7919e+05    |      100000 |   5.581e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 33.74683193645776 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.21523126089677952, dt = 0.004768739103220482 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.31 us    (1.8%)
   patch tree reduce : 1402.00 ns (0.3%)
   gen split merge   : 1102.00 ns (0.3%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.2%)
   LB compute        : 385.71 us  (94.8%)
   LB move op cnt    : 0
   LB apply          : 3.42 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.41 us    (71.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001144487121132724,-4.1567593996788595e-05,-9.364345965761066e-08)
    sum a = (-0.0005048876215611898,-0.0004837928154204137,-2.127812313368435e-06)
    sum e = 0.05030738648243619
    sum de = 0.0009329664576468403
Info: CFL hydro = 0.004939638822942556 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004939638822942556 cfl multiplier : 0.9999999976710322        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7986e+05    |      100000 |   5.560e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.877227984311652 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 50                                                      [SPH][rank=0]
Info: time since start : 1225.653937067 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000011.vtk   [VTK Dump][rank=0]
              - took 8.47 ms, bandwidth = 630.24 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000011.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.18 us    (53.0%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000011.sham [Shamrock Dump][rank=0]
              - took 13.34 ms, bandwidth = 858.52 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1755.41 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.22, dt = 0.004939638822942556 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.79 us    (1.5%)
   patch tree reduce : 1383.00 ns (0.3%)
   gen split merge   : 791.00 ns  (0.1%)
   split / merge op  : 0/0
   apply split merge : 1031.00 ns (0.2%)
   LB compute        : 515.44 us  (95.9%)
   LB move op cnt    : 0
   LB apply          : 3.36 us    (0.6%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.58 us    (70.7%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00011693084942296017,-4.400041706848475e-05,-1.043681459437834e-07)
    sum a = (-0.000499037787802607,-0.0005028066141163477,-2.2208712750084037e-06)
    sum e = 0.0503125361536966
    sum de = 0.0009443725995870475
Info: CFL hydro = 0.004770239360440624 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004770239360440624 cfl multiplier : 0.9999999984473549        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7916e+05    |      100000 |   5.582e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.85956147487819 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.22493963882294254, dt = 0.004770239360440624 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.90 us    (1.7%)
   patch tree reduce : 1443.00 ns (0.3%)
   gen split merge   : 882.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 394.45 us  (95.0%)
   LB move op cnt    : 0
   LB apply          : 3.34 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.40 us    (70.8%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00011929693108771138,-4.6445885618936984e-05,-1.1519207234416746e-07)
    sum a = (-0.0004926657396273997,-0.0005214360347042113,-2.3107015513533266e-06)
    sum e = 0.050316591012126585
    sum de = 0.0009547045476582186
Info: CFL hydro = 0.004874531507280875 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004874531507280875 cfl multiplier : 0.9999999989649032        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7957e+05    |      100000 |   5.569e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.836592276796722 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.22970987818338318, dt = 0.000290121816616834 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.26 us    (2.0%)
   patch tree reduce : 1523.00 ns (0.4%)
   gen split merge   : 691.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.3%)
   LB compute        : 346.67 us  (94.3%)
   LB move op cnt    : 0
   LB apply          : 3.21 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.28 us    (68.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00011942466606956388,-4.664159898625191e-05,-1.160767132358929e-07)
    sum a = (-0.0004922547359013618,-0.0005225767299025201,-2.316161001088324e-06)
    sum e = 0.05031007517348761
    sum de = 0.0009566873484344456
Info: CFL hydro = 0.00486691499197302 sink sink = inf                                 [SPH][rank=0]
Info: cfl dt = 0.00486691499197302 cfl multiplier : 0.9999999993099354         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8320e+05    |      100000 |   5.458e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 1.9134621513759063 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 53                                                      [SPH][rank=0]
Info: time since start : 1229.136257594 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1753.81 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.23, dt = 0.00486691499197302 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.52 us    (1.5%)
   patch tree reduce : 1413.00 ns (0.3%)
   gen split merge   : 841.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1052.00 ns (0.2%)
   LB compute        : 473.06 us  (95.5%)
   LB move op cnt    : 0
   LB apply          : 3.79 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.50 us    (72.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00012182036840301793,-4.918510097775189e-05,-1.2735006388865758e-07)
    sum a = (-0.00048493945379644653,-0.0005418368523772717,-2.4076691536748227e-06)
    sum e = 0.050321837995574464
    sum de = 0.0009646291822538901
Info: CFL hydro = 0.005072567511612904 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005072567511612904 cfl multiplier : 0.999999999539957         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7927e+05    |      100000 |   5.578e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.409868097857075 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.23486691499197304, dt = 0.005072567511612904 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.79 us    (1.8%)
   patch tree reduce : 1462.00 ns (0.4%)
   gen split merge   : 741.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.3%)
   LB compute        : 347.16 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.11 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.35 us    (70.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00012426245509337367,-5.1980473681131225e-05,-1.397858094161518e-07)
    sum a = (-0.0004764523320227156,-0.0005621289111145361,-2.5028014342402053e-06)
    sum e = 0.05032738700007817
    sum de = 0.0009731293035311213
Info: CFL hydro = 0.005060858059484368 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005060858059484368 cfl multiplier : 0.9999999996933046        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7919e+05    |      100000 |   5.581e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 32.72233497533821 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.23993948250358593, dt = 6.051749641405868e-05 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.00 us    (1.9%)
   patch tree reduce : 1513.00 ns (0.4%)
   gen split merge   : 682.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 357.15 us  (94.5%)
   LB move op cnt    : 0
   LB apply          : 3.24 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.37 us    (69.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00012426976304657938,-5.206595873443998e-05,-1.4017855515082864e-07)
    sum a = (-0.00047634567155145847,-0.0005623721419857693,-2.5039343893133213e-06)
    sum e = 0.05031970621117299
    sum de = 0.0009748566608152707
Info: CFL hydro = 0.0050603195354948696 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0050603195354948696 cfl multiplier : 0.9999999997955363       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8552e+05    |      100000 |   5.390e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 0.40418005625960685 (tsim/hr)                           [sph::Model][rank=0]
Info: iteration since start : 56                                                      [SPH][rank=0]
Info: time since start : 1232.579170293 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000012.vtk   [VTK Dump][rank=0]
              - took 8.40 ms, bandwidth = 635.44 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000012.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.63 us    (53.8%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000012.sham [Shamrock Dump][rank=0]
              - took 13.28 ms, bandwidth = 861.94 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1758.37 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.24, dt = 0.0050603195354948696 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.81 us    (1.5%)
   patch tree reduce : 1743.00 ns (0.4%)
   gen split merge   : 751.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.2%)
   LB compute        : 436.01 us  (95.3%)
   LB move op cnt    : 0
   LB apply          : 3.33 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.73 us    (67.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00012668022112656752,-5.491174883060893e-05,-1.5284929753846624e-07)
    sum a = (-0.0004669553719186044,-0.0005827963794899516,-2.598496703759751e-06)
    sum e = 0.05033237346966689
    sum de = 0.0009806990918487964
Info: CFL hydro = 0.005090958016911512 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005090958016911512 cfl multiplier : 0.9999999998636909        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7878e+05    |      100000 |   5.593e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 32.56877810596079 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.24506031953549487, dt = 0.004939680464505131 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.15 us    (1.7%)
   patch tree reduce : 1302.00 ns (0.3%)
   gen split merge   : 942.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1012.00 ns (0.2%)
   LB compute        : 390.79 us  (94.9%)
   LB move op cnt    : 0
   LB apply          : 3.47 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.65 us    (72.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001289630724966956,-5.784225330517759e-05,-1.6592429870666812e-07)
    sum a = (-0.0004568772734069038,-0.0006028617950678783,-2.6904053420994834e-06)
    sum e = 0.05033689042192584
    sum de = 0.0009869431225157247
Info: CFL hydro = 0.0048805040139376175 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0048805040139376175 cfl multiplier : 0.9999999999091272       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7890e+05    |      100000 |   5.590e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.81300421775962 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 58                                                      [SPH][rank=0]
Info: time since start : 1235.516428848 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1754.40 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.25, dt = 0.0048805040139376175 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.54 us    (1.7%)
   patch tree reduce : 1673.00 ns (0.4%)
   gen split merge   : 872.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 428.06 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.36 us    (0.7%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.77 us    (67.7%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001311679725702607,-6.083408108652977e-05,-1.792818324305621e-07)
    sum a = (-0.0004460064667504979,-0.0006227661618434107,-2.7807636995723828e-06)
    sum e = 0.05034156374265078
    sum de = 0.0009919302374818233
Info: CFL hydro = 0.004696757787427285 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004696757787427285 cfl multiplier : 0.9999999999394182        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7703e+05    |      100000 |   5.649e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.103038274184083 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.25488050401393764, dt = 0.004696757787427285 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.38 us    (1.5%)
   patch tree reduce : 1553.00 ns (0.4%)
   gen split merge   : 1112.00 ns (0.3%)
   split / merge op  : 0/0
   apply split merge : 1002.00 ns (0.2%)
   LB compute        : 403.81 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.64 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.50 us    (71.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00013323622940846497,-6.380763457788584e-05,-1.9256290315469696e-07)
    sum a = (-0.00043466366148021214,-0.0006419517757687679,-2.8672489628770702e-06)
    sum e = 0.05034571723478276
    sum de = 0.0009957389656806818
Info: CFL hydro = 0.0048609649644203974 sink sink = inf                               [SPH][rank=0]
Info: cfl dt = 0.0048609649644203974 cfl multiplier : 0.9999999999596122       [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7930e+05    |      100000 |   5.577e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.31693559982341 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.2595772618013649, dt = 0.00042273819863508644 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.45 us    (2.0%)
   patch tree reduce : 1794.00 ns (0.5%)
   gen split merge   : 751.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1021.00 ns (0.3%)
   LB compute        : 358.25 us  (94.3%)
   LB move op cnt    : 0
   LB apply          : 3.29 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.33 us    (71.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00013339334113722824,-6.412406720598971e-05,-1.9397809898325755e-07)
    sum a = (-0.00043359994634584073,-0.0006436783969457973,-2.8750083535274696e-06)
    sum e = 0.050339501731684784
    sum de = 0.0009974737328381064
Info: CFL hydro = 0.004852744156919284 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004852744156919284 cfl multiplier : 0.9999999999730749        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8397e+05    |      100000 |   5.436e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 2.7997635157326415 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 61                                                      [SPH][rank=0]
Info: time since start : 1238.9718627510001 (s)                                       [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000013.vtk   [VTK Dump][rank=0]
              - took 8.70 ms, bandwidth = 613.70 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000013.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.02 us    (52.6%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000013.sham [Shamrock Dump][rank=0]
              - took 13.27 ms, bandwidth = 863.03 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1759.39 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.26, dt = 0.004852744156919284 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.75 us    (1.7%)
   patch tree reduce : 1442.00 ns (0.3%)
   gen split merge   : 892.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1073.00 ns (0.2%)
   LB compute        : 433.93 us  (95.2%)
   LB move op cnt    : 0
   LB apply          : 3.55 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.48 us    (70.3%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00013549726590679232,-6.72480387400707e-05,-2.0793141906734957e-07)
    sum a = (-0.0004208662051369969,-0.0006634821940662463,-2.9637882860597594e-06)
    sum e = 0.05035145025681746
    sum de = 0.0009984375561984205
Info: CFL hydro = 0.004710359789939905 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004710359789939905 cfl multiplier : 0.9999999999820499        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7885e+05    |      100000 |   5.591e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.24511941845655 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.2648527441569193, dt = 0.004710359789939905 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.49 us    (2.0%)
   patch tree reduce : 1503.00 ns (0.4%)
   gen split merge   : 662.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.3%)
   LB compute        : 349.03 us  (94.2%)
   LB move op cnt    : 0
   LB apply          : 3.39 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.44 us    (70.2%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001374488003622904,-7.042132996871808e-05,-2.2210734138532356e-07)
    sum a = (-0.00040758122844273295,-0.0006826378624643208,-3.049409389328264e-06)
    sum e = 0.050355761760102286
    sum de = 0.0009998201459939385
Info: CFL hydro = 0.004782659067463012 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004782659067463012 cfl multiplier : 0.9999999999880332        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7909e+05    |      100000 |   5.584e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.36911499446493 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.2695631039468592, dt = 0.000436896053140845 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.95 us    (1.8%)
   patch tree reduce : 1373.00 ns (0.4%)
   gen split merge   : 751.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 369.37 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.39 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.81 us    (69.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00013759558248230792,-7.076468680164045e-05,-2.236412694129356e-07)
    sum a = (-0.00040630265906255246,-0.0006844093579805253,-3.057321321746986e-06)
    sum e = 0.0503494980948489
    sum de = 0.0010013747525659108
Info: CFL hydro = 0.004773907706402027 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004773907706402027 cfl multiplier : 0.999999999992022         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8360e+05    |      100000 |   5.447e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 2.8877688461402675 (tsim/hr)                            [sph::Model][rank=0]
Info: iteration since start : 64                                                      [SPH][rank=0]
Info: time since start : 1242.455266178 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1758.25 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.27, dt = 0.004773907706402027 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 33.91 us   (7.5%)
   patch tree reduce : 1443.00 ns (0.3%)
   gen split merge   : 732.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 403.49 us  (89.5%)
   LB move op cnt    : 0
   LB apply          : 3.43 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.73 us    (70.1%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001395349545765844,-7.40323808897341e-05,-2.3823836757779423e-07)
    sum a = (-0.0003918044265981054,-0.0007036925676263797,-3.1434520429190506e-06)
    sum e = 0.0503611762301031
    sum de = 0.0009998430987971717
Info: CFL hydro = 0.004632159271225786 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004632159271225786 cfl multiplier : 0.9999999999946813        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7781e+05    |      100000 |   5.624e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.5577883417698 (tsim/hr)                              [sph::Model][rank=0]
---------------- t = 0.27477390770640203, dt = 0.004632159271225786 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.30 us    (1.8%)
   patch tree reduce : 1714.00 ns (0.4%)
   gen split merge   : 681.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1062.00 ns (0.3%)
   LB compute        : 383.16 us  (94.6%)
   LB move op cnt    : 0
   LB apply          : 3.83 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.48 us    (69.9%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001413152484719129,-7.733802507253164e-05,-2.5300492815883703e-07)
    sum a = (-0.0003768129599945333,-0.000722235826623063,-3.2264253091370554e-06)
    sum e = 0.05036541669722349
    sum de = 0.000998762900282307
Info: CFL hydro = 0.005267146203385006 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005267146203385006 cfl multiplier : 0.9999999999964541        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7871e+05    |      100000 |   5.596e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 29.801381860413343 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.2794060669776278, dt = 0.0005939330223722061 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.72 us    (1.7%)
   patch tree reduce : 1473.00 ns (0.4%)
   gen split merge   : 811.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1032.00 ns (0.3%)
   LB compute        : 378.69 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.57 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.36 us    (70.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.0001415043287016056,-7.78099324444336e-05,-2.551133813863438e-07)
    sum a = (-0.00037482500028219795,-0.000724598763110912,-3.2370187114644785e-06)
    sum e = 0.05035955333228594
    sum de = 0.0009999897833653441
Info: CFL hydro = 0.005299073689785458 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005299073689785458 cfl multiplier : 0.999999999997636         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8355e+05    |      100000 |   5.448e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 3.924492956321879 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 67                                                      [SPH][rank=0]
Info: time since start : 1245.914779312 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000014.vtk   [VTK Dump][rank=0]
              - took 8.57 ms, bandwidth = 623.31 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000014.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.52 us    (55.4%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000014.sham [Shamrock Dump][rank=0]
              - took 13.28 ms, bandwidth = 861.95 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1764.42 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.28, dt = 0.005299073689785458 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 25.50 us   (5.7%)
   patch tree reduce : 1703.00 ns (0.4%)
   gen split merge   : 982.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1092.00 ns (0.2%)
   LB compute        : 409.23 us  (91.0%)
   LB move op cnt    : 0
   LB apply          : 3.84 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.43 us    (70.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00014348996364141965,-8.165033639869582e-05,-2.722697279593322e-07)
    sum a = (-0.00035641036041744126,-0.000745509467170333,-3.331088498580162e-06)
    sum e = 0.05037334537933459
    sum de = 0.0009949530378502423
Info: CFL hydro = 0.005096274369083579 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.005096274369083579 cfl multiplier : 0.999999999998424         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7792e+05    |      100000 |   5.620e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 33.94135720745144 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.28529907368978547, dt = 0.004700926310214515 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.32 us    (1.8%)
   patch tree reduce : 1593.00 ns (0.4%)
   gen split merge   : 651.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1031.00 ns (0.2%)
   LB compute        : 391.24 us  (94.7%)
   LB move op cnt    : 0
   LB apply          : 3.54 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.32 us    (69.5%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00014511663221512277,-8.521032514829183e-05,-2.8817817089091826e-07)
    sum a = (-0.00033906094127563027,-0.0007637549565429962,-3.413828796285887e-06)
    sum e = 0.050376201283145225
    sum de = 0.0009914489201992916
Info: CFL hydro = 0.00491751548516008 sink sink = inf                                 [SPH][rank=0]
Info: cfl dt = 0.00491751548516008 cfl multiplier : 0.9999999999989493         [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7874e+05    |      100000 |   5.595e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.24931945252101 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 69                                                      [SPH][rank=0]
Info: time since start : 1248.862275893 (s)                                           [SPH][rank=0]
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1760.00 ms                           [sph::CartesianRender][rank=0]
---------------- t = 0.29, dt = 0.00491751548516008 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.44 us    (1.8%)
   patch tree reduce : 1713.00 ns (0.4%)
   gen split merge   : 782.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1062.00 ns (0.3%)
   LB compute        : 390.34 us  (94.8%)
   LB move op cnt    : 0
   LB apply          : 3.21 us    (0.8%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.42 us    (72.0%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00014674319047381084,-8.900898732446719e-05,-3.0516020488153965e-07)
    sum a = (-0.000319888276890693,-0.0007824728549156159,-3.4996442011258565e-06)
    sum e = 0.05038172429742387
    sum de = 0.0009857538248504603
Info: CFL hydro = 0.004734211710643756 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004734211710643756 cfl multiplier : 0.9999999999992996        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7836e+05    |      100000 |   5.607e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.57579311144661 (tsim/hr)                             [sph::Model][rank=0]
---------------- t = 0.29491751548516004, dt = 0.004734211710643756 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.38 us    (1.5%)
   patch tree reduce : 1493.00 ns (0.4%)
   gen split merge   : 822.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1022.00 ns (0.2%)
   LB compute        : 397.61 us  (95.0%)
   LB move op cnt    : 0
   LB apply          : 3.58 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.54 us    (72.6%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.000148210468363354,-9.275940225503145e-05,-3.21939260732676e-07)
    sum a = (-0.00030044355078907213,-0.0008000783248113816,-3.5815243687001822e-06)
    sum e = 0.05038585105092304
    sum de = 0.0009792838496861312
Info: CFL hydro = 0.004817016521798478 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004817016521798478 cfl multiplier : 0.9999999999995332        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.7773e+05    |      100000 |   5.627e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 30.290871943374132 (tsim/hr)                            [sph::Model][rank=0]
---------------- t = 0.2996517271958038, dt = 0.00034827280419619244 ----------------
Info: summary :                                                               [LoadBalance][rank=0]
Info:  - strategy "psweep" : max = 100000 min = 100000                        [LoadBalance][rank=0]
Info:  - strategy "round robin" : max = 100000 min = 100000                   [LoadBalance][rank=0]
Info: Loadbalance stats :                                                     [LoadBalance][rank=0]
    npatch = 1
    min = 100000
    max = 100000
    avg = 100000
    efficiency = 100.00%
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 7.14 us    (1.9%)
   patch tree reduce : 1472.00 ns (0.4%)
   gen split merge   : 691.00 ns  (0.2%)
   split / merge op  : 0/0
   apply split merge : 1042.00 ns (0.3%)
   LB compute        : 350.50 us  (94.4%)
   LB move op cnt    : 0
   LB apply          : 3.42 us    (0.9%)
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 2.38 us    (69.4%)
Info: free boundaries skipping geometry update                            [PositionUpdated][rank=0]
Info: conservation infos :                                                     [sph::Model][rank=0]
    sum v = (-0.00014826907695628873,-9.307972178765475e-05,-3.23380427291958e-07)
    sum a = (-0.00029897537019068414,-0.0008013557377120074,-3.5875177554412657e-06)
    sum e = 0.05037930943320254
    sum de = 0.0009803007582085326
Info: CFL hydro = 0.004812157229628818 sink sink = inf                                [SPH][rank=0]
Info: cfl dt = 0.004812157229628818 cfl multiplier : 0.9999999999996888        [sph::Model][rank=0]
Info: processing rate infos :                                                  [sph::Model][rank=0]
---------------------------------------------------------------------------------------
| rank |  rate  (N.s^-1)  |     Nobj    | t compute (s) |  MPI   | alloc |  mem (max) |
---------------------------------------------------------------------------------------
| 0    |    1.8297e+05    |      100000 |   5.465e-01   |    0 % |   0 % |    1.16 GB |
---------------------------------------------------------------------------------------
Info: estimated rate : 2.294010449332591 (tsim/hr)                             [sph::Model][rank=0]
Info: iteration since start : 72                                                      [SPH][rank=0]
Info: time since start : 1252.326775214 (s)                                           [SPH][rank=0]
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000015.vtk   [VTK Dump][rank=0]
              - took 8.09 ms, bandwidth = 660.41 MB/s
Info: Dumping state to _to_trash/circular_disc_central_pot_100000/dump/dump_0000015.sham  [SPH][rank=0]
Info: Scheduler step timings :                                                  [Scheduler][rank=0]
   metadata sync     : 6.88 us    (56.8%)
Info: dump to _to_trash/circular_disc_central_pot_100000/dump/dump_0000015.sham [Shamrock Dump][rank=0]
              - took 12.69 ms, bandwidth = 902.13 MB/s
Info: compute_column_integ field_name: rho, rays count: 1048576      [sph::CartesianRender][rank=0]
Info: compute_column_integ took 1766.18 ms                           [sph::CartesianRender][rank=0]

Plot generation (make_plots.py)

Load the on-the-fly analysis after the run to make the plots (everything in this section can be in another file)

import matplotlib
import matplotlib.pyplot as plt

# Uncomment this and replace by you dump folder, here since it is just above i comment it out
# dump_folder = "my_masterpiece"
# dump_folder += "/"


def plot_rho_integ(metadata, arr_rho, iplot):

    ext = metadata["extent"]

    dpi = 200

    # Reset the figure using the same memory as the last one
    plt.figure(num=1, clear=True, dpi=dpi)
    import copy

    my_cmap = matplotlib.colormaps["gist_heat"].copy()  # copy the default cmap
    my_cmap.set_bad(color="black")

    res = plt.imshow(
        arr_rho, cmap=my_cmap, origin="lower", extent=ext, norm="log", vmin=1e-8, vmax=1e-4
    )

    plt.xlabel("x")
    plt.ylabel("y")
    plt.title(f"t = {metadata['time']:0.3f} [years]")

    cbar = plt.colorbar(res, extend="both")
    cbar.set_label(r"$\int \rho \, \mathrm{d}z$ [code unit]")

    plt.savefig(plot_folder + "rho_integ_{:04}.png".format(iplot))
    plt.close()


def get_list_dumps_id():
    import glob

    list_files = glob.glob(plot_folder + "rho_integ_*.npy")
    list_files.sort()
    list_dumps_id = []
    for f in list_files:
        list_dumps_id.append(int(f.split("_")[-1].split(".")[0]))
    return list_dumps_id


def load_rho_integ(iplot):
    with open(plot_folder + f"rho_integ_{iplot:07}.json") as fp:
        metadata = json.load(fp)
    return np.load(plot_folder + f"rho_integ_{iplot:07}.npy"), metadata


if shamrock.sys.world_rank() == 0:
    for iplot in get_list_dumps_id():
        print("Rendering rho integ plot for dump", iplot)
        arr_rho, metadata = load_rho_integ(iplot)
        plot_rho_integ(metadata, arr_rho, iplot)

Rendering rho integ plot for dump 0
Rendering rho integ plot for dump 1
Rendering rho integ plot for dump 2
Rendering rho integ plot for dump 3
Rendering rho integ plot for dump 4
Rendering rho integ plot for dump 5
Rendering rho integ plot for dump 6
Rendering rho integ plot for dump 7
Rendering rho integ plot for dump 8
Rendering rho integ plot for dump 9
Rendering rho integ plot for dump 10
Rendering rho integ plot for dump 11
Rendering rho integ plot for dump 12
Rendering rho integ plot for dump 13
Rendering rho integ plot for dump 14
Rendering rho integ plot for dump 15
Rendering rho integ plot for dump 16
Rendering rho integ plot for dump 17
Rendering rho integ plot for dump 18
Rendering rho integ plot for dump 19
Rendering rho integ plot for dump 20
Rendering rho integ plot for dump 21
Rendering rho integ plot for dump 22
Rendering rho integ plot for dump 23
Rendering rho integ plot for dump 24
Rendering rho integ plot for dump 25
Rendering rho integ plot for dump 26
Rendering rho integ plot for dump 27
Rendering rho integ plot for dump 28
Rendering rho integ plot for dump 29
Rendering rho integ plot for dump 30

Make gif for the doc (plot_to_gif.py)

Convert PNG sequence to Image sequence in mpl

import matplotlib.animation as animation


def show_image_sequence(glob_str, render_gif):

    if render_gif and shamrock.sys.world_rank() == 0:

        import glob

        files = sorted(glob.glob(glob_str))

        from PIL import Image

        image_array = []
        for my_file in files:
            image = Image.open(my_file)
            image_array.append(image)

        if not image_array:
            raise RuntimeError(f"Warning: No images found for glob pattern: {glob_str}")

        pixel_x, pixel_y = image_array[0].size

        # Create the figure and axes objects
        # Remove axes, ticks, and frame & set aspect ratio
        dpi = 200
        fig = plt.figure(dpi=dpi)
        plt.gca().set_position((0, 0, 1, 1))
        plt.gcf().set_size_inches(pixel_x / dpi, pixel_y / dpi)
        plt.axis("off")

        # Set the initial image with correct aspect ratio
        im = plt.imshow(image_array[0], animated=True, aspect="auto")

        def update(i):
            im.set_array(image_array[i])
            return (im,)

        # Create the animation object
        ani = animation.FuncAnimation(
            fig,
            update,
            frames=len(image_array),
            interval=50,
            blit=True,
            repeat_delay=10,
        )

        return ani

Do it for rho integ

render_gif = True
glob_str = os.path.join(plot_folder, "rho_integ_*.png")

# If the animation is not returned only a static image will be shown in the doc
ani = show_image_sequence(glob_str, render_gif)

if render_gif and shamrock.sys.world_rank() == 0:
    # To save the animation using Pillow as a gif
    writer = animation.PillowWriter(fps=15, metadata=dict(artist="Me"), bitrate=1800)
    ani.save(analysis_folder + "rho_integ.gif", writer=writer)

    # Show the animation
    plt.show()

Once Loop Reflect

helper function to load data from JSON files

def load_data_from_json(filename, key):
    filepath = os.path.join(analysis_folder, filename)
    with open(filepath, "r") as fp:
        data = json.load(fp)[key]
    t = [d["t"] for d in data]
    values = [d[key] for d in data]
    return t, values

load the json file for barycenter

t, barycenter = load_data_from_json("barycenter.json", "barycenter")
barycenter_x = [d[0] for d in barycenter]
barycenter_y = [d[1] for d in barycenter]
barycenter_z = [d[2] for d in barycenter]

plt.figure(figsize=(8, 5), dpi=200)

plt.plot(t, barycenter_x)
plt.plot(t, barycenter_y)
plt.plot(t, barycenter_z)
plt.xlabel("t")
plt.ylabel("barycenter")
plt.legend(["x", "y", "z"])
plt.savefig(analysis_folder + "barycenter.png")
plt.show()

load the json file for disc_mass

t, disc_mass = load_data_from_json("disc_mass.json", "disc_mass")

plt.figure(figsize=(8, 5), dpi=200)

plt.plot(t, disc_mass)
plt.xlabel("t")
plt.ylabel("disc_mass")
plt.savefig(analysis_folder + "disc_mass.png")
plt.show()

load the json file for total_momentum

t, total_momentum = load_data_from_json("total_momentum.json", "total_momentum")
total_momentum_x = [d[0] for d in total_momentum]
total_momentum_y = [d[1] for d in total_momentum]
total_momentum_z = [d[2] for d in total_momentum]

plt.figure(figsize=(8, 5), dpi=200)

plt.plot(t, total_momentum_x)
plt.plot(t, total_momentum_y)
plt.plot(t, total_momentum_z)
plt.xlabel("t")
plt.ylabel("total_momentum")
plt.legend(["x", "y", "z"])
plt.savefig(analysis_folder + "total_momentum.png")
plt.show()

load the json file for energies

t, potential_energy = load_data_from_json("potential_energy.json", "potential_energy")
_, kinetic_energy = load_data_from_json("kinetic_energy.json", "kinetic_energy")

total_energy = [p + k for p, k in zip(potential_energy, kinetic_energy)]

plt.figure(figsize=(8, 5), dpi=200)
plt.plot(t, potential_energy)
plt.plot(t, kinetic_energy)
plt.plot(t, total_energy)
plt.xlabel("t")
plt.ylabel("energy")
plt.legend(["potential_energy", "kinetic_energy", "total_energy"])
plt.savefig(analysis_folder + "energies.png")
plt.show()

load the json file for sim_time_delta

t, sim_time_delta = load_data_from_json("sim_time_delta.json", "sim_time_delta")

plt.figure(figsize=(8, 5), dpi=200)
plt.plot(t, sim_time_delta)
plt.xlabel("t")
plt.ylabel("sim_time_delta")
plt.savefig(analysis_folder + "sim_time_delta.png")
plt.show()

load the json file for sim_step_count_delta

t, sim_step_count_delta = load_data_from_json("sim_step_count_delta.json", "sim_step_count_delta")

plt.figure(figsize=(8, 5), dpi=200)
plt.plot(t, sim_step_count_delta)
plt.xlabel("t")
plt.ylabel("sim_step_count_delta")
plt.savefig(analysis_folder + "sim_step_count_delta.png")
plt.show()

Time per step

t, sim_time_delta = load_data_from_json("sim_time_delta.json", "sim_time_delta")
_, sim_step_count_delta = load_data_from_json("sim_step_count_delta.json", "sim_step_count_delta")
_, part_count = load_data_from_json("part_count.json", "part_count")

time_per_step = []

for td, sc, pc in zip(sim_time_delta, sim_step_count_delta, part_count):
    if sc > 0:
        time_per_step.append(td / sc)
    else:
        # NAN here because the step count is 0
        time_per_step.append(np.nan)

plt.figure(figsize=(8, 5), dpi=200)
plt.plot(t, time_per_step, "+-")
plt.xlabel("t")
plt.ylabel("time_per_step")
plt.savefig(analysis_folder + "time_per_step.png")
plt.show()

rate = []

for td, sc, pc in zip(sim_time_delta, sim_step_count_delta, part_count):
    if sc > 0:
        rate.append(pc / (td / sc))
    else:
        # NAN here because the step count is 0
        rate.append(np.nan)

plt.figure(figsize=(8, 5), dpi=200)
plt.plot(t, rate, "+-")
plt.xlabel("t")
plt.ylabel("Particles / second")
plt.yscale("log")
plt.savefig(analysis_folder + "rate.png")
plt.show()

• 
• 

Estimated memory usage: 546 MB