'Fill ary with normally distributed random numbers.
Arguments:
ary (Tensor): Tensor to fill with random values
mean (float): Mean value. Default 0
stdv (float): standard deviation value.  Default 1'

def fill_normal(self, ary, mean=0, stdv=1):

    self.pcg.fill_normal(p_gpuarray(ary.shape, ary.dtype, gpudata=ary.gpudata))
    if (not all([(mean == 0), (stdv == 1)])):
        ary[:] = ((ary * stdv) + mean)
