Newton's method in inversion.solvers

Implemented Newton's method in solvers. The stats of the gradient
descent are recorded in an optional stats dict.