removed legacy __main__ sections

Author: paolomorettin

src/wmipa/core/polynomial.py

@@ -199,41 +199,3 @@ def _is_integral(v: SupportsInt) -> bool:
         return v.denominator == 1
     else:
         return False
-
-
-if __name__ == "__main__":
-    import pysmt.shortcuts as smt
-
-    x = smt.Symbol("x", REAL)
-    y = smt.Symbol("y", REAL)
-    z = smt.Symbol("z", REAL)
-    vv = [x, y, z]
-
-    """
-
-    p1 = Plus(Times(Real(3), Pow(x,Real(2))), Times(Real(4), Pow(y,Real(7)), z), Real(1))
-    print("(p1)", str(Polynomial(p1, variables)), "\n\n\n")
-
-    p2 = Plus(p1, p1)
-    print("(p2 = p1 + p1)", str(Polynomial(p2, variables)), "\n\n\n")
-
-    p3 = Times(p1, p1)
-    print("(p3 = p1 * p1)", str(Polynomial(p3, variables)), "\n\n\n")
-
-    p4 = Plus(p3, p2, p1)
-    print("(p4 = p3 + p2 + p1)", str(Polynomial(p4, variables)), "\n\n\n")
-
-    """
-
-    p = Polynomial(
-        smt.Plus(smt.Times(smt.Real(3), smt.Pow(x, smt.Real(2))), smt.Times(y, z)),
-        vv,
-        smt.get_env(),
-    )
-    f = p.to_numpy()
-
-    x = np.array(list(range(12))).reshape(-1, 3)
-    print("p:", str(p))
-    print("x:", x)
-
-    print("p(x):", f(x))

src/wmipa/core/polytope.py

@@ -41,24 +41,3 @@ def to_numpy(self) -> tuple[np.ndarray, np.ndarray]:
 
     def __str__(self) -> str:
         return "\n".join(["[" + str(b) + "]" for b in self.inequalities])
-
-
-if __name__ == "__main__":
-
-    import pysmt.shortcuts as smt
-
-    x = smt.Symbol("x", smt.REAL)
-    y = smt.Symbol("y", smt.REAL)
-    variables = [x, y]
-
-    h1 = smt.LE(
-        smt.Plus(smt.Times(smt.Real(3), x), smt.Times(smt.Real(5), y), smt.Real(7)), x
-    )
-    h2 = smt.LE(x, y)
-
-    p = Polytope([h1, h2], variables, smt.get_env())
-    print("p:", p)
-    print("pysmt:", smt.serialize(p.to_pysmt()))
-    A, b = p.to_numpy()
-    print("A:", A)
-    print("b:", b)

src/wmipa/integration/rejection.py

@@ -55,31 +55,3 @@ def integrate_batch(
             volumes.append(self.integrate(polytope, integrand))
 
         return np.array(volumes)
-
-
-if __name__ == "__main__":
-    # avoid polluting the global namespace
-    import pysmt.shortcuts as smt
-    import pysmt.typing as smt_typing
-
-    import wmipa.core as wmipa_ds
-
-    x = smt.Symbol("x", smt_typing.REAL)
-    y = smt.Symbol("y", smt_typing.REAL)
-    env = smt.get_env()
-
-    variables = [x, y]
-
-    h1 = smt.LE(smt.Real(0), x)
-    h2 = smt.LE(smt.Real(0), y)
-    h3 = smt.LE(smt.Plus(x, y), smt.Real(1))
-
-    # w = smt.Minus(smt.Real(1), smt.Plus(x, y))
-    w = smt.Plus(x, y)
-
-    pt = wmipa_ds.Polytope([h1, h2, h3], variables, env=env)
-    pn = wmipa_ds.Polynomial(w, variables, env=env)
-
-    integrator = RejectionIntegrator(100000)
-
-    print("integral:", integrator.integrate(pt, pn))