Evgeni Burovski

Commented in the scipy issue. Basically, CubicSpline only guarantees the value of derivatives at the boundary, and its extrapolation mode may be not very useful. However it provides all the...

For me as well

Thanks a lot Nicolas! It did help, https://github.com/ReScience/submissions/issues/42

**Paper number**: 36 **Reference**: Evgeni Burovski, Nikolay Prokof’ev, Boris Svistunov, and Matthias Troyer Phys. Rev. Lett. 96, 160402 (2006). https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.96.160402 **Open Access version**: https://arxiv.org/abs/cond-mat/0602224 **Article age**: 14 years **Confidence in...

Not sure what the right way is, can probably move the loop to cython + call a relevant lapack function via cython_lapack. If that's what you're after, is easiest to...

Re accuracy: hard to tell in general, depends on the quantity.Offhand, 1e-14 sounds ok.

ISTM (superficially, without looking at any code) that this pbinom would better fit among multivariate distributions. They are in _multivariate.py, a simple example is e.g. in https://github.com/scipy/scipy/pull/5410

@josef-pkt Indeed. Then it might be easier to register it as having no shape parameters and set the probabilities in the constructor.

Something like this should work: ``` In [26]: class pbinom(rv_discrete): def __init__(self, prob, seed=None): self.prob = np.asarray(prob) # validate etc super(pbinom, self).__init__(seed=seed, a=0, b=len(prob)) def _pmf(self, x): ... can use...

It should be relatively easy to add a higher-order stencil to approx_derivative. IIRC approx_derivative does not broadcast for scalar functions of a scalar arguments. In retrospect, it's unfortunate that it...