wildart / computationalhomology.jl Goto Github PK

using ComputationalHomology

dim_0_simplices = [Simplex(i) for i in 0:5]
dim_1_simplices = [Simplex(0, 1), Simplex(0, 2), Simplex(0, 3), Simplex(0,4),
Simplex(0,5),Simplex(1, 2), Simplex(1,3), 
    Simplex(1, 4), Simplex(1, 5), Simplex(2, 3), 
    Simplex(2, 4), Simplex(2, 5), Simplex(3,4), 
    Simplex(3, 5), Simplex(4, 5)]
dim_2_simplices = [Simplex(0,1,3), Simplex(0, 1, 4), Simplex(0, 2, 3), 
    Simplex(0, 2, 5), Simplex(0, 3, 4), Simplex(0, 3, 5), Simplex(1, 2, 4), 
    Simplex(1, 2, 5), Simplex(1, 3, 4), Simplex(1, 4, 5), Simplex(2, 3, 5), Simplex(2, 4, 5)]
cplx = SimplicialComplex((vcat(vcat(dim_0_simplices, dim_1_simplices), dim_2_simplices)...))

h = homology(Int, cplx)

betti(h) 

>> (1,1)

using Polymake

dim_0_simplices = [Set([i]) for i in 0:5]
dim_1_simplices = [Set((0, 1)), Set((0, 2)), Set((0, 3)), Set((0,4)),
Set((0,5)),Set((1, 2)), Set((1,3)), 
    Set((1, 4)), Set((1, 5)), Set((2, 3)), 
    Set((2, 4)), Set((2, 5)), Set((3,4)), 
    Set((3, 5)), Set((4, 5))]
dim_2_simplices = [Set((0,1,3)), Set((0, 1, 4)), Set((0, 2, 3)), 
    Set((0, 2, 5)), Set((0, 3, 4)), Set((0, 3, 5)), Set((1, 2, 4)),
    Set((1, 2, 5)), Set((1, 3, 4)), Set((1, 4, 5)), Set((2, 3, 5)), Set((2, 4, 5))]

cplx = topaz.SimplicialComplex(INPUT_FACES=vcat(vcat(dim_0_simplices, dim_1_simplices), dim_2_simplices))

f = cplx.HOMOLOGY

>>
 ({} 0)
({} 1)
({} 3)