ReverseDiff fixes (#96)

Project.toml

@@ -21,7 +21,7 @@ StatsFuns = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
 ArgCheck = "1, 2.0"
 Combinatorics = "0.7"
 Compat = "3.0"
-Distributions = "0.21.11, 0.22"
+Distributions = "0.21.11, 0.22, 0.23"
 DistributionsAD = "0.4.6"
 ForwardDiff = "0.10.3"
 MappedArrays = "0.2.2"

src/compat/reversediff.jl

@@ -6,16 +6,6 @@ using ReverseDiff: SpecialInstruction
 
 _eps(::Type{<:RTR{T}}) where {T} = _eps(T)
 
-Base.:*(A::RTM{<:Any, D}, B::AbstractTriangular) where {D} = record_mul(A, B, D)
-Base.:*(A::AbstractTriangular, B::RTV{<:Any, D}) where {D} = record_mul(A, B, D)
-Base.:*(A::AbstractTriangular, B::RTM{<:Any, D}) where {D} = record_mul(A, B, D)
-function Base.:*(A::Adjoint{T, <:AbstractTriangular{T}}, B::RTM{<:Any, D}) where {T, D}
-    return record_mul(A, B, D)
-end
-function Base.:*(A::Adjoint{T, <:AbstractTriangular{T}}, B::RTV{<:Any, D}) where {T, D}
-    return record_mul(A, B, D)
-end
-
 function replace_diag(::typeof(log), X::RTM{<:Any, D}) where {D}
     tp = ReverseDiff.tape(X)
     X_value = ReverseDiff.value(X)
@@ -54,7 +44,7 @@ end
     return nothing
 end
 
-@noinline function special_forward_exec!(instruction::SpecialInstruction{typeof(replace_diag)})
+@noinline function ReverseDiff.special_forward_exec!(instruction::SpecialInstruction{typeof(replace_diag)})
     output, input = instruction.output, instruction.input
     out_value = replace_diag(ReverseDiff.value(input[1]), ReverseDiff.value(input[2]))
     ReverseDiff.value!(output, out_value)