U Zʹj< @sddlZddlTddlmZmZmZmZmZmZm Z m Z m Z m Z m Z ddZddZdd Zd d Zd d ZddZddZddZdS)N)*) SecantMuller BisectionIllinoisPegasusAndersonRidderANewtonNewtonMNewtonMDNewtonc Csdt_tddtdds"ttttdts:ttttdtdftsZttddtd d svtd d}t t t t t fD]&}t|d |d }t||tkstqt t tttttfD]&}t|d|d }t||tkstqdd}tt|dddtstttj|dddtjs.tttj|dddtsJtttj|dddtsfttt0tdtddtdddW5QRXW5QRXtttdddW5QRXtttdddW5QRXdS)NcSs d|dS)Nxrr\/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mpmath/tests/test_rootfinding.py ztest_findroot..g?rgQ @cSs ||dSNrrrrrr ry@@y?cSst|SN)cosrrrrrr@solver)?rcSs |ddSNrrrrrrrrg|=)Ztolry?r!cSs|dd|dS)Nr!iClF[rrrrrr"rz 5.278e+13z -5.278e+13)realimagcSsdS)Nrrrrrr'rrcSsd|ddS)Nrr)prrrr+rg?)mpZdpsfindrootZmpfZaeAssertionErrorsinpiZmpcr rr rr absepsrrrrr isinstanceZivfpfloatcomplexpytestZraises ValueErrorZworkprec)frrrrr test_findroots8       r5cCstddddddkstdS)NcSs |ddS)Nr!rrrrrrr/rz test_bisection..)rr!bisectrrr(r)rrrrtest_bisection-sr8cCs.dd}t|ddd}t||tks*tdS)NcSstddddg|S)Nrr)Zpolyvalrrrrr2rztest_mnewton..gZmnewtonr)r(r,r-r)r4rrrr test_mnewton1sr:cCs.dd}t|dtd}t||tks*tdS)NcSs |ddS)Nr!drrrrrr7rztest_anewton..rr)r(r r,r-r)r9rrr test_anewton6sr<cCs.dd}t|dtd}t||tks*tdS)NcSsd|ddS)Nr!rrrrrrr<rztest_muller..rr)r(rr,r-r)r9rrr test_muller;sr=csDtddD]tfdddks tq tddddks@tdS)Nrrcs |dSrrrirrrBrz#test_multiplicity..cSs|dSr rrrrrrCrr)rangeZ multiplicityr)rrr>rtest_multiplicity@srAcsdd}tt|dtddgddggdd ks4ttt|ddd d d D]\}}qJt||d dksjtdd dd dd fdd}t|d}dd|DddgkstdS)NcWsXd|ddd|ddd|ddd|d|ddd|ddgS)Nrrr!rrrrrrr4Fsz test_multidimensional..f)rrBrrCrgHz>cSs t|tSr)norminfrrrrrJrz'test_multidimensional..)verboserEr!g+=cSs| d|Sr rryrrrrOrcSs(|d||ddd||dS)Nr!rrrHrrrrPrcSst|d|dSr )sqrtrHrrrrQrcs&||}||||||fSrr)rrIZf1xf1f2Zf3rrr4Rs ) rNcSsg|]}tt|qSr)intround).0r?rrr Vsz)test_multidimensional..rr)ZmnormZjacobianmatrixr)r r'rEr()r4rerrorrrKrtest_multidimensionalEs,  rUcCs0tddddksttddddks,tdS)NcSsdS)NrrrrrrrYrztest_trivial..rcSs|SrrrrrrrZrrr7rrrr test_trivialXsrV)r2ZmpmathZmpmath.calculus.optimizationrrrrrrr r r r r r5r8r:r<r=rArUrVrrrrs4'