Merge pull request #527 from JuliaReach/schillic/expm

#121 - Fix 'expm' deprecation
JuliaReach · Aug 12, 2018 · c4768cd · c4768cd
2 parents 3bb074d + e079b0a
commit c4768cd
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/src/compat.jl b/src/compat.jl
@@ -15,11 +15,13 @@ if VERSION < v"0.7-"
     @inline function _At_mul_B(A, B)
         return At_mul_B(A, B)
     end
+    expmat = expm
 else
     using SparseArrays, LinearAlgebra
     import LinearAlgebra: norm, checksquare, LAPACKException, SingularException,
                           eye, ×
     @inline function _At_mul_B(A, B)
         return transpose(A) * B
     end
+    expmat = exp
 end
diff --git a/test/unit_ExponentialMap.jl b/test/unit_ExponentialMap.jl
@@ -68,7 +68,7 @@ for N in [Float64, Float32]
     σ(sparsevec(d), emap)
 
     # check consistency with respect to explicit computation of the matrix exponential
-    svec_explicit = σ(d, expm(Matrix(m)) * b)
+    svec_explicit = σ(d, expmat(Matrix(m)) * b)
     if N == Float64
         # precision with Float32 is not sufficient
         @test svec ≈ svec_explicit
@@ -114,7 +114,7 @@ for N in [Float64, Float32]
     σ(sparsevec(d), projmap)
 
     # check consistency with respect to explicit computation of the matrix exponential
-    P = L * expm(Matrix(m)) * R
+    P = L * expmat(Matrix(m)) * R
     svec_explicit = σ(d, P*b)
     @test svec ≈ svec_explicit
 end