Note
Go to the end to download the full example code.
Test the precision of SG methods in SPH#
Test that all methods give expected deviation 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 = 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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 3.0e+06 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.007780961 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.007480337 s
SPH setup: the setup took : 0.021973607000000003 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 9.7e+06 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0018191580000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005663593000000001 s
SPH setup: the setup took : 0.010357169000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 2.0e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002663902 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.0048852750000000006 s
SPH setup: the setup took : 0.010786834 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.8e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0017547480000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005477895 s
SPH setup: the setup took : 0.010108503000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.2e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0021742940000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.006077639 s
SPH setup: the setup took : 0.011163901 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.9e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0023130250000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.0056951530000000005 s
SPH setup: the setup took : 0.010862818000000002 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.8e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002757237 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.0055024620000000005 s
SPH setup: the setup took : 0.011026674 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.9e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002425135 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.0058442830000000005 s
SPH setup: the setup took : 0.011021244000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.7e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002484407 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005967063000000001 s
SPH setup: the setup took : 0.011379666 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.9e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002438089 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005780042 s
SPH setup: the setup took : 0.011143534 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.7e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.003344137 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.006314604000000001 s
SPH setup: the setup took : 0.013342424 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.8e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.001309124 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005871223 s
SPH setup: the setup took : 0.010098565 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.7e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.002482873 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.00572554 s
SPH setup: the setup took : 0.011163642000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.8e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0013732030000000001 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005597941 s
SPH setup: the setup took : 0.00976489 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 2.0e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.001198636 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005836107 s
SPH setup: the setup took : 0.009933896000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 2.1e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0011790700000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005597801 s
SPH setup: the setup took : 0.009734162000000001 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 2.0e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.001185913 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.005705432000000001 s
SPH setup: the setup took : 0.0098878 s
---------------- t = 0, dt = 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
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,
"dust_config": {
"ballabio_ts_limiter": false,
"drag_mode": {
"type": "none"
},
"evol_mode": {
"type": "none"
},
"mode": {
"type": "none"
}
},
"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"
},
"neigh_cache_strategy": "two_stage",
"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_cfl_detail": false,
"show_ghost_zone_graph": false,
"show_neigh_stats": false,
"smoothing_length_config": {
"type": "density_based"
},
"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
}
}
]
------------------------------------
SPH setup: generating particles ...
SPH setup: Nstep = 10087 ( 1.0e+04 ) Ntotal = 10087 ( 1.0e+04 rank min = 1.9e+07 max = 1.0e+04) rate = 1.008700e+04 N.s^-1
SPH setup: the generation step took : 0.0012778350000000002 s
SPH setup: final particle count = 10087 beginning injection ...
SPH setup: injected 10087 / 10087 => 100.0% | ranks with patchs = 1 / 1 <- global loop -> (msg count : 0)
SPH setup: the injection step took : 0.006007158 s
SPH setup: the setup took : 0.010273443 s
---------------- t = 0, dt = 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
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 34.078 seconds)
Estimated memory usage: 160 MB