mathai 0.3.9__py3-none-any.whl → 0.4.0__py3-none-any.whl

mathai/__init__.py

@@ -36,6 +36,8 @@ from .console import console
 
 from .limit import limit
 
+from .search import dfs_simplify as search0
+
 from .univariate_inequality import wavycurvy, absolute, domain, handle_sqrt
 
 from .base import *

mathai/apart.py

@@ -1,18 +1,14 @@
 from .linear import linear_solve
 from .expand import expand
 from .simplify import simplify
-
 from .diff import diff
 from .inverse import inverse
 from .base import *
 import math
 from .tool import poly, enclose_const
 
-def _apart(eq, v=None):
-    if v is None:
-        if len(vlist(eq)) == 0:
-            return eq
-        v = vlist(eq)[0]
+def _apart(eq, v="v_0"):
+    
     origv = vlist(eq)
     eq = simplify(eq)
     if eq.name != "f_mul":
@@ -46,9 +42,8 @@ def _apart(eq, v=None):
     s = []
     facd = [simplify(x) for x in factor_generation(simplify(d))]
     
-    
-    facd2 = remove_duplicates_custom(facd, lambda m, n: simplify(expand(simplify(m-n))) == tree_form("d_0"))
-
+     
+    facd2 = remove_duplicates_custom(facd, lambda m, n: simplify(m-n) == tree_form("d_0"))
     if len(facd2) == 1:
         return eq
     x = tree_form(v)
@@ -63,7 +58,6 @@ def _apart(eq, v=None):
             if n > 2:
                 return eq
             n = tree_form("d_"+str(n))
-            
             l = len(poly(item, v))
             if l == 3:
                 a = alloclst.pop(0)
@@ -99,15 +93,14 @@ def _apart(eq, v=None):
     final = summation(final2)
    
     s = simplify(TreeNode("f_eq", [final-eq2, tree_form("d_0")]))
-
+    
     lst = poly(s.children[0], v)
-
+    
     lst = [TreeNode("f_eq", [item, tree_form("d_0")]) for item in lst if "v_" in str_form(item)]
     lst2 = []
     for item in lst:
         lst2+=vlist(item)
     origv = list(set(lst2)-set(origv))
-    
     out = linear_solve(TreeNode("f_and", lst), [tree_form(item) for item in origv])
     for item in out.children:
         
@@ -115,9 +108,6 @@ def _apart(eq, v=None):
     return simplify(final3)
 def apart(eq):
     eq, fx = enclose_const(eq)
-    def helper(eq):
-        eq2 = _apart(eq)
-        if eq != eq2:
-            return eq2
-        return TreeNode(eq.name, [helper(child) for child in eq.children])
-    return fx(helper(eq))
+    
+    eq = _apart(eq)
+    return fx(eq)

mathai/base.py

@@ -131,8 +131,6 @@ def frac(eq):
         return Fraction(int(eq.name[2:]))
     if eq.name == "f_add":
         p = frac(eq.children[0])
-        if p is None:
-            return None
         for child in eq.children[1:]:
             tmp = frac(child)
             if isinstance(tmp, Fraction):
@@ -142,8 +140,6 @@ def frac(eq):
         return p
     if eq.name == "f_mul":
         p = frac(eq.children[0])
-        if p is None:
-            return None
         for child in eq.children[1:]:
             tmp = frac(child)
             if isinstance(tmp, Fraction):

mathai/integrate.py

@@ -104,6 +104,8 @@ def handle_try(eq):
     else:
         return TreeNode(eq.name, [handle_try(child) for child in eq.children])
 def inteq(eq):
+    if "f_ref" not in str_form(eq):
+        return eq
     if eq.name == "f_try":
         eq2 = None
         output = []
@@ -111,19 +113,13 @@ def inteq(eq):
             if child.name == "f_ref":
                 eq2 = child.children[0]
                 break
-        if eq2 is None:
-            return eq
         for child in eq.children:
             if child.name == "f_ref":
                 output.append(child)
             else:
                 eq3 = simplify(expand(simplify(eq2 - child)))
                 if contain(eq3, eq2):
-                    out = inverse(eq3, str_form(eq2))
-                    if out is None:
-                        output.append(child)
-                    else:
-                        output.append(out)
+                    output.append(inverse(eq3, str_form(eq2)))
                 else:
                     output.append(child)
         return TreeNode("f_try", output)
