There are couple of problems with the Laplace approximation:
- For Laplace approximation we use Optim.jl package to evaluate derivatives, however in current version of Optim.jl derivative syntax for 1D functions has changed leading to failures in ReactiveMP.jl.
- Current implementation of Laplace is not numerically stable.
There are couple of problems with the Laplace approximation: