Allow PyMC3 to work with GW likelihood
This MR is to allow the PyMC3 sampler to be used with the GW likelihood functions.
It creates a Theano Op to wrap the log likelihood function in, and another Theano Op to numerically calculate the gradient of the log likelihood function (just using finite difference).
Edited by Matthew Pitkin