BoundsError in CartesianIndices on 32-bit Windows


I experienced Integer overflow on Windows due to the fact that Int32 is the default one there:

length(CartesianIndices((1000000, 1000000)))

Succeeds on linux64 but fails on Windows (something like A[-72232] is out of bounds).

I copied randsubseq! implementation, where the overflow happens, and adjusted it so index vars are Int64 now:

randsubseq(r::AbstractRNG, A::AbstractArray{T}, p::Real) where {T} =
    randsubseq!(r, T[], A, p)

# Fill S (resized as needed) with a random subsequence of A, where
# each element of A is included in S with independent probability p.
function randsubseq!(r::AbstractRNG, S::AbstractArray, A::AbstractArray, p::Real)
    Base.require_one_based_indexing(S, A)
    0 <= p <= 1 || throw(ArgumentError("probability $p not in [0,1]"))
    n = prod(size(A), init=Int64(1))  # length - there the overflow happened originally
    p == 1 && return copyto!(resize!(S, n), A)
    p == 0 && return S
    nexpected = p * n
    sizehint!(S, round(Int64, nexpected + 5*sqrt(nexpected)))
    if p > 0.15
        for i = 1:n
            rand(r) <= p && push!(S, A[i])
        L = -1 / log1p(-p)
        i = Int64(0)
        while true
            s = randexp(r) * L
            s >= n - i && return S
            push!(S, A[i += ceil(Int64, s)])
    return S

Surprisingly I get:

  BoundsError: attempt to access 1000000×1000000 CartesianIndices{2, Tuple{Base.OneTo{Int32}, Base.OneTo{Int32}}} at index [10504]

on Windows, which doesn’t really make sense. How come I can’t access something that is in bounds?

Welcome! So, the reason that CartesianIndices exists is to be used as an index — and on a 32-bit system, you cannot create a thing that has more than the ~2 billion addressable elements. A 1000000x1000000 array has a trillion elements.

What are you trying to do here? I suspect there’s a much better solution than a trillion iterations in the first place.

So I work with a sparse tensor library Finch.jl.

There’s a method there fsprand, fiber-sparse-rand that returns a COO tensor, given dimensions and density p:


fsprand(1_000_000, 1_000_000, 1e-06)

uses CartesianIndices to describe only non-zeros, and calls:

randsubseq(r, CartesianIndices(M), p)

where M=(1_000_000, 1_000_000) and p=1e-06, which fails on Windows.

We have an open PR with some discussion: Wma/fix 488 by willow-ahrens · Pull Request #489 · willow-ahrens/Finch.jl · GitHub

1 Like

Ah, of course, the non-sensical bounds error is because, when you use linear indexing, it’ll check the index against the length. And the length overflows to negative, so everything is going to be OOB.

The crux of it is that arrays (even sparse ones) with more than 2^31 elements are going to be troublesome on 32-bit systems.

1 Like