Sampler initiation stuck

You can wrap arbitrary scipy code in a pytensor Op for use in PyMC models. The wrinkle will be to compute gradients :slight_smile: