a 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?)mpZdpsfindrootZmpfZaesinpiZmpcr rr rr absepsrrrrr isinstanceZivfpfloatcomplexpytestZraises ValueErrorZworkprec)frrrrr test_findroots8     D . r4cCstddddddksJdS)NcSs |ddS)Nr!rrrrrrr/rz test_bisection..)rr!bisectrrr(rrrrtest_bisection-sr7cCs.dd}t|ddd}t||tks*JdS)NcSstgd|S)N)rrrr)Zpolyvalrrrrr2rztest_mnewton..gZmnewtonr)r(r+r,r3rrrr test_mnewton1sr9cCs.dd}t|dtd}t||tks*JdS)NcSs |ddS)Nr!drrrrrr7rztest_anewton..rr)r(r r+r,r8rrr test_anewton6sr;cCs.dd}t|dtd}t||tks*JdS)NcSsd|ddS)Nr!rrrrrrr<rztest_muller..rr)r(rr+r,r8rrr test_muller;sr<csDtddD]tfdddks Jq tddddks@JdS)Nrrcs |dSrrrirrrBrz#test_multiplicity..cSs|dSr rrrrrrCrr)rangeZ multiplicityrrr=rtest_multiplicity@sr@csdd}tt|dtddgddggdd ks4Jtt|ddd d d D]\}}qJt||d dksjJdd dd dd fdd}t|d}dd|DddgksJdS)NcWsXd|ddd|ddd|ddd|d|ddd|ddgS)Nrrr!rrrrrrr3Fsz test_multidimensional..f)rrArrBrgHz>cSs t|tSr)norminfrrrrrJrz'test_multidimensional..)verboserDr!g+=cSs| d|Sr rryrrrrOrcSs(|d||ddd||dS)Nr!rrrGrrrrPrcSst|d|dSr )sqrtrGrrrrQrcs&||}||||||fSrr)rrHZf1xf1f2Zf3rrr3Rs ) rMcSsg|]}tt|qSr)intround).0r>rrr Vrz)test_multidimensional..rr)ZmnormZjacobianZmatrixr r'rDr()r3rerrorrrJrtest_multidimensionalEs,  rScCs0tddddksJtddddks,JdS)NcSsdS)NrrrrrrrYrztest_trivial..rcSs|SrrrrrrrZrrr6rrrr test_trivialXsrT)r1ZmpmathZmpmath.calculus.optimizationrrrrrrr r r r r r4r7r9r;r<r@rSrTrrrrs4'