Update nd.py

estimator/nd.py

@@ -301,6 +301,34 @@ def CenteredBinomial(eta, n=None):
             bounds=(-eta, eta),
             tag="CenteredBinomial",
         )
+
+    @staticmethod
+    def TUniform(b, n=None):
+        """
+        TUniform distribution ∈ ``[-2^b,2^b]``, endpoints inclusive.
+        This distribution samples the two endpoints with a probability of
+        1/2^(b+2) and the intermediate points with a probability of
+        1/2^(b+1).
+
+        e.g. TUniform(0) samples -1,1 each with a probability 1/4 and 0
+        with a probability 1/2.
+
+        EXAMPLE::
+
+            >>> from estimator.nd import NoiseDistribution as ND
+            >>> ND.TUniform(0)
+            >>> ND.Uniform(10)
+
+        """
+        if b < 0:
+            raise ValueError("TUniform parameter b needs to be larger than zero")
+        mean = 0
+        stddev = sqrt(2**(2*b + 1) / 6)
+        density = 2**(b+1) / (2**(b+1) + 1)
+
+        return NoiseDistribution(
+            n=n, stddev=stddev, mean=mean, bounds=(-2**b, 2**b), density=density, tag="TUniform"
+        )
 
     @staticmethod
     def Uniform(a, b, n=None):