Memory Required for Large Matrices

This table shows how much memory is required for four matrices of the given size and type. (We need four matrices: A, B, C1, and C2.)

Generating These Tables

mem_req(s, ::Type{T}) where {T} = 4s^2*sizeof(T) / (1 << 30)

function print_table(types::Vector{DataType}, Ns = nothing)
    println("| Matrix Size | $(join(types, " Memory | ")) Memory |")
    println("| ----------- | $(repeat(" ------ |", length(types)))")
    if Ns isa Nothing
        _Ns = sort(unique(vcat(collect(1:1:20), collect(20:10:100))))
    else
        _Ns = Ns
    end
    for N in _Ns
        mem = mem_req.(N * 1_000, types)
        m = round.(mem; digits = 2)
        println("| $(N)k by $(N)k | $(join(m, " GiB | ")) GiB |")
    end
    return nothing
end

julia> print_table([Float64, Float32, Int64, Int32])