U Zʹj @s,ddlTddZddZddZdd Zd S) )*cCsfdt_ttdddtds"tttdtd s ztest_diff..g@quadmethodcSs d|dSNr r r r r r rrr)rricSs d|dSrr r r r rrrstep) mpdpsdifflogaeAssertionErrorcossinabsexpejZdiffunr r r r test_diffs"  r#cCsdt_ddttddDddgks(tddttddddDddgksNtd dttdd Ddddgksrtd dttdd ddDdddgkstdS) NrcSsg|] }t|qSr Zchop.0dr r r sztest_diffs..rrcSsg|] }t|qSr r$r%r r rr(srrcSsg|] }t|qSr r$r%r r rr(srcSsg|] }t|qSr r$r%r r rr(s)rrZdiffsrrr r r r test_diffss "&$r)cCs(dt_ttdddddddgks$tdS)Nrrg?gg?g)rrZtaylorsqrtrr r r r test_taylorsr,cCsdt_d\}}}}dd}t||dds2tt||ddsHtt||d d s^tt||d d sttt||d dstt||ddstt||ddstt||ddstt||ddstt||ddstt||ddstt||dd s(tt||dd s@tt||ddsXtt||ddsptt||ddstt||ddstt||ddstt||d d!stt||d"d#stt||d$d%stt||d&d'stt||d(d)s0tt||d*d+sHtt||d,d-s`tt||d.d/sxtt||d0d1stdS)2Nr)rr cSs d|d|dd|dS)Nr rr r )r yzr r rr"rz#test_diff_partial..)rrri)rrri)rrriP)rrri)rrriݤ)rrri05^)rrriR\)rrriA)rrri%)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrriBW)rrrib)rrriN)rrrigs)rrrinR)rrri/)rrril).)rrri )rrrip)rrrrr)r r.r/Zxyzfr r rtest_diff_partials<r1N)Zmpmathr#r)r,r1r r r rs