radioactiveshrimp 0.1.1__cp39-abi3-manylinux_2_28_x86_64.whl → 0.1.2__cp39-abi3-manylinux_2_28_x86_64.whl

radioactiveshrimp/distributions/__init.py__

@@ -0,0 +1,5 @@
+"""
+Distribution subpackage containing distrobution functions.
+"""
+from .cvdistributions import uniform, exponentialdist, poissondist
+__all__=['uniform,exponentialdist, poissondist']

radioactiveshrimp/distributions/cvdistributions.py

@@ -0,0 +1,111 @@
+import secrets
+import numpy as np
+# import matplotlib.pyplot as plt
+import math as m
+
+def uniform(a:float = 0.0, b:float = 1.0)->float:
+    """ 
+    Generate cryptographically Secure uniform sample
+
+    Args:
+        a (float): lower bound value (default 0.0)
+        b (float): upper bound value (default 1.0)
+
+    Returns:
+        float representing uniform sample
+    
+    Example:
+        >>> uniform()
+        0.3546927444895387
+        >>> uniform()
+        0.21098189931876055
+    """
+    # 53 random bits gives 53-bit precision double
+    u = secrets.randbits(53)/(1<<53) # in [0,1)
+    return a+(b-a)*u
+
+def exponentialdist(lamb):
+    """ 
+    Generate exponentially distributed sample
+
+    Args:
+        lamb (float or int >0): lambda value of exponential distrobution
+
+    Returns:
+        x: float representing the exponentially distributed sample
+
+    Raises:
+        ValueError: if lambda is not greater than 0
+    
+    Example:
+        >>> uexponentialdist(1)
+        0.6839328417240588
+        >>> uexponentialdist(1)
+        2.901353101723147
+        >>> exponentialdist(.5)
+        2.451134119865936
+    """
+    if lamb<=0:
+        raise ValueError('Lambda must be greater than 0')
+    
+    y = uniform()
+    x = -(1/lamb)*np.log(y)
+    return x
+
+
+def poissondist(lamb):
+    """ 
+    Generate poisson distributed sample
+
+    Args:
+        lamb (int or float > 0): lambda value of poisson distribution
+
+    Returns:
+        i: int representing poisson distributed sample
+    
+    Example:
+        >>> poissondist(1)
+        1
+        >>> poissondist(1)
+        3
+        >>> poissondist(10)
+        9
+    """
+    if lamb<=0:
+        raise ValueError('Lambda must be greater than 0')
+    
+    elamb = np.exp(-lamb)
+    i = 0
+    y = uniform()
+    prob = elamb
+    while y>prob:
+        i = i+1
+        factorial = m.factorial(i)
+        power = lamb**i
+        # print(y,i,elamb, power, factorial)
+        prob = prob + (power/factorial)*elamb
+
+    return i
+
+
+# # draw histogram of 100,000 randomly distributed samples using exponential distribution
+# samplexs = []
+# for _ in range(100000):
+#     samplexs.append(exponentialdist(1))
+
+# plt.figure()
+# hist = plt.hist(samplexs, bins='auto', orientation='horizontal', histtype='bar')
+# plt.xlabel('count')
+# plt.ylabel('X')
+# plt.show()
+
+# # draw histogram of 100,000 randomly distributed samples using exponential distribution
+# samplexs = []
+# for _ in range(100000):
+#     samplexs.append(poissondist(10))
+
+# plt.figure()
+# hist = plt.hist(samplexs, bins='auto', histtype='bar')
+# plt.xlabel('X')
+# plt.ylabel('count')
+# plt.show()

{radioactiveshrimp-0.1.1.dist-info → radioactiveshrimp-0.1.2.dist-info}/METADATA

@@ -1,6 +1,7 @@
 Metadata-Version: 2.4
 Name: radioactiveshrimp
-Version: 0.1.1
+Version: 0.1.2
+Requires-Dist: matplotlib>=3.10.6
 Requires-Dist: numpy>=2.3.3
 Requires-Dist: torch>=2.8.0
 Requires-Dist: twine>=6.1.0

{radioactiveshrimp-0.1.1.dist-info → radioactiveshrimp-0.1.2.dist-info}/RECORD

@@ -1,11 +1,13 @@
-radioactiveshrimp-0.1.1.dist-info/METADATA,sha256=4miZ8FtVUImI28hoh1lyXSbo-IuHwLBUmaK838D378Y,378
-radioactiveshrimp-0.1.1.dist-info/WHEEL,sha256=VCoUqQjnl5wExkC-2PG3I2s7-FBXYum-vf1jh8qiYg8,97
+radioactiveshrimp-0.1.2.dist-info/METADATA,sha256=eIFvbRoU_QV2ADRG--hldQQI256vfOcPQ6leK0jYZXI,412
+radioactiveshrimp-0.1.2.dist-info/WHEEL,sha256=VCoUqQjnl5wExkC-2PG3I2s7-FBXYum-vf1jh8qiYg8,97
 radioactiveshrimp/__init__.py,sha256=QY4ait-Lo2WQdyx7MU8N6a1v2ZYyAmYRqzYqdNbiavI,101
 radioactiveshrimp/_core.abi3.so,sha256=3IMwsP7UrgNWdKdw2Pl11_VUQF8NqPy1xyyPuzkREdA,489672
 radioactiveshrimp/_core.pyi,sha256=b6oJaUXUzEzqUE5rpqefV06hl8o_JCU8pgKgIIzQgmc,33
 radioactiveshrimp/differential/__init.py__,sha256=jOhQKr6xM_7TVUt0GWRddtkgBvNn5I873WZcArT-GtE,101
 radioactiveshrimp/differential/discrete.py,sha256=9-w6Avz4MZfu8gFhJU8kTvryn9oT2xS1TvvYQ_8pmPU,1371
+radioactiveshrimp/distributions/__init.py__,sha256=hZX2UFt54wJG4zDLuTTTonDSA-anWvfjhDmR-akYcn4,182
+radioactiveshrimp/distributions/cvdistributions.py,sha256=UEiO_RYFQYqHybI_wcOdpOrxHToXRpLFsB58dwvI4vU,2599
 radioactiveshrimp/matrix/__init.py__,sha256=EoerrXlfmd2lI9grb3jxqIhHn-bzLtYhRKyxI0RVgaE,166
 radioactiveshrimp/matrix/elementary.py,sha256=zPiBG24F5zIncrrVWbiXEp1oFnHw2xNKy5quF0GBOuc,6333
 radioactiveshrimp/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-radioactiveshrimp-0.1.1.dist-info/RECORD,,
+radioactiveshrimp-0.1.2.dist-info/RECORD,,