B <Z @s,ddlTddZddZddZdd Zd S) )*cCsfdt_ttdddtds"tttdtd s sztest_diff..g@quad)methodcSs d|dS)Nr r r )r r r r rs)rricSs d|dS)Nr r r )r r r r rsstep) mpdpsdifflogaeAssertionErrorZcossinabsZexpejZdiffunr r r r test_diffs"  rcCsdt_ddttddDddgks(tddttddddDddgksNtd dttdd Ddddgksrtd dttdd ddDdddgkstdS) NrcSsg|] }t|qSr )chop).0dr r r sztest_diffs..rrcSsg|] }t|qSr )r)rr r r r r!sr)rcSsg|] }t|qSr )r)rr r r r r!srcSsg|] }t|qSr )r)rr r r r r!s)rrZdiffsrrr r r r test_diffss "&$r"cCs(dt_ttdddddddgks$tdS)Nrrg?gg?g)rrZtaylorZsqrtrr 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 r r"sz#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 r test_diff_partials<r)N)Zmpmathrr"r$r)r r r r s