o o h@sbdZddlZddlmZmZmZddlmZddlm Z ddl m Z m Z m Z mZddlmZmZddlmZmZmZmZdd lmZdd lmZdd lmZdd lmZdd lm Z m!Z!m"Z"ddl#m$Z$ddl%m&Z&m'Z'ddl(m)Z)m*Z*ddl+m,Z,ddZ-ddZ.ddZ/d*ddZ0ddZ1ddZ2dej3dfd d!Z4d"d#Z5d+d$d%Z6d&d'Z7gfd(d)Z8dS),z>> from sympy import solve_poly_inequality, Poly >>> from sympy.abc import x >>> solve_poly_inequality(Poly(x, x, domain='ZZ'), '==') [{0}] >>> solve_poly_inequality(Poly(x**2 - 1, x, domain='ZZ'), '!=') [Interval.open(-oo, -1), Interval.open(-1, 1), Interval.open(1, oo)] >>> solve_poly_inequality(Poly(x**2 - 1, x, domain='ZZ'), '==') [{-1}, {1}] See Also ======== solve_poly_inequalities z8For efficiency reasons, `poly` should be a Poly instancer%could not determine truth value of %sF)multiple==!=T)NF><>=)r$T<=)r%Tz'%s' is not a valid relation) isinstancer ValueErroras_expr is_numberrrtrueRealsfalseEmptySetNotImplementedError real_rootsrappendNegativeInfinityInfinityLCreversedinsert)polyreltreals intervalsroot_intervalleftrightsigneq_signequal right_open multiplicityrJn/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/sympy/solvers/inequalities.pysolve_poly_inequalitysx       0 )         rLcCstdd|DS)aSolve polynomial inequalities with rational coefficients. Examples ======== >>> from sympy import Poly >>> from sympy.solvers.inequalities import solve_poly_inequalities >>> from sympy.abc import x >>> solve_poly_inequalities((( ... Poly(x**2 - 3), ">"), ( ... Poly(-x**2 + 1), ">"))) Union(Interval.open(-oo, -sqrt(3)), Interval.open(-1, 1), Interval.open(sqrt(3), oo)) cSsg|] }t|D]}|qqSrJ)rL).0psrJrJrK sz+solve_poly_inequalities..)r)polysrJrJrKsolve_poly_inequalitiesqsrRcCstj}|D]j}|s qttjtjg}|D]P\\}}}t|||}t|d}g} t||D]\} } | | } | tjurA| | q.| }g} |D]} |D]} | | 8} qL| tjur]| | qH| }|sdnq|D]} | | }qgq|S)a3Solve a system of rational inequalities with rational coefficients. Examples ======== >>> from sympy.abc import x >>> from sympy import solve_rational_inequalities, Poly >>> solve_rational_inequalities([[ ... ((Poly(-x + 1), Poly(1, x)), '>='), ... ((Poly(-x + 1), Poly(1, x)), '<=')]]) {1} >>> solve_rational_inequalities([[ ... ((Poly(x), Poly(1, x)), '!='), ... ((Poly(-x + 1), Poly(1, x)), '>=')]]) Union(Interval.open(-oo, 0), Interval.Lopen(0, 1)) See Also ======== solve_poly_inequality r") rr2rr6r7rL itertoolsproduct intersectr5union)eqsresult_eqsglobal_intervalsnumerdenomr<numer_intervalsdenom_intervalsr?numer_intervalglobal_intervalrBdenom_intervalrJrJrKsolve_rational_inequalitiess@        rbTc sd}g}tj}|D]}|sq g}tj}|D]} t| tr!| \} } n| jr/| j| j| j} } n| d} } | tj urDtj tj d} } } n| tj urTtj tj d} } } n| \} } z t| | f\\} } } Wn tyvttdw| jjs| | d} } }| j}|js|js| | } t| d| } |t| ddM}q|| | f| fq|r|t|gM}tfdd|Dg}||8}||O}q |s|r|}|r|}|S) a8Reduce a system of rational inequalities with rational coefficients. Examples ======== >>> from sympy import Symbol >>> from sympy.solvers.inequalities import reduce_rational_inequalities >>> x = Symbol('x', real=True) >>> reduce_rational_inequalities([[x**2 <= 0]], x) Eq(x, 0) >>> reduce_rational_inequalities([[x + 2 > 0]], x) -2 < x >>> reduce_rational_inequalities([[(x + 2, ">")]], x) -2 < x >>> reduce_rational_inequalities([[x + 2]], x) Eq(x, -2) This function find the non-infinite solution set so if the unknown symbol is declared as extended real rather than real then the result may include finiteness conditions: >>> y = Symbol('y', extended_real=True) >>> reduce_rational_inequalities([[y + 2 > 0]], y) (-2 < y) & (y < oo) Tr"z only polynomials and rational functions are supported in this context. Fr) relationalcs6g|]}|D]\\}}}|r||jfdfqqS)r")hasone)rMindrAgenrJrKrPs z0reduce_rational_inequalities..)rr2r0r+tuple is_Relationallhsrhsrel_opr/ZeroOner1togetheras_numer_denomrrrdomainis_Exactto_exact get_exactis_ZZis_QQrsolve_univariate_inequalityr5rbevalf as_relational)exprsrjrcexactrWsolution_exprsrY_solexprr<r[r\optrtexcluderJrirKreduce_rational_inequalitiess\              rcs|jdur ttdfddddd}g}|D]"\}}||vr-t|d|}n t| d||}||g|qt||S) aReduce an inequality with nested absolute values. Examples ======== >>> from sympy import reduce_abs_inequality, Abs, Symbol >>> x = Symbol('x', real=True) >>> reduce_abs_inequality(Abs(x - 5) - 3, '<', x) (2 < x) & (x < 8) >>> reduce_abs_inequality(Abs(x + 2)*3 - 13, '<', x) (-19/3 < x) & (x < 7/3) See Also ======== reduce_abs_inequalities Fzs Cannot solve inequalities with absolute values containing non-real variables. csg}|js|jr)|j|jD]}|}|s|}qfddt||D}q|S|jrG|jjs6t d| fdd|j D|St |t rw|jd}|D]\}}|||t|dgf|| |t|dgfqU|S|gfg}|S)Ncs*g|]\\}}\}}||||fqSrJrJ)rMrconds_expr_conds)oprJrKrPEs*zBreduce_abs_inequality.._bottom_up_scan..z'Only Integer Powers are allowed on Abs.c3s |] \}}||fVqdSNrJ)rMrr)rgrJrK LszAreduce_abs_inequality.._bottom_up_scan..r)is_Addis_MulfuncargsrSrTis_Powexp is_Integerr,extendbaser+rr5r r )rr}argrr_bottom_up_scan)rgrrKr9s2       z.reduce_abs_inequality.._bottom_up_scanr&r(r'r)r)is_extended_real TypeErrorrkeysrr5r)rr<rjmapping inequalitiesrrJrrKreduce_abs_inequalitys      rcstfdd|DS)aReduce a system of inequalities with nested absolute values. Examples ======== >>> from sympy import reduce_abs_inequalities, Abs, Symbol >>> x = Symbol('x', extended_real=True) >>> reduce_abs_inequalities([(Abs(3*x - 5) - 7, '<'), ... (Abs(x + 25) - 13, '>')], x) (-2/3 < x) & (x < 4) & (((-oo < x) & (x < -38)) | ((-12 < x) & (x < oo))) >>> reduce_abs_inequalities([(Abs(x - 4) + Abs(3*x - 5) - 7, '<')], x) (1/2 < x) & (x < 4) See Also ======== reduce_abs_inequality csg|] \}}t||qSrJ)r)rMrr<rirJrKrP{sz+reduce_abs_inequalities..r)r}rjrJrirKreduce_abs_inequalitiesfs rFc) sddlm}|tjdurttd|tjur-td|d|}|r+| }|S }|}j durCtj }|s>|S| |Sj duret ddd z  |iWn tydttd wd}tjurp|}ntjurztj }njj} t| } | tjkrt| } | d} | tjur|}na| tjurtj }nX| durt| |} j} | d vrˈ| jdr|}n#| jdstj }n| d vr| jdr|}n | jdstj }|j|j}}||tjurtd| dd|}|}|dur]| \}}z|jvrt | jd krt!t"| |}|dur&t!Wnt!tfy?ttd#t$dwt| fdd}g}|D] }|%t"||qS|sit&|}djvosjdk}zet'|j(t)|j|j}t)||t*|t|j|j|j|v|j|v}t+dd|Drt,|ddd}n*t-|dd}|drtz|d}t |d krt.|}Wn tytwWn tytdwtj}/tj0}tjkrd}t)}zt1||}t2|ts"|D]}||vr||r|j r|t)|7}q no|j|j} }!t,|t)|!D]S}|| }"| |!kr||}#t3| |}$|$|vr|$j r||$r|"ra|#ra|t| |7}n |"rm|t4| |7}n|#ry|t5| |7}n|t6| |7}|} q1|D] }%|t)|%8}qWntytj}d}Ynw|tj urt!td#||f||}tj g}&|j} | |vr|| r| j7r|&8t)| |D]?}'|'}!|t3| |!r|&8t| |!dd|'|vr|9|'n|'|vr |9|'||'}(n|}(|(r|&8t)|'|!} q|j}!|!|vr1||!r1|!j7r1|&8t)|!|t3| |!rB|&8t6| |!|tjkrQ|rQ||}n t:t;|&||#|}|sb|S| |S)aTSolves a real univariate inequality. Parameters ========== expr : Relational The target inequality gen : Symbol The variable for which the inequality is solved relational : bool A Relational type output is expected or not domain : Set The domain over which the equation is solved continuous: bool True if expr is known to be continuous over the given domain (and so continuous_domain() does not need to be called on it) Raises ====== NotImplementedError The solution of the inequality cannot be determined due to limitation in :func:`sympy.solvers.solveset.solvify`. Notes ===== Currently, we cannot solve all the inequalities due to limitations in :func:`sympy.solvers.solveset.solvify`. Also, the solution returned for trigonometric inequalities are restricted in its periodic interval. See Also ======== sympy.solvers.solveset.solvify: solver returning solveset solutions with solve's output API Examples ======== >>> from sympy import solve_univariate_inequality, Symbol, sin, Interval, S >>> x = Symbol('x') >>> solve_univariate_inequality(x**2 >= 4, x) ((2 <= x) & (x < oo)) | ((-oo < x) & (x <= -2)) >>> solve_univariate_inequality(x**2 >= 4, x, relational=False) Union(Interval(-oo, -2), Interval(2, oo)) >>> domain = Interval(0, S.Infinity) >>> solve_univariate_inequality(x**2 >= 4, x, False, domain) Interval(2, oo) >>> solve_univariate_inequality(sin(x) > 0, x, relational=False) Interval.open(0, pi) rdenomsFz| Inequalities in the complex domain are not supported. Try the real domain by setting domain=S.Reals)rc continuousNrjT extended_realz When gen is real, the relational has a complex part which leads to an invalid comparison like I < 0. r)r&r(r$z The inequality, %s, cannot be solved using solve_univariate_inequality. xcst|}z|d}Wn tytj}Ynw|tjtjfvr'|S|jdur/tjS|d}|j r=|dSt d|)NrFr*z!relationship did not evaluate: %s) subsrrrrr1r/rrg is_comparabler3)rvr expanded_errjrJrKvalids       z*solve_univariate_inequality..valid=r#css|]}|jVqdSr)r.)rMrrJrJrKr@sz.solve_univariate_inequality..) separatedcSs|jSrr)rrJrJrKCsz-solve_univariate_inequality..z'sorting of these roots is not supportedz %s contains imaginary parts which cannot be made 0 for any value of %s satisfying the inequality, leading to relations like I < 0. )<sympy.solvers.solversr is_subsetrr0r3rrz intersectionr|rr2r xreplacerr/r1rmrnrrprrrrosupinfr7rrUrs free_symbolslenr,rrr rrsetboundaryrlistallrrsortedcoeff ImaginaryUnitrr+_ptRopenLopenopen is_finiter5removerr))rrjrcrtrrrv_gen_domaineperiodconstfranger<rrrgrhsolnsr singularities include_xdiscontinuitiescritical_pointsr>sifted make_realcoeffIcheckim_solazstartend valid_startvalid_zptrOsol_setsr_validrJrrKrzsf 9                                          rzcCs|js|js||d}|S|jr|jrtj}|S|jr!|jdus)|jr-|jdur-td|jr3|js9|jr>|jr>||}}|jrZ|jrJ|d}|S|jrT|tj}|S|d}|S|jrt|jrg|tj}|S|jrp|d}|S|d}|S)z$Return a point between start and endr*Nz,cannot proceed with unsigned infinite valuesr$) is_infiniterrpis_extended_positiver,is_extended_negativeHalf)rrrrJrJrKrsF          rc Cshddlm}||jvr |S|j|kr|j}|j|kr"||jjvr"|Sdd}d}tj}|j|j}z t||}| dkrF| | d}n |sP| dkrPt Wnt t fy|sz t|gg|}Wnt yst||}Ynw||||} | tjur||||tjur|||kd}|||| } | tjur|||| tjur|| |kd}||| kd}|tjur| tjur||kn||k}| tjurt| |k|}nt|}Ynwg} |dur| } d} | j|dd\}}| |8} | |8} t| }|j|d d\}} |jd ks&|j|jkrdur+nn |jd vr+|} tj}| |} |jr:| | | }n|j | | }||j||jB}||}||D]*}tt|d||d }t|tr||j|kr|||||jtjur|| |qS| |fD]'}||||tjur||||tjur| ||ur||kn||kq| |t| S) aReturn the inequality with s isolated on the left, if possible. If the relationship is non-linear, a solution involving And or Or may be returned. False or True are returned if the relationship is never True or always True, respectively. If `linear` is True (default is False) an `s`-dependent expression will be isolated on the left, if possible but it will not be solved for `s` unless the expression is linear in `s`. Furthermore, only "safe" operations which do not change the sense of the relationship are applied: no division by an unsigned value is attempted unless the relationship involves Eq or Ne and no division by a value not known to be nonzero is ever attempted. Examples ======== >>> from sympy import Eq, Symbol >>> from sympy.solvers.inequalities import _solve_inequality as f >>> from sympy.abc import x, y For linear expressions, the symbol can be isolated: >>> f(x - 2 < 0, x) x < 2 >>> f(-x - 6 < x, x) x > -3 Sometimes nonlinear relationships will be False >>> f(x**2 + 4 < 0, x) False Or they may involve more than one region of values: >>> f(x**2 - 4 < 0, x) (-2 < x) & (x < 2) To restrict the solution to a relational, set linear=True and only the x-dependent portion will be isolated on the left: >>> f(x**2 - 4 < 0, x, linear=True) x**2 < 4 Division of only nonzero quantities is allowed, so x cannot be isolated by dividing by y: >>> y.is_nonzero is None # it is unknown whether it is 0 or not True >>> f(x*y < 1, x) x*y < 1 And while an equality (or inequality) still holds after dividing by a non-zero quantity >>> nz = Symbol('nz', nonzero=True) >>> f(Eq(x*nz, 1), x) Eq(x, 1/nz) the sign must be known for other inequalities involving > or <: >>> f(x*nz <= 1, x) nz*x <= 1 >>> p = Symbol('p', positive=True) >>> f(x*p <= 1, x) x <= 1/p When there are denominators in the original expression that are removed by expansion, conditions for them will be returned as part of the result: >>> f(x < x*(2/x - 1), x) (x < 1) & Ne(x, 0) rrcSsJz|||}|tjur|WS|dvrWdS|WSty$tjYSw)NTF)rrNaNr)ierOrfrrJrJrKclassifys    z#_solve_inequality..classifyNr$T)as_AddF)r#r")linear)rrrrnr9rmrr7rdegreerr-r3rrrzr/r1rras_independentris_zero is_negative is_positiverorq_solve_inequalityr r+r5)rrOrrrroorrNokoooknoorrrnbaxefrbeginning_denomscurrent_denomsrhcrfrJrJrKrs J                  rc s8ii}}g}|D]x}|j|j}}|t}t|dkr"|n"|j|@} t| dkr>| |tt |d|q t t d| rU| g||fq |fdd} | rutdd| Dru| g||fq |tt |d|q dd |D} d d |D} t| | |S) Nr$rzZ inequality has more than one symbol of interest. cs |o|jp|jo|jj Sr)rd is_Functionrrr)urirJrKrs z&_reduce_inequalities..css|]}t|tVqdSr)r+rrMrfrJrJrKrz'_reduce_inequalities..cSsg|] \}}t|g|qSrJ)rrMrjr}rJrJrKrPsz(_reduce_inequalities..cSsg|] \}}t||qSrJ)rrrJrJrKrPs)rmroatomsr rpoprr5rrr3r is_polynomial setdefaultfindritemsr) rsymbols poly_partabs_partother inequalityrr<genscommon components poly_reduced abs_reducedrJrirK_reduce_inequalitiests,        rcsJt|s|g}dd|D}tjdd|D}t|s |g}t|p%||@}tdd|Dr7ttddd|Dfd d|D}fd d |D}g}|D]<}t|trj||j |j d }n |d vrst |d }|dkrxqT|dkrt jS|j jrtd|||qT|}~t||}|ddDS)aEReduce a system of inequalities with rational coefficients. Examples ======== >>> from sympy.abc import x, y >>> from sympy import reduce_inequalities >>> reduce_inequalities(0 <= x + 3, []) (-3 <= x) & (x < oo) >>> reduce_inequalities(0 <= x + y*2 - 1, [x]) (x < oo) & (x >= 1 - 2*y) cSsg|]}t|qSrJrrrJrJrKrPsz'reduce_inequalities..cSsg|]}|jqSrJ)rrrJrJrKrPscss|]}|jduVqdS)FNrrrJrJrKrrz&reduce_inequalities..zP inequalities cannot contain symbols that are not real. cSs&i|]}|jdur|t|jddqS)NTr)rr namerrJrJrK s z'reduce_inequalities..csg|]}|qSrJrrrecastrJrKrPcsh|]}|qSrJrrrrJrK rz&reduce_inequalities..rrTFr cSsi|]\}}||qSrJrJ)rMkrrJrJrKrr)rrrVanyrrr+rrrmr-rnr rr1r.r3r5rrr)rrr keeprfrrJrrKreduce_inequalitiessB      r)T)F)9__doc__rSsympy.calculus.utilrrr sympy.corersympy.core.exprtoolsrsympy.core.relationalrr r r sympy.core.symbolr r sympy.sets.setsrrrrsympy.core.singletonrsympy.core.functionr$sympy.functions.elementary.complexesr sympy.logicr sympy.polysrrrsympy.polys.polyutilsrsympy.solvers.solvesetrrsympy.utilities.iterablesrrsympy.utilities.miscrrLrRrbrrrr0rzrrrrrJrJrJrKs>        [ B[G) !.-