Crash with BigFloats #10
Description
I have a crash which I've narrowed down to the following... no idea whether it's this package's fault though:
using LinearAlgebra, KissMCMC
f(θs, ts) = [-t*exp(θ-exp(θ)*t)/length(θs) for t in ts, θ in θs]
g(v, ts) = begin
j = f(v, ts)
x = j' * j
# s, ld = LinearAlgebra.logabsdet(x)
# ifelse(s>0, Float64(ld), -Inf)
d = det(x)
d>0 ? Float64(log(d)) : -Inf
end
KissMCMC.emcee(v -> g(v, 1:3), [rand(3) for _ in 1:10]; niter=10^4) # fine!
versioninfo()
KissMCMC.emcee(v -> g(v, big.(1:3.0)), [big.(rand(3)) for _ in 1:10]; niter=10^4)julia> versioninfo()
Julia Version 1.5.0
Commit 96786e22cc (2020-08-01 23:44 UTC)
Platform Info:
OS: macOS (x86_64-apple-darwin18.7.0)
CPU: Intel(R) Core(TM) i5-7360U CPU @ 2.30GHz
WORD_SIZE: 64
LIBM: libopenlibm
LLVM: libLLVM-9.0.1 (ORCJIT, skylake)
Environment:
JULIA_NUM_THREADS = 4
JULIA_PKG_SERVER = https://geo.pkg.julialang.org
julia> KissMCMC.emcee(v -> g(v, big.(1:3.0)), [big.(rand(3)) for _ in 1:10]; niter=10^4)
julia(303,0x70000e9e8000) malloc: *** error for object 0x7fb6d3ea8d70: pointer being freed was not allocated
julia(303,0x108ff35c0) malloc: *** error for object 0x7fb6d306b460: pointer being freed was not allocated
julia(303,0x70000e9e8000) malloc: *** set a breakpoint in malloc_error_break to debug
julia(303,0x108ff35c0) malloc: *** set a breakpoint in malloc_error_break to debug
signal (6): Abort trap: 6
in expression starting at REPL[5]:1
in expression starting at Rin expression starting at REPL[5]:1
__pthread_kill at /usr/lib/system/libsystem_kernel.dylib (unknown line)
Allocations: 34375251 (Pool: 34366235; Big: 9016); GC: 38
Abort trap: 6
(And the same on Julia 1.4.2)