Note
Go to the end to download the full example code.
Test the precision of SG methods in SPH#
Test that all methodes give expected deviantion compare to the reference mode
8 import matplotlib.pyplot as plt
9 import numpy as np
10
11 import shamrock
Shamrock init
16 # Self-gravity is still an experimental feature
17 shamrock.enable_experimental_features()
18
19 # If we use the shamrock executable to run this script instead of the python interpreter,
20 # we should not initialize the system as the shamrock executable needs to handle specific MPI logic
21 if not shamrock.sys.is_initialized():
22 shamrock.change_loglevel(1)
23 shamrock.sys.init("0:0")
Parameters of the test
28 # tolerances for the test for each quantity [min, max] outside = fail
29 # here the challenge is to find the right tolerances for many combinations of hardware
30 TOL_HPC_CUBE = {
31 "direct": {
32 "max_rel_delta": [0.0, 1e-15],
33 "avg_rel_delta": [0.0, 1e-15],
34 "min_rel_delta": [0.0, 1e-16],
35 "std_rel_delta": [0.0, 1e-16],
36 },
37 "mm1": {
38 "max_rel_delta": [0.10237644408204995 - 2e-6, 0.10237644408204995 + 2e-6],
39 "avg_rel_delta": [0.02854797753451957 - 1e-5, 0.02854797753451957 + 1e-5],
40 "min_rel_delta": [0.0011855528704246662 - 2e-6, 0.0011855528704246662 + 2e-6],
41 "std_rel_delta": [0.012057071282465973 - 1e-5, 0.012057071282465973 + 1e-5],
42 },
43 "mm2": {
44 "max_rel_delta": [0.09870604133669476 - 2e-7, 0.09870604133669476 + 2e-7],
45 "avg_rel_delta": [0.02809153086497972 - 4e-6, 0.02809153086497972 + 4e-6],
46 "min_rel_delta": [0.0007660948531784951 - 2e-7, 0.0007660948531784951 + 2e-7],
47 "std_rel_delta": [0.011728053659573516 - 3e-6, 0.011728053659573516 + 3e-6],
48 },
49 "mm3": {
50 "max_rel_delta": [0.053440736724222254 - 1e-2, 0.053440736724222254 + 1e-2],
51 "avg_rel_delta": [0.016082633182446748 - 1e-5, 0.016082633182446748 + 1e-5],
52 "min_rel_delta": [0.0006196291223064161 - 1e-5, 0.0006196291223064161 + 1e-5],
53 "std_rel_delta": [0.007676767861132326 - 1e-5, 0.007676767861132326 + 1e-5],
54 },
55 "mm4": {
56 "max_rel_delta": [0.05833850788436308 - 4e-7, 0.05833850788436308 + 4e-7],
57 "avg_rel_delta": [0.014909667170457708 - 4e-6, 0.014909667170457708 + 4e-6],
58 "min_rel_delta": [0.000900026203337091 - 4e-7, 0.000900026203337091 + 4e-7],
59 "std_rel_delta": [0.007051020680900702 - 1e-5, 0.007051020680900702 + 1e-5],
60 },
61 "mm5": {
62 "max_rel_delta": [0.03130228236510434 - 4e-8, 0.03130228236510434 + 4e-8],
63 "avg_rel_delta": [0.005541631318650798 - 2e-6, 0.005541631318650798 + 2e-6],
64 "min_rel_delta": [0.00011105937841131173 - 4e-8, 0.00011105937841131173 + 4e-8],
65 "std_rel_delta": [0.003575474772831352 - 4e-7, 0.003575474772831352 + 4e-7],
66 },
67 "fmm1": {
68 "max_rel_delta": [0.09457251908559332 - 1e-15, 0.09457251908559332 + 1e-15],
69 "avg_rel_delta": [0.025849750275017017 - 1.5e-7, 0.025849750275017017 + 1.5e-7],
70 "min_rel_delta": [0.0009107123147511316 - 1e-15, 0.0009107123147511316 + 1e-15],
71 "std_rel_delta": [0.014771871930712405 - 5e-7, 0.014771871930712405 + 5e-7],
72 },
73 "fmm2": {
74 "max_rel_delta": [0.0637884157532678 - 2e-15, 0.0637884157532678 + 2e-15],
75 "avg_rel_delta": [0.011974487570609172 - 1e-6, 0.011974487570609172 + 1e-6],
76 "min_rel_delta": [0.0003606563787779815 - 1e-15, 0.0003606563787779815 + 1e-15],
77 "std_rel_delta": [0.005685736194766207 - 1e-6, 0.005685736194766207 + 1e-6],
78 },
79 "fmm3": {
80 "max_rel_delta": [0.041451600220479765 - 1e-14, 0.041451600220479765 + 1e-14],
81 "avg_rel_delta": [0.00590471908748555 - 1e-7, 0.00590471908748555 + 1e-7],
82 "min_rel_delta": [0.00011353264639102378 - 1e-14, 0.00011353264639102378 + 1e-14],
83 "std_rel_delta": [0.003551982999913884 - 1e-6, 0.003551982999913884 + 1e-6],
84 },
85 "fmm4": {
86 "max_rel_delta": [0.027041008234275913 - 1e-14, 0.027041008234275913 + 1e-14],
87 "avg_rel_delta": [0.004870267930794465 - 5e-7, 0.004870267930794465 + 5e-7],
88 "min_rel_delta": [9.263728961868266e-05 - 1e-14, 9.263728961868266e-05 + 1e-14],
89 "std_rel_delta": [0.0028526970191818392 - 1e-7, 0.0028526970191818392 + 1e-7],
90 },
91 "fmm5": {
92 "max_rel_delta": [0.013154949137016961 - 1e-14, 0.013154949137016961 + 1e-14],
93 "avg_rel_delta": [0.001933834336767182 - 5e-7, 0.001933834336767182 + 5e-7],
94 "min_rel_delta": [3.7696982571244474e-05 - 1e-14, 3.7696982571244474e-05 + 1e-14],
95 "std_rel_delta": [0.0014635310344450105 - 5e-7, 0.0014635310344450105 + 5e-7],
96 },
97 "sfmm1": {
98 "max_rel_delta": [0.21892380131546568 - 1e-15, 0.21892380131546568 + 1e-15],
99 "avg_rel_delta": [0.05782934468995147 - 5e-7, 0.05782934468995147 + 5e-7],
100 "min_rel_delta": [0.0013277688426936505 - 1e-14, 0.0013277688426936505 + 1e-14],
101 "std_rel_delta": [0.03078712426421105 - 5e-7, 0.03078712426421105 + 5e-7],
102 },
103 "sfmm2": {
104 "max_rel_delta": [0.03804965294061554 - 5e-15, 0.03804965294061554 + 5e-15],
105 "avg_rel_delta": [0.011068506783910114 - 5e-7, 0.011068506783910114 + 5e-7],
106 "min_rel_delta": [0.0005623871943121567 - 1e-15, 0.0005623871943121567 + 1e-15],
107 "std_rel_delta": [0.005375924880946539 - 5e-7, 0.005375924880946539 + 5e-7],
108 },
109 "sfmm3": {
110 "max_rel_delta": [0.01693824120108178 - 5e-15, 0.01693824120108178 + 5e-15],
111 "avg_rel_delta": [0.004724856483131658 - 1e-6, 0.004724856483131658 + 1e-6],
112 "min_rel_delta": [0.00011119860215853415 - 1e-15, 0.00011119860215853415 + 1e-15],
113 "std_rel_delta": [0.0026448823012749684 - 5e-7, 0.0026448823012749684 + 5e-7],
114 },
115 "sfmm4": {
116 "max_rel_delta": [0.011356406701671626 - 1e-15, 0.011356406701671626 + 1e-15],
117 "avg_rel_delta": [0.0022876264690043753 - 5e-8, 0.0022876264690043753 + 5e-8],
118 "min_rel_delta": [5.847968796399966e-05 - 5e-15, 5.847968796399966e-05 + 5e-15],
119 "std_rel_delta": [0.0015169740765526084 - 5e-9, 0.0015169740765526084 + 5e-9],
120 },
121 "sfmm5": {
122 "max_rel_delta": [0.010203458622545465 - 1e-15, 0.010203458622545465 + 1e-15],
123 "avg_rel_delta": [0.0014135757954817333 - 5e-7, 0.0014135757954817333 + 5e-7],
124 "min_rel_delta": [3.503562542652889e-05 - 1e-15, 3.503562542652889e-05 + 1e-15],
125 "std_rel_delta": [0.001217910806477318 - 5e-7, 0.001217910806477318 + 5e-7],
126 },
127 }
Helper functions for this test
134 # helper to run one case (SG config & setup)
135 def run_case(setup_func, setup_name, sg_setup_func):
136 ctx = shamrock.Context()
137 ctx.pdata_layout_new()
138
139 model = shamrock.get_Model_SPH(context=ctx, vector_type="f64_3", sph_kernel="M4")
140
141 cfg = model.gen_default_config()
142
143 sg_setup_func(cfg)
144
145 setup_func(model, cfg)
146
147 model.timestep()
148
149 data = ctx.collect_data()
150
151 return data
152
153
154 # Compare the SG method to the reference and the one without SG and return error metrics
155 def compare_sg_methods_data(no_sg_data, reference_data, data_to_comp, sat_relative_error=1):
156 a_sg = data_to_comp["axyz"] - no_sg_data["axyz"]
157
158 a_sg_ref = reference_data["axyz"] - no_sg_data["axyz"]
159 delta_sg = a_sg - a_sg_ref
160
161 delta_sg_norm = np.linalg.norm(delta_sg, axis=1)
162 rel_delta_norm = delta_sg_norm / (np.max(np.linalg.norm(a_sg_ref, axis=1)))
163
164 return delta_sg, rel_delta_norm, data_to_comp["xyz"]
165
166
167 to_raise = []
168
169
170 # Compute error related quantities and check if they are within the tolerances
171 def check_print_errors(rel_delta, setup_name, method_name, tols):
172 global to_raise
173 max_rel_delta = np.max(np.abs(rel_delta))
174 if shamrock.sys.world_rank() == 0:
175 print(f"max relative error {method_name}: {max_rel_delta} for {setup_name}")
176 avg_rel_delta = np.mean(np.abs(rel_delta))
177 if shamrock.sys.world_rank() == 0:
178 print(f"avg relative error {method_name}: {avg_rel_delta} for {setup_name}")
179 min_rel_delta = np.min(np.abs(rel_delta))
180 if shamrock.sys.world_rank() == 0:
181 print(f"min relative error {method_name}: {min_rel_delta} for {setup_name}")
182 std_rel_delta = np.std(np.abs(rel_delta))
183 if shamrock.sys.world_rank() == 0:
184 print(f"std relative error {method_name}: {std_rel_delta} for {setup_name}")
185
186 delta_max_tol = (tols["max_rel_delta"][1] + tols["max_rel_delta"][0]) / 2 - max_rel_delta
187 delta_avg_tol = (tols["avg_rel_delta"][1] + tols["avg_rel_delta"][0]) / 2 - avg_rel_delta
188 delta_min_tol = (tols["min_rel_delta"][1] + tols["min_rel_delta"][0]) / 2 - min_rel_delta
189 delta_std_tol = (tols["std_rel_delta"][1] + tols["std_rel_delta"][0]) / 2 - std_rel_delta
190
191 if max_rel_delta > tols["max_rel_delta"][1] or max_rel_delta < tols["max_rel_delta"][0]:
192 to_raise.append(
193 f"max relative error {method_name} is out of tolerance for {setup_name}: {max_rel_delta} not in [{tols['max_rel_delta'][0]}, {tols['max_rel_delta'][1]}], delta = {delta_max_tol}"
194 )
195 if avg_rel_delta > tols["avg_rel_delta"][1] or avg_rel_delta < tols["avg_rel_delta"][0]:
196 to_raise.append(
197 f"avg relative error {method_name} is out of tolerance for {setup_name}: {avg_rel_delta} not in [{tols['avg_rel_delta'][0]}, {tols['avg_rel_delta'][1]}], delta = {delta_avg_tol}"
198 )
199 if min_rel_delta > tols["min_rel_delta"][1] or min_rel_delta < tols["min_rel_delta"][0]:
200 to_raise.append(
201 f"min relative error {method_name} is out of tolerance for {setup_name}: {min_rel_delta} not in [{tols['min_rel_delta'][0]}, {tols['min_rel_delta'][1]}], delta = {delta_min_tol}"
202 )
203 if std_rel_delta > tols["std_rel_delta"][1] or std_rel_delta < tols["std_rel_delta"][0]:
204 to_raise.append(
205 f"std relative error {method_name} is out of tolerance for {setup_name}: {std_rel_delta} not in [{tols['std_rel_delta'][0]}, {tols['std_rel_delta'][1]}], delta = {delta_std_tol}"
206 )
207
208
209 # Compare the SG method to the reference and the one without SG and return error metrics
210 def compare_sg_methods(setup_func, setup_name, tols):
211 def sg_case_none(cfg):
212 cfg.set_self_gravity_none()
213
214 def sg_case_reference(cfg):
215 cfg.set_self_gravity_direct(reference_mode=True)
216 cfg.set_softening_plummer(epsilon=1e-9)
217
218 def sg_case_direct(cfg):
219 cfg.set_self_gravity_direct(reference_mode=False)
220 cfg.set_softening_plummer(epsilon=1e-9)
221
222 def sg_case_mm1(cfg):
223 cfg.set_self_gravity_mm(order=1, opening_angle=0.5, reduction_level=3)
224 cfg.set_softening_plummer(epsilon=1e-9)
225
226 def sg_case_mm2(cfg):
227 cfg.set_self_gravity_mm(order=2, opening_angle=0.5, reduction_level=3)
228 cfg.set_softening_plummer(epsilon=1e-9)
229
230 def sg_case_mm3(cfg):
231 cfg.set_self_gravity_mm(order=3, opening_angle=0.5, reduction_level=3)
232 cfg.set_softening_plummer(epsilon=1e-9)
233
234 def sg_case_mm4(cfg):
235 cfg.set_self_gravity_mm(order=4, opening_angle=0.5, reduction_level=3)
236 cfg.set_softening_plummer(epsilon=1e-9)
237
238 def sg_case_mm5(cfg):
239 cfg.set_self_gravity_mm(order=5, opening_angle=0.5, reduction_level=3)
240 cfg.set_softening_plummer(epsilon=1e-9)
241
242 def sg_case_fmm1(cfg):
243 cfg.set_self_gravity_fmm(order=1, opening_angle=0.5, reduction_level=3)
244 cfg.set_softening_plummer(epsilon=1e-9)
245
246 def sg_case_fmm2(cfg):
247 cfg.set_self_gravity_fmm(order=2, opening_angle=0.5, reduction_level=3)
248 cfg.set_softening_plummer(epsilon=1e-9)
249
250 def sg_case_fmm3(cfg):
251 cfg.set_self_gravity_fmm(order=3, opening_angle=0.5, reduction_level=3)
252 cfg.set_softening_plummer(epsilon=1e-9)
253
254 def sg_case_fmm4(cfg):
255 cfg.set_self_gravity_fmm(order=4, opening_angle=0.5, reduction_level=3)
256 cfg.set_softening_plummer(epsilon=1e-9)
257
258 def sg_case_fmm5(cfg):
259 cfg.set_self_gravity_fmm(order=5, opening_angle=0.5, reduction_level=3)
260 cfg.set_softening_plummer(epsilon=1e-9)
261
262 def sg_case_sfmm1(cfg):
263 cfg.set_self_gravity_sfmm(order=1, opening_angle=0.5, reduction_level=3)
264 cfg.set_softening_plummer(epsilon=1e-9)
265
266 def sg_case_sfmm2(cfg):
267 cfg.set_self_gravity_sfmm(order=2, opening_angle=0.5, reduction_level=3)
268 cfg.set_softening_plummer(epsilon=1e-9)
269
270 def sg_case_sfmm3(cfg):
271 cfg.set_self_gravity_sfmm(order=3, opening_angle=0.5, reduction_level=3)
272 cfg.set_softening_plummer(epsilon=1e-9)
273
274 def sg_case_sfmm4(cfg):
275 cfg.set_self_gravity_sfmm(order=4, opening_angle=0.5, reduction_level=3)
276 cfg.set_softening_plummer(epsilon=1e-9)
277
278 def sg_case_sfmm5(cfg):
279 cfg.set_self_gravity_sfmm(order=5, opening_angle=0.5, reduction_level=3)
280 cfg.set_softening_plummer(epsilon=1e-9)
281
282 no_sg_data = run_case(setup_func, setup_name, sg_case_none)
283 reference_data = run_case(setup_func, setup_name, sg_case_reference)
284
285 direct_data = run_case(setup_func, setup_name, sg_case_direct)
286 mm1_data = run_case(setup_func, setup_name, sg_case_mm1)
287 mm2_data = run_case(setup_func, setup_name, sg_case_mm2)
288 mm3_data = run_case(setup_func, setup_name, sg_case_mm3)
289 mm4_data = run_case(setup_func, setup_name, sg_case_mm4)
290 mm5_data = run_case(setup_func, setup_name, sg_case_mm5)
291 fmm1_data = run_case(setup_func, setup_name, sg_case_fmm1)
292 fmm2_data = run_case(setup_func, setup_name, sg_case_fmm2)
293 fmm3_data = run_case(setup_func, setup_name, sg_case_fmm3)
294 fmm4_data = run_case(setup_func, setup_name, sg_case_fmm4)
295 fmm5_data = run_case(setup_func, setup_name, sg_case_fmm5)
296 sfmm1_data = run_case(setup_func, setup_name, sg_case_sfmm1)
297 sfmm2_data = run_case(setup_func, setup_name, sg_case_sfmm2)
298 sfmm3_data = run_case(setup_func, setup_name, sg_case_sfmm3)
299 sfmm4_data = run_case(setup_func, setup_name, sg_case_sfmm4)
300 sfmm5_data = run_case(setup_func, setup_name, sg_case_sfmm5)
301
302 delta_sg_direct, rel_delta_direct, xyz_direct = compare_sg_methods_data(
303 no_sg_data, reference_data, direct_data
304 )
305
306 delta_sg_mm1, rel_delta_mm1, xyz_mm1 = compare_sg_methods_data(
307 no_sg_data, reference_data, mm1_data
308 )
309 delta_sg_mm2, rel_delta_mm2, xyz_mm2 = compare_sg_methods_data(
310 no_sg_data, reference_data, mm2_data
311 )
312 delta_sg_mm3, rel_delta_mm3, xyz_mm3 = compare_sg_methods_data(
313 no_sg_data, reference_data, mm3_data
314 )
315 delta_sg_mm4, rel_delta_mm4, xyz_mm4 = compare_sg_methods_data(
316 no_sg_data, reference_data, mm4_data
317 )
318 delta_sg_mm5, rel_delta_mm5, xyz_mm5 = compare_sg_methods_data(
319 no_sg_data, reference_data, mm5_data
320 )
321
322 delta_sg_fmm1, rel_delta_fmm1, xyz_fmm1 = compare_sg_methods_data(
323 no_sg_data, reference_data, fmm1_data
324 )
325 delta_sg_fmm2, rel_delta_fmm2, xyz_fmm2 = compare_sg_methods_data(
326 no_sg_data, reference_data, fmm2_data
327 )
328 delta_sg_fmm3, rel_delta_fmm3, xyz_fmm3 = compare_sg_methods_data(
329 no_sg_data, reference_data, fmm3_data
330 )
331 delta_sg_fmm4, rel_delta_fmm4, xyz_fmm4 = compare_sg_methods_data(
332 no_sg_data, reference_data, fmm4_data
333 )
334 delta_sg_fmm5, rel_delta_fmm5, xyz_fmm5 = compare_sg_methods_data(
335 no_sg_data, reference_data, fmm5_data
336 )
337
338 delta_sg_sfmm1, rel_delta_sfmm1, xyz_sfmm1 = compare_sg_methods_data(
339 no_sg_data, reference_data, sfmm1_data
340 )
341 delta_sg_sfmm2, rel_delta_sfmm2, xyz_sfmm2 = compare_sg_methods_data(
342 no_sg_data, reference_data, sfmm2_data
343 )
344 delta_sg_sfmm3, rel_delta_sfmm3, xyz_sfmm3 = compare_sg_methods_data(
345 no_sg_data, reference_data, sfmm3_data
346 )
347 delta_sg_sfmm4, rel_delta_sfmm4, xyz_sfmm4 = compare_sg_methods_data(
348 no_sg_data, reference_data, sfmm4_data
349 )
350 delta_sg_sfmm5, rel_delta_sfmm5, xyz_sfmm5 = compare_sg_methods_data(
351 no_sg_data, reference_data, sfmm5_data
352 )
353
354 check_print_errors(rel_delta_direct, setup_name, "direct", tols["direct"])
355 check_print_errors(rel_delta_mm1, setup_name, "mm1", tols["mm1"])
356 check_print_errors(rel_delta_mm2, setup_name, "mm2", tols["mm2"])
357 check_print_errors(rel_delta_mm3, setup_name, "mm3", tols["mm3"])
358 check_print_errors(rel_delta_mm4, setup_name, "mm4", tols["mm4"])
359 check_print_errors(rel_delta_mm5, setup_name, "mm5", tols["mm5"])
360 check_print_errors(rel_delta_fmm1, setup_name, "fmm1", tols["fmm1"])
361 check_print_errors(rel_delta_fmm2, setup_name, "fmm2", tols["fmm2"])
362 check_print_errors(rel_delta_fmm3, setup_name, "fmm3", tols["fmm3"])
363 check_print_errors(rel_delta_fmm4, setup_name, "fmm4", tols["fmm4"])
364 check_print_errors(rel_delta_fmm5, setup_name, "fmm5", tols["fmm5"])
365 check_print_errors(rel_delta_sfmm1, setup_name, "sfmm1", tols["sfmm1"])
366 check_print_errors(rel_delta_sfmm2, setup_name, "sfmm2", tols["sfmm2"])
367 check_print_errors(rel_delta_sfmm3, setup_name, "sfmm3", tols["sfmm3"])
368 check_print_errors(rel_delta_sfmm4, setup_name, "sfmm4", tols["sfmm4"])
369 check_print_errors(rel_delta_sfmm5, setup_name, "sfmm5", tols["sfmm5"])
370
371 return (
372 {
373 "direct": delta_sg_direct,
374 "mm1": delta_sg_mm1,
375 "mm2": delta_sg_mm2,
376 "mm3": delta_sg_mm3,
377 "mm4": delta_sg_mm4,
378 "mm5": delta_sg_mm5,
379 "fmm1": delta_sg_fmm1,
380 "fmm2": delta_sg_fmm2,
381 "fmm3": delta_sg_fmm3,
382 "fmm4": delta_sg_fmm4,
383 "fmm5": delta_sg_fmm5,
384 "sfmm1": delta_sg_sfmm1,
385 "sfmm2": delta_sg_sfmm2,
386 "sfmm3": delta_sg_sfmm3,
387 "sfmm4": delta_sg_sfmm4,
388 "sfmm5": delta_sg_sfmm5,
389 },
390 {
391 "direct": rel_delta_direct,
392 "mm1": rel_delta_mm1,
393 "mm2": rel_delta_mm2,
394 "mm3": rel_delta_mm3,
395 "mm4": rel_delta_mm4,
396 "mm5": rel_delta_mm5,
397 "fmm1": rel_delta_fmm1,
398 "fmm2": rel_delta_fmm2,
399 "fmm3": rel_delta_fmm3,
400 "fmm4": rel_delta_fmm4,
401 "fmm5": rel_delta_fmm5,
402 "sfmm1": rel_delta_sfmm1,
403 "sfmm2": rel_delta_sfmm2,
404 "sfmm3": rel_delta_sfmm3,
405 "sfmm4": rel_delta_sfmm4,
406 "sfmm5": rel_delta_sfmm5,
407 },
408 {
409 "direct": xyz_direct,
410 "mm1": xyz_mm1,
411 "mm2": xyz_mm2,
412 "mm3": xyz_mm3,
413 "mm4": xyz_mm4,
414 "mm5": xyz_mm5,
415 "fmm1": xyz_fmm1,
416 "fmm2": xyz_fmm2,
417 "fmm3": xyz_fmm3,
418 "fmm4": xyz_fmm4,
419 "fmm5": xyz_fmm5,
420 "sfmm1": xyz_sfmm1,
421 "sfmm2": xyz_sfmm2,
422 "sfmm3": xyz_sfmm3,
423 "sfmm4": xyz_sfmm4,
424 "sfmm5": xyz_sfmm5,
425 },
426 )
427
428
429 # Plot the 3D delta of the SG method
430 def plot3d_delta_sg(delta_sg_norm, xyz, case_name, method_name):
431 fig = plt.figure()
432 ax = fig.add_subplot(111, projection="3d")
433 dat = ax.scatter(xyz[:, 0], xyz[:, 1], xyz[:, 2], c=delta_sg_norm, s=1, cmap="viridis")
434 ax.set_xlabel("X")
435 ax.set_ylabel("Y")
436 ax.set_zlabel("Z")
437 ax.set_title(f"{method_name} - {case_name} relative error")
438 ax.set_aspect("equal")
439 fig.colorbar(dat)
440
441 return fig
Setup for the test
448 def setup_cube_hcp(model, cfg):
449 si = shamrock.UnitSystem()
450 sicte = shamrock.Constants(si)
451 codeu = shamrock.UnitSystem(
452 unit_time=sicte.year(),
453 unit_length=sicte.au(),
454 unit_mass=sicte.sol_mass(),
455 )
456 ucte = shamrock.Constants(codeu)
457
458 gamma = 5.0 / 3.0
459 rho_g = 100
460 initial_u = 10
461
462 sphere_radius = 0.1
463 sim_radius = 0.5
464
465 Npart = 1e4
466
467 bmin = (-sim_radius, -sim_radius, -sim_radius)
468 bmax = (sim_radius, sim_radius, sim_radius)
469
470 init_part_bmin = (-sphere_radius, -sphere_radius, -sphere_radius)
471 init_part_bmax = (sphere_radius, sphere_radius, sphere_radius)
472
473 scheduler_split_val = int(2e7)
474 scheduler_merge_val = int(1)
475
476 N_target = Npart
477 xm, ym, zm = init_part_bmin
478 xM, yM, zM = init_part_bmax
479 vol_b = (xM - xm) * (yM - ym) * (zM - zm)
480
481 if shamrock.sys.world_rank() == 0:
482 print("Npart", Npart)
483 print("scheduler_split_val", scheduler_split_val)
484 print("scheduler_merge_val", scheduler_merge_val)
485 print("N_target", N_target)
486 print("vol_b", vol_b)
487
488 part_vol = vol_b / N_target
489
490 # lattice volume
491 part_vol_lattice = 0.74 * part_vol
492
493 dr = (part_vol_lattice / ((4.0 / 3.0) * 3.1416)) ** (1.0 / 3.0)
494
495 cfg.set_artif_viscosity_VaryingCD10(
496 alpha_min=0.0, alpha_max=1, sigma_decay=0.1, alpha_u=1, beta_AV=2
497 )
498 cfg.set_boundary_periodic()
499 cfg.set_eos_adiabatic(gamma)
500
501 cfg.set_units(codeu)
502 cfg.print_status()
503 model.set_solver_config(cfg)
504 model.init_scheduler(scheduler_split_val, scheduler_merge_val)
505
506 model.resize_simulation_box(bmin, bmax)
507
508 setup = model.get_setup()
509 gen = setup.make_generator_lattice_hcp(dr, init_part_bmin, init_part_bmax)
510
511 # On aurora /2 was correct to avoid out of memory
512 setup.apply_setup(gen, insert_step=int(scheduler_split_val / 2))
513
514 vol_b = (xM - xm) * (yM - ym) * (zM - zm)
515 totmass = rho_g * vol_b
516 pmass = model.total_mass_to_part_mass(totmass)
517 model.set_particle_mass(pmass)
518
519 model.set_value_in_a_box("uint", "f64", initial_u, bmin, bmax)
520
521 model.set_cfl_cour(0.1)
522 model.set_cfl_force(0.1)
Run the tests#
529 delta_sg_dict, rel_delta_dict, xyz_dict = compare_sg_methods(
530 setup_cube_hcp, "cube_hcp", TOL_HPC_CUBE
531 )
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "none"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 2.421620e+06 N.s^-1
SPH setup: the generation step took : 0.009570488 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 6.45 us (69.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1583.00 ns (0.3%)
patch tree reduce : 901.00 ns (0.2%)
gen split merge : 852.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 891.00 ns (0.2%)
LB compute : 560.43 us (96.4%)
LB move op cnt : 0
LB apply : 12.80 us (2.2%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.75 us (65.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1793.00 ns (0.3%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.1%)
LB compute : 535.76 us (98.5%)
LB move op cnt : 0
LB apply : 2.02 us (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.60 us (64.8%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1353.00 ns (0.3%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 331.00 ns (0.1%)
LB compute : 419.23 us (98.2%)
LB move op cnt : 0
LB apply : 1904.00 ns (0.4%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.009974183000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.7% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.024593368 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.16 us (0.6%)
patch tree reduce : 541.00 ns (0.1%)
gen split merge : 431.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 381.00 ns (0.1%)
LB compute : 501.46 us (98.0%)
LB move op cnt : 0
LB apply : 1984.00 ns (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.62 us (69.6%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (1.8465799040655867e-14,-2.5967529900922313e-15,-1.4426405500512397e-15)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 3.2268505544103883e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 3.2268505544103883e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7172e+04 | 10087 | 1 | 3.712e-01 | 0.0% | 0.8% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"reference_mode": true,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "direct"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.674005e+07 N.s^-1
SPH setup: the generation step took : 0.001200597 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.01 us (69.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1543.00 ns (0.4%)
patch tree reduce : 941.00 ns (0.2%)
gen split merge : 692.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1242.00 ns (0.3%)
LB compute : 372.15 us (97.0%)
LB move op cnt : 0
LB apply : 3.56 us (0.9%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.81 us (63.3%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1443.00 ns (0.4%)
patch tree reduce : 401.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.1%)
LB compute : 387.44 us (98.0%)
LB move op cnt : 0
LB apply : 2.04 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 (66.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1953.00 ns (0.5%)
patch tree reduce : 471.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 340.00 ns (0.1%)
LB compute : 424.47 us (98.1%)
LB move op cnt : 0
LB apply : 2.00 us (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006995309000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.4% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.01081699 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.06 us (0.7%)
patch tree reduce : 491.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 360.00 ns (0.1%)
LB compute : 402.38 us (97.6%)
LB move op cnt : 0
LB apply : 1984.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.50 us (70.1%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-2.5967529900922314e-13,1.788874282063537e-13,2.2505192580799338e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.458107303830286e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.458107303830286e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 1.0830e+04 | 10087 | 1 | 9.314e-01 | 0.0% | 0.3% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"reference_mode": false,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "direct"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 9.738761e+06 N.s^-1
SPH setup: the generation step took : 0.001508222 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.00 us (68.3%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1423.00 ns (0.4%)
patch tree reduce : 882.00 ns (0.2%)
gen split merge : 812.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 952.00 ns (0.2%)
LB compute : 392.45 us (97.3%)
LB move op cnt : 0
LB apply : 3.17 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.83 us (65.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1443.00 ns (0.4%)
patch tree reduce : 441.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 310.00 ns (0.1%)
LB compute : 349.32 us (97.9%)
LB move op cnt : 0
LB apply : 1963.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.32 us (68.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1202.00 ns (0.4%)
patch tree reduce : 410.00 ns (0.1%)
gen split merge : 410.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 310.00 ns (0.1%)
LB compute : 318.56 us (97.8%)
LB move op cnt : 0
LB apply : 1773.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.00687246 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.5% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.010991045000000001 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.97 us (0.9%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 336.51 us (97.2%)
LB move op cnt : 0
LB apply : 1934.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.44 us (71.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-2.5967529900922314e-13,1.788874282063537e-13,2.2505192580799338e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.458107303830286e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.458107303830286e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 1.6407e+04 | 10087 | 1 | 6.148e-01 | 0.0% | 0.4% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 1,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "mm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 8.177358e+06 N.s^-1
SPH setup: the generation step took : 0.0017808210000000001 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.89 us (68.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1362.00 ns (0.4%)
patch tree reduce : 751.00 ns (0.2%)
gen split merge : 811.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 902.00 ns (0.3%)
LB compute : 322.63 us (96.9%)
LB move op cnt : 0
LB apply : 3.18 us (1.0%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.31 us (68.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.87 us (1.3%)
patch tree reduce : 1934.00 ns (0.6%)
gen split merge : 390.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1132.00 ns (0.4%)
LB compute : 291.63 us (96.0%)
LB move op cnt : 0
LB apply : 1643.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.14 us (66.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1392.00 ns (0.4%)
patch tree reduce : 321.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.1%)
LB compute : 303.29 us (97.7%)
LB move op cnt : 0
LB apply : 1723.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006838935000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.5% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.011365353 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.04 us (0.8%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 431.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 360.00 ns (0.1%)
LB compute : 384.25 us (97.5%)
LB move op cnt : 0
LB apply : 1974.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.62 us (73.2%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-0.03597005929433212,0.752719449000333,0.11263499641180232)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4384953098798583e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4384953098798583e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7853e+04 | 10087 | 1 | 3.621e-01 | 0.0% | 0.7% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 2,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "mm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.460546e+07 N.s^-1
SPH setup: the generation step took : 0.00119232 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.67 us (67.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1643.00 ns (0.4%)
patch tree reduce : 811.00 ns (0.2%)
gen split merge : 622.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 971.00 ns (0.3%)
LB compute : 373.21 us (97.1%)
LB move op cnt : 0
LB apply : 3.49 us (0.9%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.62 us (64.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1292.00 ns (0.4%)
patch tree reduce : 491.00 ns (0.2%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 301.00 ns (0.1%)
LB compute : 313.65 us (97.8%)
LB move op cnt : 0
LB apply : 1764.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.17 us (64.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1253.00 ns (0.4%)
patch tree reduce : 411.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.1%)
LB compute : 309.06 us (97.7%)
LB move op cnt : 0
LB apply : 1814.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006640634 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.010439862000000001 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.96 us (0.9%)
patch tree reduce : 481.00 ns (0.2%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 350.00 ns (0.1%)
LB compute : 310.70 us (97.1%)
LB move op cnt : 0
LB apply : 1823.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.40 us (69.7%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-0.03192867068634885,0.46157601590495817,0.03227743556291892)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4374467718013695e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4374467718013695e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.8024e+04 | 10087 | 1 | 3.599e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 3,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "mm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.267772e+07 N.s^-1
SPH setup: the generation step took : 0.001233327 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.26 us (63.8%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1562.00 ns (0.5%)
patch tree reduce : 822.00 ns (0.2%)
gen split merge : 881.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 841.00 ns (0.2%)
LB compute : 335.43 us (96.8%)
LB move op cnt : 0
LB apply : 2.69 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.26 us (67.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1282.00 ns (0.4%)
patch tree reduce : 401.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 : 316.03 us (97.8%)
LB move op cnt : 0
LB apply : 1804.00 ns (0.6%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.36 us (68.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1222.00 ns (0.4%)
patch tree reduce : 451.00 ns (0.2%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 291.00 ns (0.1%)
LB compute : 290.20 us (97.7%)
LB move op cnt : 0
LB apply : 1643.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006268177000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.01028908 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.03 us (0.9%)
patch tree reduce : 521.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 371.00 ns (0.1%)
LB compute : 341.54 us (97.2%)
LB move op cnt : 0
LB apply : 1853.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.38 us (70.6%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.04492773368219582,0.38065217229475995,-0.03294442418733314)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.450328726610016e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.450328726610016e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7575e+04 | 10087 | 1 | 3.658e-01 | 0.0% | 0.7% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 4,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "mm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.484494e+07 N.s^-1
SPH setup: the generation step took : 0.0018557320000000001 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.93 us (69.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1422.00 ns (0.4%)
patch tree reduce : 791.00 ns (0.2%)
gen split merge : 682.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 922.00 ns (0.3%)
LB compute : 334.75 us (97.0%)
LB move op cnt : 0
LB apply : 3.29 us (1.0%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.27 us (67.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1293.00 ns (0.4%)
patch tree reduce : 381.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 290.00 ns (0.1%)
LB compute : 313.46 us (97.9%)
LB move op cnt : 0
LB apply : 1644.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.22 us (66.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1292.00 ns (0.4%)
patch tree reduce : 360.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 291.00 ns (0.1%)
LB compute : 317.64 us (97.8%)
LB move op cnt : 0
LB apply : 1763.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.00678309 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.011265637 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.88 us (0.9%)
patch tree reduce : 521.00 ns (0.2%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 318.60 us (97.1%)
LB move op cnt : 0
LB apply : 1963.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.53 us (68.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.062141306734179426,0.07186267864905517,-0.11227436639894799)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4486430391159233e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4486430391159233e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7342e+04 | 10087 | 1 | 3.689e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 5,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "mm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.648485e+07 N.s^-1
SPH setup: the generation step took : 0.0010672770000000002 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.54 us (71.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1703.00 ns (0.5%)
patch tree reduce : 861.00 ns (0.3%)
gen split merge : 672.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 1032.00 ns (0.3%)
LB compute : 319.28 us (96.7%)
LB move op cnt : 0
LB apply : 3.31 us (1.0%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.17 us (67.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1242.00 ns (0.4%)
patch tree reduce : 451.00 ns (0.2%)
gen split merge : 380.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 290.63 us (97.7%)
LB move op cnt : 0
LB apply : 1643.00 ns (0.6%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.31 us (66.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1233.00 ns (0.4%)
patch tree reduce : 431.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.1%)
LB compute : 281.98 us (97.6%)
LB move op cnt : 0
LB apply : 1714.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006612073 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.010490429 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.06 us (0.9%)
patch tree reduce : 511.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 371.00 ns (0.1%)
LB compute : 333.19 us (97.2%)
LB move op cnt : 0
LB apply : 1793.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.34 us (70.1%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.014077502276772519,-0.1098138243846681,0.006886515733428269)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4508458811090976e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4508458811090976e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.6691e+04 | 10087 | 1 | 3.779e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 1,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 9.566732e+06 N.s^-1
SPH setup: the generation step took : 0.004063051000000001 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.27 us (65.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.01 us (0.6%)
patch tree reduce : 762.00 ns (0.2%)
gen split merge : 641.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 812.00 ns (0.2%)
LB compute : 328.07 us (96.8%)
LB move op cnt : 0
LB apply : 3.11 us (0.9%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.41 us (69.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1252.00 ns (0.4%)
patch tree reduce : 380.00 ns (0.1%)
gen split merge : 361.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.1%)
LB compute : 307.31 us (97.8%)
LB move op cnt : 0
LB apply : 1583.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.34 us (63.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1253.00 ns (0.4%)
patch tree reduce : 401.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.1%)
LB compute : 282.64 us (97.6%)
LB move op cnt : 0
LB apply : 1733.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006726003 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.2% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.013385974 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.04 us (0.9%)
patch tree reduce : 541.00 ns (0.2%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 336.44 us (97.3%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.34 us (70.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.30467386170117133,0.2920960474805367,-0.9049382257560593)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4273068424046396e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4273068424046396e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7845e+04 | 10087 | 1 | 3.623e-01 | 0.0% | 0.7% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 2,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.506280e+07 N.s^-1
SPH setup: the generation step took : 0.0028389720000000004 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.74 us (68.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1483.00 ns (0.4%)
patch tree reduce : 1202.00 ns (0.3%)
gen split merge : 651.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 912.00 ns (0.3%)
LB compute : 352.03 us (97.0%)
LB move op cnt : 0
LB apply : 3.01 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.37 us (68.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1252.00 ns (0.4%)
patch tree reduce : 391.00 ns (0.1%)
gen split merge : 370.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.1%)
LB compute : 296.40 us (97.7%)
LB move op cnt : 0
LB apply : 1733.00 ns (0.6%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.25 us (67.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1353.00 ns (0.4%)
patch tree reduce : 410.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 301.00 ns (0.1%)
LB compute : 328.36 us (97.8%)
LB move op cnt : 0
LB apply : 1783.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006603745 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.012011623 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.90 us (0.9%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 371.00 ns (0.1%)
LB compute : 330.80 us (97.2%)
LB move op cnt : 0
LB apply : 1963.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.39 us (70.9%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.25573720161271046,0.2762846937078224,-0.9841687710942018)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4470077486491645e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4470077486491645e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7979e+04 | 10087 | 1 | 3.605e-01 | 0.0% | 0.7% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 3,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.587057e+07 N.s^-1
SPH setup: the generation step took : 0.0025714220000000003 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.12 us (66.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1423.00 ns (0.2%)
patch tree reduce : 912.00 ns (0.2%)
gen split merge : 671.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 821.00 ns (0.1%)
LB compute : 583.72 us (98.2%)
LB move op cnt : 0
LB apply : 3.18 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.69 us (66.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1242.00 ns (0.4%)
patch tree reduce : 501.00 ns (0.2%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 320.00 ns (0.1%)
LB compute : 293.75 us (97.5%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.6%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.13 us (66.2%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1152.00 ns (0.4%)
patch tree reduce : 391.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 302.41 us (97.7%)
LB move op cnt : 0
LB apply : 1743.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.0067852550000000005 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.2% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.011893252 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.97 us (0.9%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 370.00 ns (0.1%)
LB compute : 324.98 us (97.2%)
LB move op cnt : 0
LB apply : 1794.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.58 us (69.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-0.16115349478998361,-0.11515753811449056,-0.23452529682196394)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4547816995790715e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4547816995790715e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7913e+04 | 10087 | 1 | 3.614e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 4,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.770500e+07 N.s^-1
SPH setup: the generation step took : 0.00098328 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.56 us (67.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1863.00 ns (0.5%)
patch tree reduce : 992.00 ns (0.3%)
gen split merge : 681.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 912.00 ns (0.2%)
LB compute : 378.15 us (97.1%)
LB move op cnt : 0
LB apply : 3.25 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.60 us (65.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1413.00 ns (0.4%)
patch tree reduce : 390.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.1%)
LB compute : 358.96 us (97.9%)
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.48 us (65.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1272.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 : 311.00 ns (0.1%)
LB compute : 305.34 us (97.7%)
LB move op cnt : 0
LB apply : 1683.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.007019984000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.2% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.010603659000000001 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.25 us (1.0%)
patch tree reduce : 561.00 ns (0.2%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 380.00 ns (0.1%)
LB compute : 316.60 us (97.0%)
LB move op cnt : 0
LB apply : 1774.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.65 us (70.2%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-0.17542776055539155,-0.12132397085809142,-0.22366739669452143)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4544547001222264e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4544547001222264e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.7506e+04 | 10087 | 1 | 3.667e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"opening_angle": 0.5,
"order": 5,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "fmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.836745e+07 N.s^-1
SPH setup: the generation step took : 0.002622146 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.27 us (61.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1403.00 ns (0.4%)
patch tree reduce : 812.00 ns (0.2%)
gen split merge : 892.00 ns (0.3%)
split / merge op : 0/0
apply split merge : 932.00 ns (0.3%)
LB compute : 325.64 us (96.9%)
LB move op cnt : 0
LB apply : 3.11 us (0.9%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.25 us (67.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1262.00 ns (0.4%)
patch tree reduce : 351.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.1%)
LB compute : 316.90 us (97.8%)
LB move op cnt : 0
LB apply : 1783.00 ns (0.6%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.19 us (67.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1713.00 ns (0.6%)
patch tree reduce : 401.00 ns (0.1%)
gen split merge : 361.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.1%)
LB compute : 297.00 us (97.6%)
LB move op cnt : 0
LB apply : 1553.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.006632379000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.3% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.011893752 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.13 us (1.0%)
patch tree reduce : 511.00 ns (0.2%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 307.81 us (97.0%)
LB move op cnt : 0
LB apply : 1833.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (69.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (0.020950059647177183,-0.005001869815935811,-0.008611166154152295)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.458022683047013e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.458022683047013e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.6848e+04 | 10087 | 1 | 3.757e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"leaf_lowering": true,
"opening_angle": 0.5,
"order": 1,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "sfmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.264285e+07 N.s^-1
SPH setup: the generation step took : 0.0029929590000000003 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.03 us (67.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1833.00 ns (0.4%)
patch tree reduce : 922.00 ns (0.2%)
gen split merge : 671.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 922.00 ns (0.2%)
LB compute : 424.37 us (97.4%)
LB move op cnt : 0
LB apply : 3.29 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.52 us (65.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1292.00 ns (0.3%)
patch tree reduce : 500.00 ns (0.1%)
gen split merge : 410.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 300.00 ns (0.1%)
LB compute : 362.65 us (98.0%)
LB move op cnt : 0
LB apply : 1904.00 ns (0.5%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.40 us (65.3%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
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 : 390.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 291.00 ns (0.1%)
LB compute : 379.56 us (98.0%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.5%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.0071283660000000006 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.2% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.012664615 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.10 us (1.0%)
patch tree reduce : 561.00 ns (0.2%)
gen split merge : 430.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 390.00 ns (0.1%)
LB compute : 309.95 us (97.0%)
LB move op cnt : 0
LB apply : 1813.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.35 us (69.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
SPH M2L interact count: 128994
SPH P2P interact count: 62733
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (8.540432056303338e-14,4.616449760163967e-14,-9.232899520327933e-15)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.417582508264753e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.417582508264753e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.3186e+04 | 10087 | 1 | 4.350e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"leaf_lowering": true,
"opening_angle": 0.5,
"order": 2,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "sfmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 1.280460e+07 N.s^-1
SPH setup: the generation step took : 0.003129034 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.17 us (67.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1523.00 ns (0.4%)
patch tree reduce : 831.00 ns (0.2%)
gen split merge : 681.00 ns (0.2%)
split / merge op : 0/0
apply split merge : 942.00 ns (0.3%)
LB compute : 344.63 us (97.1%)
LB move op cnt : 0
LB apply : 2.98 us (0.8%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (65.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1472.00 ns (0.4%)
patch tree reduce : 391.00 ns (0.1%)
gen split merge : 370.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 290.00 ns (0.1%)
LB compute : 334.29 us (97.8%)
LB move op cnt : 0
LB apply : 1814.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.44 us (67.6%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1363.00 ns (0.4%)
patch tree reduce : 461.00 ns (0.1%)
gen split merge : 371.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 310.00 ns (0.1%)
LB compute : 310.51 us (97.7%)
LB move op cnt : 0
LB apply : 1794.00 ns (0.6%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.007987985000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.2% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.015564991 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.25 us (0.9%)
patch tree reduce : 511.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 336.81 us (97.2%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.44 us (67.8%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
SPH M2L interact count: 128994
SPH P2P interact count: 62733
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (1.9619911480696857e-14,4.616449760163967e-14,1.2118180620430413e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4405154618237877e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4405154618237877e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.1869e+04 | 10087 | 1 | 4.612e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"leaf_lowering": true,
"opening_angle": 0.5,
"order": 3,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "sfmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 3.547900e+06 N.s^-1
SPH setup: the generation step took : 0.005822142000000001 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.44 us (61.9%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1813.00 ns (0.2%)
patch tree reduce : 841.00 ns (0.1%)
gen split merge : 751.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1052.00 ns (0.1%)
LB compute : 754.03 us (98.5%)
LB move op cnt : 0
LB apply : 3.14 us (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.79 us (68.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1452.00 ns (0.2%)
patch tree reduce : 431.00 ns (0.1%)
gen split merge : 410.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 320.00 ns (0.0%)
LB compute : 829.05 us (99.1%)
LB move op cnt : 0
LB apply : 1763.00 ns (0.2%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.31 us (64.7%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1253.00 ns (0.1%)
patch tree reduce : 521.00 ns (0.1%)
gen split merge : 401.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.0%)
LB compute : 831.16 us (99.1%)
LB move op cnt : 0
LB apply : 1944.00 ns (0.2%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.011150963 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 0.7% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.020816457 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.12 us (0.9%)
patch tree reduce : 530.00 ns (0.2%)
gen split merge : 420.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.1%)
LB compute : 329.63 us (97.2%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.5%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.46 us (68.5%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
SPH M2L interact count: 128994
SPH P2P interact count: 62733
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (3.577748564127074e-14,-4.616449760163967e-14,8.021081458284892e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.4550924945738166e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.4550924945738166e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.2012e+04 | 10087 | 1 | 4.583e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"leaf_lowering": true,
"opening_angle": 0.5,
"order": 4,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "sfmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 8.353374e+06 N.s^-1
SPH setup: the generation step took : 0.0022298400000000003 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 5.15 us (69.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1443.00 ns (0.2%)
patch tree reduce : 792.00 ns (0.1%)
gen split merge : 611.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 801.00 ns (0.1%)
LB compute : 724.37 us (98.5%)
LB move op cnt : 0
LB apply : 3.42 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.38 us (69.4%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1473.00 ns (0.2%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 390.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.0%)
LB compute : 795.79 us (99.1%)
LB move op cnt : 0
LB apply : 1844.00 ns (0.2%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.22 us (67.5%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1252.00 ns (0.2%)
patch tree reduce : 490.00 ns (0.1%)
gen split merge : 390.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 301.00 ns (0.0%)
LB compute : 777.48 us (99.1%)
LB move op cnt : 0
LB apply : 1804.00 ns (0.2%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.011018505000000001 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 1.9% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.01758068 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.00 us (0.6%)
patch tree reduce : 500.00 ns (0.1%)
gen split merge : 420.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 500.72 us (98.1%)
LB move op cnt : 0
LB apply : 1894.00 ns (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.58 us (68.3%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
SPH M2L interact count: 128994
SPH P2P interact count: 62733
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-2.077402392073785e-14,-2.654458612094281e-14,-3.750865430133223e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.457081216599208e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.457081216599208e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.1577e+04 | 10087 | 1 | 4.675e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
Npart 10000.0
scheduler_split_val 20000000
scheduler_merge_val 1
N_target 10000.0
vol_b 0.008000000000000002
----- SPH Solver configuration -----
[
{
"artif_viscosity": {
"alpha_max": 1.0,
"alpha_min": 0.0,
"alpha_u": 1.0,
"av_type": "varying_cd10",
"beta_AV": 2.0,
"sigma_decay": 0.1
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"enable_particle_reordering": false,
"eos_config": {
"Tvec": "f64_3",
"eos_type": "adiabatic",
"gamma": 1.6666666666666667
},
"epsilon_h": 1e-06,
"ext_force_config": {
"force_list": []
},
"gpart_mass": 0.0,
"h_iter_per_subcycles": 50,
"h_max_subcycles_count": 100,
"htol_up_coarse_cycle": 1.1,
"htol_up_fine_cycle": 1.1,
"kernel_id": "M4<f64>",
"mhd_config": {
"mhd_type": "none"
},
"particle_killing": [],
"particle_reordering_step_freq": 1000,
"save_dt_to_fields": false,
"self_grav_config": {
"leaf_lowering": true,
"opening_angle": 0.5,
"order": 5,
"reduction_level": 3,
"softening_length": 1e-09,
"softening_mode": "plummer",
"type": "sfmm"
},
"show_ghost_zone_graph": false,
"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
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 ) rate = 8.517009e+06 N.s^-1
SPH setup: the generation step took : 0.004038206 s
SPH setup: final particle count = 10087 begining injection ...
Info: --------------------------------------------- [DataInserterUtility][rank=0]
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 4.78 us (65.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1492.00 ns (0.2%)
patch tree reduce : 862.00 ns (0.1%)
gen split merge : 662.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 922.00 ns (0.1%)
LB compute : 821.26 us (98.7%)
LB move op cnt : 0
LB apply : 3.58 us (0.4%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.50 us (69.0%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1342.00 ns (0.2%)
patch tree reduce : 441.00 ns (0.1%)
gen split merge : 371.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 290.00 ns (0.0%)
LB compute : 845.48 us (99.1%)
LB move op cnt : 0
LB apply : 1753.00 ns (0.2%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.22 us (67.1%)
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1243.00 ns (0.2%)
patch tree reduce : 410.00 ns (0.1%)
gen split merge : 391.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 261.00 ns (0.0%)
LB compute : 789.69 us (99.1%)
LB move op cnt : 0
LB apply : 1753.00 ns (0.2%)
Info: --------------------------------------------- [DataInserterUtility][rank=0]
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop ->
SPH setup: the injection step took : 0.011035948 s
Info: injection perf report: [SPH setup][rank=0]
+======+====================+=======+=============+=============+=============+
| rank | rank get (sum/max) | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+====================+=======+=============+=============+=============+
| 0 | 0.00s / 0.00s | 0.00s | 0.8% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.019819341 s
---------------- t = 0, dt = 0 ----------------
Info: summary : [LoadBalance][rank=0]
Info: - strategy "psweep" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: - strategy "round robin" : max = 10087 min = 10087 [LoadBalance][rank=0]
Info: Loadbalance stats : [LoadBalance][rank=0]
npatch = 1
min = 10087
max = 10087
avg = 10087
efficiency = 100.00%
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 3.03 us (0.6%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 361.00 ns (0.1%)
LB compute : 527.42 us (98.2%)
LB move op cnt : 0
LB apply : 1954.00 ns (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.57 us (70.4%)
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.005729763709954786 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006302740080950265 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.006933014089045292 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.007626315497949822 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.008388947047744805 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.009227841752519287 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.010150625927771216 unconverged cnt = 10087
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.011165688520548339 unconverged cnt = 3970
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.012282257372603174 unconverged cnt = 1758
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.013510483109863492 unconverged cnt = 182
Warning: smoothing length is not converged, rerunning the iterator ... [Smoothinglength][rank=0]
largest h = 0.014861531420849842 unconverged cnt = 52
SPH M2L interact count: 128994
SPH P2P interact count: 62733
Info: conservation infos : [sph::Model][rank=0]
sum v = (0,0,0)
sum a = (-8.194198324291041e-14,2.4236361240860825e-14,1.3272293060471405e-14)
sum e = 8.000000000000002
sum de = 0
Info: CFL hydro = 2.45816497524108e-06 sink sink = inf [SPH][rank=0]
Info: cfl dt = 2.45816497524108e-06 cfl multiplier : 0.01 [sph::Model][rank=0]
Info: Timestep perf report: [sph::Model][rank=0]
+======+============+=======+========+===========+======+=============+=============+=============+
| rank | rate (N/s) | Nobj | Npatch | tstep | MPI | alloc d% h% | mem (max) d | mem (max) h |
+======+============+=======+========+===========+======+=============+=============+=============+
| 0 | 2.0264e+04 | 10087 | 1 | 4.978e-01 | 0.0% | 0.6% 0.0% | 1.16 GB | 5.04 MB |
+------+------------+-------+--------+-----------+------+-------------+-------------+-------------+
Info: estimated rate : 0 (tsim/hr) [sph::Model][rank=0]
Info: collected : 1 patches [PatchScheduler][rank=0]
max relative error direct: 0.0 for cube_hcp
avg relative error direct: 0.0 for cube_hcp
min relative error direct: 0.0 for cube_hcp
std relative error direct: 0.0 for cube_hcp
max relative error mm1: 0.10237644408204856 for cube_hcp
avg relative error mm1: 0.028547050794827368 for cube_hcp
min relative error mm1: 0.001185552870420316 for cube_hcp
std relative error mm1: 0.012055180964609495 for cube_hcp
max relative error mm2: 0.09870604133669349 for cube_hcp
avg relative error mm2: 0.02809018282580241 for cube_hcp
min relative error mm2: 0.0007660948531796988 for cube_hcp
std relative error mm2: 0.011725871267882113 for cube_hcp
max relative error mm3: 0.060028282345646375 for cube_hcp
avg relative error mm3: 0.01608490739010445 for cube_hcp
min relative error mm3: 0.0006196291223047618 for cube_hcp
std relative error mm3: 0.00768369412614662 for cube_hcp
max relative error mm4: 0.05833850788436179 for cube_hcp
avg relative error mm4: 0.014911420251122888 for cube_hcp
min relative error mm4: 0.0009000262033407927 for cube_hcp
std relative error mm4: 0.007058594416839622 for cube_hcp
max relative error mm5: 0.03130228236510673 for cube_hcp
avg relative error mm5: 0.005542488773510675 for cube_hcp
min relative error mm5: 0.00011105937841265464 for cube_hcp
std relative error mm5: 0.0035758340149414244 for cube_hcp
max relative error fmm1: 0.09457251908559332 for cube_hcp
avg relative error fmm1: 0.025849631024874815 for cube_hcp
min relative error fmm1: 0.0009107123147506598 for cube_hcp
std relative error fmm1: 0.014771379187304317 for cube_hcp
max relative error fmm2: 0.06378841575326645 for cube_hcp
avg relative error fmm2: 0.01197538537799122 for cube_hcp
min relative error fmm2: 0.0003606563787785033 for cube_hcp
std relative error fmm2: 0.0056847862791565874 for cube_hcp
max relative error fmm3: 0.041451600220482575 for cube_hcp
avg relative error fmm3: 0.005904750341965699 for cube_hcp
min relative error fmm3: 0.00011353264639237627 for cube_hcp
std relative error fmm3: 0.003552441879243517 for cube_hcp
max relative error fmm4: 0.02704100823427953 for cube_hcp
avg relative error fmm4: 0.0048700473223711035 for cube_hcp
min relative error fmm4: 9.263728962019831e-05 for cube_hcp
std relative error fmm4: 0.002852652151208191 for cube_hcp
max relative error fmm5: 0.013154949137015834 for cube_hcp
avg relative error fmm5: 0.0019333759110011861 for cube_hcp
min relative error fmm5: 3.7696982567621314e-05 for cube_hcp
std relative error fmm5: 0.0014632873063487252 for cube_hcp
max relative error sfmm1: 0.21892380131546468 for cube_hcp
avg relative error sfmm1: 0.057829039225082396 for cube_hcp
min relative error sfmm1: 0.0013277688426966254 for cube_hcp
std relative error sfmm1: 0.030787601102585056 for cube_hcp
max relative error sfmm2: 0.03804965294061314 for cube_hcp
avg relative error sfmm2: 0.011068208855935383 for cube_hcp
min relative error sfmm2: 0.0005623871943127493 for cube_hcp
std relative error sfmm2: 0.005376106489515813 for cube_hcp
max relative error sfmm3: 0.016938241201078053 for cube_hcp
avg relative error sfmm3: 0.004725763495042012 for cube_hcp
min relative error sfmm3: 0.00011119860215946636 for cube_hcp
std relative error sfmm3: 0.0026446593252427626 for cube_hcp
max relative error sfmm4: 0.011356406701672126 for cube_hcp
avg relative error sfmm4: 0.002287613747357328 for cube_hcp
min relative error sfmm4: 5.847968796672419e-05 for cube_hcp
std relative error sfmm4: 0.0015169692087367923 for cube_hcp
max relative error sfmm5: 0.010203458622545838 for cube_hcp
avg relative error sfmm5: 0.001413447866169747 for cube_hcp
min relative error sfmm5: 3.503562542590484e-05 for cube_hcp
std relative error sfmm5: 0.0012180716522390488 for cube_hcp
534 fig = plot3d_delta_sg(rel_delta_dict["direct"], xyz_dict["direct"], "cube_hcp", "direct")
535 plt.show()
538 fig = plot3d_delta_sg(rel_delta_dict["mm1"], xyz_dict["mm1"], "cube_hcp", "mm1")
539 plt.show()
542 fig = plot3d_delta_sg(rel_delta_dict["mm2"], xyz_dict["mm2"], "cube_hcp", "mm2")
543 plt.show()
546 fig = plot3d_delta_sg(rel_delta_dict["mm3"], xyz_dict["mm3"], "cube_hcp", "mm3")
547 plt.show()
550 fig = plot3d_delta_sg(rel_delta_dict["mm4"], xyz_dict["mm4"], "cube_hcp", "mm4")
551 plt.show()
554 fig = plot3d_delta_sg(rel_delta_dict["mm5"], xyz_dict["mm5"], "cube_hcp", "mm5")
555 plt.show()
558 fig = plot3d_delta_sg(rel_delta_dict["fmm1"], xyz_dict["fmm1"], "cube_hcp", "fmm1")
559 plt.show()
562 fig = plot3d_delta_sg(rel_delta_dict["fmm2"], xyz_dict["fmm2"], "cube_hcp", "fmm2")
563 plt.show()
566 fig = plot3d_delta_sg(rel_delta_dict["fmm3"], xyz_dict["fmm3"], "cube_hcp", "fmm3")
567 plt.show()
570 fig = plot3d_delta_sg(rel_delta_dict["fmm4"], xyz_dict["fmm4"], "cube_hcp", "fmm4")
571 plt.show()
574 fig = plot3d_delta_sg(rel_delta_dict["fmm5"], xyz_dict["fmm5"], "cube_hcp", "fmm5")
575 plt.show()
579 fig = plot3d_delta_sg(rel_delta_dict["sfmm1"], xyz_dict["sfmm1"], "cube_hcp", "sfmm1")
580 plt.show()
583 fig = plot3d_delta_sg(rel_delta_dict["sfmm2"], xyz_dict["sfmm2"], "cube_hcp", "sfmm2")
584 plt.show()
587 fig = plot3d_delta_sg(rel_delta_dict["sfmm3"], xyz_dict["sfmm3"], "cube_hcp", "sfmm3")
588 plt.show()
591 fig = plot3d_delta_sg(rel_delta_dict["sfmm4"], xyz_dict["sfmm4"], "cube_hcp", "sfmm4")
592 plt.show()
595 fig = plot3d_delta_sg(rel_delta_dict["sfmm5"], xyz_dict["sfmm5"], "cube_hcp", "sfmm5")
596 plt.show()
Report errors
Total running time of the script: (0 minutes 39.097 seconds)
Estimated memory usage: 160 MB