g< zddlZddlTddlmZmZmZmZmZmZm Z m Z m Z m Z m Z dZdZdZdZdZd Zd Zd ZdS) N)*) SecantMuller BisectionIllinoisPegasusAndersonRidderANewtonNewtonMNewtonMDNewtonc dt_tdtddsJtt tdt sJtt tdtdft sJtdtdd sJd }tttttfD]7}t|d | }t||tksJ8tttt t"t$t&fD]7}t|d | }t||tksJ8d}t)t|ddtsJt)t*|ddt*jsJt)t,|ddt.sJt)t,|ddt0sJt3jt65t9d5tdtdddddn #1swxYwYdddn #1swxYwYt3jt65tdddddn #1swxYwYt3jt65tddddddS#1swxYwYdS)Ncd|zdz S)Nxs Y/var/www/html/ai-engine/env/lib/python3.11/site-packages/mpmath/tests/test_rootfinding.pyztest_findroot.. ac!eg?rgQ @c||zdzSNrrs rrztest_findroot.. rry@@y?c t|SN)cosrs rrztest_findroot..s #a&&r@solver)?r"c|dz dzSNrrs rrztest_findroot..s1q51*rrg|=)tolr%y?r(c|dzd|zz dzS)Nr(iClF[rrs rrztest_findroot.."sq!tgai/.@rz 5.278e+13z -5.278e+13)realimagcdS)Nrrs rrztest_findroot..'s2rrcd|z dzdz S)Nrr)ps rrztest_findroot..+sAEB;?rg?)mpdpsfindrootmpfaesinpimpcr rr rr absepsrrrr r isinstanceivfpfloatcomplexpytestraises ValueErrorworkprec)fr$rs r test_findrootrFs BF OOSVV , , / / 5 5555 CQ # #B ' '''' C#a&&#d)), - - 0 0 4 4444 OOSYY / / 2 22 6 6666A67FG< Q6 * * *11Q44yy369h Q 0 0 011Q44yy3A hq!/// 5 5555 bkk!RUk33RV < <<<< bkk!QEk22E : :::: bkk!Tuk55w ? ???? z " "?? a[[ ? ? @@k === ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ???????????????? z " """q!!!""""""""""""""" z " "11**C000111111111111111111sZK+'!K K+K K+K K++K/2K/L--L14L1M00M47M4c8tddddksJdS)Nc|dzdz S)Nr(rrrs rrz test_bisection../sad1fr)rr(bisectr#rr4rrrtest_bisectionrK-s, $$U( ; ; ;q @ @ @ @ @ @rcvd}t|dd}t||tksJdS)Nc&tgd|S)N)rrrr)polyvalrs rrztest_mnewton..2s')))A&&rgmnewtonr#)r4r:r;rErs r test_mnewtonrQ1s=&&AD+++A qqtt99s??????rcd}t|dt}t||tksJdS)Nc|dz dzS)Nr(drrs rrztest_anewton..7s1q53,rr%r#)r4r r:r;rPs r test_anewtonrU6s=ABw'''A qqtt99s??????rcd}t|dt}t||tksJdS)Ncd|zdzdzS)Nr(rrrs rrztest_muller..<s1q51*q.rr%r#)r4rr:r;rPs r test_mullerrX;s=  ABv&&&A qqtt99s??????rctddD]tfddksJtdddksJdS)Nrrc|dz zSrr)ris rrz#test_multiplicity..Bsq1uqjrc |dzSr'rrs rrz#test_multiplicity..Cs !Q$rr)range multiplicity)r[s@rtest_multiplicityr_@se 1a[[::0000!4499999  * *a / / / / / /rc`d}tt|dtddgddggz ddksJtt|ddd D]\}}t ||d d ksJd ddfd}t |d}d|DddgksJdS)Ncd|ddzzd|ddzzz dz |ddzd|dzz |ddzzd|dzzdz gS)Nrrr(rrrs rrEz test_multidimensional..fFsb!A$' !AaD!G)#A%qtQwq1v~adAg'=a!f'DQ'FGGr)rrbrrcrgHz>c,t|tSr )norminfrs rrz'test_multidimensional..JsDCLLr)verboserfr(g+=c| d|zzSr'rrys rrz'test_multidimensional..Osqb1Q3hrc<|dz||dzdz zzd|zz |dzz S)Nr(rrrjs rrz'test_multidimensional..Ps-q!taAl*QqS0q1u=rc2t|dz|dzzSr')sqrtrjs rrz'test_multidimensional..Qsd1a4!Q$;''rcX||}|||z |||z fSr r)rrkf1xf1f2f3s rrEz test_multidimensional..fRs9bAhh1a31a3//r) rtcFg|]}tt|Sr)intround).0r[s r z)test_multidimensional..Vs$ % % %aCaMM % % %rrr)mnormjacobianmatrixrr2rfr4)rErerrorrqrrrss @@@rtest_multidimensionalr~EsHHH !V$$v!uaVn'='==a @ @5 H H H HRFA"8"8:::  5 1q>>E ! ! ! !  B = =B ' 'B0000000 HA % %1 % % %!Q / / / / / /rcbtdddksJtdddksJdS)NcdS)Nrrrs rrztest_trivial..Yarrc|Sr rrs rrztest_trivial..ZrrrrJrrr test_trivialrXsB KK # #q ( ( ( ( KK # #q ( ( ( ( ( (r)rAmpmathmpmath.calculus.optimizationrrrrrr r r r r rrFrKrQrUrXr_r~rrrrrs< BBBBBBBBBBBBBBBBBBBBBBBBBB%1%1%1NAAA   000 000&)))))r