search

bosscha / gaia-shock

Code to extract stellar clusters and tidal tails using Gaia data. The main stellar cluster is extracted by optimizing the DBSCAN parameters and the isolated stars ("tails") by minimizing the CMD distance.
MIT License
2 stars 1 forks source link

voronoi tesselation called from scipy seems to crash #19

Open bosscha opened 2 years ago

bosscha commented 2 years ago

It seems that with a lot of points (epsilon large) the voronoi tesselation called by PyCall crashes.

munmap_chunk(): invalid pointer

signal (6): Aborted in expression starting at /home/stephane/Science/GAIA/dev/scripts/extra.jl:116 gsignal at /lib/x86_64-linux-gnu/libc.so.6 (unknown line) abort at /lib/x86_64-linux-gnu/libc.so.6 (unknown line) unknown function (ip: 0x7f292a388735) unknown function (ip: 0x7f292a39108b) unknown function (ip: 0x7f292a39145b) cfree at /lib/x86_64-linux-gnu/libc.so.6 (unknown line) jl_gc_free_array at /buildworker/worker/package_linux64/build/src/gc.c:1111 [inlined] sweep_malloced_arrays at /buildworker/worker/package_linux64/build/src/gc.c:1132 [inlined] gc_sweep_other at /buildworker/worker/package_linux64/build/src/gc.c:1467 [inlined] _jl_gc_collect at /buildworker/worker/package_linux64/build/src/gc.c:3144 jl_gc_collect at /buildworker/worker/package_linux64/build/src/gc.c:3240 maybe_collect at /buildworker/worker/package_linux64/build/src/gc.c:879 [inlined] jl_gc_pool_alloc at /buildworker/worker/package_linux64/build/src/gc.c:1203 jl_gcalloc at /buildworker/worker/package_linux64/build/src/julia_internal.h:315 [inlined] jl_new_uninitialized_datatype at /buildworker/worker/package_linux64/build/src/datatype.c:91 inst_datatype_inner at /buildworker/worker/package_linux64/build/src/jltypes.c:1367 inst_typew at /buildworker/worker/package_linux64/build/src/jltypes.c:1693 inst_datatype_inner at /buildworker/worker/package_linux64/build/src/jltypes.c:1435 inst_typew at /buildworker/worker/package_linux64/build/src/jltypes.c:1693 jl_instantiate_unionall at /buildworker/worker/package_linux64/build/src/jltypes.c:898 rename_unionall at /buildworker/worker/package_linux64/build/src/subtype.c:469 [inlined] unalias_unionall at /buildworker/worker/package_linux64/build/src/subtype.c:808 subtype_unionall at /buildworker/worker/package_linux64/build/src/subtype.c:836 [inlined] subtype at /buildworker/worker/package_linux64/build/src/subtype.c:1285 with_tvar at /buildworker/worker/package_linux64/build/src/subtype.c:745 subtype_unionall at /buildworker/worker/package_linux64/build/src/subtype.c:839 [inlined] subtype at /buildworker/worker/package_linux64/build/src/subtype.c:1282 with_tvar at /buildworker/worker/package_linux64/build/src/subtype.c:745 subtype_unionall at /buildworker/worker/package_linux64/build/src/subtype.c:839 [inlined] subtype at /buildworker/worker/package_linux64/build/src/subtype.c:1282 exists_subtype at /buildworker/worker/package_linux64/build/src/subtype.c:1426 [inlined] forall_exists_subtype at /buildworker/worker/package_linux64/build/src/subtype.c:1454 jl_types_equal at /buildworker/worker/package_linux64/build/src/subtype.c:1980 sig_match_simple at /buildworker/worker/package_linux64/build/src/typemap.c:216 [inlined] jl_typemap_entry_assoc_exact at /buildworker/worker/package_linux64/build/src/typemap.c:967 jl_typemap_assoc_exact at /buildworker/worker/package_linux64/build/src/julia_internal.h:1218 [inlined] jl_typemap_level_assoc_exact at /buildworker/worker/package_linux64/build/src/typemap.c:1009 jl_typemap_assoc_exact at /buildworker/worker/package_linux64/build/src/julia_internal.h:1222 [inlined] jl_typemap_level_assoc_exact at /buildworker/worker/package_linux64/build/src/typemap.c:1022 jl_typemap_assoc_exact at /buildworker/worker/package_linux64/build/src/julia_internal.h:1222 [inlined] jl_lookupgeneric at /buildworker/worker/package_linux64/build/src/gf.c:2359 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2415 convert at /home/stephane/.julia/packages/PyCall/BD546/src/conversions.jl:835 _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 py2array at /home/stephane/.julia/packages/PyCall/BD546/src/conversions.jl:350 unknown function (ip: 0x7f28c7ff4af6) _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 py2array at /home/stephane/.julia/packages/PyCall/BD546/src/conversions.jl:407 convert at /home/stephane/.julia/packages/PyCall/BD546/src/conversions.jl:421 _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 convert at /home/stephane/.julia/packages/PyCall/BD546/src/conversions.jl:835 getproperty at /home/stephane/.julia/packages/PyCall/BD546/src/PyCall.jl:312 _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 voronoi at /home/stephane/Science/GAIA/run/src/geometry.jl:80 voronoi at /home/stephane/Science/GAIA/run/src/geometry.jl:59 [inlined] metricHRD at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:24 metric2 at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:480 metric2 at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:465 [inlined] find_clusters2 at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:217 find_clusters2 at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:204 [inlined] abc_mcmc_dbscan_full2 at /home/stephane/Science/GAIA/run/src/mcmc.jl:269 cycle_extraction_optim at /home/stephane/Science/GAIA/run/src/stellarcluster.jl:848 main at /home/stephane/Science/GAIA/dev/scripts/extra.jl:101 unknown function (ip: 0x7f28c7fbeef6) _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 jl_apply at /buildworker/worker/package_linux64/build/src/julia.h:1703 [inlined] do_call at /buildworker/worker/package_linux64/build/src/interpreter.c:115 eval_value at /buildworker/worker/package_linux64/build/src/interpreter.c:204 eval_stmt_value at /buildworker/worker/package_linux64/build/src/interpreter.c:155 [inlined] eval_body at /buildworker/worker/package_linux64/build/src/interpreter.c:562 jl_interpret_toplevel_thunk at /buildworker/worker/package_linux64/build/src/interpreter.c:670 jl_toplevel_eval_flex at /buildworker/worker/package_linux64/build/src/toplevel.c:877 jl_toplevel_eval_flex at /buildworker/worker/package_linux64/build/src/toplevel.c:825 jl_toplevel_eval_in at /buildworker/worker/package_linux64/build/src/toplevel.c:929 eval at ./boot.jl:360 [inlined] include_string at ./loading.jl:1116 _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 _include at ./loading.jl:1170 include at ./Base.jl:386 _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 exec_options at ./client.jl:285 _start at ./client.jl:485 jfptrstart_43689.clone_1 at /home/stephane/soft/julia-1.6.3/lib/julia/sys.so (unknown line) _jl_invoke at /buildworker/worker/package_linux64/build/src/gf.c:2237 [inlined] jl_apply_generic at /buildworker/worker/package_linux64/build/src/gf.c:2419 jl_apply at /buildworker/worker/package_linux64/build/src/julia.h:1703 [inlined] true_main at /buildworker/worker/package_linux64/build/src/jlapi.c:560 repl_entrypoint at /buildworker/worker/package_linux64/build/src/jlapi.c:702 main at julia (unknown line) libc_start_main at /lib/x86_64-linux-gnu/libc.so.6 (unknown line) unknown function (ip: 0x4007d8) Allocations: 797508797 (Pool: 783951798; Big: 13556999); GC: 3282 Aborted (core dumped)

bosscha commented 2 years ago

A pure julia solution should be implemented