mathai 0.2.7__py3-none-any.whl → 0.2.9__py3-none-any.whl

mathai/univariate_inequality.py

@@ -1,414 +1,410 @@
-import itertools
-
-from .base import *
-from .inverse import inverse
-from collections import Counter
-from .factor import factor2
-from .simplify import simplify, solve
-from .expand import expand
-from .fraction import fraction
-import copy
-from .diff import diff
-from .logic import logic0
-def intersection2(domain, lst):
-    domain = copy.deepcopy(domain)
-    if domain == [True]:
-        return lst
-    elif domain == [True]:
-        return []
-    lst = [item for item in lst if item not in domain]
-    out = []
-    for item2 in lst:
-        for index in range(len(domain)):
-            
-            if isinstance(domain[index], bool) and domain[index]:
-                
-                if index == 0 and compute(item2) < compute(domain[index+1]):
-                    
-                    out.append(item2)
-                    break
-                elif index == len(domain)-1 and compute(domain[index-1]) < compute(item2):
-                    out.append(item2)
-                    break
-                elif index != 0 and index != len(domain)-1 and compute(domain[index-1]) < compute(item2) and compute(item2) < compute(domain[index+1]):
-                    
-                    out.append(item2)
-                    break
-                
-    return list(set(out))
-def flip_less_than(inter):
-    inter = copy.deepcopy(inter)
-    return [not item if isinstance(item, bool) else item for item in inter]
-def intersection(domain_1, domain_2):
-    domain_1, domain_2 = copy.deepcopy(domain_1), copy.deepcopy(domain_2)
-    if domain_1 == [True]:
-        return domain_2
-    if domain_2 == [True]:
-        return domain_1
-    if domain_1 == [False] or domain_2 == [False]:
-        return [False]
-    def simplify_ranges(ranges):
-        simplified_ranges = []
-        i = 0
-        while i < len(ranges):
-            if i + 2 < len(ranges) and ranges[i] is True and ranges[i + 2] is True:
-                simplified_ranges.append(True)
-                i += 3
-            elif i + 2 < len(ranges) and ranges[i] is False and ranges[i + 2] is False:
-                simplified_ranges.append(False)
-                i += 3
-            else:
-                simplified_ranges.append(ranges[i])
-                i += 1
-        return simplified_ranges
-    result = domain_1 + domain_2
-    result = [item for item in result if not isinstance(item, bool)]
-    result = list(set(result))
-    result = sorted(result, key=lambda x: compute(x))
-    i = len(result)
-    while i>=0:
-        result.insert(i, True)
-        i = i - 1
-    result[0] = domain_1[0] and domain_2[0]
-    result[-1] = domain_1[-1] and domain_2[-1]
-    def find_fraction_in_list(fraction_list, target_fraction):
-        for i in range(1, len(fraction_list)-1, 2):
-            if fraction_list[i] == target_fraction:
-                return i
-        return -1
-    for i in range(2, len(result)-1, 2):
-        if result[i+1] in domain_1:
-            result[i] = result[i] and domain_1[find_fraction_in_list(domain_1, result[i+1])-1]
-        if result[i+1] in domain_2:
-            result[i] = result[i] and domain_2[find_fraction_in_list(domain_2, result[i+1])-1]
-        if result[i-1] in domain_1:
-            result[i] = result[i] and domain_1[find_fraction_in_list(domain_1, result[i-1])+1]
-        if result[i-1] in domain_2:
-            result[i] = result[i] and domain_2[find_fraction_in_list(domain_2, result[i-1])+1]
-
-    result = simplify_ranges(result)
-    return result
-class Range:
-    def __init__(self, r=[True], p=[], z=[]):
-        self.r = r
-        self.p = p
-        self.z = z
-        self.do = True
-        
-    def unfix(self):
-        self.do = False
-        return self
-    def fix(self):
-        if not self.do:
-            return
-        def simplify_ranges(ranges):
-            simplified_ranges = []
-            i = 0
-            while i < len(ranges):
-                if i + 2 < len(ranges) and ranges[i] is True and ranges[i + 2] is True:
-                    simplified_ranges.append(True)
-                    i += 3
-                elif i + 2 < len(ranges) and ranges[i] is False and ranges[i + 2] is False:
-                    simplified_ranges.append(False)
-                    i += 3
-                else:
-                    simplified_ranges.append(ranges[i])
-                    i += 1
-            return simplified_ranges
-        
-        self.r = simplify_ranges(self.r)
-        
-        common = set(self.p) & set(self.z)
-        self.z = list(set(self.z) - common)
-        self.p = list(set(self.p) - common)
-        
-        self.p = list(set(self.p) - set(intersection2(self.r, self.p)))
-        self.z = list(set(intersection2(self.r, self.z)))
-        return self
-        
-    def __or__(self, other):
-        return (self.unfix().__invert__().unfix() & other.unfix().__invert__().unfix()).unfix().__invert__().fix()
-    def __invert__(self):
-        tmp =  Range(flip_less_than(self.r), self.z, list(set(self.p)-set(self.z)))
-
-        return tmp
-    def __and__(self, other):
-        a = intersection(self.r, other.r)
-        b = intersection2(self.r, other.p)
-        c = intersection2(other.r, self.p)
-        tmp = Range(a, list(set(b)|set(c)|(set(self.p)&set(other.p))), list(set(self.z)|set(other.z)))
-        return tmp
-    def __str__(self):
-
-        if self.r == [False] and self.p == [] and self.z == []:
-            return "{}"
-        out = []
-        out2 = ""
-        if self.r != [False]:
-            for i in range(0, len(self.r), 2):
-                string = ""
-                if self.r[i]:
-                    if i == 0:
-                        string += "(-inf,"
-                        if len(self.r)==1:
-                            string += "+inf)"
-                        else:
-                            string += str(self.r[i+1])+")"
-                    elif i == len(self.r)-1 and len(self.r)!=1:
-                        string += "("+str(self.r[i-1])+",+inf)"
-                    else:
-                        string += "("+str(self.r[i-1])+","+str(self.r[i+1])+")"
-                    out.append(string)
-        if self.p != []:
-            out.append("{"+",".join([str(item) for item in self.p])+"}")
-        if self.z != []:
-            out2 = "{"+",".join([str(item) for item in self.z])+"}"
-        if out2 == "":
-            return "U".join(out)
-        else:
-            return "U".join(out)+"-"+out2
-
-dic_table = {}
-def helper(eq, var="v_0"):
-    global dic_table
-    
-    eq2 = copy.deepcopy(eq)
-    
-    if eq2 in dic_table.keys():
-        return dic_table[eq2]
-    
-    if eq.children[0].name == "f_add":
-        
-        eq.children[0] = simplify(expand(eq.children[0]))
-        eq = simplify(factor2(eq))
-    
-        
-    equ = False
-    sign= True
-    if eq.name in ["f_gt", "f_ge"]:
-        sign = True
-    elif eq.name in ["f_lt", "f_le"]:
-        sign = False
-    if eq.name in ["f_ge", "f_le"]:
-        equ = True
-    if eq.name == "f_eq":
-        equ= True
-    critical = []
-    equal = []
-    more = []
-    
-    _, d = num_dem(eq.children[0])
-    d = factor2(d)
-    
-    for item in factor_generation(d):
-        
-        item = simplify(expand(item))
-        if len(vlist(item)) != 0:
-            v = vlist(item)[0]
-            if diff(diff(item, v), v) != tree_form("d_0"):
-                continue
-            out = inverse(item, vlist(item)[0])
-            more.append(simplify(out))
-    
-    eq.children[0] = factor2(eq.children[0])
-    
-    
-    for item in factor_generation(eq.children[0]):
-        item = simplify(expand(item))
-        
-        if len(vlist(item)) == 0:
-            if compute(item) <0:
-                sign = not sign
-            continue
-        v = vlist(item)[0]
-       
-        if item.name == "f_pow" and item.children[1].name== "d_-1":
-            
-            item = item.children[0]
-
-            tmp2 = diff(copy.deepcopy(item))
-            if "v_" in str_form(tmp2):
-                return None
-
-            
-            if diff(diff(item, v), v) != tree_form("d_0"):
-                
-                a = replace(diff(diff(item, v), v), tree_form(v), tree_form("d_0"))/tree_form("d_2")
-                if compute(a) < 0:
-                    sign = not sign
-                continue
-            
-            
-            if compute(tmp2)<0:
-                sign = not sign
-                item = simplify(item * tree_form("d_-1"))
-            out = inverse(item, vlist(item)[0])
-            critical.append(out)
-        else:
-            tmp2 = diff(copy.deepcopy(item))
-            if "v_" in str_form(tmp2):
-                return None
-            
-            if diff(diff(item, v), v) != tree_form("d_0"):
-                a = replace(diff(diff(item, v), v), tree_form(v), tree_form("d_0"))/tree_form("d_2")
-                if compute(a) < 0:
-                    sign = not sign
-                continue
-            
-            
-            if compute(tmp2)<0:
-                sign = not sign
-                item = simplify(item * tree_form("d_-1"))
-            out = inverse(item, vlist(item)[0])
-            critical.append(out)
-            if equ:
-                equal.append(out)
-    equal = list(set([simplify(item) for item in equal]))
-    more = list(set([simplify(item) for item in more]))
-    critical = [simplify(item) for item in critical]
-    critical = Counter(critical)
-    
-    critical = sorted(critical.items(), key=lambda x: compute(x[0]))
-
-    i = len(critical)
-    element = sign
-    while i>=0:
-        critical.insert(i, element)
-        if i>0 and critical[i-1][1] % 2 != 0:
-            element = not element
-        i = i - 1
-    for i in range(1, len(critical), 2):
-        critical[i] = critical[i][0]
-
-    
-    
-    if eq.name == "f_eq":
-        final = Range([False], equal, more)
-        dic_table[eq2] = final
-        return final
-  
-    final = Range(critical, equal, more)
-    dic_table[eq2] = final
-    return final
-def wavycurvy(eq):
-    if eq.name == "s_true":
-        return Range([True])
-    if eq.name == "s_false":
-        return Range([False])
-    if eq.name not in ["f_and", "f_or", "f_not"]:
-        
-        out = helper(eq)
-        if out is None:
-            return None
-        return out
-    lst= [wavycurvy(child) for child in eq.children]
-    if None in lst:
-        return None
-    ra = lst[0]
-    if eq.name == "f_and":
-        for child in lst[1:]:
-            ra = ra & child
-    elif eq.name == "f_or":
-        for child in lst[1:]:
-            ra = ra | child
-    elif eq.name == "f_not":
-        ra = ~ra
-    return ra
-
-def absolute(equation):
-    if equation.name in ["f_and", "f_or", "f_not"]:
-        tmp = TreeNode(equation.name, [absolute(child) for child in equation.children])
-        if len(tmp.children)==1:
-            tmp =tmp.children[0]
-        return tmp
-    def mul_abs(eq):
-        if eq.name == "f_abs" and eq.children[0].name == "f_mul":
-            return simplify(product([item.fx("abs") for item in factor_generation(eq.children[0])]))
-        return TreeNode(eq.name, [mul_abs(child) for child in eq.children])
-    equation = mul_abs(equation)
-            
-    def collectabs(eq):
-        out = []
-        if eq.name == "f_abs":
-            out.append(eq.children[0])
-            return out
-        for child in eq.children:
-            out += collectabs(child)
-        return out
-    def abc(eq, arr):
-        def trans(eq):
-            nonlocal arr
-            out = {}
-            if eq.name == "f_abs":
-                x = arr.pop(0)
-                if x == 0 or x==2:
-                    return eq.children[0]
-                else:
-                    return -eq.children[0]
-            else:  
-                return TreeNode(eq.name, [trans(child) for child in eq.children])
-        return trans(eq)
-    out = list(set(collectabs(equation)))
-    
-    out2 = []
-    for item in itertools.product([0,1,2], repeat=len(out)):
-        out3 = []
-        for i in range(len(item)):
-            out3.append(copy.deepcopy(TreeNode({0:"f_gt", 1:"f_lt", 2:"f_eq"}[item[i]], [out[i], tree_form("d_0")])))
-        out3 = TreeNode("f_and", out3+[abc(copy.deepcopy(equation), list(item))])
-        if len(out3.children) == 1:
-            out3 = out3.children[0]
-        out2.append(out3)
-    if len(out2) == 1:
-        return out2[0]
-    else:
-        equation = TreeNode("f_or", out2)
-    return equation
-def handle_sqrt(eq):
-    def helper2(eq):
-        if eq.name in ["f_lt", "f_gt", "f_le", "f_ge","f_eq"]:
-            out = []
-            def helper(eq):
-                nonlocal out
-                if eq.name == "f_pow" and frac(eq.children[1]) == Fraction(1,2):
-                    out.append(simplify(eq))
-                x = [helper(child) for child in eq.children]
-            helper(eq)
-            for item in out:
-                n = tree_form("d_1")
-                if eq.name == "f_eq":
-                    eq2 = inverse(simplify(eq.children[0]), str_form(item))
-                else:
-                    eq2, sgn = inverse(simplify(eq.children[0]), str_form(item), True)
-                    if sgn == False:
-                        n = tree_form("d_-1")
-                eq3 = simplify(expand(simplify(eq2**2)))
-                
-                return simplify(TreeNode(eq.name, [simplify(n*item.children[0]-eq3*n), tree_form("d_0")]))
-        return TreeNode(eq.name, [helper2(child) for child in eq.children])
-    return helper2(eq)
-def domain(eq):
-    eq = solve(eq, True)
-    out = []
-    def helper2(eq):
-        nonlocal out
-        if eq.name == "f_pow" and frac(eq.children[1]) is not None and frac(eq.children[1]).denominator == 2:
-            if "v_" in str_form(eq.children[0]):
-                out.append(TreeNode("f_ge", [eq.children[0], tree_form("d_0")]))
-                out.append(TreeNode("f_ge", [eq, tree_form("d_0")]))
-        if eq.name == "f_pow" and frac(eq.children[1]) is not None and frac(eq.children[1]) <0:
-            tmp = TreeNode("f_eq", [eq.children[0], tree_form("d_0")]).fx("not")
-            if "v_" in str_form(tmp):
-                out.append(tmp)
-        x = [helper2(child) for child in eq.children]
-    helper2(eq)
-    out = list(set([simplify(item) for item in out]))
-    if out == []:
-        return eq
-    if len(out)==1:
-        out = out[0]
-    else:
-        out = TreeNode("f_and", list(out))
-    if eq.name in ["f_lt", "f_gt", "f_le", "f_ge", "f_eq"]:
-        return eq & out
-    return out
+import itertools
+
+from .base import *
+from .inverse import inverse
+from collections import Counter
+from .factor import factor2
+from .simplify import simplify, solve
+from .expand import expand
+from .fraction import fraction
+import copy
+from .diff import diff
+from .logic import logic0
+def intersection2(domain, lst):
+    domain = copy.deepcopy(domain)
+    if domain == [True]:
+        return lst
+    elif domain == [True]:
+        return []
+    lst = [item for item in lst if item not in domain]
+    out = []
+    for item2 in lst:
+        for index in range(len(domain)):
+            
+            if isinstance(domain[index], bool) and domain[index]:
+                
+                if index == 0 and compute(item2) < compute(domain[index+1]):
+                    
+                    out.append(item2)
+                    break
+                elif index == len(domain)-1 and compute(domain[index-1]) < compute(item2):
+                    out.append(item2)
+                    break
+                elif index != 0 and index != len(domain)-1 and compute(domain[index-1]) < compute(item2) and compute(item2) < compute(domain[index+1]):
+                    
+                    out.append(item2)
+                    break
+                
+    return list(set(out))
+def flip_less_than(inter):
+    inter = copy.deepcopy(inter)
+    return [not item if isinstance(item, bool) else item for item in inter]
+def intersection(domain_1, domain_2):
+    domain_1, domain_2 = copy.deepcopy(domain_1), copy.deepcopy(domain_2)
+    if domain_1 == [True]:
+        return domain_2
+    if domain_2 == [True]:
+        return domain_1
+    if domain_1 == [False] or domain_2 == [False]:
+        return [False]
+    def simplify_ranges(ranges):
+        simplified_ranges = []
+        i = 0
+        while i < len(ranges):
+            if i + 2 < len(ranges) and ranges[i] is True and ranges[i + 2] is True:
+                simplified_ranges.append(True)
+                i += 3
+            elif i + 2 < len(ranges) and ranges[i] is False and ranges[i + 2] is False:
+                simplified_ranges.append(False)
+                i += 3
+            else:
+                simplified_ranges.append(ranges[i])
+                i += 1
+        return simplified_ranges
+    result = domain_1 + domain_2
+    result = [item for item in result if not isinstance(item, bool)]
+    result = list(set(result))
+    result = sorted(result, key=lambda x: compute(x))
+    i = len(result)
+    while i>=0:
+        result.insert(i, True)
+        i = i - 1
+    result[0] = domain_1[0] and domain_2[0]
+    result[-1] = domain_1[-1] and domain_2[-1]
+    def find_fraction_in_list(fraction_list, target_fraction):
+        for i in range(1, len(fraction_list)-1, 2):
+            if fraction_list[i] == target_fraction:
+                return i
+        return -1
+    for i in range(2, len(result)-1, 2):
+        if result[i+1] in domain_1:
+            result[i] = result[i] and domain_1[find_fraction_in_list(domain_1, result[i+1])-1]
+        if result[i+1] in domain_2:
+            result[i] = result[i] and domain_2[find_fraction_in_list(domain_2, result[i+1])-1]
+        if result[i-1] in domain_1:
+            result[i] = result[i] and domain_1[find_fraction_in_list(domain_1, result[i-1])+1]
+        if result[i-1] in domain_2:
+            result[i] = result[i] and domain_2[find_fraction_in_list(domain_2, result[i-1])+1]
+
+    result = simplify_ranges(result)
+    return result
+class Range:
+    def __init__(self, r=[True], p=[], z=[]):
+        self.r = r
+        self.p = p
+        self.z = z
+        self.do = True
+        
+    def unfix(self):
+        self.do = False
+        return self
+    def fix(self):
+        if not self.do:
+            return
+        def simplify_ranges(ranges):
+            simplified_ranges = []
+            i = 0
+            while i < len(ranges):
+                if i + 2 < len(ranges) and ranges[i] is True and ranges[i + 2] is True:
+                    simplified_ranges.append(True)
+                    i += 3
+                elif i + 2 < len(ranges) and ranges[i] is False and ranges[i + 2] is False:
+                    simplified_ranges.append(False)
+                    i += 3
+                else:
+                    simplified_ranges.append(ranges[i])
+                    i += 1
+            return simplified_ranges
+        
+        self.r = simplify_ranges(self.r)
+        
+        common = set(self.p) & set(self.z)
+        self.z = list(set(self.z) - common)
+        self.p = list(set(self.p) - common)
+        
+        self.p = list(set(self.p) - set(intersection2(self.r, self.p)))
+        self.z = list(set(intersection2(self.r, self.z)))
+        return self
+        
+    def __or__(self, other):
+        return (self.unfix().__invert__().unfix() & other.unfix().__invert__().unfix()).unfix().__invert__().fix()
+    def __invert__(self):
+        tmp =  Range(flip_less_than(self.r), self.z, list(set(self.p)-set(self.z)))
+
+        return tmp
+    def __and__(self, other):
+        a = intersection(self.r, other.r)
+        b = intersection2(self.r, other.p)
+        c = intersection2(other.r, self.p)
+        tmp = Range(a, list(set(b)|set(c)|(set(self.p)&set(other.p))), list(set(self.z)|set(other.z)))
+        return tmp
+    def __str__(self):
+
+        if self.r == [False] and self.p == [] and self.z == []:
+            return "{}"
+        out = []
+        out2 = ""
+        if self.r != [False]:
+            for i in range(0, len(self.r), 2):
+                string = ""
+                if self.r[i]:
+                    if i == 0:
+                        string += "(-inf,"
+                        if len(self.r)==1:
+                            string += "+inf)"
+                        else:
+                            string += str(self.r[i+1])+")"
+                    elif i == len(self.r)-1 and len(self.r)!=1:
+                        string += "("+str(self.r[i-1])+",+inf)"
+                    else:
+                        string += "("+str(self.r[i-1])+","+str(self.r[i+1])+")"
+                    out.append(string)
+        if self.p != []:
+            out.append("{"+",".join([str(item) for item in self.p])+"}")
+        if self.z != []:
+            out2 = "{"+",".join([str(item) for item in self.z])+"}"
+        if out2 == "":
+            return "U".join(out)
+        else:
+            return "U".join(out)+"-"+out2
+
+dic_table = {}
+def helper(eq, var="v_0"):
+    global dic_table
+    
+    eq2 = copy.deepcopy(eq)
+    
+    if eq2 in dic_table.keys():
+        return dic_table[eq2]
+    
+    if eq.children[0].name == "f_add":
+        
+        eq.children[0] = simplify(expand(eq.children[0]))
+        eq = simplify(factor2(eq))
+    
+        
+    equ = False
+    sign= True
+    if eq.name in ["f_gt", "f_ge"]:
+        sign = True
+    elif eq.name in ["f_lt", "f_le"]:
+        sign = False
+    if eq.name in ["f_ge", "f_le"]:
+        equ = True
+    if eq.name == "f_eq":
+        equ= True
+    critical = []
+    equal = []
+    more = []
+    
+    _, d = num_dem(eq.children[0])
+    d = factor2(d)
+    
+    for item in factor_generation(d):
+        
+        item = simplify(expand(item))
+        if len(vlist(item)) != 0:
+            v = vlist(item)[0]
+            if diff(diff(item, v), v) != tree_form("d_0"):
+                continue
+            out = inverse(item, vlist(item)[0])
+            more.append(simplify(out))
+    
+    eq.children[0] = factor2(eq.children[0])
+    
+    for item in factor_generation(eq.children[0]):
+        item = simplify(expand(item))
+        
+        if len(vlist(item)) == 0:
+            if compute(item) <0:
+                sign = not sign
+            continue
+        v = vlist(item)[0]
+        
+        if item.name == "f_pow" and item.children[1].name== "d_-1":
+            
+            item = item.children[0]
+            
+            if diff(diff(item, v), v) != tree_form("d_0"):
+                
+                a = replace(diff(diff(item, v), v), tree_form(v), tree_form("d_0"))/tree_form("d_2")
+                if "v_" in str_form(a):
+                    return None
+                if compute(a) < 0:
+                    sign = not sign
+                continue
+            else:
+                tmp2 = diff(copy.deepcopy(item))
+                if compute(tmp2)<0:
+                    sign = not sign
+                    item = simplify(item * tree_form("d_-1"))
+                out = inverse(item, vlist(item)[0])
+                critical.append(out)
+        else:
+            
+            
+            if diff(diff(item, v), v) != tree_form("d_0"):
+                a = replace(diff(diff(item, v), v), tree_form(v), tree_form("d_0"))/tree_form("d_2")
+                if "v_" in str_form(a):
+                    return None
+                if compute(a) < 0:
+                    sign = not sign
+                continue
+            else:
+                tmp2 = diff(copy.deepcopy(item))
+                if compute(tmp2)<0:
+                    sign = not sign
+                    item = simplify(item * tree_form("d_-1"))
+                out = inverse(item, vlist(item)[0])
+                critical.append(out)
+                if equ:
+                    equal.append(out)
+    equal = list(set([simplify(item) for item in equal]))
+    more = list(set([simplify(item) for item in more]))
+    critical = [simplify(item) for item in critical]
+    critical = Counter(critical)
+    
+    critical = sorted(critical.items(), key=lambda x: compute(x[0]))
+
+    i = len(critical)
+    element = sign
+    while i>=0:
+        critical.insert(i, element)
+        if i>0 and critical[i-1][1] % 2 != 0:
+            element = not element
+        i = i - 1
+    for i in range(1, len(critical), 2):
+        critical[i] = critical[i][0]
+
+    
+    
+    if eq.name == "f_eq":
+        final = Range([False], equal, more)
+        dic_table[eq2] = final
+        return final
+  
+    final = Range(critical, equal, more)
+    dic_table[eq2] = final
+    return final
+def wavycurvy(eq):
+    if eq.name == "s_true":
+        return Range([True])
+    if eq.name == "s_false":
+        return Range([False])
+    if eq.name not in ["f_and", "f_or", "f_not"]:
+        
+        out = helper(eq)
+        if out is None:
+            return None
+        return out
+    lst= [wavycurvy(child) for child in eq.children]
+    if None in lst:
+        return None
+    ra = lst[0]
+    if eq.name == "f_and":
+        for child in lst[1:]:
+            ra = ra & child
+    elif eq.name == "f_or":
+        for child in lst[1:]:
+            ra = ra | child
+    elif eq.name == "f_not":
+        ra = ~ra
+    return ra
+
+def absolute(equation):
+    if equation.name in ["f_and", "f_or", "f_not"]:
+        tmp = TreeNode(equation.name, [absolute(child) for child in equation.children])
+        if len(tmp.children)==1:
+            tmp =tmp.children[0]
+        return tmp
+    def mul_abs(eq):
+        if eq.name == "f_abs" and eq.children[0].name == "f_mul":
+            return simplify(product([item.fx("abs") for item in factor_generation(eq.children[0])]))
+        return TreeNode(eq.name, [mul_abs(child) for child in eq.children])
+    equation = mul_abs(equation)
+            
+    def collectabs(eq):
+        out = []
+        if eq.name == "f_abs":
+            out.append(eq.children[0])
+            return out
+        for child in eq.children:
+            out += collectabs(child)
+        return out
+    def abc(eq, arr):
+        def trans(eq):
+            nonlocal arr
+            out = {}
+            if eq.name == "f_abs":
+                x = arr.pop(0)
+                if x == 0 or x==2:
+                    return eq.children[0]
+                else:
+                    return -eq.children[0]
+            else:  
+                return TreeNode(eq.name, [trans(child) for child in eq.children])
+        return trans(eq)
+    out = list(set(collectabs(equation)))
+    
+    out2 = []
+    for item in itertools.product([0,1,2], repeat=len(out)):
+        out3 = []
+        for i in range(len(item)):
+            out3.append(copy.deepcopy(TreeNode({0:"f_gt", 1:"f_lt", 2:"f_eq"}[item[i]], [out[i], tree_form("d_0")])))
+        out3 = TreeNode("f_and", out3+[abc(copy.deepcopy(equation), list(item))])
+        if len(out3.children) == 1:
+            out3 = out3.children[0]
+        out2.append(out3)
+    if len(out2) == 1:
+        return out2[0]
+    else:
+        equation = TreeNode("f_or", out2)
+    return equation
+def handle_sqrt(eq):
+    def helper2(eq):
+        if eq.name in ["f_lt", "f_gt", "f_le", "f_ge","f_eq"]:
+            out = []
+            def helper(eq):
+                nonlocal out
+                if eq.name == "f_pow" and frac(eq.children[1]) == Fraction(1,2):
+                    out.append(simplify(eq))
+                x = [helper(child) for child in eq.children]
+            helper(eq)
+            for item in out:
+                n = tree_form("d_1")
+                if eq.name == "f_eq":
+                    eq2 = inverse(simplify(eq.children[0]), str_form(item))
+                else:
+                    eq2, sgn = inverse(simplify(eq.children[0]), str_form(item), True)
+                    if sgn == False:
+                        n = tree_form("d_-1")
+                eq3 = simplify(expand(simplify(eq2**2)))
+                
+                return simplify(TreeNode(eq.name, [simplify(n*item.children[0]-eq3*n), tree_form("d_0")]))
+        return TreeNode(eq.name, [helper2(child) for child in eq.children])
+    return helper2(eq)
+def domain(eq):
+    eq = solve(eq, True)
+    out = []
+    def helper2(eq):
+        nonlocal out
+        if eq.name == "f_pow" and frac(eq.children[1]) is not None and frac(eq.children[1]).denominator == 2:
+            if "v_" in str_form(eq.children[0]):
+                out.append(TreeNode("f_ge", [eq.children[0], tree_form("d_0")]))
+                out.append(TreeNode("f_ge", [eq, tree_form("d_0")]))
+        if eq.name == "f_pow" and frac(eq.children[1]) is not None and frac(eq.children[1]) <0:
+            tmp = TreeNode("f_eq", [eq.children[0], tree_form("d_0")]).fx("not")
+            if "v_" in str_form(tmp):
+                out.append(tmp)
+        x = [helper2(child) for child in eq.children]
+    helper2(eq)
+    out = list(set([simplify(item) for item in out]))
+    if out == []:
+        return eq
+    if len(out)==1:
+        out = out[0]
+    else:
+        out = TreeNode("f_and", list(out))
+    if eq.name in ["f_lt", "f_gt", "f_le", "f_ge", "f_eq"]:
+        return eq & out
+    return out