@@ -138,8 +134,7 @@ def solve_integrate(eq):
     eq2 = dowhile(eq, _solve_integrate)
     eq2 = dowhile(eq2, handle_try)
     eq2 = rm(eq2)
-    if eq2.name == "f_try":
-        eq2.children = list(set(eq2.children))
+    eq2.children = list(set(eq2.children))
     return eq2
 def integrate_subs(equation, term, v1, v2):
     output = []
@@ -170,7 +165,7 @@ def integrate_subs(equation, term, v1, v2):
         
         return none
 
-    return TreeNode("f_subs", [TreeNode("f_integrate", [simplify(fraction(expand(simplify(equation)))), tree_form(origv2)]),tree_form(origv2) ,g])
+    return TreeNode("f_subs", [TreeNode("f_integrate", [simplify(expand(simplify(equation))), tree_form(origv2)]),tree_form(origv2) ,g])
 
 def integrate_subs_main(equation):
     if equation.name == "f_ref":
@@ -381,7 +376,6 @@ def rm_const(equation):
                 return rm_const(TreeNode("f_integrate",[equation, wrt])) *const
         equation = eq2
     return TreeNode(equation.name, [rm_const(child)  for child in equation.children])
-
 def integrate_formula(equation):
     if equation.name == "f_ref":
         return equation.copy_tree()
@@ -390,18 +384,17 @@ def integrate_formula(equation):
         integrand = eq2.children[0]
         wrt = eq2.children[1]
         if integrand == wrt:
-            return wrt**2/2  # x^2/2
+            return TreeNode("f_add", [TreeNode("f_power", [wrt.copy_tree(), TreeNode("2")]), TreeNode("f_div", [TreeNode("1"), TreeNode("2")])])  # x^2/2
         if not contain(integrand, wrt):
-            return integrand*wrt
+            return TreeNode("f_mul", [wrt.copy_tree(), integrand.copy_tree()])  # constant * dx
         out = transform_formula(simplify(trig0(integrand)), wrt.name, formula_gen[0], formula_gen[1], formula_gen[2])
         if out is not None:
-            
             return out
         expr_str = str_form(integrand)
         if expr_str.count("f_sin") + expr_str.count("f_cos") > 2:
             out = transform_formula(integrand, wrt.name, formula_gen4[0], formula_gen4[1], formula_gen4[2])
             if out is not None:
-                
+                print(out, 111)
                 return out
     return TreeNode(eq2.name, [integrate_formula(child) for child in eq2.children])
 

mathai/linear.py

@@ -3,9 +3,6 @@ from .simplify import simplify, solve
 from .fraction import fraction
 from .expand import expand
 from .base import *
-from .factor import factorconst
-def ss(eq):
-    return dowhile(eq, lambda x: fraction(expand(simplify(x))))
 def rref(matrix):
     rows, cols = len(matrix), len(matrix[0])
     lead = 0
@@ -13,7 +10,7 @@ def rref(matrix):
         if lead >= cols:
             return matrix
         i = r
-        while ss(matrix[i][lead]) == tree_form("d_0"):
+        while fraction(simplify(matrix[i][lead])) == tree_form("d_0"):
             i += 1
             if i == rows:
                 i = r
@@ -22,11 +19,11 @@ def rref(matrix):
                     return matrix
         matrix[i], matrix[r] = matrix[r], matrix[i]
         lv = matrix[r][lead]
-        matrix[r] = [ss(m / lv) for m in matrix[r]]
+        matrix[r] = [fraction(simplify(m / lv)) for m in matrix[r]]
         for i in range(rows):
             if i != r:
                 lv = matrix[i][lead]
-                matrix[i] = [ss(m - lv * n) for m, n in zip(matrix[i], matrix[r])]
+                matrix[i] = [fraction(simplify(m - lv * n)) for m, n in zip(matrix[i], matrix[r])]
         lead += 1
     return matrix
 def islinear(eq, fxconst):
@@ -94,7 +91,7 @@ def linear(eqlist, fxconst):
     for i in range(len(m)):
         for j in range(len(m[i])):
             m[i][j] = fraction(m[i][j])
-
+    #print(m)
     for item in m:
         if all(item2==tree_form("d_0") for item2 in item[:-1]) and item[-1] != tree_form("d_0"):
             return tree_form("s_false")

mathai/search.py

