a z$test_approximation..rT)errorgh㈵> g@)mpdpsZchebyfitrangeabspolyval)fperrirrrr test_approximations   rcCsBdt_tdddtdds&Jtddtts>JdS)NrcSs|t||dS)N)sinrrrr r r ztest_limits..rcSsdd||SNrr)nrrr r r )rrlimitaempfinferrrr test_limits s r%cCsntgddksJtdgddks&Jtdgddks:Jgd}t|ddksTJt|dddd ksjJdS) Nrr)r rr&r T)Z derivative)r()rrrrr test_polyvals r+cCstddg}|ddsJtgd\}}|dtdtsHJ|dtdtsbJtdggkstJttdd dS) Nrrr )rrrrcSs tdgS)Nr) polyrootsrrrr r #r z test_polyroots..)r.r!sqrtjpytestraises ValueError)rqrrr test_polyrootss r5c Csd}gd}tdttj&t|ddd|ddWdn1sP0Yt|d dd|dd}d d |D}|gd ksJWdn1s0YdS) N@Al #!yUMDQj~U4 rl`< && xXD`Drl #~O6dIW6o b9rl N$UVOF >`}Drl h a0r3(| Wo9Erlcs9lM hPc;/,3rl p~ #]lsrl J={|L m~;`rl`|fl?'h_Krl eo.!03a7# rl /$+-kMsf]`:rl ;bfLHQL Rrl / v\QB8%rl e]r+a ;M~DMLrl`|]B`3*0Vwk _rl J5]c<OJg1rl x,a6F~E8Prl V'6@h>)rl~&-ME A"rl<-AGn,a;0$SrlW?n05%Zurlpbz^.qrlc dqTTrlh 9N$$).lF4rl Ns,1j rl#vsKn?b;rl`<x!url#!Oerr&TFrmaxstepscleanupr extraprec2cSsg|] }t|qSr)str).0rrrr Pr z+test_polyroots_legendre..)@-0.999-0.996-0.991-0.983-0.973-0.961-0.946-0.93-0.911-0.889-0.866-0.841-0.813-0.784-0.753-0.72-0.685-0.649-0.611-0.572-0.531-0.489-0.446-0.402-0.357-0.311-0.265-0.217-0.17-0.121-0.073-0.02430.02430.0730.1210.170.2170.265z0.3110.3570.4020.4460.4890.5310.5720.6110.6490.6850.720.7530.7840.8130.8410.8660.8890.9110.930.9460.9610.9730.9830.9910.9960.999)rworkdpsr1r2 NoConvergencer.)rcoeffsrootsrrr test_polyroots_legendre%s"  $ rc Csd}gd}tgd}tdtj&t|dddd|d}Wdn1sT0Yttj"t|d dd|dWdn1s0Yt|d dd||d \}}tt|t| t |ksJt|d dd||dd d \}}tt|t| t |ksJdS) Ndr7)ArArBrCrDrErFrGrHrIrJrKrLrMrNrOrPrQrRrSrTrUrVrWrXrYrZr[r\r]r^r_r`rarbrcrdrerfz 0.311rgrhrirjrkrlrmrnrorprqrrrsrtrurvrwrxryrzr{r|r}r~rz1.0rr<TFr8r&)r9r:rr; roots_init<) Zmatrixrrrr.r1r2rmaxapplyr)Z extra_precrrZ roots_exactrrZroots1Zerr1rrr test_polyroots_legendre_init]s&"  $ $  "  rc Cstd}dt_d}|g}d}td|dD]}||9}|||q*t||d|d\}}tdddD]@}t|ddd|t|ddd|}|t |dsjJqjd t_dS) Nrrrrg?r-g|=r) r"rrrappendZpadeZarangerr!exp) ZoneNakrrr4rr?rrr test_pades(rcCsdt_tddddgd\}}|dds0J|dd td dtsRJ|dd td d tstJ|ddksJ|dd dtsJ|dd d tsJt||fddgdd sJdS) NrcSs|dSrrrrrr r r ztest_fourier..r-rrg?rrr g?)rrZfourierr!r/piZ fourierval)csrrr test_fouriers""rcCs4dt_tdddddtdtds0JdS)NrcSs|SNrtrrr r r z test_differint..rgr r)rrZ differintr!r/rrrrr test_differintsrcCsdt_d}dd}dd}||}t||dd|s:Jt||dd|sRJt||d d|sjJd }||}t||dd|sJt||dd|sJt||d d|sJd}d d}d d}||}t||dd|sJt||dd|s Jt||d d|s&Jd }||}t||dd|sLJt||dd|sfJt||d d|sJdS) Nrg{Gz?cSsd|ddS)Nrrrr*rrr r r ztest_invlap..cSs|t| Sr)rrrrr r r Ztalbot)methodZstehfestZdehoogg?cSs t||Sr)logr*rrr r r cSst t|Sr)Zeulerrrrrr r r )rrZ invertlaplacer!)rfpftZfttrrr test_invlaps2r) r1Zmpmathrr%r+r5rrrrrrrrrr s   8>