added division by scalar

Manifest.toml

@@ -1,6 +1,6 @@
 # This file is machine-generated - editing it directly is not advised
 
-julia_version = "1.7.0"
+julia_version = "1.7.1"
 manifest_format = "2.0"
 
 [[deps.Artifacts]]

Project.toml

@@ -1,7 +1,7 @@
 name = "WoodburyFactorizations"
 uuid = "9f1bac23-581c-4ebc-bd36-df60d764636d"
 authors = ["Sebastian Ament <sebastianeament@gmail.com>"]
-version = "1.0.1"
+version = "1.0.2"
 
 [deps]
 LazyInverses = "9f18896c-49a9-43cc-8eef-c455a8a119c6"

src/WoodburyFactorizations.jl

@@ -192,8 +192,12 @@ end
 
 ######################## Linear Algebra primitives #############################
 # IDEA: have efficient methods of hermitian WoodburyFactorization
-Base.:*(a::Number, W::Woodbury) = Woodbury(a*W.A, W.U, a*W.C, W.V, W.α) # IDEA: could take over logabsdet efficiently
-Base.:*(W::Woodbury, a::Number) = a*W
+Base.:*(a::Number, W::Woodbury) = Woodbury(a * W.A, W.U, a * W.C, W.V, W.α) # IDEA: could take over logabsdet efficiently
+Base.:*(W::Woodbury, a::Number) = a * W
+
+Base.:\(a::Number, W::Woodbury) = W / a # IDEA: could take over logabsdet efficiently
+Base.:/(W::Woodbury, a::Number) = Woodbury(W.A / a, W.U, W.C / a, W.V, W.α)
+
 
 Base.:*(W::Woodbury, x::AbstractVecOrMat) = mul!(zero(x), W, x)
 Base.:*(B::AbstractMatrix, W::Woodbury) = adjoint(W'*B')

test/WoodburyFactorizations.jl

@@ -48,7 +48,10 @@ Base.AbstractMatrix(F::BunchKaufman) = Symmetric(Matrix(F))
     @test Y ≈ B
 
     @test α*W isa Woodbury
-    @test Matrix(α*W) ≈ α*Matrix(W) # with scalar
+    @test Matrix(α * W) ≈ α * Matrix(W) # with scalar
+    @test Matrix(W * α) ≈ Matrix(W) * α # with scalar
+    @test Matrix(W / α) ≈ Matrix(W) / α # with scalar
+    @test Matrix(α \ W) ≈ α \ Matrix(W) # with scalar
     @test logdet(abs(α)*W) ≈ logdet(abs(α)*Matrix(W)) # checking that logdet works correctly
 
     @test eltype(W) == Float64