@@ -0,0 +1,117 @@
+from mathai import *
+import copy
+from concurrent.futures import ThreadPoolExecutor, TimeoutError
+
+def dfs_simplify(equation, functions, true_expr, false_expr,
+                 max_timeout=25, max_small=4,
+                 base_timeout=1, time_per_char=0.1, timeout_increase=0.5):
+    """
+    Perform DFS simplification on a given equation using provided functions.
+
+    Args:
+        equation: The starting expression (TreeNode or parsed equation)
+        functions: List of simplification functions
+        true_expr: Expression representing True (immediate termination)
+        false_expr: Expression representing False (immediate termination)
+        max_timeout: Maximum timeout allowed for any function
+        max_small: Number of smallest expressions to track
+        base_timeout: Base timeout in seconds
+        time_per_char: Additional timeout per character of expression
+        timeout_increase: Factor to increase timeout for consecutive timeouts
+
+    Returns:
+        tuple(found_boolean, boolean_path, smallest_expressions)
+    """
+    original_eq = simplify(equation)
+    smallest_four = []
+
+    stack = [(copy.deepcopy(original_eq), [copy.deepcopy(original_eq)])]
+    visited = set()
+
+    found_boolean = False
+    boolean_path = None
+    boolean_expr = None
+
+    executor = ThreadPoolExecutor(max_workers=3)
+    consecutive_timeouts = 0
+
+    while stack and not found_boolean:
+        current_eq, path = stack.pop()
+        expr_str = str(current_eq)
+
+        if expr_str in visited:
+            continue
+        visited.add(expr_str)
+
+        # Thinking message
+        printeq(current_eq)
+
+        # Immediate termination using predicate functions
+        if true_expr(current_eq):
+            found_boolean = True
+            boolean_path = path
+            boolean_expr = current_eq
+            break
+        if false_expr(current_eq):
+            found_boolean = True
+            boolean_path = path
+            boolean_expr = current_eq
+            break
+
+
+        # Insert into smallest_four if qualifies
+        inserted = False
+        for j in range(len(smallest_four)):
+            if len(expr_str) < len(str(smallest_four[j][0])):
+                smallest_four.insert(j, (copy.deepcopy(current_eq), copy.deepcopy(path)))
+                inserted = True
+                break
+        if not inserted and len(smallest_four) < max_small:
+            smallest_four.append((copy.deepcopy(current_eq), copy.deepcopy(path)))
+        if len(smallest_four) > max_small:
+            smallest_four = smallest_four[:max_small]
+
+        # Calculate adaptive timeout with cap
+        timeout = (base_timeout + time_per_char * len(expr_str)) * (1 + timeout_increase * consecutive_timeouts)
+        if timeout > max_timeout:
+            timeout = max_timeout
+
+        # Try functions that reduce length first
+        reduced_any = False
+        for fx in functions:
+            print(f"[Thinking] Executing {fx.__name__} on current expression (timeout={timeout:.2f}s):")
+            printeq(current_eq)
+            future = executor.submit(fx, current_eq)
+            try:
+                new_expr = future.result(timeout=timeout)
+                new_expr_str = str(new_expr)
+                if len(new_expr_str) <= len(expr_str) and new_expr_str != expr_str:
+                    reduced_any = True
+                    stack.append((new_expr, path + [copy.deepcopy(new_expr)]))
+                consecutive_timeouts = 0  # reset after success
+            except TimeoutError:
+                print(f"[Thinking] {fx.__name__} timed out, skipping.")
+                consecutive_timeouts += 1
+                continue
+
+        # If no reducing function worked, try growing functions
+        if not reduced_any:
+            for fx in functions:
+                print(f"[Thinking] Trying growing {fx.__name__} on current expression (timeout={timeout:.2f}s):")
+                printeq(current_eq)
+                future = executor.submit(fx, current_eq)
+                try:
+                    new_expr = future.result(timeout=timeout)
+                    new_expr_str = str(new_expr)
+                    if new_expr_str != expr_str:
+                        stack.append((new_expr, path + [copy.deepcopy(new_expr)]))
+                        consecutive_timeouts = 0
+                        break  # only take one growing function
+                except TimeoutError:
+                    print(f"[Thinking] {fx.__name__} (growing) timed out, skipping.")
+                    consecutive_timeouts += 1
+                    continue
+
+    executor.shutdown(wait=True)
+
+    return found_boolean, boolean_path, smallest_four

mathai/simplify.py

@@ -286,23 +286,6 @@ def simplify(eq):
                 error = True
             else:
                 eq = tree_form("d_0")
