Merge pull request #35 from bcbi/dpa/pseudoinverse

Add the private/internal `_pseudoinverse` function

Author: DilumAluthge

Project.toml

@@ -1,7 +1,7 @@
 name = "MaximumLikelihoodProblems"
 uuid = "e149cee4-a9c8-4a03-a7b6-d91720412337"
 authors = ["Dilum Aluthge"]
-version = "0.1.5"
+version = "0.1.6"
 
 [deps]
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"

src/MaximumLikelihoodProblems.jl

@@ -5,6 +5,7 @@ include("public.jl")
 include("types.jl")
 
 include("fit.jl")
+include("linear_algebra.jl")
 include("loglikelihood.jl")
 
 include("Internal.jl")

src/linear_algebra.jl

@@ -0,0 +1,10 @@
+import LinearAlgebra
+
+function _pseudoinverse(m::AbstractMatrix)
+    try
+        return LinearAlgebra.inv(m)
+    catch ex
+        @debug("", exception = (ex, catch_backtrace()))
+    end
+    return LinearAlgebra.pinv(m)
+end

test/runtests.jl

@@ -1,5 +1,6 @@
 import MaximumLikelihoodProblems
 
+import LinearAlgebra
 import Random
 import Statistics
 import Test
@@ -11,5 +12,14 @@ examples_directory = joinpath(root_directory, "examples")
 test_directory = joinpath(root_directory, "test")
 
 Test.@testset "MaximumLikelihoodProblems.jl" begin
-    include("examples.jl")
+    Test.@testset "linear_algebra.jl" begin
+        Test.@testset "_pseudoinverse" begin
+            Test.@test MaximumLikelihoodProblems._pseudoinverse([1 2; 3 4]) == LinearAlgebra.inv([1 2; 3 4])
+            Test.@test MaximumLikelihoodProblems._pseudoinverse([0 0; 0 0]) == [0 0; 0 0]
+        end
+    end
+    
+    Test.@testset "examples" begin
+        include("examples.jl")
+    end
 end