test: update the test to check distances

Author: tpoisot

src/lib/overallmean.jl

@@ -18,10 +18,10 @@ function Statistics.mean(w::Vector{T}) where {T <: Womble}
         di = filter(!isnan, [deg2rad(womble.θ[_idx]) for womble in w])
         if !isempty(ch)
             m[_idx] = mean(ch)
-            α[_idx] = atan(mean(sin.(di)), mean(cos.(di)))
+            α[_idx] = atan(sum(sin.(di)), sum(cos.(di)))
         end
     end
-    average_direction = rad2deg.(α) .+ 180.0
+    average_direction = rad2deg.(α .- π)
     if w[1] isa LatticeWomble
         return LatticeWomble(m, average_direction, w[1].x, w[1].y)
     end

test/runtests.jl

@@ -1,12 +1,14 @@
 global anyerrors = false
 
 tests = [
-    "mock test" => "00_allgood.jl",
+    "Mock test" => "00_allgood.jl",
     "Linear gradient" => "01_gradient.jl",
     "Boundaries detection" => "02_boundaries.jl",
     "Triangulation wombling" => "03_delaunay.jl",
     "Direction of change" => "04_direction.jl",
     "Overall mean wombling" => "05_mean.jl",
+    "Correct values are returned" => "06_values.jl",
+    "Overall mean correct values are returned" => "07_overallmean.jl"
 ]
 
 for test in tests

test/units/06_values.jl

@@ -0,0 +1,21 @@
+module SBTestGradient
+
+using SpatialBoundaries
+using NeutralLandscapes
+using Test
+
+angles = rand(100) .* 360
+
+for angle in angles
+    landscape = rand(PlanarGradient(angle), 20, 20)
+    x = collect(LinRange(0.2, 1.8, size(landscape, 1)))
+    y = collect(LinRange(0.2, 1.8, size(landscape, 2)))
+
+    W = wombling(x, y, landscape)
+
+    for θ in W.θ
+        @test deg2rad(abs(θ  - angle)) ≈ π atol = 0.01
+    end
+end
+
+end

test/units/07_overallmean.jl

@@ -0,0 +1,32 @@
+module SBTestGradient
+
+using SpatialBoundaries
+using NeutralLandscapes
+using Test
+
+angles = rand(100) .* 360
+
+for replicate in Base.OneTo(10)
+    θ1 = rand(angles)
+    θ2 = rand(angles)
+    L1 = rand(PlanarGradient(θ1), 20, 20)
+    L2 = rand(PlanarGradient(θ2), size(L1)...)
+    x = collect(LinRange(0.2, 1.8, size(L1, 1)))
+    y = collect(LinRange(0.2, 1.8, size(L1, 2)))
+
+    W = mean([wombling(x, y, L) for L in [L1, L2]])
+
+    α1 = deg2rad(θ1)
+    α2 = deg2rad(θ2)
+    target = atan(sin(α1)+sin(α2), cos(α1)+cos(α2))
+    observed = deg2rad(rand(W.θ))
+
+    diff = (target-observed) - 2π
+    if abs(diff) > 0.01
+        @test diff ≈ -2π atol = 0.001
+    else
+        @test diff ≈ 0.0 atol = 0.001
+    end
+end
+
+end