Define with_logabsdet_jacobian for unary plus and minus

Author: oschulz

src/with_ladj.jl

@@ -118,6 +118,8 @@ with_logabsdet_jacobian(::typeof(inv), A::AbstractMatrix) = (inv(A), -2 * _ndof(
 with_logabsdet_jacobian(::typeof(adjoint), x) = (adjoint(x), zero(eltype(x)))
 with_logabsdet_jacobian(::typeof(transpose), x) = (transpose(x), zero(eltype(x)))
 
+with_logabsdet_jacobian(::typeof(+), x) = (+(x), zero(eltype(x)))
+with_logabsdet_jacobian(::typeof(-), x) = (-(x), zero(eltype(x)))
 
 with_logabsdet_jacobian(::typeof(exp), x) = (exp(x), x)
 with_logabsdet_jacobian(::typeof(exp2), x) = (exp2(x), log(2) * x + log(log(oftype(x, 2))))

test/test_with_ladj.jl

@@ -57,7 +57,7 @@ include("getjacobian.jl")
     end
 
     @testset "with_logabsdet_jacobian on log and exp functions" begin
-        for f in (exp, log, exp2, log2, exp10, log10, expm1, log1p)
+        for f in (+, -, exp, log, exp2, log2, exp10, log10, expm1, log1p)
             test_with_logabsdet_jacobian(f, x, getjacobian)
         end
     end