(x, sigma, N=10, PI=np.pi)
| 19 | |
| 20 | |
| 21 | def grad(x, sigma, N=10, PI=np.pi): |
| 22 | p_ = 0 |
| 23 | for i in tqdm.trange(-N, N + 1): |
| 24 | p_ += (x + 2 * PI * i) / sigma ** 2 * np.exp(-(x + 2 * PI * i) ** 2 / 2 / sigma ** 2) |
| 25 | return p_ |
| 26 | |
| 27 | |
| 28 | def sample(sigma, PI=np.pi): |