@testset "Radiality" begin
    gint = loadgraph(joinpath(testdir, "testdata", "graph-50-500.jgz"), "graph-50-500")

    c = vec(readcsv(joinpath(testdir, "testdata", "graph-50-500-rc.txt")))
    for g in testdigraphs(gint)
        z  = @inferred(radiality_centrality(g))
        zp = @inferred(parallel_radiality_centrality(g))
        @test z == zp == c
    end

    g1 = CycleGraph(4)
    add_vertex!(g1)
    add_edge!(g1, 4, 5)
    for g in testgraphs(g1)
        z  = @inferred(radiality_centrality(g))
        zp = @inferred(parallel_radiality_centrality(g))
        @test z ≈ zp ≈ [5//6, 3//4, 5//6, 11//12, 2//3]
    end
end
