ó ¡¼™\c@sddlmZddlmZmZmZmZmZmZm Z ddlm Z m Z m Z m Z mZmZmZmZddlmZddlmZddlmZmZmZddlmZmZedƒ\ZZZZd „Zd „Zd „Z d „Z!d „Z"dS(iÿÿÿÿ(t Quaternion(tsymbolstretimtAddtMultItAbs(tcostsintsqrtt conjugatetlogtacostEtpi(traises(tMatrix(tdifft integratettrigsimp(tStRationalsx y z wc Cs‰tttttƒ}||tdtdtdtdtƒksJt‚tjtdƒdtdƒdtdƒdfdtdƒ}|tt ddƒt ddƒt ddƒt ddƒƒksÐt‚t t tƒt tƒ dgt tƒt tƒdgdddggƒ}t tj|ƒƒ}|ttdƒtt tƒdƒdddtdt tƒdƒdƒks…t‚dS(Niiiiiþÿÿÿ(RtxtytztwtAssertionErrortfrom_axis_angleR RRRRR Rtfrom_rotation_matrix(tqtq2tMtq3((sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyttest_quaternion_construction s5A'Ic Cstddtƒ}tddtƒ}tddtƒ}tttttƒ}||ttt|ƒtt |ƒttƒks„t ‚d|tdttttƒks­t ‚t |ttdtttƒksÖt ‚||tt|tttƒksÿt ‚||t |tttttƒdtƒks3t ‚||t tttttƒ|dtƒksgt ‚||t |tttttƒdtƒks›t ‚| tt t t t ƒksÁt ‚td d t d d t d ddt dtƒ}tdd ddƒ}|d d t td dt d d t d ddt ƒksNt ‚|d d t td dddƒks{t ‚|d d t td d t d d t d d t d d t d d d t ddt ƒksät ‚|d d t tddddƒkst ‚dS(Ntatcomplextbtrealtct commutativeitevaluateiiiiiiit real_fieldiöÿÿÿi i&iûÿÿÿ(RtTruetFalseRRRRRRRRRRR(R#R%R'Rtq1R((sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyt%test_quaternion_complex_real_additions&9)))444&3E-Zc sótttttƒ}tddddƒ}tddddƒ‰t|ƒttt t t ƒksit‚|jƒttdtdtdtdƒks£t‚|j ƒtttttƒttdtdtdtdƒksðt‚|j ƒttt t t ƒtdtdtdtdks:t‚t t ‡fd†ƒ|j dƒttd tdtdtddttdttdttƒks¯t‚|j dƒttdƒ d tdƒ d tdƒ d tdƒ d ƒkst‚|j d ƒtks t‚|jƒttttd ƒƒdtd ƒtttd ƒƒd dtd ƒtttd ƒƒd dtd ƒtttd ƒƒd ƒks½t‚|jƒtttd ƒƒdtd ƒttd ƒd ƒd dtd ƒttd ƒd ƒd dtd ƒttd ƒd ƒd ƒksVt‚|jdƒtd tdttd ƒd ƒƒdtd ƒtdttd ƒd ƒƒd dtd ƒtdttd ƒd ƒƒd dtd ƒtdttd ƒd ƒƒd ƒks"t‚ttttttƒtƒtddddƒks[t‚ttttttƒtƒttddtddtddtddƒks´t‚tjd|ƒtdƒddtdƒdfksït‚dS(Niiiiics ˆjƒS(N(tinverse((tq0(sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyt<tiþÿÿÿiiái–gà¿iii<iZixi(iii(RRRRRR RtnormR t normalizeR/Rt ValueErrortpowRtNotImplementedtexpRRR t_lnR R t pow_cos_sinRRt rotate_point(RR-((R0sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyttest_quaternion_functions3s<*:MJ_V %%4 %%4 $//>9>c Cs¢tddddƒ}|jƒdtdƒddtdƒddtdƒdfdttdƒdƒfkszt‚|jƒttdƒ dtdƒdtdƒdgtdƒdtdƒ dtdƒdgtdƒdtd ƒdtdƒdggƒkst‚|jdƒttdƒ dtdƒdtdƒdtdƒd gtdƒdtdƒ dtdƒdtd ƒgtdƒdtd ƒdtdƒdtdƒ d gtd ƒtd ƒtd ƒtdƒggƒkst‚td d t ƒ}tt |dƒd d t |dƒƒ}t |jƒƒtt |ƒt |ƒ d gt |ƒt |ƒd gd d dggƒks t‚|jƒd d t |dƒt t |dƒƒfdtt |dƒƒfksùt‚t |jdƒƒtt |ƒt |ƒ d t |ƒt |ƒdgt |ƒt |ƒd t |ƒ t |ƒdgd d dd gd d d dggƒksžt‚dS(Niiiiiiii iiitthetaR&(iii(iii(Rt to_axis_angleR R Rtto_rotation_matrixRRRR+RR RR(R-R=R((sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyttest_quaternion_conversionsZs2&7+<G489)3&00cCsŽtttdƒttdƒddƒ}t|jƒƒtdddgdttƒttƒ gdttƒttƒggƒksŠt‚dS(s3 There was a sign mistake in the definition, of the rotation matrix. This tests that particular sign mistake. See issue 1593 for reference. See wikipedia https://en.wikipedia.org/wiki/Quaternions_and_spatial_rotation#Quaternion-derived_rotation_matrix for the correct definition iiiN(RRRR RR?RR(R((sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyt test_quaternion_rotation_iss1593}s ) N(#tsympy.algebras.quaternionRtsympyRRRRRRRRR R R R R RRtsympy.utilities.pytestRRRRRRRRRRRR"R.R<R@RA(((sClib/python2.7/site-packages/sympy/algebras/tests/test_quaternion.pyts4:  ' #