a ztest_diff..g@quadmethodcSs d|dSNr r r r r r rrr)rricSs d|dSrr r r r rrrstep) mpdpsdifflogaecossinabsexpejZdiffunr r r r test_diffs"  r"cCsdt_ddttddDddgks(JddttddddDddgksNJd dttdd Dgd kspJd dttdd ddDgd ksJdS) NrcSsg|] }t|qSr Zchop.0dr r r rztest_diffs..rrcSsg|] }t|qSr r#r$r r rr'rrrcSsg|] }t|qSr r#r$r r rr'rrrrrcSsg|] }t|qSr r#r$r r rr'r)rrZdiffsrr r r r test_diffss "&"r)cCs"dt_ttddgdksJdS)Nrr)rg?gg?g)rrZtaylorsqrtr r r r test_taylorsr,cCsdt_d\}}}}dd}t||dds2Jt||ddsHJt||d d s^Jt||d d stJt||d dsJt||ddsJt||ddsJt||ddsJt||ddsJt||ddsJt||ddsJt||dd s(Jt||dd s@Jt||ddsXJt||ddspJt||ddsJt||ddsJt||ddsJt||d d!sJt||d"d#sJt||d$d%sJt||d&d'sJt||d(d)s0Jt||d*d+sHJt||d,d-s`Jt||d.d/sxJt||d0d1sJdS)2Nr)rr cSs d|d|dd|dS)Nr rr r )r yzr r rr"rz#test_diff_partial..)rrri)rrri)rrriPr(i)rrriݤ)rrri05^)rrriR\)rrriA)rrri%)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrr)rrriBW)rrrib)rrriN)rrrigs)rrrinR)rrri/)rrril).)rrri )rrrip)rrrr)r r.r/Zxyzfr r rtest_diff_partials<r1N)Zmpmathr"r)r,r1r r r rs