Solved by googlebot in post #2 
              
              
                distributions/normal.py has it: 
        return 0.5 * (1 + torch.erf((value - self.loc) * self.scale.reciprocal() / math.sqrt(2)))