ó ~9­\c@sYdZddlmZmZmZmZmZmZddlmZm Z m Z m Z edƒZ ede e ƒZ ede dƒZede e ƒZed e e ƒZed e e ƒZed e e ƒZed e dƒZed e dƒZede dƒZede dƒZd„Zd„Zd„Zd„Zd„Zd„Zd„ZdS(sa Some examples have been taken from: http://www.math.uwaterloo.ca/~hwolkowi//matrixcookbook.pdf iÿÿÿÿ(t MatrixSymboltInversetsymbolst DeterminanttTracet Derivative(tMatAddtIdentitytMatMult ZeroMatrixtktXtxitAtBtCtDtatbtctdcCsàtjtƒtttƒks$t‚tjjtƒttjtƒksNt‚dtjtƒtdttƒkszt‚tttƒjtƒttttƒtƒks°t‚ttjtƒttttƒksÜt‚dS(Ni(R tdiffRtAssertionErrortTRR(((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt(test_matrix_derivative_non_matrix_results $*,6cCs(tjtƒtttƒks$t‚dS(N(R RR RR(((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt$test_matrix_derivative_trivial_cases"scCsútjttƒt}|jtƒttƒj ttjttƒjksTt‚tttƒƒ}tt ttƒt ƒ}|jtƒtdt t tdj ks³t‚tttt ƒƒ}|jtƒttt ƒjd ksöt‚dS(Niÿÿÿÿi( RRRR RRRRRR R(texpr((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt#test_matrix_derivative_with_inverse(s=3cCs!tjtƒttƒks!t‚tjjtƒttƒksEt‚tjt}|jtƒtksmt‚tjt}|jtƒtks•t‚tjtt}|jtƒttjksÈt‚tjtjt}|jtƒttjksþt‚tjtt}|jtƒttjks1t‚tjtjt}|jtƒttjksgt‚tjtjtt }|jtƒttt jtt tjks´t‚t ttjt t tt }|jtƒt jt t tt t jt jt ttkst‚tjt t}|jtƒt tt jtksXt‚tjtjt tt }|jtƒt jttt jt tt tjks´t‚ttt jt ttt }|jtƒt ttt tjt jttt tjkst‚dS(N(R RRR RRRR RRRRRR(R((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt*test_matrix_derivative_vectors_and_scalars<s0!$  """"5!H*@!cCs£ ttƒt}|jtƒtttƒttƒks>t‚ttƒ}|jtƒttƒkskt‚tttƒ}|jtƒtjks™t‚tttt ƒ}|jtƒtjt jksÒt‚tttjt ƒ}|jtƒt tkst‚ttjtƒ}|jtƒtks6t‚tttjƒ}|jtƒtksdt‚ttdƒ}|jtƒdtjks–t‚ttdt ƒ}|jtƒtt t tjksÔt‚tt ttt ƒƒ}|jtƒtt t tjkst‚ttjt tƒ}|jtƒt tt jtksWt‚tt ttjƒ}|jtƒt tt jtks˜t‚tttjt ƒ}|jtƒt tt jtksÙt‚ttt tjƒ}|jtƒtt jtt kst‚tt tjtƒ}|jtƒtt jtt ks[t‚ttjtt ƒ}|jtƒtt jtt ksœt‚tttt tƒ}|jtƒtjtjt jt jtjtjksõt‚ttjtƒ}|jtƒdtks't‚tttjƒ}|jtƒdtksYt‚tt jtjt tt ƒ}|jtƒt jtt t jt tt t jks»t‚ttjt tt ƒ}|jtƒt tt t jtt jks t‚tttt tjt ƒ}|jtƒtjt jtt jt ttt ksjt‚tttt t ttt t jƒ}|jtƒdtjttt t t jksÎt‚tttƒ}ttttƒ}tt jtjt ttjt tt ƒ}|jtƒt ttjt tt t jt jtt t jtjt jtt tt t jtjt tt jttjt jtt t jksÎt‚tttdt ƒ}|jtƒt tƒj tjt jt tƒjks&t‚tt tjt tƒtƒ}|jtƒtj ƒj tjtj ƒt j ƒjtj ƒjtj ƒjttj ƒt j ƒtj ƒjksÉt‚ttjt tj ƒtjt tƒ}|jtƒdt ttjt tj ƒtjt ttjt tj ƒdt ttjt tj ƒksrt‚tttjt tj ƒtjt tƒ}|jtƒt tt ttjt tƒt tt ttjt tƒtjt tt ttjt tƒt jtt tjt tjtƒtjt jtt tjt tjtƒt jtt tjt tjtƒksŸ t‚dS(Niiÿÿÿÿiþÿÿÿ(RR RRRR RR RRRRRtinv(R((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt!test_matrix_derivatives_of_tracesjst. !% "********A"@5@+9 1«@!‚,}0c Csmtjtttjttttjtjttjtttjttjttttttjtjtttttjttttjdtttdtttjtjttjt}tttttjdttttjttjttttjttttttjtjtjtttttjtjtdttjtjtjttjttdtttjtjttjtttdttjtjtjttjttttjttdttjttjtttjttjttttjtjdttjttjttttjtjdtjtjttjtjdttjtttjtjdtjttjttjtttjtjdttttjtjttjttjtjtjtjttjtjttjttttjtjdtjtjttttjtjtdtjttttjtjtjtjttdttjtjtjttjtjtjttjtj}|jtƒ|ksit ‚dS(Nii*i( RRR R RRRRRR(Rtresult((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt+test_derivatives_of_complicated_matrix_exprõs ÜÿÿÿucCs‰dttƒ}|jtƒdttƒks5t‚t}|jtƒ}t|tƒs_t‚|tttƒkszt‚ttƒd}|jtƒdttƒttƒks¹t‚ttƒt}|jtƒt|tƒksít‚tttƒtƒ}|jtƒdttƒttƒks2t‚ttttƒtƒtƒ}|jtƒdttƒdttƒks…t‚dS(Nii( RR RRR Rt isinstanceR R(Rtderiv((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyt"test_mixed_deriv_mixed_expressionsûs%/$/ N(t__doc__tsympyRRRRRRRRRR R R R R RRRRRRRRRRRRR R#(((sPlib/python2.7/site-packages/sympy/matrices/expressions/tests/test_derivatives.pyts(."    . ‹