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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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.342761e+06 N.s^-1
SPH setup: the generation step took : 0.010730926 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.84 us (65.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 : 1804.00 ns (0.0%)
patch tree reduce : 1573.00 ns (0.0%)
gen split merge : 782.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 992.00 ns (0.0%)
LB compute : 9.31 ms (99.9%)
LB move op cnt : 0
LB apply : 3.60 us (0.0%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.96 us (66.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 : 1593.00 ns (0.1%)
patch tree reduce : 491.00 ns (0.0%)
gen split merge : 430.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 320.00 ns (0.0%)
LB compute : 1341.93 us (99.4%)
LB move op cnt : 0
LB apply : 1843.00 ns (0.1%)
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 (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 : 1714.00 ns (0.1%)
patch tree reduce : 441.00 ns (0.0%)
gen split merge : 351.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.0%)
LB compute : 1271.70 us (99.4%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.1%)
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.021563116 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.043964909 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.33 us (0.4%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.0%)
LB compute : 761.66 us (98.7%)
LB move op cnt : 0
LB apply : 1824.00 ns (0.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.58 us (65.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 = (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 | 1.4647e+04 | 10087 | 1 | 6.887e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 3.662704e+06 N.s^-1
SPH setup: the generation step took : 0.004607754 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.27 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 : 1653.00 ns (0.2%)
patch tree reduce : 1022.00 ns (0.1%)
gen split merge : 671.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 922.00 ns (0.1%)
LB compute : 930.46 us (98.8%)
LB move op cnt : 0
LB apply : 3.11 us (0.3%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.46 us (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 : 1533.00 ns (0.2%)
patch tree reduce : 530.00 ns (0.1%)
gen split merge : 380.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.0%)
LB compute : 988.80 us (99.3%)
LB move op cnt : 0
LB apply : 1663.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.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 : 1574.00 ns (0.2%)
patch tree reduce : 441.00 ns (0.0%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.0%)
LB compute : 882.12 us (99.2%)
LB move op cnt : 0
LB apply : 1613.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.010765779000000001 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.020351436 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.17 us (0.5%)
patch tree reduce : 521.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.1%)
LB compute : 636.10 us (98.5%)
LB move op cnt : 0
LB apply : 1803.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.22 us (67.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
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 | 8.3245e+03 | 10087 | 1 | 1.212e+00 | 0.0% | 0.2% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 1.711898e+07 N.s^-1
SPH setup: the generation step took : 0.002660061 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.41 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 : 1643.00 ns (0.2%)
patch tree reduce : 921.00 ns (0.1%)
gen split merge : 681.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 952.00 ns (0.1%)
LB compute : 722.77 us (98.5%)
LB move op cnt : 0
LB apply : 3.45 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.20 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 : 1192.00 ns (0.1%)
patch tree reduce : 430.00 ns (0.0%)
gen split merge : 391.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 271.00 ns (0.0%)
LB compute : 880.33 us (99.3%)
LB move op cnt : 0
LB apply : 1552.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.28 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 : 1213.00 ns (0.1%)
patch tree reduce : 461.00 ns (0.1%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 261.00 ns (0.0%)
LB compute : 862.49 us (99.3%)
LB move op cnt : 0
LB apply : 1533.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.008758515 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.016612062 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.26 us (0.6%)
patch tree reduce : 510.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 : 542.14 us (98.3%)
LB move op cnt : 0
LB apply : 1643.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1943.00 ns (66.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 = (-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.0984e+04 | 10087 | 1 | 9.183e-01 | 0.0% | 0.2% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 1.899373e+07 N.s^-1
SPH setup: the generation step took : 0.002376081 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.94 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 : 1653.00 ns (0.4%)
patch tree reduce : 1052.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 : 386.44 us (97.2%)
LB move op cnt : 0
LB apply : 3.15 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.50 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 : 1302.00 ns (0.4%)
patch tree reduce : 561.00 ns (0.2%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 330.00 ns (0.1%)
LB compute : 300.21 us (97.8%)
LB move op cnt : 0
LB apply : 1513.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.16 us (64.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 : 1854.00 ns (0.6%)
patch tree reduce : 451.00 ns (0.1%)
gen split merge : 381.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.1%)
LB compute : 306.27 us (97.7%)
LB move op cnt : 0
LB apply : 1612.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.013574767000000001 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.01s / 0.01s | 0.00s | 0.4% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.020030698 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.14 us (0.5%)
patch tree reduce : 531.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.0%)
LB compute : 682.50 us (98.6%)
LB move op cnt : 0
LB apply : 1844.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.56 us (67.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
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 | 1.5482e+04 | 10087 | 1 | 6.515e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 3.754460e+06 N.s^-1
SPH setup: the generation step took : 0.004708511 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.51 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 : 1623.00 ns (0.2%)
patch tree reduce : 1062.00 ns (0.1%)
gen split merge : 691.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1302.00 ns (0.1%)
LB compute : 927.68 us (98.8%)
LB move op cnt : 0
LB apply : 3.08 us (0.3%)
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 (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 : 2.08 us (0.2%)
patch tree reduce : 431.00 ns (0.0%)
gen split merge : 390.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 331.00 ns (0.0%)
LB compute : 1021.09 us (99.2%)
LB move op cnt : 0
LB apply : 1693.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 (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 : 1594.00 ns (0.2%)
patch tree reduce : 481.00 ns (0.0%)
gen split merge : 371.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 1047.85 us (99.3%)
LB move op cnt : 0
LB apply : 1664.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.013631974000000002 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.5% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.023170083 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.98 us (0.5%)
patch tree reduce : 511.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 : 631.59 us (98.5%)
LB move op cnt : 0
LB apply : 1884.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.40 us (68.0%)
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 | 1.5442e+04 | 10087 | 1 | 6.532e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 5.238253e+06 N.s^-1
SPH setup: the generation step took : 0.003065978 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.57 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 : 1433.00 ns (0.2%)
patch tree reduce : 1313.00 ns (0.2%)
gen split merge : 691.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 832.00 ns (0.1%)
LB compute : 848.45 us (98.7%)
LB move op cnt : 0
LB apply : 2.94 us (0.3%)
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 (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 : 2.09 us (0.2%)
patch tree reduce : 491.00 ns (0.0%)
gen split merge : 361.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 301.00 ns (0.0%)
LB compute : 1075.91 us (99.3%)
LB move op cnt : 0
LB apply : 1633.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.38 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 : 1724.00 ns (0.2%)
patch tree reduce : 531.00 ns (0.1%)
gen split merge : 351.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.0%)
LB compute : 1023.38 us (99.2%)
LB move op cnt : 0
LB apply : 1903.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.012056285 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.020157134 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 (0.7%)
patch tree reduce : 511.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 : 408.58 us (97.7%)
LB move op cnt : 0
LB apply : 1783.00 ns (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.27 us (68.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.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 | 1.5442e+04 | 10087 | 1 | 6.532e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 3.990957e+06 N.s^-1
SPH setup: the generation step took : 0.003922676 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.64 us (67.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 : 1693.00 ns (0.2%)
patch tree reduce : 1092.00 ns (0.1%)
gen split merge : 681.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1242.00 ns (0.2%)
LB compute : 800.07 us (98.6%)
LB move op cnt : 0
LB apply : 2.85 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.18 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 : 1713.00 ns (0.2%)
patch tree reduce : 461.00 ns (0.1%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 271.00 ns (0.0%)
LB compute : 891.54 us (99.2%)
LB move op cnt : 0
LB apply : 1594.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.54 us (64.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 : 1723.00 ns (0.2%)
patch tree reduce : 410.00 ns (0.0%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.0%)
LB compute : 913.21 us (99.2%)
LB move op cnt : 0
LB apply : 1573.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.013679703000000001 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.022446935 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.08 us (0.5%)
patch tree reduce : 521.00 ns (0.1%)
gen split merge : 411.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 331.00 ns (0.1%)
LB compute : 556.10 us (98.3%)
LB move op cnt : 0
LB apply : 1793.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.24 us (66.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 | 1.4085e+04 | 10087 | 1 | 7.162e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 5.648730e+06 N.s^-1
SPH setup: the generation step took : 0.0027818580000000003 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.90 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 : 1994.00 ns (0.2%)
patch tree reduce : 1062.00 ns (0.1%)
gen split merge : 812.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 952.00 ns (0.1%)
LB compute : 810.63 us (98.6%)
LB move op cnt : 0
LB apply : 2.93 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.27 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 : 2.21 us (0.2%)
patch tree reduce : 431.00 ns (0.0%)
gen split merge : 380.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.0%)
LB compute : 956.52 us (99.2%)
LB move op cnt : 0
LB apply : 1663.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.20 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 : 1263.00 ns (0.1%)
patch tree reduce : 451.00 ns (0.1%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 240.00 ns (0.0%)
LB compute : 867.07 us (99.2%)
LB move op cnt : 0
LB apply : 1663.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.011918238000000001 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.019607128 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 (0.5%)
patch tree reduce : 531.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 320.00 ns (0.1%)
LB compute : 591.61 us (98.4%)
LB move op cnt : 0
LB apply : 1703.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (66.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
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 | 1.4316e+04 | 10087 | 1 | 7.046e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 5.071495e+06 N.s^-1
SPH setup: the generation step took : 0.010539156000000001 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.42 us (66.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 : 1503.00 ns (0.2%)
patch tree reduce : 912.00 ns (0.1%)
gen split merge : 692.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 952.00 ns (0.1%)
LB compute : 789.81 us (98.7%)
LB move op cnt : 0
LB apply : 2.85 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.23 us (9.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.02 us (0.2%)
patch tree reduce : 461.00 ns (0.1%)
gen split merge : 351.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 300.00 ns (0.0%)
LB compute : 886.54 us (99.2%)
LB move op cnt : 0
LB apply : 1603.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.17 us (64.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 : 1252.00 ns (0.1%)
patch tree reduce : 451.00 ns (0.0%)
gen split merge : 370.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 271.00 ns (0.0%)
LB compute : 902.29 us (99.2%)
LB move op cnt : 0
LB apply : 1723.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.012269784 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.027895379 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.19 us (0.4%)
patch tree reduce : 491.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 321.00 ns (0.0%)
LB compute : 753.82 us (98.8%)
LB move op cnt : 0
LB apply : 1633.00 ns (0.2%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (65.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 | 1.4587e+04 | 10087 | 1 | 6.915e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 4.358787e+06 N.s^-1
SPH setup: the generation step took : 0.0032225500000000002 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.86 us (62.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 : 1533.00 ns (0.2%)
patch tree reduce : 921.00 ns (0.1%)
gen split merge : 711.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 831.00 ns (0.1%)
LB compute : 791.72 us (98.7%)
LB move op cnt : 0
LB apply : 2.99 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.30 us (66.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 : 1273.00 ns (0.1%)
patch tree reduce : 431.00 ns (0.0%)
gen split merge : 371.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.0%)
LB compute : 869.68 us (99.2%)
LB move op cnt : 0
LB apply : 1683.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.47 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 : 1322.00 ns (0.1%)
patch tree reduce : 380.00 ns (0.0%)
gen split merge : 341.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 886.79 us (99.2%)
LB move op cnt : 0
LB apply : 1673.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.014357858000000001 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.023017698000000003 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.15 us (0.5%)
patch tree reduce : 491.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.1%)
LB compute : 580.49 us (98.4%)
LB move op cnt : 0
LB apply : 1854.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.27 us (65.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.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 | 1.4539e+04 | 10087 | 1 | 6.938e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 5.983744e+06 N.s^-1
SPH setup: the generation step took : 0.004322774000000001 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.82 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 : 1723.00 ns (0.2%)
patch tree reduce : 1062.00 ns (0.1%)
gen split merge : 711.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1102.00 ns (0.1%)
LB compute : 868.39 us (98.7%)
LB move op cnt : 0
LB apply : 3.00 us (0.3%)
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 (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 : 1333.00 ns (0.1%)
patch tree reduce : 381.00 ns (0.0%)
gen split merge : 391.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 290.00 ns (0.0%)
LB compute : 906.34 us (99.3%)
LB move op cnt : 0
LB apply : 1643.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.12 us (65.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 : 1192.00 ns (0.1%)
patch tree reduce : 391.00 ns (0.0%)
gen split merge : 351.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 261.00 ns (0.0%)
LB compute : 862.09 us (99.3%)
LB move op cnt : 0
LB apply : 1572.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.012046092000000001 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.021146398 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.01 us (0.5%)
patch tree reduce : 470.00 ns (0.1%)
gen split merge : 421.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 618.03 us (98.5%)
LB move op cnt : 0
LB apply : 1753.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.38 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 | 1.4114e+04 | 10087 | 1 | 7.147e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 4.653461e+06 N.s^-1
SPH setup: the generation step took : 0.006112012 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.86 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 : 1523.00 ns (0.2%)
patch tree reduce : 922.00 ns (0.1%)
gen split merge : 702.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1202.00 ns (0.1%)
LB compute : 963.77 us (98.9%)
LB move op cnt : 0
LB apply : 3.19 us (0.3%)
Info: Compute load ... [DataInserterUtility][rank=0]
Info: run scheduler step ... [DataInserterUtility][rank=0]
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.63 us (64.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.1%)
patch tree reduce : 480.00 ns (0.0%)
gen split merge : 401.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 311.00 ns (0.0%)
LB compute : 1033.33 us (99.3%)
LB move op cnt : 0
LB apply : 1743.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.25 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 : 1223.00 ns (0.1%)
patch tree reduce : 431.00 ns (0.0%)
gen split merge : 361.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 864.70 us (99.2%)
LB move op cnt : 0
LB apply : 1533.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.011822946 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.02374145 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.19 us (0.5%)
patch tree reduce : 501.00 ns (0.1%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 340.00 ns (0.1%)
LB compute : 610.48 us (98.5%)
LB move op cnt : 0
LB apply : 1793.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.35 us (68.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.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 | 1.4535e+04 | 10087 | 1 | 6.940e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 7.360861e+06 N.s^-1
SPH setup: the generation step took : 0.003000748 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.41 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 : 1523.00 ns (0.2%)
patch tree reduce : 1152.00 ns (0.1%)
gen split merge : 691.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 931.00 ns (0.1%)
LB compute : 788.09 us (98.7%)
LB move op cnt : 0
LB apply : 2.92 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.27 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 : 1282.00 ns (0.1%)
patch tree reduce : 401.00 ns (0.0%)
gen split merge : 370.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 281.00 ns (0.0%)
LB compute : 883.68 us (99.2%)
LB move op cnt : 0
LB apply : 1563.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 : 21.65 us (94.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 : 1703.00 ns (0.2%)
patch tree reduce : 431.00 ns (0.0%)
gen split merge : 340.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 885.23 us (99.2%)
LB move op cnt : 0
LB apply : 1644.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.011643401000000001 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.019426080000000002 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.36 us (0.6%)
patch tree reduce : 511.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 : 586.61 us (98.4%)
LB move op cnt : 0
LB apply : 1824.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.29 us (67.0%)
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 | 1.3542e+04 | 10087 | 1 | 7.449e-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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 6.354092e+06 N.s^-1
SPH setup: the generation step took : 0.005401747 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.39 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 : 1542.00 ns (0.2%)
patch tree reduce : 1212.00 ns (0.1%)
gen split merge : 792.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 952.00 ns (0.1%)
LB compute : 891.52 us (98.7%)
LB move op cnt : 0
LB apply : 3.68 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.42 us (64.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 : 1673.00 ns (0.2%)
patch tree reduce : 441.00 ns (0.0%)
gen split merge : 321.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 972.50 us (99.3%)
LB move op cnt : 0
LB apply : 1733.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.19 us (62.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.1%)
patch tree reduce : 401.00 ns (0.0%)
gen split merge : 341.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 251.00 ns (0.0%)
LB compute : 874.39 us (99.2%)
LB move op cnt : 0
LB apply : 1623.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.012759513 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.6% 0.0% | 1.16 GB | 5.04 MB |
+------+--------------------+-------+-------------+-------------+-------------+
SPH setup: the setup took : 0.023067052 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.32 us (0.6%)
patch tree reduce : 501.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 : 563.62 us (98.3%)
LB move op cnt : 0
LB apply : 1864.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.54 us (64.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
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 | 1.3768e+04 | 10087 | 1 | 7.326e-01 | 0.0% | 0.5% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 6.233901e+06 N.s^-1
SPH setup: the generation step took : 0.002366154 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 : 8.02 us (63.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 : 1222.00 ns (0.1%)
patch tree reduce : 2.06 us (0.2%)
gen split merge : 1012.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1162.00 ns (0.1%)
LB compute : 889.23 us (98.5%)
LB move op cnt : 0
LB apply : 3.97 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.37 us (60.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 : 1192.00 ns (0.1%)
patch tree reduce : 351.00 ns (0.0%)
gen split merge : 361.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 260.00 ns (0.0%)
LB compute : 814.83 us (99.2%)
LB move op cnt : 0
LB apply : 1583.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.02 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 : 1113.00 ns (0.1%)
patch tree reduce : 311.00 ns (0.0%)
gen split merge : 330.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 261.00 ns (0.0%)
LB compute : 799.43 us (99.2%)
LB move op cnt : 0
LB apply : 1493.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.010953013000000001 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.017563345 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.82 us (0.6%)
patch tree reduce : 611.00 ns (0.1%)
gen split merge : 401.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 351.00 ns (0.1%)
LB compute : 601.91 us (98.2%)
LB move op cnt : 0
LB apply : 1914.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.33 us (66.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 = (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 | 1.2919e+04 | 10087 | 1 | 7.808e-01 | 0.0% | 0.5% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 6.858993e+06 N.s^-1
SPH setup: the generation step took : 0.0020123140000000003 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 : 7.26 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 : 1393.00 ns (0.2%)
patch tree reduce : 1372.00 ns (0.2%)
gen split merge : 1042.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 942.00 ns (0.1%)
LB compute : 870.57 us (98.5%)
LB move op cnt : 0
LB apply : 4.19 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.16 us (67.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 : 1313.00 ns (0.1%)
patch tree reduce : 531.00 ns (0.1%)
gen split merge : 331.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 270.00 ns (0.0%)
LB compute : 881.35 us (99.2%)
LB move op cnt : 0
LB apply : 1734.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.05 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 : 1152.00 ns (0.1%)
patch tree reduce : 431.00 ns (0.1%)
gen split merge : 361.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.0%)
LB compute : 833.26 us (99.2%)
LB move op cnt : 0
LB apply : 1433.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.011003597 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.01802196 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.97 us (0.7%)
patch tree reduce : 1001.00 ns (0.2%)
gen split merge : 431.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 391.00 ns (0.1%)
LB compute : 576.91 us (98.0%)
LB move op cnt : 0
LB apply : 2.24 us (0.4%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.01 us (67.0%)
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 | 1.4819e+04 | 10087 | 1 | 6.807e-01 | 0.0% | 0.5% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 6.135357e+06 N.s^-1
SPH setup: the generation step took : 0.0021450510000000002 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.27 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 : 1272.00 ns (0.2%)
patch tree reduce : 1292.00 ns (0.2%)
gen split merge : 1142.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 1032.00 ns (0.1%)
LB compute : 796.51 us (98.5%)
LB move op cnt : 0
LB apply : 4.01 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.09 us (64.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 : 1132.00 ns (0.1%)
patch tree reduce : 381.00 ns (0.0%)
gen split merge : 351.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 291.00 ns (0.0%)
LB compute : 869.41 us (99.3%)
LB move op cnt : 0
LB apply : 1493.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.09 us (67.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 : 1132.00 ns (0.1%)
patch tree reduce : 371.00 ns (0.0%)
gen split merge : 361.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 250.00 ns (0.0%)
LB compute : 812.30 us (99.2%)
LB move op cnt : 0
LB apply : 1633.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.010885487000000001 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.017909538000000003 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 : 5.84 us (1.0%)
patch tree reduce : 1342.00 ns (0.2%)
gen split merge : 671.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 771.00 ns (0.1%)
LB compute : 574.03 us (97.5%)
LB move op cnt : 0
LB apply : 1984.00 ns (0.3%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 2.12 us (68.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
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 | 1.4439e+04 | 10087 | 1 | 6.986e-01 | 0.0% | 0.5% 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,
"beta_AV": 2.0,
"sigma_decay": 0.1,
"type": "varying_cd10"
},
"boundary_config": {
"bc_type": "periodic"
},
"cfl_config": {
"cfl_cour": 0.0,
"cfl_force": 0.0,
"cfl_multiplier_stiffness": 2.0,
"eta_sink": 0.05
},
"combined_dtdiv_divcurlv_compute": false,
"debug_dump_filename": "",
"do_debug_dump": false,
"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,
"scheduler_config": {
"merge_load_value": 0,
"split_load_value": 0
},
"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 = 3.255109e+06 N.s^-1
SPH setup: the generation step took : 0.0037509780000000003 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.02 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 : 1282.00 ns (0.2%)
patch tree reduce : 981.00 ns (0.1%)
gen split merge : 811.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 962.00 ns (0.1%)
LB compute : 754.23 us (98.5%)
LB move op cnt : 0
LB apply : 3.69 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.23 us (63.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 : 1612.00 ns (0.2%)
patch tree reduce : 461.00 ns (0.1%)
gen split merge : 371.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 280.00 ns (0.0%)
LB compute : 890.48 us (99.2%)
LB move op cnt : 0
LB apply : 1673.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.13 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 : 1222.00 ns (0.1%)
patch tree reduce : 380.00 ns (0.0%)
gen split merge : 381.00 ns (0.0%)
split / merge op : 0/0
apply split merge : 240.00 ns (0.0%)
LB compute : 856.87 us (99.2%)
LB move op cnt : 0
LB apply : 1643.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.010943726 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.018915588 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.55 us (1.2%)
patch tree reduce : 501.00 ns (0.2%)
gen split merge : 391.00 ns (0.1%)
split / merge op : 0/0
apply split merge : 341.00 ns (0.1%)
LB compute : 288.67 us (96.8%)
LB move op cnt : 0
LB apply : 1754.00 ns (0.6%)
Info: Scheduler step timings : [Scheduler][rank=0]
metadata sync : 1964.00 ns (66.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
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 | 1.3894e+04 | 10087 | 1 | 7.260e-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]
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.011725871267882111 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.0076836941261466205 for cube_hcp
max relative error mm4: 0.05833850788436179 for cube_hcp
avg relative error mm4: 0.01491142025112289 for cube_hcp
min relative error mm4: 0.0009000262033407927 for cube_hcp
std relative error mm4: 0.007058594416839623 for cube_hcp
max relative error mm5: 0.03130228236510673 for cube_hcp
avg relative error mm5: 0.005542488773510676 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.02584963102487481 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.011068208855935385 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 37.524 seconds)
Estimated memory usage: 155 MB