gh XdZddlmZmZmZddlmZddlmZm Z m Z ddl m Z m Z mZddlmZddlmZmZddlmZdd lmZmZmZdd lmZdd lmZdd lmZdd l m!Z!m"Z"m#Z#m$Z$m%Z%ddl&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0ddl1m2Z2ddl3m4Z4m5Z5ddl6m7Z7ddl8m9Z9m:Z:ddl;mm?Z?e d\ Z@ZAZBZCZDZEZFZGZHZIdZJdZKdZLdZMdZNdZOdZPdZQdZRd ZSd!ZTd"ZUd#ZVd$ZWd%ZXd&ZYd'ZZd(Z[d)Z\d*Z]d+Z^d,Z_d-Z`d.Zad/Zbd0Zce:d1Zdd2Zed3Zfd4Zgd5Zhd6Zid7ZjdS)8zBTests for algorithms for computing symbolic roots of polynomials. )IRationalpi)S)SymbolWildsymbols) conjugateimre)exp)rootsqrt) Piecewise)acoscossin)ZZInterval)powsimp)Polycyclotomic_poly intervalsnrootsrootof) root_factors roots_linearroots_quadratic roots_cubic roots_quartic roots_quinticroots_cyclotomicroots_binomialpreprocess_rootsroots) legendre_poly)PolynomialErrorUnsolvableFactorError)_nsort)raisesslow)verify_numericallyN)productza,b,c,d,e,q,t,x,y,zcRtd|D}t|d|t|d|ksJt|t |dD]<}||dz||ks%||dzt ||ksdS=dS)Nc3,K|]}|jrdndVdS)rNis_real.0is \/var/www/html/ai-engine/env/lib/python3.11/site-packages/sympy/polys/tests/test_polyroots.py z_check..(s-55!QY%A555555r1FT)sumsortedlistrangelenr )r&nrealixs r7_checkrB$s 55u555 5 5E %- Dvv$7$7 7 7 7 7E3u::q))b1f r**b1f 59!5!55555 4r9cttdtzdzttddgksJdS)Nr:r1)rrxrr9r7test_roots_linearrG2s> QqS1Wa(( ) )hr1oo-> > > > > > >r9c  ttdtdzztddgksJttdtdzzdtzzttdddgksJttdtdzzdztt t dzdz tt dzdz gksJttdtdzzdtzzdztdtt dzdz z dtt dzdz zgksJt tdtdzzdtzzdzttdzdtztzdtztzzttz z tzzttz ttdzzzttdzzz ttz z z}tt|tt ttzzttz z t ttzttzzttzz ttzz dtztztdzzzttz z z t ttzzttz z t ttzttzzttzz ttzz dtztztdzzzttz z zgksJtttdzzdtzz dtzz t}t|t dtdzzdz tz dtz zt dtdzzdztz dtz zgksJttdztdz dz tzztdzzdzt}t|dtdzdzgksJtt dtdzzdz t}t|}|t|ksJtd tdzzd tzz d z}d |d Dd|dDksJtdddD]R\}}}t|tdzz|tzz|z}t|}|t|ksJSdS)Nr:rrDr1iic8g|]}|dSr:nr5ws r7 z(test_roots_quadratic..Ps" 7 7 7qACCFF 7 7 7r9T)radicalsc8g|]}|dSrOrPrRs r7rTz(test_roots_quadratic..Qs" 8 8 8qACCFF 8 8 8r9F)r:)rrD)rrrErrrrB all_rootsaecdbyr*r.)fr_a_b_cr&s r7test_roots_quadraticrd6s 4!Q$?? + +1v 5 5 5 5 4!Q$1 a00 1 1hr1ooq5I I I I I 4!Q$ A.. / /QBtAwwJqL!DGG)A+3N N N N N 4!Q$1 q 0!44 5 5"qay{:JBQRSWXYSZSZQZ[\Q\L\9] ] ] ] ] 4!Q$1 q ! $ $ . . 0 0111 1!A!A A&q((AEAadFNQq!tV,Ca!e+LLA 41:: & & QUQU d1Q319qs?QqS#81Q3q5A:#EFFAN N QUQU d1Q319qs?QqS#81Q3q5A:#EFFAN N P P P P P Qq!tVac\AaC ##A 1   qAvz   1 qs "D1a4!$4$4Q$61$<= > > > > QTadUQYM !AqD (1 ,a00A 1   1a4!8      T!WWQT\A q!!AA q >>>> SAXA  #$$A 7 7AKKK66 7 7 7 8 8AKKK77 8 8 8 9 9 9 9gw88&& B AqD2a4"$ % %""u %%%%%&&r9c 8ttdztztdzztzt}t|t tdz t dtzdz zdt t d tzdz tdz zdt tdz t dtzdz z dt t d tzdz tdz z diksJdS)NrLr:r1)rrErr&reqs r7test_issue_7724riXs ad1fq!tma # #B 99 QqS47719Q;  d1ggXaZ\AaC   ! acDGGAIaK  ! tAwwhqjlQqS ! !!1 & & & & & & &r9c tdtddgtt ttt}t ddt dzdz z fd|D}t|tt tttksJtfd|DsJdS)Nr1rWrJr:cFg|]}|tSrF)subsr^)r5r`zs r7rTz#test_issue_8438..es% ( ( (QAFF1aLL ( ( (r9c3K|]@}tt|idddkVAdS)r:TchoprN)rmr^rErQ)r5r6prns r7r8z"test_issue_8438..isOGGQqvvq!Ql##%%ad%33q8GGGGGGr9) rr^rEas_exprr rrsetrmall)r&postrrrns @@r7test_issue_8438rwas aB^Q''))A Q A & &EQ!A#a%A ( ( ( (% ( ( (D t99 D1q))**++ , , , , GGGGG$GGG G GGGGGGr9ctdtdzzdz tttdzdzz}t|sJttdzdtdzzzdztfdt dD}ttdzdtdzzzdzt}||ksJt|sJtdtdzzdz }t|sJt tdtdzzdz dksJdS) NrLr1r:rfrKc0g|]}t|SrFr)r5r6r_s r7rTz#test_issue_8285..ps! ) ) )1&A,, ) ) )r9 )rrErXrBr>r?)r&ror_s @r7test_issue_8285r~lsY !AqD&1*a adQh / : : < .s! * * *A6!Q<< * * *r9)rrErXrrr>rB)ansrrs @r7test_issue_14291rs !a%!a1q51*q.11q59  )++!QUAE1tAwwqy=!d1ggai-H I I I I 1r!Q$wA * * * *q * * *C 77    # % % % % 3KKKKKr9cttdztttzztztd}t |}t |dksJdS)NrIEXdomain)rr^r rEr&r?)rhroots_ds r7test_issue_13340rsU adSVVAXo!1T 2 2 2BBiiG w<<1      r9cFttdztdzddtzzzztdzddtzzzztdd tzz zzd z d tzz tt}t fd |DsJdS) NrLrI r:iiiii ic36K|]}|dkVdSrNrFr5r`rhs r7r8z#test_issue_14522..s/,,arr!uuz,,,,,,r9)rrErr&ru)roots_eqrhs @r7test_issue_14522rs adQT219%%1dSUl(;;aQ>OORVVY]^_Y__ab c cBRyyH ,,,,8,,, , ,,,,,,r9ctttdzdtdzzz ttz zdz t}|dtsJdS)NrLrKr:rIr)r!rtrEhas)sols r7test_issue_15076rsW QTAadF]QqS014a88 9 9C q6::a==r9cttdzdtdztdzzz dtdzzzdtdztdzzz dtzzt}t|}d|vsJdS)NrLryr:rI@ir)rrErr&)rhrs r7test_issue_16589rsu adQtAwwYq!t^#a1f,r$q''z!Q$>aG K KBRyyH ======r9c ttdtdzztgdksJtttdzdtdzzz dtzzdz tgdksJttd}tttdztz t||t j tdtzdz zz|t j tdtzdz z zgksJtttdzdz tdt jttdzdz z t jttdzdz zgksJttdtdzzdtdzzz dtzz dz tdt jdtddzdz zdtddzdz zksJtdz dtdzzzdtzzdz }t|dd tt|td cxkr&tdddtd zttd td zd z dz zdz zdtd zttd td zd z  dz tdz zzdz tddzdtd zttd td zd z  dz tdz zzdz tddzgksnJdS)Nr:rI)rrrr1r1r1r1rDrT)trigmultiple)r ryrWrK)r rrErr^rHalfrrrr&rrrr)r`rhs r7test_roots_cubicrs` tAadFA ' '999 4 4 4 4 tAqD1QT6MAaC/!3Q77 8 8III E E E E Q A tAqD1Ha(( ) )a AF7T!WWQYq[ ! AF7T!WWQYq[ !.# # # # # tAqD2Iq)) * * QVaQik !16Ad1ggIaK#78 9 9 9 9 tAadFQq!tVOac1A5q99 : :1 = !Xa^^#A% %8Aq>>(9!(; ; < < < < Q$1a4!A#  !B $ . . . tB{{ . . .    A488CQtBxxZ^(<($?$??AA 488 CaRjn---a/"Q$67779HQNNJ 488 CaRjn---a/"Q$67779HQNNJ3         r9c tttdztgdksJtttdztdzztgdgdgdgdfvsJtttdztdzz tgdgd gd gd fvsJtttdztzt}tjt t dzd z ztjt t dzd z z tjtj g}t|t t|t ksJtgdD]\}\}}}}|d kr-| |d ztdz |td z z z}nH|dkrB|||d ztddz|tdz z z|tdz zz}tdz|tdzzz|td zzz|tzz|z tt t}t fd|DsJtttztdz ztdzztdzzd td zzztddz t t }t fd|DsJtdd tdzd tdzzzdtd zzzt dzzzdz tt t} t fd| DsJtt tdzztdzztz zt t }td|DsJtddtdddtddtdd dtddddf} | D]b tt  t}t fdt%||DsJcdS) NrLrrrrrI)rDrrr)rrDrr)rrrDr)rrrrD)r1rrr)rr1rr)rrr1r)rrrr1r:key))r1r:rIr)rIi r1r:rIrLr)rrJrI)rJrfr)rKirJryrc3~K|]7}t|ddkV8dS)TrprNrmrErQ)r5airhs r7r8z%test_roots_quartic..sEBB2771b>>###..!3BBBBBBr9r1c3jK|]-}tt|dV.dSrr-rmrE)r5r6rhs r7r8z%test_roots_quartic..s:AA!"''!Q--33AAAAAAr9rnT)negative rfc3tK|]2}tt|fdffdV3dS)rDrNr)r5r6rhrns r7r8z%test_roots_quartic..sGOOQ!"''Aq6Ar7*;"<".s-11tAww)#111111r9rD)r^rnrWc3jK|]-\}}t||z dV.dSr)r-rm)r5rSsreps r7r8z%test_roots_quartic..s@TTda%affSkkAoq99TTTTTTr9)r!rrErrrrZeroOner<hash enumeraterruqr r^rmzip)lhsrhsr6rYr]r[r\rrzansrepsrhrrns @@@r7test_roots_quarticrs adA ' '<<< 7 7 7 7 adQTk1-- . .    3     adQTk1-- . .    3     QTAXq)) * *C 6Ad1ggIaK !DGG)A+!5qvv FC #4 F3D$9$9$9 9 9 9 9%&8&8&899 C C.sDBB:a++BBBBBBr9rLr:rIrfrDr1)anyr!rrErrFr9r7test_issue_21287rsBB] QTAqD!ac'N "QqS"q&\ 1A 5! ;Q??6A6ABBBBBBBBBBr9cvtdzdz tdz dzdzdtdz dzzz dztdzdtdzzz dtdzzz dtzzd zf}|D]Ptt}t|dksJt fd |DsJQdS) Nrfr:r1 nrI7i ic3K|]J}t|dddkVKdS)r|gh㈵>rprNrrs r7r8z%test_roots_quintic..sQJJq2771acc"gg&&(((55:JJJJJJr9)rEr"rr?ru)eqsr&rhs @r7test_roots_quinticrs a4!8 qS1WqL1acAg; & + qD3q!t8ObAg %Q . 4 6CKKd2hh''5zzQJJJJEJJJJJJJJJKKr9cN ttdtddgksJttdtddgksJttdtdtddtt dzdz z tddtt dzdz zgksJttdtdt tgksJttdtdt jtt dzdz z t jtt dzdz zgksJttd tdttd z  tttd z zz ttd z  tttd z zztttdd z ttttdd zzz tttdd z ttttdd zzztttdd zttttdd zzz tttdd zttttdd zzzgksJttd tdt d dz tt dzdz z t d dz tt dzdz zt ddz tt dzdz z t ddz tt dzdz zgksJttd tdt d dz tdz z t d dz tdz zt ddz tdz z t ddz tdz zgksJttdtdd dgksJttdtdd dgksJttdtdd tdd dtddzgksJttdtdd t tgksJttd tdd tdd  tdd dz tdd dzdtdd dzz tdd ztdd dzzgksJttdtdd dtddz tddgksJdS)Nr1Tpolysr:rDrIrLrKrkryr)factorrf) r#rrErrrrrrrrrrFr9r7test_roots_cyclotomicrs* OAq=== > >1# E E E E OAq=== > >2$ F F F F O 1D  'A4771[\ ] ] ] ] OAq=== > >A2q' I I I I O 1D   v$q'' ! 3QVaQik5IJ K K K K OAq=== > > RT Qs2a4yy[ RT Qs2a4yy[ RA   !C8Aq>>(9$:$:":: RA   !C8Aq>>(9$:$:":: Bx1~~ 3r(1a..'8#9#9!99 Bx1~~ 3r(1a..'8#9#9!99 C     OAq=== > > a QtAwwYq[ a QtAwwYq[  Q Ad1ggIaK Q Ad1ggIaK C     OB>>> ? ? a QqS a QqS Q AaC Q AaC D     1D)))$ 8 8 8<=3 ? ? ? ? 1D)))$ 8 8 8<>4 @ @ @ @ OAq===d K K K r1++rDQKK'( ) ) ) ) OAq===d K K K Q     OAq===d K K K r1++R QR QT"a[[!^8KdSUWXkk8Y\`acef\g\gij\j8jk l l l l OAq===d K K K T"a[[$r1++& ' ' ' ' ' 'r9c ttdtztdgksJttdtdzztgdksJttdtzdzttddgksJdtddzdz }ttdtdzzdzt| |tzz | |tzz||tzz ||tzzgksJt tttd zdz t d d }t d d }tt|tdzz|zt}tt|tdzz|zt}t|dt|dksJt|dt|dksJtddddD]X\}}}}||kr|dkrt|t|zz||zz} t| } | t| ksJYttdtdzzdtdzzz tdtztddtdtzdz z zddtzddtztddtdtzdz zzdiksJdS)NrfrrLrr:rWr|rIrya1T) nonnegativeb1r1)r1r:)rDr1)r:rIrLrfrrD)r$rrErrrBrrrr.r*r&r^r) Arrr0r1rYr]rrQrrrs r7test_roots_binomialr7s $qsA,, ' 'A3 . . . . $qAvq// * *lll : : : : $qsQw** + +Q/@ @ @ @ @ HQNN2A $qAvz1-- . . acA2!8Q1Wa!A#g. / / / / >$q!tax.. ) )*** $ ' ' 'B $ ' ' 'B bAglA.. / /B R1Wr\1-- . .B 2a5>>WRU^^ + + + + 2a5>>WRU^^ + + + +ffg|DD"" 1a 66a1ff  1a4!A#  QfSkk!!!!! a1fr!Q$w&** + + !Xb!__tAwwqy{ *+Q !Q !Xb!__tAwwqy{ *+Q00 0 0 0 0 0 0r9c  ttztdzztztz }t t |t\}}|dksJ|t ttztdzztztz tksJt dztdzzt dztdzzzt tzztz}t t |t\}}|dt z ksJ|t tdztdzztztztksJt dztdzzt dztdzzztz}t t |t\}}|dt z ksJ|t tdztdzztztksJt dztdzzt tzztz}t t |t\}}|dt z ksJ|t tdztztztksJt dztdzztz}t t |t\}}|dt z ksJ|t tdztztksJtd\}}}}d|dzz|dzz|dzz d|ztz|d zz|d zz|d zz zd |d zz|dzz|d zztdzz|d zz z d|dzz|dzz|dzztdzz|dzz zd|dzz|dzz|dzztdzz|dzz z d|dzz|dzz|dzztdzz|d zz zd|dzz|d zz|dzztd zzd|dzzz zd|z|z|d zztd zzd|dzzz z d|dzztdzzdz z}t t |t\}}|d|z|z||dzzz ksJ|dtdzzdtd zzz dtd zzzdtdzzzdtdzzz dtdzzzdtdzzz dtzzd z ksJt tdz tdzttzztttttf!}t tdz ttztttt!}t |t|fksJdS)"Nr:r1rIE,F,J,L[r|)s2^ryr83ll8rk]$=,yrKr/@:rfr|D!^rLq Pyqidi Bi߰lul/1g&l glFW&kr) rYr^rEr]r%rr[r r r)r_coeffpolyEFJLgs r7test_roots_preprocessingrXsp !AqD1 qA"41::..KE4 A:::: 4!AqD1 q(!,, , , , , 1QT AqDAI!#a'A"41::..KE4 AaC<<<< 41q!t a!+Q// / / / / 1QT AqDAI!A"41::..KE4 AaC<<<< 41q!t a++ + + + + 1QT AaC!A"41::..KE4 AaC<<<< 41q1 a(( ( ( ( ( 1QT A A"41::..KE4 AaC<<<< 41q!$$ $ $ $ $##JAq!Qad"1a4'2-1Qq!t#AqD(B. /Aad"1a4'1,QU2 3 q!tAqD A%ad*1b50 1 AqDAad"1a4'1,  - A ad1a41$QT)  * ad 1a41QT!1QT6*  + QqA adBq!tG$ % AqDA f A#41::..KE4 BqDFAadFO # # # # 3q!t8fQTk)GAqDL89QT>IKXY[\X\L\\QT*1a4/02@2BCETU U U U U adUQT#a&&[ !Bq#a&&yM:::A adUSVV^Qr#a&&z222A A  1a& ( ( ( ( ( (r9c|5tdtiksJttttjdiksJttdzttjdiksJttdz tdzztdz zttd dtddtddiksJtdtzdztt dddiksJtdtzdzdztt dddiksJtdtzdzdztt dddiksJtdtzdzdztt dddiksJttdzdz tt dtjdtj dt diksJttdzdz dztt dtjdtj dt diksJtdtzdz dztt dddiksJtdtzdzdztt d ddiksJtdtzdz dztt dddiksJtdtzdzdztt d ddiksJtdtzdz dztt dddiksJtdtzdzdztt d ddiksJtttztz dztttz diksJtttztzdztt tz diksJttdzt dz tzztzttdtjdiksJttdzdtdzzz dzttjdtj diksJttd zdtdzzz dtdzzztdzz ttjddtdz dtjddtdzdiksJttd zdz ttddz t tdzdz zdtddz t tdzdz z dtd dz t tdzdz zdtd dz t tdzdz z dtjdtj dt dt diksJd tdzzd tdzzz dtdzzz dtdzzzdtd zzzdtdzzz dtd zzz dtdzzz }t|tjdtd dtddt dddt d ddt dddt dddiksJtttztztzttztztzz tttztztzttztzz diksJttdztdzztz dztdiksJttdz tdzztdz ztdtd dtddtddiksJttdz tdzztdz ztdz ztdtd dtddtddtddiksJttdzdtdzzzdtzzd zttd ddt zddt zdiksJttdzdtdzzzdtzzd ztddt zddt zdtd diksJttdztz tdzdtdzzzdtzzd zzttjdtjdtd ddt zddt zdiksJtjt dd}}dtdzdz|z}d|z dz }|dz }tdt zdz }||z }| |z }ttdztdzztz dztd| |z |z d| |z ||zz |z d| |z ||zz |z diksJtdzdtzzdztdtdzzdtzztdtzdz }ddD\} } td} t|t| | tdd t z| zz zzd| | tdd t z| zz zz d| | tdd t z| zzzzd| | tdd t z| zzzz diksJtdztdzztdzztzdz}ddD\} } }t|t| | d|zzzt || d|zzz|zzzd| | d|zzzt || d|zzz|zzz d| | d|zzz t || d|zzz |zzzd| | d|zzz t || d|zzz |zzz diksJt dzdt"z t dzzzddt"zzdtdzzz t zzt"z dtdzzz}t|t tjdtjtjt"zztjtddt"zz t"dzzd tdzzzzzdtjtjt"zztjtddt"zz t"dzzd tdzzzzz diksJtttztztdzzdtdzzzdtzzd ztdiksJtttztztdzzdtdzzzdtzzd ztdiksJttdzdz td !tjdtj diksJttdzdz td"!t dt diksJttdz tdzzttjdtj diksJttdz tdzztd#$tjdiksJttdzdz td d%tj tjgksJttdzdz td"d%t t gksJt%d&d'\}}tdz||z dzzdtz||z zzdz}t|td(!d||z z diksJttdztd)tjtjtjgksJtd*td)gksJtd ztdzz tdzztdzz tdzztz dz}t|t t't(dz zt+t(dz zdt t't(t ddzzt+t(t ddzz dt t't(t ddzzt+t(t ddzzdt t't(dz zt+t(dz zdt t't(t ddzzt+t(t ddzz dt t't(t ddzzt+t(t ddzzdiksJtdzdz|dzz}t|t t't(dz zt+t(dz zdt t't(t ddzzt+t(t ddzz dt t't(t ddzzt+t(t ddzzdt t't(dz zt+t(dz zdt t't(t ddzzt+t(t ddzz dt t't(t ddzzt+t(t ddzzdt dt diksJttdzd+tzzdz}d,|Dd-}d.d td/zz}|dt-|dzdz d0t-|dz zksJt/ddtdzztdzzd1dtdzz tdzzztd dz tzzdztd23}t|dtdzdddtdzzdtd dzdiksJt/d4tdzzd5tdztdzzzd6tdzzz d.tdztdzzz d7tdzzzd8tdztdzzzd9tdzzz d:tdztdzzz d;tdztzz d;tzz dtzzd?z d?tdzztd23}t|dtdzdzddtdzdzddtdzdz diksJtt/ttdzdzdz td23td t-ddtj tdt zdz z zzdtd t-ddtj tdt zdz zzzdtd t-ddzdiksJdS)@Nr1rr:rIrLrDrfr|rJrKryi iii irkirrF)cubicsrWT!c g|] }t| SrFrr5r`s r7rTztest_roots0..s,''' l'''r9))rr)r1rc g|] }t| SrFrrs r7rTztest_roots0..sIII!1IIIr9))r1rL)r1ry)rfryZfilterrc|jSr) is_positive)r`s r7ztest_roots0..s r9) predicate)rrza, b)realRri(c g|] }|j | SrFr2)r5rxs r7rTztest_roots0..*s...2:....r9rli2rrrr)}-rr)r&rErrexpandrrrrYr] NegativeOnerr[r\rrmrnr^r rrrrr)r_r1_2r1_3x0x1x2x3x4x5r13_20r1_20s2r1_4r1_8r5_8arbrrrrr` real_rootcrrhs r7 test_roots0r-si A;;"     A;;161+ % % % % Aq>>afa[ ( ( ( ( 1q51q5/1q5)1133Q 7 7QqTTE1addAqQRttUV;W W W W W 1q!  "a! 4 4 4 4 4 !A#'Aq ! !hr1ooq%9 9 9 9 9 !A#'Aq ! !hr1ooq%9 9 9 9 9 !A#'B " "xA&; ; ; ; ; A1  !Qq15&!aR!C C C C C !Q$(Q " "q!QUAvq1"a&H H H H H 1Q37Q,&&((! , ,1a!0D D D D D 1Q37Q,&&((! , ,"a!0D D D D D 1Q37Q,&&((! , ,1a!0D D D D D 1Q37Q,&&((! , ,"a!0D D D D D 1Q37Q,&&((! , ,1a!0D D D D D 1Q37Q,&&((! , ,"a!0D D D D D 1Q37Q,&&((! , ,!A#q 9 9 9 9 1Q37Q,&&((! , ,!Aq 9 9 9 9 A!a "Q& * *q!QUA.> > > > > A!Q$"A & &15!Q]A*F F F F F A!Q$1a4'!Q$. 2 2 2Q<AFArDGG|Q? @ @ @ @ A1   Q Ad1ggIaK Q Ad1ggIaK a QtAwwYq[ ! a QtAwwYq[ ! q15&!QA2q "     ad T!Q$Yad*T!Q$Y6adB AqD q!t8  Ag &A 88AaDD5!QqTT1hr1ooqQHROOQAK K K K K !a%!)Q!a%!)a-0! 4 4!a%!)a-!a%RS)9TVW8X X X X X A1q1$a 6 6 6" < < < < 1q5 AA!E 3 3 389!ua1q!A$$PQ6R S S S S 1q51q5/1q5)1q5199;;Qu M M M A$$1Q44AaDD!QqTT1- . . . . A!Q$1$q(! , ,!A$$2a4AaC0K K K K K A!Q$1$q(!D 9 9 9 Aq!A#q1Q44%# $ $ $ $ !Q$(QTAadF]QqS0145q : : 161qtteQ1a1a8 9 9 9 9!Q$D DHH*r/D B 2aB AB a1B dB tB A1q1$a 5 5 5 b4 BBq BBq:     A!aa1a4!A#..33AqsQw??A''$'''MFE aB A;;tA!BJ''''tA!BJ''''tA!BJ''''tA!BJ''''      1q!t adQ"AII/GIIID$ A;; QW q$ag"5!<<>(9$:$::A 3r(1a.. ! !!C8Aq>>(9$:$::A #bd)) c"Q$ii #b!Q  3r(1a..'8#9#991 #b!Q  3r(1a..'8#9#991      Q$(AqD  ""A 88 3r!t99 s2a4yy ! 3r(1a.. ! !!C8Aq>>(9$:$::A 3r(1a.. ! !!C8Aq>>(9$:$::A #bd)) c"Q$ii #b!Q  3r(1a..'8#9#991 #b!Q  3r(1a..'8#9#991 Aq!     adRTkBA..a...q1I qd| B 4A;;q(2d2qkk>9 9 9 9 9 q1T!WW9}ad"b1T!WW9nad%::tAwwhlA=MMPQQST]a b b bB 99d1ggq"qay.!d1ggX\1M M M M M bAg477 1a4'#ad(2Sa[A5EE1H 477{1a4 "$QT'*,.tAwwJq!tO<>@ajlKqD    B 99$q''Aq"qay.!R$q''\1aQiB477NA? ? ? ? ? adQq!tVmaQi1n,qay{:RTABF tAwwJ$ ( ( (B 99DGGaBtAwwJNAr$q''zA~qQ Q Q Q Q q477{Q&*Ad;;; < < q''DAJJ$q''!)A+ 56 6 q''DAJJ$q''!)A+ 56 6 q''DAJJ  # # # # # # #r9ctd\}}}}dz|z||z z|z|zz|z }dz||||z z|zzz|z|z|z||z|z z z||z||z z z}tt|ddgksJtt|ddgksJtd\}}}} } } } || z || z z|| z z| | z| zz| z| z | z } || z | z| z| || z z| zz| | z|| z zz}tt| |z | gdksJdzddzzzdztt}t dfd |DDrJd S) z6Just test that calculating these roots does not hang. z a,b,c,d,xr:r1zzz,yy,xx,zy,zx,yx,krrIryc3K|]}|VdSrrFr4s r7r8z"test_roots_slow..Ts"GG1GGGGGGr9cdg|],}|d-S)Trp)rmrQ)r5rir_rEs r7rTz#test_roots_slow..Ts3FFFbqvva}}D99FFFr9N)r r=r&valueskeysr)rYr]r[r\f1f2zzyyxxzyzxyxke1e2r r_rEs @@r7test_roots_slowr?@s%K((MAq!Q Aa1Q3!A#a% ! #B Aq1a!e9Q; !A#a%'1Q37"3 3qsQqSy AB b! ##%% & &1a& 0 0 0 0 b! ##%% & &1a& 0 0 0 0")*?"@"@RRRQ q&26 BF #beBh . 3b 82 =B q&"R"b1f+b. (2b5"q&> 9B b2gq!!((** + +yyy 8 8 8 8 1qAv A U1XX]]__AGGFFFFFAFFFGGGGGGGGGGr9cPttdztzdztd}ttdztzdztd}t||D]\}}t||z dksJtdzdt d ztd zzzd d t dzztdzzz d dt dzdzzz}t|d}d}t||D]\}}t||z dksJdS)Nr:r1Tr g?g-q=rLg@g@rIgS@g8@gb@g"@)g؋~)g)c g3! @g @绽|=)r&rErabsr)R1R2rr2r_s r7test_roots_inexactrFWs5 q!tax!|Q . . .B q!tax#~q4 0 0 0Bb"++$$B27||e##### 1s499}QT!!TDcN%:AqD$@@ 499s" # $A q4 B .Bb"++$$B27||e#####$$r9c6td\}}}}d|dzz|dzz|dzz d|ztz|dzz|dzz|dzz zd|d zz|d zz|d zztd zz|d zz z d |d zz|dzz|d zztdzz|dzz zd|dzz|dzz|dzztdzz|dzz z d|dzz|d zz|dzztd zz|d zz zd|d zz|d zz|d zztd zzd |dzzz zd|z|z|dzztdzzd|d zzz z d|dzztdzzdz z}t|tiksJt|td}gd}t d}||z|z||d zzz }t |t |ksJt ||D]:\} } | |} | t| || z dksJ;dS)NrrryrrrkrrrKr:rrrfrIr|rrLrrrrrrTr )g?|YcgOXLzGX@g'I^g@g |۰n@gdZw@g5`e@gS$Ȫ@gJm@rSrA) r rEr&evalfrr?rmatchrB) rrrrr_rCrDrSrrrrErIs r7test_roots_preprocessedrJis##JAq!Qad"1a4'2-1Qq!t#AqD(B. /Aad"1a4'1,QU2 3 q!tAqD A%ad*1b50 1 AqDAad"1a4'1,  - A ad1a41$QT)  * ad 1a41QT!1QT6*  + QqA adBq!tG$ % AqDA f A A;;"     qwwyy!d + + +B R R RB S A !AqAvA r77c"gg    b"++@@B  SqB%7%7%%?%?%?%?%?@@r9c~ttdzdtzz dzdddiksJttdzdtzz dzdddiksJttdzdtdzzz tdzz dtzzdz dddiksJttd dS) Nr:r1FstrictTrKrfrIcttdzdtdzzz tdzz dtzzdz dS)NrKr:rfrITrL)r&rErFr9r7rz#test_roots_strict..s;%1qAv 10Dqs0JQ0NW[*\*\*\r9)r&rEr+r)rFr9r7test_roots_strictrOs A!a . . .1a& 8 8 8 8 A!a - - -!Q 7 7 7 7 A!Q$A%!+a/ > > >1a& H H H H "\"\]]]]]r9cddtzz dtdzzz dtdzzzdtdzzz }t|d \}}t|}t|d }d |D}d |D}||z}d|D}t |d}t |d}||fD]H}t ||D]5\}}|jr||vsJt|t|f|vsJ6IdS)Nipiir:i rI1rLT)rur c:g|]\\}}}t||SrFrr5rYr]_s r7rTz$test_roots_mixed..s( 1 1 1yv1qHQNN 1 1 1r9c g|]]\\}}}tt|t|tt|t|z^SrF)rr r rSs r7rTz$test_roots_mixed..sd   ?fq!  RUUBqEE " "8BqEE2a55#9#9 9   r9c6g|]}|SrF)rHrs r7rTz$test_roots_mixed..s ,,,a ,,,r9c*|Srsort_keyrEs r7rz"test_roots_mixed..AJJLLr9rc*|SrrXrZs r7rz"test_roots_mixed..r[r9) rErrr&r<rr3r r ) r__re_im_nroots_sroots _intervals_rootsr6r`s r7test_roots_mixedrcsa Qad"T!Q$Y.AqD8A%%%HCQiiGA%%%G 1 1C 1 1 1C     CsJ,,7,,,GW"8"8999GW"8"8999GG$++ F++ + +DAqy +Avvvvv1r!uu~*****  +++r9c XttdttdtgksJtttttttgksJttdzdz ttdztdz gksJttdztz ttt tz tt tzgksJttdzdz dztdztdztdz tdz tt z tt z tt ztt zgksJtttdzdz tdttdztttdz tttdzdztgksJtdtdzzdtdzzzdtdzzztdzztd tttdzdtdzzzdtdzzzdzgksJdS) Nr1r:rLrrryrrKQ)rrrEr^rrrFr9r7test_root_factorsrfs Q # #Q | 3 3 3 3 Q # #Q | 3 3 3 3 1q! $ $QA 6 6 6 6 1q! $ $T!WWa$q''k(B B B B B AA & & QAq1ua!eQUAE1q5!a%@ A A A A QTAXq))# 6 6 6 a!eQa!eQadQh):):; < < < < !Q$AqD(1QT61AqD8!C H H H Aq!ta1f}r!Q$w&*+ , , , , , ,r9cd}t|tdttjjfdd}d|DgdksJdS) NrTrc2ddS)NrIrf)rQmaxsteps)r)rrsr7rztest_nroots1..sAHHq1H,E,Er9rIrPc,g|]}t|SrFstrrs r7rTz test_nroots1.. " " "qCFF " " "r9)@z-0.999z-0.996z-0.991z-0.983z-0.973z-0.961z-0.946z-0.930z-0.911z-0.889z-0.866z-0.841z-0.813z-0.784z-0.753z-0.720z-0.685z-0.649z-0.611z-0.572z-0.531z-0.489z-0.446z-0.402z-0.357z-0.311z-0.265z-0.217z-0.170z-0.121z-0.0730z-0.0243z0.0243z0.0730z0.121z0.170z0.217z0.265z0.311z0.357z0.402z0.446z0.489z0.531z0.572z0.611z0.649z0.685z0.720z0.753z0.784z0.813z0.841z0.866z0.889z0.911z0.930z0.946z0.961z0.973z0.983z0.991z0.996z0.999)r'rEr+mpmathmp NoConvergencer)rQr&rrs @r7 test_nroots1rqs Aa$'''A 69 "$E$E$E$EFFF HHqHMME # "E " " " K K K K K K K K Kr9cttdzdtzzdzt}|d}d|DgdksJ|d}d|DgdksJdS) NrfrIr1rPc,g|]}t|SrFrkrs r7rTz test_nroots2..rmr9)z-0.332z-0.839 - 0.944*Iz-0.839 + 0.944*Iz1.01 - 0.937*Iz1.01 + 0.937*Ic,g|]}t|SrFrkrs r7rTz test_nroots2..rmr9)z-0.33199z-0.83907 - 0.94385*Iz-0.83907 + 0.94385*Iz1.0051 - 0.93726*Iz1.0051 + 0.93726*I)rrEr)rrr&s r7 test_nroots2rus QTAaCZ!^QA HHqHMME # "E " " " 4 4 4 4 4 4 4 HHqHMME " "E " " " : : : : : : : : :r9c ttttdztdzt t zztzt ztddksJdS)NrIr:T) composite)r?r&rr^rrErFr9r7test_roots_compositerxsY uT!Q$Ad1gg-1A5qDIIIJJ K Kq P P P P P Pr9cttdzttz dzttfddS)NrIr1c"tSrrrgsr7rz"test_issue_19113..s E"IIr9)rrEr+r(rgs@r7test_issue_19113r{sB QSVV a B ?----.....r9c|tddddtzz dzzdz dtzzdgd ddgksJdS) Nr1rJrrLr:ryrrTr )r&rrFr9r7test_issue_17454r}sP !RacA %a'"Q$.2T B B Bq!f L L L L L Lr9c ttdztdgksJttdzdztdtddz z gksJdS)Nl l rIrLr:)rrE real_rootsrrFr9r7test_issue_20913rsz --q 1 1 < < > >C]B^ ^ ^ ^ ^ 1q!   ' ' ) )q1Q446{l^ ; ; ; ; ; ;r9c ttdd}dtz |ztdzd|zzz}ttttdzztdzdzzt|d| dt dt zzzdz d|dt dt zzzdz diksJdS)Nr1rIrDrfr:)rrYr&rrErr)rZr`s r7test_issue_22768rsAA A 1q5AaC. A a1fAz)1-- . . 1 AQ M1a 2Q >1a3! ! ! ! ! ! !r9)k__doc__sympy.core.numbersrrrsympy.core.singletonrsympy.core.symbolrrr $sympy.functions.elementary.complexesr r r &sympy.functions.elementary.exponentialr (sympy.functions.elementary.miscellaneousrr$sympy.functions.elementary.piecewiser(sympy.functions.elementary.trigonometricrrrsympy.polys.domains.integerringrsympy.sets.setsrsympy.simplify.powsimpr sympy.polysrrrrrsympy.polys.polyrootsrrrr r!r"r#r$r%r&sympy.polys.orthopolysr'sympy.polys.polyerrorsr(r)sympy.polys.polyutilsr*sympy.testing.pytestr+r,sympy.core.randomr-rn itertoolsr.rYr]r[r\rZrrrEr^rnrBrGrdrirwr~rrrrrrrrrrrrrr-r?rFrJrOrcrfrqrurxr{r}rrrFr9r7rsHH0000000000""""""5555555555DDDDDDDDDD666666AAAAAAAA::::::EEEEEEEEEE......$$$$$$******HHHHHHHHHHHHHH????????????????????????100000((((((--------000000  'w'<==1aAq!Q1   ???&&&D&&&HHH 4 4 4    ---      .7U7U7UtBBB KKK-'-'-'`000B9)9)9)xj#j#j#XHHH.$$$$@@@:^^^+++2 , , , KKK*:::"QQQ/// MMM<<< !!!!!r9