Add better optimizer (e.g. l-bfgs)

[meissnereric]

Optimizing GP models with Adam is less than ideal, would be great if we could bring in a better optimizer. It's not clear how exactly that should happen (or if it should even be in the MXFusion library vs. MXNet itself).

The simplest solution is probably just to wrap the scipy optimizer.