-        
-        if eq.name == "f_mul" and str_form(eq).count("f_pow")>1:
-            dic = {}
-            for child in eq.children:
-                head = child
-                tail = None
-                if child.name == "f_pow":
-                    head = child.children[0]
-                    tail = child.children[1]
-                if tail is None:
-                    tail = tree_form("d_1")
-                if head not in dic.keys():
-                    dic[head] = tail
-                else:
-                    dic[head] += tail
-            if len(eq.children) != len(dic.keys()):
-                eq = product([key if dic[key] == 1 else key**dic[key] for key in dic.keys()])
         if eq.name == "f_pow" and eq.children[0].name == "f_pow" and eq.children[0].children[1] == tree_form("d_2")**-1 and eq.children[1] == tree_form("d_2"):
             eq = eq.children[0].children[0]
         if (eq.name == "f_sin" and eq.children[0].name == "f_arcsin") or (eq.name == "f_cos" and eq.children[0].name == "f_arccos") or (eq.name == "f_tan" and eq.children[0].name == "f_arctan"):

mathai/trig.py

@@ -199,37 +199,16 @@ def trig4(eq, numer=True):
             return tree_form("d_1")/(1+a**2)**(tree_form("d_2")**-1)
     
     return TreeNode(eq.name, [trig4(child, False) if not numer or (eq.name == "f_pow" and frac(eq.children[1]) is not None and frac(eq.children[1]) < 0) else trig4(child, True) for child in eq.children])
-
 def trig2(eq):
     if eq.name == "f_add":
         for item in itertools.combinations(range(len(eq.children)), 2):
-            child1, child2 = eq.children[item[0]], eq.children[item[1]]
-
-            # Check if both are sin or cos
-            if child1.name in ["f_sin", "f_cos"] and child2.name in ["f_sin", "f_cos"]:
-                a, b = child1.children[0], child2.children[0]
-                
-                # Compute the rest of the sum
-                rest = [eq.children[i] for i in range(len(eq.children)) if i not in item]
-                if len(rest) == 0:
-                    rest_tree = tree_form("d_0")
-                else:
-                    rest_tree = summation(rest)
-
-                # Now handle the sin/cos combination formula
-                if child1.name == "f_sin" and child2.name == "f_sin":
-                    # sin A + sin B = 2 sin((A+B)/2) cos((A-B)/2)
-                    two = tree_form("d_2")
-                    combined = two * ((a + b) / two).fx("sin") * ((a - b) / two).fx("cos")
-                elif child1.name == "f_cos" and child2.name == "f_cos":
-                    # cos A + cos B = 2 cos((A+B)/2) cos((A-B)/2)
-                    two = tree_form("d_2")
-                    combined = two * ((a + b) / two).fx("cos") * ((a - b) / two).fx("cos")
+            if all(eq.children[item2].name == "f_sin" for item2 in item):
+                a, b = eq.children[item[0]].children[0], eq.children[item[1]].children[0]
+                rest = [item2 for index, item2 in enumerate(eq.children) if index not in item]
+                if len(rest)==0:
+                    rest = tree_form("d_0")
                 else:
-                    # sin A + cos B = sin A + cos B (leave unchanged, or implement formula if desired)
-                    continue  # skip for now, keep original
-
-                return rest_tree + combined
-
-    # Recurse for other nodes
+                    rest = summation(rest)
+                two = tree_form("d_2")
+                return rest + two*((a+b)/two).fx("sin")*((a-b)/two).fx("cos")
     return TreeNode(eq.name, [trig2(child) for child in eq.children])

