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

mathai/__init__.py

@@ -36,8 +36,6 @@ 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,14 +1,18 @@
 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="v_0"):
-    
+def _apart(eq, v=None):
+    if v is None:
+        if len(vlist(eq)) == 0:
+            return eq
+        v = vlist(eq)[0]
     origv = vlist(eq)
     eq = simplify(eq)
     if eq.name != "f_mul":
@@ -42,8 +46,9 @@ def _apart(eq, v="v_0"):
     s = []
     facd = [simplify(x) for x in factor_generation(simplify(d))]
     
-     
-    facd2 = remove_duplicates_custom(facd, lambda m, n: simplify(m-n) == tree_form("d_0"))
+    
+    facd2 = remove_duplicates_custom(facd, lambda m, n: simplify(expand(simplify(m-n))) == tree_form("d_0"))
+
     if len(facd2) == 1:
         return eq
     x = tree_form(v)
@@ -58,6 +63,7 @@ def _apart(eq, v="v_0"):
             if n > 2:
                 return eq
             n = tree_form("d_"+str(n))
+            
             l = len(poly(item, v))
             if l == 3:
                 a = alloclst.pop(0)
@@ -93,14 +99,15 @@ def _apart(eq, v="v_0"):
     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:
         
@@ -108,6 +115,9 @@ def _apart(eq, v="v_0"):
     return simplify(final3)
 def apart(eq):
     eq, fx = enclose_const(eq)
-    
-    eq = _apart(eq)
-    return fx(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))

mathai/base.py

@@ -131,6 +131,8 @@ 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):
@@ -140,6 +142,8 @@ 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):
@@ -286,6 +290,8 @@ def flatten_tree(node):
         node.children = [flatten_tree(child) for child in node.children]
         return node
 def dowhile(eq, fx):
+    if eq is None:
+        return None
     while True:
         orig = eq.copy_tree()
         eq2 = fx(eq)

mathai/integrate.py

@@ -104,8 +104,6 @@ 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 = []
@@ -113,13 +111,19 @@ 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):
-                    output.append(inverse(eq3, str_form(eq2)))
+                    out = inverse(eq3, str_form(eq2))
+                    if out is None:
+                        output.append(child)
+                    else:
+                        output.append(out)
                 else:
                     output.append(child)
         return TreeNode("f_try", output)
@@ -128,13 +132,14 @@ def inteq(eq):
 def rm(eq):
     if eq.name == "f_try":
         eq = TreeNode(eq.name, list(set(eq.children)))
-    return TreeNode(eq.name, [rm(child) for child in eq.children])
+    return TreeNode(eq.name, [rm(child) for child in eq.children if child is not None])
 def solve_integrate(eq):
     
     eq2 = dowhile(eq, _solve_integrate)
     eq2 = dowhile(eq2, handle_try)
     eq2 = rm(eq2)
-    eq2.children = list(set(eq2.children))
+    if eq2.name == "f_try":
+        eq2.children = list(set(eq2.children))
     return eq2
 def integrate_subs(equation, term, v1, v2):
     output = []
@@ -165,7 +170,7 @@ def integrate_subs(equation, term, v1, v2):
         
         return none
 
-    return TreeNode("f_subs", [TreeNode("f_integrate", [simplify(expand(simplify(equation))), tree_form(origv2)]),tree_form(origv2) ,g])
+    return TreeNode("f_subs", [TreeNode("f_integrate", [simplify(fraction(expand(simplify(equation)))), tree_form(origv2)]),tree_form(origv2) ,g])
 
 def integrate_subs_main(equation):
     if equation.name == "f_ref":
@@ -291,7 +296,7 @@ def _sqint(equation):
     return coll
 
 def sqint(eq):
-    out = _sqint(eq)
+    out = simplify(_sqint(eq))
     if out is None:
         return eq
     return out
@@ -376,6 +381,7 @@ 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()
@@ -384,17 +390,18 @@ def integrate_formula(equation):
         integrand = eq2.children[0]
         wrt = eq2.children[1]
         if integrand == wrt:
-            return TreeNode("f_add", [TreeNode("f_power", [wrt.copy_tree(), TreeNode("2")]), TreeNode("f_div", [TreeNode("1"), TreeNode("2")])])  # x^2/2
+            return wrt**2/2  # x^2/2
         if not contain(integrand, wrt):
-            return TreeNode("f_mul", [wrt.copy_tree(), integrand.copy_tree()])  # constant * dx
+            return integrand*wrt
         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,6 +3,9 @@ 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
@@ -10,7 +13,7 @@ def rref(matrix):
         if lead >= cols:
             return matrix
         i = r
-        while fraction(simplify(matrix[i][lead])) == tree_form("d_0"):
+        while ss(matrix[i][lead]) == tree_form("d_0"):
             i += 1
             if i == rows:
                 i = r
@@ -19,11 +22,11 @@ def rref(matrix):
                     return matrix
         matrix[i], matrix[r] = matrix[r], matrix[i]
         lv = matrix[r][lead]
