bugs in faem()

[githubtomtom]

calling fit(FactorAnalysis, x; method = :em) sometimes get:

ERROR: DomainError with -5.4988886678356523e-17:
log will only return a complex result if called with a complex argument. Try log(Complex(x)).
Stacktrace:
 [1] throw_complex_domainerror(f::Symbol, x::Float64)
   @ Base.Math ./math.jl:33
 [2] log(x::Float64)
   @ Base.Math ./special/log.jl:285
 [3] faem(S::Matrix{Float64}, mv::Vector{Float64}, n::Int64; maxoutdim::Int64, tol::Float64, maxiter::Int64)
   @ MultivariateStats ~/.julia/packages/MultivariateStats/HTpHt/src/fa.jl:72
 [4] fit(::Type{MultivariateStats.FactorAnalysis}, X::Matrix{Float64}; method::Symbol, maxoutdim::Int64, mean::Int64, tol::Float64, η::Float64, maxiter::Int64)
   @ MultivariateStats ~/.julia/packages/MultivariateStats/HTpHt/src/fa.jl:171

it is due to the following line inside faem() where detΣ can sometimes be slightly negative:

        L = (-n/2)*(d*log(2π) + log(detΣ) + tr(Σ⁻¹*S))