{mathai-0.3.9.dist-info → mathai-0.4.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mathai
-Version: 0.3.9
+Version: 0.4.0
 Summary: Mathematics solving Ai tailored to NCERT
 Home-page: https://github.com/infinity390/mathai4
 Author: educated indians are having a low iq and are good for nothing

{mathai-0.3.9.dist-info → mathai-0.4.0.dist-info}/RECORD

@@ -1,24 +1,25 @@
-mathai/__init__.py,sha256=hxMCrvwj98yE5ChWSzvkn0gMxqAQnh8oze3ZlImWyD0,1170
-mathai/apart.py,sha256=8IlJ8X6SKAjPunUPHVLmgNB5GuEi4_XpY0oLkQqrHKc,3770
-mathai/base.py,sha256=zJqozbSeesbxw8YJ_eNXmaoy8K5VmkPHECnuAFGRvW4,12683
+mathai/__init__.py,sha256=GsDxDy6iWjoJ7CJty2nzE4Jx03mZkHFQW6hkmPJHkPE,1217
+mathai/apart.py,sha256=Fhl-OHZzPGIOj8YFyGQddjuh-2HV7uVDhkIzQ-YoMO0,3492
+mathai/base.py,sha256=x4Rjgz5mVblclsycKtDJjSudFpyIBWSEUX83_120BZY,12587
 mathai/console.py,sha256=Sn58iwYE79MLEh67s8X3q6vZjw6g7f9XM1T8_dBBR2o,3048
 mathai/diff.py,sha256=YUBpRsz0qmBkq5vGxeGnvR4nMKjdOQiIXlNMxpij2ns,3051
 mathai/expand.py,sha256=SnBltkpIENMGkP0AYmbMlSc4H-CF5RslO2PcBEkn1BQ,3359
 mathai/factor.py,sha256=NPXxET52TacNExuvw6p1jbC6g3wY6_VOCdlGlexXZio,5916
 mathai/fraction.py,sha256=Q2ztsh5Bpz6YhML2QU0tfufbAs0Q6J319AhlzKephIY,4396
-mathai/integrate.py,sha256=OH_DqnbgoCx7e3gq4u1zEz9ABcaiUvH5t75ap67t4GI,15034
+mathai/integrate.py,sha256=UKR_GH-ifplq_VyJKZukWjVDsQKkXXIuKNhpZ4w1Mhk,15033
 mathai/inverse.py,sha256=QCvDrzKquWsZv-BDAzZd9HnU0c3gZvcc44UztHVO5LQ,2919
 mathai/limit.py,sha256=RA8YAehgYCGVWv9qBc8uQ34BQ9mFthWl2OrVTwcHl2g,4920
-mathai/linear.py,sha256=BzSnm941Zlod_l6hON4Rs6J4pdAA3MGpRVqr6-66ZBk,5524
+mathai/linear.py,sha256=wyiLpIxRDmD96xXktkgvc5gegIB3zblLB1EuV3TFdmU,5474
 mathai/logic.py,sha256=UvHzRmKcO9AD51tRzHmpNSEhgW5gmaf4XPaQKFjGfC4,9653
 mathai/parser.py,sha256=f7bemieFmp0sbup1NlraMLvZDVFvqKGFknEVtlFRMVk,6979
 mathai/printeq.py,sha256=gIes-pstFOa6FcnpVIVvkjVKuWdsVdo11LlEnmHhakU,1303
-mathai/simplify.py,sha256=4gJgEh7fmwJ3aKMliB4d1NalGAx9duIYA7bGAIYUX2Y,15858
+mathai/search.py,sha256=6S_V0PrC2ZnbhkWyT2ZTtExRD_mDlo-C3p0TDn6Dk3w,4743
+mathai/simplify.py,sha256=F37h-Z_rW35uVgx0G86vQErU2Ac6ZiTBN9KcC3RzDkg,15156
 mathai/structure.py,sha256=4Ww2IAx62RcQSO7_17TZES-DjMWBpcFQtL939FBIHwY,4103
 mathai/tool.py,sha256=UyccamiJy_CkFPakfufyPzdhtlEO6v2D7qwbXQ9V7Rg,2000
-mathai/trig.py,sha256=Zadw-uTJaMfoXwxuiio7QYHWgPxezEP2xZ2fr__PFM0,10290
+mathai/trig.py,sha256=5P0RNS4eetNds2l-wyA4BAKqJdFIUws_8RGS_dH7gp8,9348
 mathai/univariate_inequality.py,sha256=_r-kkiS4Hr-jRN7f-EL_E4svAMFWJP1Ea50HJKKbjfk,14778
-mathai-0.3.9.dist-info/METADATA,sha256=sWapzjSv4xLp5xqwXkpt0_lqv32BcPNlQvpYmwi1tiU,7021
-mathai-0.3.9.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-mathai-0.3.9.dist-info/top_level.txt,sha256=ROP4l3OhGYw3ihkQGASr18xM9GsK4z3_6whV5AyXLwE,7
-mathai-0.3.9.dist-info/RECORD,,
+mathai-0.4.0.dist-info/METADATA,sha256=SB2hfiau23qNhjHplLEs5NzT7g3vBlW63YLJnh96Iq0,7021
+mathai-0.4.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+mathai-0.4.0.dist-info/top_level.txt,sha256=ROP4l3OhGYw3ihkQGASr18xM9GsK4z3_6whV5AyXLwE,7
+mathai-0.4.0.dist-info/RECORD,,