Note
Go to the end to download the full example code.
Shearing box in SPH#
This simple example shows how to run an unstratified shearing box simulaiton
9 import shamrock
10
11 # If we use the shamrock executable to run this script instead of the python interpreter,
12 # we should not initialize the system as the shamrock executable needs to handle specific MPI logic
13 if not shamrock.sys.is_initialized():
14 shamrock.change_loglevel(1)
15 shamrock.sys.init("0:0")
16
17 from math import exp
18
19 import matplotlib.pyplot as plt
20 import numpy as np
Initialize context & attach a SPH model to it
25 ctx = shamrock.Context()
26 ctx.pdata_layout_new()
27
28 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
Setup parameters
33 gamma = 5.0 / 3.0
34 rho = 1
35 uint = 1
36
37 dr = 0.02
38 bmin = (-0.6, -0.6, -0.1)
39 bmax = (0.6, 0.6, 0.1)
40 pmass = -1
41
42 bmin, bmax = model.get_ideal_fcc_box(dr, bmin, bmax)
43 xm, ym, zm = bmin
44 xM, yM, zM = bmax
45
46 Omega_0 = 1
47 eta = 0.00
48 q = 3.0 / 2.0
49
50 shear_speed = -q * Omega_0 * (xM - xm)
Generate the config & init the scheduler
55 cfg = model.gen_default_config()
56 # cfg.set_artif_viscosity_Constant(alpha_u = 1, alpha_AV = 1, beta_AV = 2)
57 # cfg.set_artif_viscosity_VaryingMM97(alpha_min = 0.1,alpha_max = 1,sigma_decay = 0.1, alpha_u = 1, beta_AV = 2)
58 cfg.set_artif_viscosity_VaryingCD10(
59 alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
60 )
61 cfg.set_boundary_shearing_periodic((1, 0, 0), (0, 1, 0), shear_speed)
62 cfg.set_eos_adiabatic(gamma)
63 cfg.add_ext_force_shearing_box(Omega_0=Omega_0, eta=eta, q=q)
64 cfg.set_units(shamrock.UnitSystem())
65 cfg.print_status()
66 model.set_solver_config(cfg)
67
68 model.init_scheduler(int(1e7), 1)
69
70 model.resize_simulation_box(bmin, bmax)
----- SPH Solver configuration -----
units :
unit_length : 1
unit_mass : 1
unit_current : 1
unit_temperature : 1
unit_qte : 1
unit_lumint : 1
part mass 0 ( can be changed using .set_part_mass() )
cfl force 0
cfl courant 0
--- artificial viscosity config
Config Type : VaryingCD10 (Cullen & Dehnen 2010)
alpha_min = 0
alpha_max = 1
sigma_decay = 0.1
alpha_u = 1
beta_AV = 2
--- artificial viscosity config (deduced)
-------------
EOS config f64_3 :
adiabatic :
gamma 1.6666666666666667
--- Bondaries config
Config Type : ShearingPeriodic (Stone 2010)
shear_base = (1,0,0)
shear_dir = (0,1,0)
shear_speed = -1.8000000000000003
--- Bondaries config config (deduced)
-------------
------------------------------------
Add the particles & set fields values Note that every field that are not mentionned are set to zero
76 model.add_cube_fcc_3d(dr, bmin, bmax)
77
78 vol_b = (xM - xm) * (yM - ym) * (zM - zm)
79
80 totmass = rho * vol_b
81 # print("Total mass :", totmass)
82
83 pmass = model.total_mass_to_part_mass(totmass)
84
85 model.set_value_in_a_box("uint", "f64", 1, bmin, bmax)
86 # model.set_value_in_a_box("vxyz","f64_3", (-10,0,0) , bmin,bmax)
87
88 pen_sz = 0.1
89
90 mm = 1
91 MM = 0
92
93
94 def vel_func(r):
95 global mm, MM
96 x, y, z = r
97
98 s = (x - (xM + xm) / 2) / (xM - xm)
99 vel = (shear_speed) * s
100
101 mm = min(mm, vel)
102 MM = max(MM, vel)
103
104 return (0, vel, 0.0)
105 # return (1,0,0)
106
107
108 model.set_field_value_lambda_f64_3("vxyz", vel_func)
109 # print("Current part mass :", pmass)
110 model.set_particle_mass(pmass)
111
112
113 tot_u = pmass * model.get_sum("uint", "f64")
114 # print("total u :",tot_u)
115
116 print(f"v_shear = {shear_speed} | dv = {MM-mm}")
117
118
119 model.set_cfl_cour(0.3)
120 model.set_cfl_force(0.25)
Info: Add fcc lattice size : ( 31 34 6 ) [SPH][rank=0]
Info: Push particles : [Model][rank=0]
rank = 0 patch id=0, add N=6120 particles, coords = (-0.6,-0.6,-0.1) (0.6000000000000002,0.5777945491468365,0.09595917942265422)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 17.01 us (73.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.15 us (0.6%)
patch tree reduce : 8.88 us (1.7%)
gen split merge : 671.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 962.00 ns (0.2%)
LB compute : 496.98 us (95.7%)
LB move op cnt : 0
LB apply : 3.48 us (0.7%)
Info: current particle counts : min = 6120 max = 6120 [Model][rank=0]
v_shear = -1.8000000000000003 | dv = 1.7700000000000002
123 dump_folder = "_to_trash"
124 import os
125
126 os.system("mkdir -p " + dump_folder)
0
Perform the plot
131 def plot():
132 dic = ctx.collect_data()
133 fig, axs = plt.subplots(2, 1, figsize=(5, 8), sharex=True)
134 fig.suptitle("t = {:.2f}".format(model.get_time()))
135 axs[0].scatter(dic["xyz"][:, 0], dic["xyz"][:, 1], s=1)
136 axs[1].scatter(dic["xyz"][:, 0], dic["vxyz"][:, 1], s=1)
137
138 axs[0].set_ylabel("y")
139 axs[1].set_ylabel("vy")
140 axs[1].set_xlabel("x")
141
142 plt.tight_layout()
143 plt.show()
Performing the timestep loop
148 model.timestep()
149
150 dt_stop = 0.1
151 for i in range(2):
152
153 t_target = i * dt_stop
154 # skip if the model is already past the target
155 if model.get_time() > t_target:
156 continue
157
158 model.evolve_until(i * dt_stop)
159
160 # Dump name is "dump_xxxx.sham" where xxxx is the timestep
161 model.do_vtk_dump(dump_folder + "/dump_{:04}.vtk".format(i), True)
162 plot()
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 13.45 us (1.5%)
patch tree reduce : 1923.00 ns (0.2%)
gen split merge : 591.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 942.00 ns (0.1%)
LB compute : 853.48 us (96.0%)
LB move op cnt : 0
LB apply : 4.24 us (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 24.52 us (94.3%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 0.741503267973856
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.022000000000000002 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 0.7684640522875816
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.024200000000000003 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 0.7684640522875816
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.026620000000000005 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 0.794607843137255
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.029282000000000006 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 0.8504901960784312
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.03221020000000001 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.320424836601307
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.03543122000000001 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.320424836601307
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.038974342000000016 unconverged cnt = 6120
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4184640522875813
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 2.8314338416863407e-07, max = 2.831433850588215e-07
iterations = 3
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0.004154393760827157,0)
sum a = (5.6948709491889e-17,-8.766053124257185e-18,0.005082303061895428)
sum e = 0.3143698998713959
sum de = 9.604356448640433e-05
Info: cfl dt = 0.00011089173088629873 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 | 2.5309e+04 | 6120 | 2.418e-01 | 0 % | 1 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 1 [SPH][rank=0]
Info: time since start : 17.114109471000003 (s) [SPH][rank=0]
Info: dump to _to_trash/dump_0000.vtk [VTK Dump][rank=0]
- took 5.47 ms, bandwidth = 89.67 MB/s
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
---------------- t = 0, dt = 0.00011089173088629873 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 10.07 us (1.4%)
patch tree reduce : 1262.00 ns (0.2%)
gen split merge : 781.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1052.00 ns (0.2%)
LB compute : 675.22 us (96.6%)
LB move op cnt : 0
LB apply : 3.70 us (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.52 us (71.5%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4470588235294115
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 1.5839830969886152e-09, max = 3.408253654510689e-09
iterations = 0
Info: conservation infos : [sph::Model][rank=0]
sum v = (6.315145049214027e-21,-4.3007972628265963e-17,5.635853834223198e-07)
sum a = (5.19318507104892e-17,6.088338327277574e-18,0.00508230303064696)
sum e = 0.31436991054897867
sum de = 9.74466885364498e-05
Info: cfl dt = 0.0037703106533756884 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 | 6.0753e+04 | 6120 | 1.007e-01 | 0 % | 1 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 3.962930231370381 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.00011089173088629873, dt = 0.0037703106533756884 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.22 us (1.2%)
patch tree reduce : 932.00 ns (0.3%)
gen split merge : 511.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 982.00 ns (0.3%)
LB compute : 347.73 us (95.7%)
LB move op cnt : 0
LB apply : 3.20 us (0.9%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.65 us (67.4%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4388888888888887
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 7.172685023217074e-12, max = 3.3201651540132534e-11
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (2.0183625283288808e-19,3.778270492576636e-18,1.9725446641821508e-05)
sum a = (6.111891063456037e-17,-5.003473743189519e-18,0.005082264782670164)
sum e = 0.3143702833343977
sum de = 0.00013824059786473694
Info: cfl dt = 0.006209449394346196 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.1060e+05 | 6120 | 5.533e-02 | 0 % | 1 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 245.30215768507045 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.003881202384261987, dt = 0.006209449394346196 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.98 us (0.8%)
patch tree reduce : 1002.00 ns (0.3%)
gen split merge : 551.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 741.00 ns (0.2%)
LB compute : 363.77 us (96.5%)
LB move op cnt : 0
LB apply : 2.88 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.28 us (69.9%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4367647058823527
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 2.3294096404573677e-10, max = 1.0880275964821347e-09
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (5.986590816098565e-19,-8.32023395705706e-18,5.128344051510254e-05)
sum a = (5.850965248209674e-17,5.758902876894751e-19,0.005082044319848336)
sum e = 0.3143712315922659
sum de = 0.00020394600495906925
Info: cfl dt = 0.007834596847674341 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.0838e+05 | 6120 | 5.647e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 395.8594456000552 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.010090651778608184, dt = 0.007834596847674341 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.40 us (1.0%)
patch tree reduce : 1052.00 ns (0.3%)
gen split merge : 511.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1142.00 ns (0.3%)
LB compute : 330.72 us (95.9%)
LB move op cnt : 0
LB apply : 2.55 us (0.7%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.14 us (70.6%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4455882352941174
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 1.3550983973723722e-09, max = 6.545525037921152e-09
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (1.0489100801015626e-18,-1.9011935563710095e-17,9.109852444675989e-05)
sum a = (3.836615128341325e-17,2.6808543685796116e-18,0.005081486569882731)
sum e = 0.3143730498982785
sum de = 0.00028163943820289555
Info: cfl dt = 0.00891668243491238 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.1017e+05 | 6120 | 5.555e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 507.73192692718993 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.017925248626282524, dt = 0.00891668243491238 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.87 us (0.7%)
patch tree reduce : 511.00 ns (0.1%)
gen split merge : 541.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 802.00 ns (0.2%)
LB compute : 382.61 us (96.7%)
LB move op cnt : 0
LB apply : 3.10 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.17 us (70.0%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4475490196078429
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 3.615263896248438e-09, max = 1.8821916322165866e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (1.3122342115785653e-18,-3.336132456211285e-18,0.00013640634161461525)
sum a = (7.399280524491505e-17,-3.84197699408773e-18,0.005080472285302088)
sum e = 0.31437589351354017
sum de = 0.0003608568219550676
Info: cfl dt = 0.009636396780122314 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.0726e+05 | 6120 | 5.706e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 562.5900331425333 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.026841931061194903, dt = 0.009636396780122314 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.67 us (1.3%)
patch tree reduce : 962.00 ns (0.3%)
gen split merge : 541.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1182.00 ns (0.3%)
LB compute : 337.58 us (95.3%)
LB move op cnt : 0
LB apply : 3.02 us (0.9%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.22 us (70.9%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4460784313725488
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 5.752376507485498e-09, max = 3.4459765836674773e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (2.1844205723774023e-18,-1.6077746776921855e-18,0.0001853592663594488)
sum a = (4.141401476421019e-17,-7.070126341453038e-19,0.005078921976555799)
sum e = 0.3143797778012691
sum de = 0.0004329654653391086
Info: cfl dt = 0.010114315298959963 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.0749e+05 | 6120 | 5.693e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 609.3198051211052 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.036478327841317215, dt = 0.010114315298959963 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.30 us (1.0%)
patch tree reduce : 621.00 ns (0.2%)
gen split merge : 501.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 902.00 ns (0.3%)
LB compute : 331.58 us (96.1%)
LB move op cnt : 0
LB apply : 2.69 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.56 us (69.4%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4460784313725488
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 5.818251260610762e-09, max = 4.486423362805652e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (2.4457526490474497e-18,-3.472793303815121e-17,0.00023672161491404577)
sum a = (3.43322697385629e-17,4.401440189428226e-18,0.005076787484773814)
sum e = 0.31438460229587495
sum de = 0.0004912658954562416
Info: cfl dt = 0.010430975181113323 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.1014e+05 | 6120 | 5.556e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 655.311739605749 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.04659264314027718, dt = 0.010430975181113323 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.27 us (0.9%)
patch tree reduce : 1032.00 ns (0.3%)
gen split merge : 541.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1113.00 ns (0.3%)
LB compute : 363.55 us (96.3%)
LB move op cnt : 0
LB apply : 2.60 us (0.7%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.69 us (66.8%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4483660130718952
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 3.405816618081674e-09, max = 4.2910854860000134e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (2.7687402939544358e-18,-2.893994419845934e-18,0.00028966666470606504)
sum a = (5.2166007921276686e-17,-7.135285178488415e-18,0.0050740421695320065)
sum e = 0.3143901775107707
sum de = 0.0005310720353439191
Info: cfl dt = 0.01064024299013669 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.1174e+05 | 6120 | 5.477e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 685.6316996030756 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.0570236183213905, dt = 0.01064024299013669 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.24 us (1.0%)
patch tree reduce : 521.00 ns (0.2%)
gen split merge : 521.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 992.00 ns (0.3%)
LB compute : 308.68 us (95.9%)
LB move op cnt : 0
LB apply : 2.50 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.42 us (67.2%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4483660130718952
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 6.124547504410648e-10, max = 2.9116102549298282e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (3.417738011751374e-18,-6.431098710768742e-19,0.0003436413881745103)
sum a = (6.22987566126383e-17,-1.771660928529598e-17,0.0050706729699212345)
sum e = 0.3143962517375111
sum de = 0.0005501730733550375
Info: cfl dt = 0.010583003685891922 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.1094e+05 | 6120 | 5.517e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 694.3609920095762 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.0676638613115272, dt = 0.010583003685891922 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.46 us (1.0%)
patch tree reduce : 631.00 ns (0.2%)
gen split merge : 541.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 931.00 ns (0.3%)
LB compute : 328.82 us (96.0%)
LB move op cnt : 0
LB apply : 2.79 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.24 us (71.1%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4460784313725488
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 2.582426396378962e-10, max = 1.762618766251557e-08
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (4.1282441016172494e-18,-1.680124538188334e-17,0.0003972864143538687)
sum a = (4.909490707209122e-17,-2.156129468887494e-17,0.005066752443936444)
sum e = 0.31440241760719223
sum de = 0.0005615715602514454
Info: cfl dt = 0.010635095941600751 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.1090e+05 | 6120 | 5.518e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 690.3886886022775 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.07824686499741912, dt = 0.010635095941600751 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.77 us (1.1%)
patch tree reduce : 891.00 ns (0.3%)
gen split merge : 550.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 992.00 ns (0.3%)
LB compute : 314.89 us (90.4%)
LB move op cnt : 0
LB apply : 20.84 us (6.0%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.26 us (69.5%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4477124183006533
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 1.06301328742834e-09, max = 6.683179663726466e-09
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (4.579881197357637e-18,-2.5402839907536532e-17,0.0004511510672369989)
sum a = (4.0979413171646526e-17,-6.481341669446623e-19,0.005062240947242894)
sum e = 0.3144087119423163
sum de = 0.0005739543750557067
Info: cfl dt = 0.010652394617585905 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.1022e+05 | 6120 | 5.553e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 689.5060664453872 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.08888196093901987, dt = 0.010652394617585905 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.50 us (1.0%)
patch tree reduce : 471.00 ns (0.1%)
gen split merge : 501.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1082.00 ns (0.3%)
LB compute : 328.65 us (96.0%)
LB move op cnt : 0
LB apply : 2.63 us (0.8%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.43 us (68.6%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4452614379084965
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 2.3014687575276836e-10, max = 6.215082235314744e-09
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (4.970598705701095e-18,-3.890814720015089e-17,0.0005050520653562442)
sum a = (7.083001099614246e-17,-1.4940685818343064e-17,0.0050571481484674215)
sum e = 0.3144151462463685
sum de = 0.0005908429305553166
Info: cfl dt = 0.010647033142926733 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.1205e+05 | 6120 | 5.462e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 702.1035747366844 (tsim/hr) [sph::Model][rank=0]
---------------- t = 0.09953435555660578, dt = 0.0004656444433942253 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 6120 min = 6120 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 6120
max = 6120
avg = 6120
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.53 us (1.0%)
patch tree reduce : 641.00 ns (0.2%)
gen split merge : 561.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 991.00 ns (0.3%)
LB compute : 342.86 us (96.1%)
LB move op cnt : 0
LB apply : 2.60 us (0.7%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.28 us (70.1%)
Warning: High interface/patch volume ratio. [InterfaceGen][rank=0]
This can lead to high mpi overhead, try to increase the patch split crit
patch 0 high interf/patch volume: 1.4452614379084965
Info: smoothing length iteration converged [Smoothinglength][rank=0]
eps min = 1.0704970741195936e-11, max = 7.460773315354855e-07
iterations = 1
Info: conservation infos : [sph::Model][rank=0]
sum v = (5.164505124348541e-18,-2.1343208846363764e-17,0.0005073797730398672)
sum a = (5.699812447212187e-17,-3.6187490987743646e-18,0.0050569124381528545)
sum e = 0.3144144614858832
sum de = 0.0005914932822185097
Info: cfl dt = 0.010671519169061247 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.1320e+05 | 6120 | 5.407e-02 | 0 % | 0 % | 10.01 MB |
---------------------------------------------------------------------------------------
Info: estimated rate : 31.005638420253746 (tsim/hr) [sph::Model][rank=0]
Info: iteration since start : 14 [SPH][rank=0]
Info: time since start : 17.987789611 (s) [SPH][rank=0]
Info: dump to _to_trash/dump_0001.vtk [VTK Dump][rank=0]
- took 2.79 ms, bandwidth = 175.62 MB/s
adding -> xyz
adding -> vxyz
adding -> axyz
adding -> axyz_ext
adding -> hpart
adding -> uint
adding -> duint
adding -> alpha_AV
adding -> divv
adding -> dtdivv
adding -> curlv
adding -> soundspeed
Total running time of the script: (0 minutes 2.119 seconds)
Estimated memory usage: 66 MB