Remove TestImage dependency. (#38)

remove TestImage dependency and stabilize random_point on SPD.

Author: kellertuer

Project.toml

@@ -1,7 +1,7 @@
 name = "Manopt"
 uuid = "0fc0a36d-df90-57f3-8f93-d78a9fc72bb5"
 authors = ["Ronny Bergmann <manopt@ronnybergmann.net>"]
-version = "0.2.7"
+version = "0.2.8"
 
 [deps]
 ColorSchemes = "35d6a980-a343-548e-a6ea-1d62b119f2f4"
@@ -17,7 +17,6 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
-TestImages = "5e47fb64-e119-507b-a336-dd2b206d9990"
 
 [compat]
 ColorSchemes = "3.5.0"
@@ -26,7 +25,6 @@ Colors = "0.11.2, 0.12"
 Manifolds = "0.4.1"
 ManifoldsBase = "0.9"
 StaticArrays = "0.12"
-TestImages = "1.0.0"
 julia = "1.0"
 
 [extras]

src/data/artificialDataFunctions.jl

@@ -324,7 +324,7 @@ function artificial_SPD_image(pts::Int=64, stepsize=1.5)
                 1 + v3[row + col] - stepsize * (col > pts / 2 ? 1 : 0)
                 4 - v3[row + col] + stepsize * (row > pts / 2 ? 1 : 0)
             ]
-            data[row, col] = A * B * C * Diagonal(scale) * C' * B' * A'
+            data[row, col] = Matrix(Symmetric(A * B * C * Diagonal(scale) * C' * B' * A'))
         end
     end
     return data
@@ -342,7 +342,7 @@ function artificial_SPD_image2(pts=64, fraction=0.66)
     β = π / 3
     B = [1.0 0.0 0.0; 0.0 cos(β) -sin(β); 0.0 sin(β) cos(β)]
     A = [cos(α) -sin(α) 0.0; sin(α) cos(α) 0.0; 0.0 0.0 1.0]
-    Zo = A * B * Diagonal([2.0, 4.0, 8.0]) * B' * A'
+    Zo = Matrix(Symmetric(A * B * Diagonal([2.0, 4.0, 8.0]) * B' * A'))
     # create a second matrix
     α = -4.0 * π / 3
     β = -π / 3
@@ -367,7 +367,9 @@ function artificial_SPD_image2(pts=64, fraction=0.66)
                 C = exp(
                     M,
                     C,
-                    vector_transport_to(M, Zo, log(M, Zo, Zl), C, ParallelTransport()),
+                    vector_transport_to(
+                        M, Symmetric(Zo), log(M, Zo, Zl), Symmetric(C), ParallelTransport()
+                    ),
                     (col - 1.0) / (pts - 1),
                 )
             end

src/random.jl

@@ -118,7 +118,7 @@ function random_point(
 ) where {N}
     D = Diagonal(1 .+ rand(N)) # random diagonal matrix
     s = qr(σ * randn(N, N)) # random q
-    return s.Q * D * transpose(s.Q)
+    return Matrix(Symmetric(s.Q * D * transpose(s.Q)))
 end
 
 @doc raw"""