mplusa 0.0.1__py3-none-any.whl → 0.0.2__py3-none-any.whl

mplusa/maxplus.py

@@ -4,10 +4,18 @@ from numbers import Real
 
 
 def add(*args) -> Real:
+    if math.inf in args:
+        raise ValueError(
+            'Maxplus.add: value out of domain.'
+        )
     return max(args)
 
 
 def mult(*args) -> Real:
+    if math.inf in args:
+        raise ValueError(
+            'Maxplus.add: value out of domain.'
+        )
     return sum(args) if -math.inf not in args else -math.inf
 
 
@@ -44,6 +52,10 @@ def mult_matrices(A : np.ndarray,
 
 def modulo(a : Real,
            t : int) -> Real:
+    if a < 0 or t < 0:
+        raise ValueError(
+            'Maxplus.modulo: modulo operation is only defined for positive numbers.'
+        )
     if a == -math.inf:
         return -math.inf
     if a == 0:
@@ -58,7 +70,7 @@ def modulo_matrices(A : np.ndarray,
     if b.shape[1] != 1:
         raise ValueError(
             'Maxplus.modulo_matrices: given matrix b ' +\
-            'is not a properly formated vector (has shape of {}).'.format(
+            'is not a vertical vector of shape Mx1 (has shape of {}).'.format(
                 b.shape
             )
         )
@@ -89,7 +101,7 @@ def power_matrix(A : np.ndarray,
                  k : int) -> np.ndarray:
     if np.any(np.diagonal(A) != 0):
         raise ValueError(
-            'Maxplus.power_matrix: matrix contains non-zero values on diagonal.'
+            'Maxplus.power_matrix: matrix contains non-zero values on the diagonal.'
         )
     if k == 0:
         result = unit_matrix(A.shape[0], A.shape[1])
@@ -150,12 +162,18 @@ class Polynomial:
     def __call__(self, x : float) -> float:
         return add(*[mult(coefficient, power(x, i)) for i, coefficient in enumerate(self.coefficients)])
 
+    def get_lines(self) -> list[tuple[float]]:
+        """ Returns the a and b values of standard linear functions building the polynomial in form of y = ax + b. """
+        return [(a, b) for a, b in enumerate(self.coefficients) if b > -math.inf]
+
     def get_hypersurface(self) -> list[float]:
+        lines = self.get_lines()
         result = []
-        candidates = [c2 - c1 for c1, c2 in zip(self.coefficients[1:], self.coefficients[:-1])]
-        for candidate in candidates:
-            if not isinstance(candidate, Real) or candidate == math.inf:
-                continue
-            if abs(self(candidate) - self(candidate - 1)) != abs(self(candidate) - self(candidate + 1)):
-                result.append(candidate)
-        return result
+        for (a, c) in lines:
+            for (b, d) in lines:
+                if a == b or c == d:
+                    continue
+                x = (d - c) / (a - b)
+                if a * x + c == self(x):
+                    result.append(x)
+        return list(set(result))

mplusa/minplus.py

@@ -4,10 +4,18 @@ from numbers import Real
 
 
 def add(*args) -> Real:
+    if -math.inf in args:
+        raise ValueError(
+            'Minplus.add: value out of domain.'
+        )
     return min(args)
 
 
 def mult(*args) -> Real:
+    if -math.inf in args:
+        raise ValueError(
+            'Minplus.mult: value out of domain.'
+        )
     return sum(args) if math.inf not in args else math.inf
 
 
@@ -44,6 +52,10 @@ def mult_matrices(A : np.ndarray,
 
 def modulo(a : Real,
            t : int) -> Real:
+    if a < 0 or t < 0:
+        raise ValueError(
+            'Minplus.modulo: modulo operation is only defined for positive numbers.'
+        )
     if a == math.inf:
         return math.inf
     if a == 0:
@@ -58,7 +70,7 @@ def modulo_matrices(A : np.ndarray,
     if b.shape[1] != 1:
         raise ValueError(
             'Minplus.modulo_matrices: given matrix b ' +\
-            'is not a properly formated vector (has shape of {}).'.format(
+            'is not a vertical vector of shape Mx1 (has shape of {}).'.format(
                 b.shape
             )
         )
@@ -89,7 +101,7 @@ def power_matrix(A : np.ndarray,
                  k : int) -> np.ndarray:
     if np.any(np.diagonal(A) != 0):
         raise ValueError(
-            'Minplus.power_matrix: matrix contains non-zero values on diagonal.'
+            'Minplus.power_matrix: matrix contains non-zero values on the diagonal.'
         )
     if k == 0:
         result = unit_matrix(A.shape[0], A.shape[1])
@@ -150,12 +162,18 @@ class Polynomial:
     def __call__(self, x : float) -> float:
         return add(*[mult(coefficient, power(x, i)) for i, coefficient in enumerate(self.coefficients)])
 
+    def get_lines(self) -> list[tuple[float]]:
+        """ Returns the a and b values of standard linear functions building the polynomial in form of y = ax + b. """
+        return [(a, b) for a, b in enumerate(self.coefficients) if b < math.inf]
+
     def get_hypersurface(self) -> list[float]:
+        lines = self.get_lines()
         result = []
-        candidates = [c2 - c1 for c1, c2 in zip(self.coefficients[1:], self.coefficients[:-1])]
-        for candidate in candidates:
-            if not isinstance(candidate, Real) or candidate == -math.inf:
-                continue
-            if abs(self(candidate) - self(candidate - 1)) != abs(self(candidate) - self(candidate + 1)):
-                result.append(candidate)
-        return result
+        for (a, c) in lines:
+            for (b, d) in lines:
+                if a == b or c == d:
+                    continue
+                x = (d - c) / (a - b)
+                if a * x + c == self(x):
+                    result.append(x)
+        return list(set(result))

{mplusa-0.0.1.dist-info → mplusa-0.0.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: mplusa
-Version: 0.0.1
+Version: 0.0.2
 Summary: A library for calculations in tropical and arctic semirings.
 Author-email: "Maksymilian W." <maksymilian3563@gmail.com>
 License: Copyright (c) 2025 Maksymilian Wiekiera
@@ -41,24 +41,34 @@ After having installed the library one can import one of the two modules the pac
 
 **`add(*args) -> Real`**
 	Tropical addition. Essentially an alias for Python's `min` function.
+
 **`mult(*args) -> Real`**
 	Tropical multiplication. Essentially an alias for Python's `sum` function.
+
 **`add_matrices(A : np.ndarray, B : np.ndarray) -> np.ndarray`**
 	Tropical addition of NumPy arrays. The summed matrices have to be of the same shape.
+
 **`mult_matrices(A : np.ndarray, B : np.ndarray) -> np.ndarray`**
 	Tropical multiplication of NumPy arrays. The multiplied matrices have to be of sizes MxN and NxP and their order matters. The result is of shape MxP.
+
 **`modulo(a : Real, t : int) -> Real`**
 	Tropical modulo operator. It can be understood as the difference between the number $a$ and $t^k$ where $k$ is the largest integer that satisfies $a \geq t^k$.
+
 **`modulo_matrices(A : np.ndarray, b : np.ndarray) -> np.ndarray`**
 	Tropical modulo operator for NumPy arrays. The input matrices should be of size MxN and Mx1. The result is an MxN matrix.
+
 **`power(a : real, k : int) -> Real`**
 	Tropical power operator.  Applies the multiplication k times.
+
 **`power_matrix(A : np.ndarray, k : int) -> np.ndarray`**
 	Tropical power operator for NumPy arrays. It multiplies the matrix k times.
+
 **`unit_matrix(width : int, height : int) -> np.ndarray`**
 	Creates a tropical unit matrix of given width and height.
+
 **`star(A : np.ndarray) -> np.ndarray`**
 	Definition of a unique operator of tropical algebra, usually denoted as $\mathbf{A}^*$. It returns the value to which an infinite recursive sum of matrices converges. The input matrix has to be square and the series created in the process of calculating the value needs to be convergent.
+
 **`Polynomial(*coefficients)`**
 	This is a class that implements basic single-variable tropical polynomials. Calling an object of this class allows to take a value the polynomial takes at a given point, it also implements function `get_hypersurface` which returns a list of its roots.
 

mplusa-0.0.2.dist-info/RECORD

@@ -0,0 +1,8 @@
+mplusa/__init__.py,sha256=8D1JvAZfjHkQ1DOqrxLyXWMEfROgGwf8Y41-rSX51cY,44
+mplusa/maxplus.py,sha256=UnFh0lfbvarUfDHH7JQpQY5KC0_IP9uzgYYnbOkd3CY,5461
+mplusa/minplus.py,sha256=uNsvTNjjlMO4_zuFI5S2z2HP07AKWXnCWOItw2fE5Js,5460
+mplusa-0.0.2.dist-info/LICENSE,sha256=x1S-x_tM1tAmndGsdQKT4DU9LnstfQRgS4befak4XdA,1063
+mplusa-0.0.2.dist-info/METADATA,sha256=2JhnKWBWTaO2dDThFq7DdlUfbCsLMObmVr2O49qqlsE,5100
+mplusa-0.0.2.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
+mplusa-0.0.2.dist-info/top_level.txt,sha256=W5b7P8CkZ-DB3-2K0Rcf0T0tpEuWLCbZBaVgrd7FEcM,7
+mplusa-0.0.2.dist-info/RECORD,,

{mplusa-0.0.1.dist-info → mplusa-0.0.2.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.9.1)
+Generator: setuptools (76.0.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

mplusa-0.0.1.dist-info/RECORD

@@ -1,8 +0,0 @@
-mplusa/__init__.py,sha256=8D1JvAZfjHkQ1DOqrxLyXWMEfROgGwf8Y41-rSX51cY,44
-mplusa/maxplus.py,sha256=q2GOM-VDGNCXEL-Gy7QXZPgU8ZyPUh89XEB8o2QNUdA,4931
-mplusa/minplus.py,sha256=QvYJ-Zg3pC134k9Qiq6qQ1SnXdE5CgjA9odcbo55uY8,4929
-mplusa-0.0.1.dist-info/LICENSE,sha256=x1S-x_tM1tAmndGsdQKT4DU9LnstfQRgS4befak4XdA,1063
-mplusa-0.0.1.dist-info/METADATA,sha256=zS3zcKjEyYylkN6VbcTf_awLCIDrxMrNYcAU31yt29c,5090
-mplusa-0.0.1.dist-info/WHEEL,sha256=EaM1zKIUYa7rQnxGiOCGhzJABRwy4WO57rWMR3_tj4I,91
-mplusa-0.0.1.dist-info/top_level.txt,sha256=W5b7P8CkZ-DB3-2K0Rcf0T0tpEuWLCbZBaVgrd7FEcM,7
-mplusa-0.0.1.dist-info/RECORD,,