Another edge case with reshaped distributions:
julia> insupport(reshape(Normal(), 1), [0.5])
0-dimensional BitArray{0}:
1
This comes from here, as size(d.dist) == () and hence reshape(x, size(d.dist)) is a 0-dim array:
|
function insupport(d::ReshapedDistribution{N}, x::AbstractArray{<:Real,N}) where {N} |
|
return size(d) == size(x) && insupport(d.dist, reshape(x, size(d.dist))) |
|
end |