Newton solver to find the mode of the Gaussian (for INLA)

[theorashid]

Part of INLA roadmap #340.

We want to do the Laplace approximation around $\mu$. We need to find it first.

Notes

• jax implementation
• GPy implementation (maybe somewhere in here)