Simple generation of locally constrained values in `numpyro`

Simple generation of locally constrained values in numpyro Just a simple experiment using numpyro to simulate values from a function with locally constrained values. This is often achieved using Gaussian processed but I thought it would be interesting to try something a bit more intuitive, plus priors on the points themselves and their interdependence can be anything, not just normal…