I noticed the default definition for the length of an abstract array is prod(size(A))
. Shouldn’t it be something more like prod(length.(axes(x)))
to account for non-zero indexing?
julia> y = OffsetArray{Float64}(undef, -1:1, -7:7, -128:512, -5:5, -1:1, -3:3, -2:2, -1:1);
julia> prod(length.(axes(y)))
99947925
julia> prod(size(y))
99947925