Add SparseRips constructor with point cloud inputs (#79)

* Add SparseRips constructor with points

* bump version

* add tests

Project.toml

@@ -1,7 +1,7 @@
 name = "Ripserer"
 uuid = "aa79e827-bd0b-42a8-9f10-2b302677a641"
 authors = ["mtsch <matijacufar@gmail.com>"]
-version = "0.14.2"
+version = "0.14.3"
 
 [deps]
 Compat = "34da2185-b29b-5c13-b0c7-acf172513d20"

src/ripsfiltration.jl

@@ -229,6 +229,9 @@ end
 function SparseRips(dist; threshold=nothing)
     return SparseRips{Int}(dist; threshold=threshold)
 end
+function SparseRips{I}(points::AbstractVector; metric=Euclidean(), kwargs...) where I
+    return SparseRips{I}(distances(metric, points); kwargs...)
+end
 
 @propagate_inbounds function dist(rips::SparseRips{<:Any, T}, i, j) where T
     res = rips.dist[i, j]

test/ripsfiltration.jl

@@ -75,13 +75,24 @@ for Filtration in (Rips, SparseRips)
 end
 
 @testset "Rips points constructor" begin
-    filtration = Rips([(sin(x), cos(x)) for x in range(0, 2π, length=100)])
-    @test all(x -> x > 0, dist(filtration, i, j) for i in 1:10 for j in i+1:9)
+    filtration = Rips(
+        [(sin(x), cos(x)) for x in range(0, 2π, length=101)[1:end-1]]
+    )
+    @test all(x -> x > 0, dist(filtration, i, j) for i in 1:100 for j in i+1:100)
     @test eltype(edges(filtration)) === Simplex{1, Float64, Int}
 
-    filtration = Rips{Int32}([(sin(x), cos(x)) for x in range(0, 2π, length=100)])
-    @test all(x -> x > 0, dist(filtration, i, j) for i in 1:10 for j in i+1:9)
-    @test eltype(edges(filtration)) === Simplex{1, Float64, Int32}
+    filtration = Rips{Int32}(
+        [(sin(x), cos(x)) for x in range(0f0, 2f0π, length=101)[1:end-1]]
+    )
+    @test all(x -> x > 0, dist(filtration, i, j) for i in 1:100 for j in i+1:100)
+    @test eltype(edges(filtration)) === Simplex{1, Float32, Int32}
+end
+
+@testset "Rips points constructor" begin
+    filtration = SparseRips(
+        [(sin(x), cos(x)) for x in range(0, 2π, length=101)[1:end-1]], threshold=0.1
+    )
+    @test maximum(dist(filtration)) ≤ 0.1
 end
 
 @testset "Errors" begin