-        matrix[r] = [fraction(simplify(m / lv)) for m in matrix[r]]
+        matrix[r] = [ss(m / lv) for m in matrix[r]]
         for i in range(rows):
             if i != r:
                 lv = matrix[i][lead]
-                matrix[i] = [fraction(simplify(m - lv * n)) for m, n in zip(matrix[i], matrix[r])]
+                matrix[i] = [ss(m - lv * n) for m, n in zip(matrix[i], matrix[r])]
         lead += 1
     return matrix
 def islinear(eq, fxconst):
@@ -91,7 +94,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/simplify.py

@@ -286,6 +286,23 @@ def simplify(eq):
                 error = True
             else:
                 eq = tree_form("d_0")
+        
+        if eq.name == "f_mul":
+            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,16 +199,37 @@ 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):
-            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")
+            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")
                 else:
-                    rest = summation(rest)
-                two = tree_form("d_2")
-                return rest + two*((a+b)/two).fx("sin")*((a-b)/two).fx("cos")
+                    # 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
     return TreeNode(eq.name, [trig2(child) for child in eq.children])

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mathai
-Version: 0.4.0
+Version: 0.4.1
 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.4.0.dist-info → mathai-0.4.1.dist-info}/RECORD

@@ -1,25 +1,25 @@
-mathai/__init__.py,sha256=GsDxDy6iWjoJ7CJty2nzE4Jx03mZkHFQW6hkmPJHkPE,1217
-mathai/apart.py,sha256=Fhl-OHZzPGIOj8YFyGQddjuh-2HV7uVDhkIzQ-YoMO0,3492
-mathai/base.py,sha256=x4Rjgz5mVblclsycKtDJjSudFpyIBWSEUX83_120BZY,12587
+mathai/__init__.py,sha256=hxMCrvwj98yE5ChWSzvkn0gMxqAQnh8oze3ZlImWyD0,1170
+mathai/apart.py,sha256=8IlJ8X6SKAjPunUPHVLmgNB5GuEi4_XpY0oLkQqrHKc,3770
+mathai/base.py,sha256=DBFHDBRSPeb8O2cH_YzTNG3h7vxrQzrEPSUpDIPCtPI,12724
 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=UKR_GH-ifplq_VyJKZukWjVDsQKkXXIuKNhpZ4w1Mhk,15033
+mathai/integrate.py,sha256=mKw_L8zccFifJINKLCFiYBrkgQM2qO_Pk8b2ShgxTa8,15065
 mathai/inverse.py,sha256=QCvDrzKquWsZv-BDAzZd9HnU0c3gZvcc44UztHVO5LQ,2919
 mathai/limit.py,sha256=RA8YAehgYCGVWv9qBc8uQ34BQ9mFthWl2OrVTwcHl2g,4920
-mathai/linear.py,sha256=wyiLpIxRDmD96xXktkgvc5gegIB3zblLB1EuV3TFdmU,5474
+mathai/linear.py,sha256=BzSnm941Zlod_l6hON4Rs6J4pdAA3MGpRVqr6-66ZBk,5524
 mathai/logic.py,sha256=UvHzRmKcO9AD51tRzHmpNSEhgW5gmaf4XPaQKFjGfC4,9653
 mathai/parser.py,sha256=f7bemieFmp0sbup1NlraMLvZDVFvqKGFknEVtlFRMVk,6979
 mathai/printeq.py,sha256=gIes-pstFOa6FcnpVIVvkjVKuWdsVdo11LlEnmHhakU,1303
 mathai/search.py,sha256=6S_V0PrC2ZnbhkWyT2ZTtExRD_mDlo-C3p0TDn6Dk3w,4743
-mathai/simplify.py,sha256=F37h-Z_rW35uVgx0G86vQErU2Ac6ZiTBN9KcC3RzDkg,15156
+mathai/simplify.py,sha256=qG23MRyVLjJL_8Xk5n4Z596bjihFvZg_zKtIvOqrLIk,15824
 mathai/structure.py,sha256=4Ww2IAx62RcQSO7_17TZES-DjMWBpcFQtL939FBIHwY,4103
 mathai/tool.py,sha256=UyccamiJy_CkFPakfufyPzdhtlEO6v2D7qwbXQ9V7Rg,2000
-mathai/trig.py,sha256=5P0RNS4eetNds2l-wyA4BAKqJdFIUws_8RGS_dH7gp8,9348
+mathai/trig.py,sha256=Zadw-uTJaMfoXwxuiio7QYHWgPxezEP2xZ2fr__PFM0,10290
 mathai/univariate_inequality.py,sha256=_r-kkiS4Hr-jRN7f-EL_E4svAMFWJP1Ea50HJKKbjfk,14778
-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,,
+mathai-0.4.1.dist-info/METADATA,sha256=m5N3sDp8R2_3ffySHyzvmo4YDpxYq6dT3FXUfkerCco,7021
+mathai-0.4.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+mathai-0.4.1.dist-info/top_level.txt,sha256=ROP4l3OhGYw3ihkQGASr18xM9GsK4z3_6whV5AyXLwE,7
+mathai-0.4.1.dist-info/RECORD,,