Hi,
The current icdf in torch.distributions.Normal 36 seems to be numerically unstable compared to scipy.stats.norm 2.
I implemented a torch version by translating the scipy source code to python and it works the same as the scipy version. I wonder if we could integrate the new code to the current torch.distributions.Normal. Thank you!