+
    i|}             	       p	   ^ RI t ^ RIt^ RIHt ^ RIHt ^ RIHtH	t	H
t
 ^ RIHtHtHtHtHtHtHtHt ^ RIHt ^ RIHt ^ RIHtHt ^ R	IHt ^ R
IHt ^ RIH t H!t! ^ RI"H#t#H$t$H%t% ^ RI&H't'H(t(H)t) ^ RI*H+t+ ^ RI,H-t- ^ RI.H/t/H0t0H1t1 ^ RI2H3t3H4t4 ^ RI5H6t6 ^ RI7H8t8 ^ RI9H:t: ^ RI;H<t<H=t= ^ RI>H?t? ^ RI@HAtAHBtBHCtC ^ RIDHEtE ^ RIFHGtG ^ RIHHItI ^ RIJHKtKHLtLHMtMHNtNHOtOHPtPHQtQHRtRHStSHTtTHUtUHVtVHWtWHXtXHYtYHZtZH[t[H\t\H]t]H^t^H_t_H`t`HataHbtb ^ RIHHctc ^ RIdHeteHftf ^ RIgHhth ^ RIiHjtjHktkHltl ^ R ImHntn ^ R!IoHptpHqtqHrtrHstsHtttHutu ^ R"IvHwtw ^ R#IxHyty ^ R$IzH{t{ ^ R%I|H}t} ^ R&I~Ht ^ R'IHt ^ R(IHtHtHtHtHtHtHtHt ]N]O]L]M3t^ R)IHtHt ^ R*IHtHt ]u! 4       ;_uu_ 4         ! R+ R,]]]4      tRRR4       ]u! 4       ;_uu_ 4         ! R- R.]]]4      tRRR4       R/ tR0 tR1 tR2 tR3 tR4 tR5 tR6 tR7 tR8 tR9 tR: tR; tR< tR= tR> tR? tR@ tRA tRB tRC tRD tRE tRF tRG tRH tRI tRJ tRK tRL tRM tRN tRO tRP tRQ tRR tRS tRT tRU tRV tRW tRX tRY tRZ tR[ tR\ tR] tR^ tR_ tR` tRa tRb tRc tRd tRe tRf tRg tRh tRi tRj tRk tRl tRm tRn tRo tRp tRq tRr tRs tRt tRu tRv tRw tRx tRy tRz tR{ tR| tR} tR~ tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR tR Et R EtR EtR EtR EtR EtR EtR EtR EtR Et	R Et
R EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR EtR Et R Et!R Et"R Et#R Et$R Et%R Et&R Et'R Et(R Et)]qR 4       Et*R Et+R Et,R Et-R Et.R Et/R Et0]rR 4       Et1R Et2R Et3R Et4R Et5R Et6R Et7R Et8R Et9R Et:R Et;R Et<R Et=R Et>R Et?R Et@R EtAR EtBR EtCR EtDR EtER EtFR EtGR EtHR EtIR EtJ]t! R4      R 4       EtKR EtLR EtMR EtNR EtOR EtPR EtQR EtRR EtSR EtTR EtUR EtVR EtWR EtXR EtYR EtZR#   + '       g   i     EL; i  + '       g   i     EL; i)    N)Hashable)Add)Functiondiffexpand)EFloatIIntegerRationalnanoopi)Pow)S)Symbolsymbols)sympify)Abs)explog)MaxMinsqrt)cossintan)	integrate)	transpose)HermitianOperatorOperatorDagger)PolyPurePoly)RootOf)sstr)	FiniteSet)signsimpsimplify)trigsimp)
ShapeErrorMatrixErrorNonSquareMatrixErrorDeferredVector)_find_reasonable_pivot_naive)	_simplify)GramSchmidtImmutableMatrixImmutableSparseMatrixMatrixSparseMatrix
casoratiandiageyehessianmatrix_multiply_elementwiseones
randMatrix	rot_axis1	rot_axis2	rot_axis3	wronskianzerosMutableDenseMatrixImmutableDenseMatrixMatrixSymboldotprodsimprot_ccw_axis1rot_ccw_axis2rot_ccw_axis3)_dotprodsimp_state)TupleWild)KroneckerDelta)flattencaptureiterable)ignore_warnings)raisesXFAILslowskipskip_under_pyodidewarns_deprecated_sympyQ)Array)ArrayDerivative)MatPow)
Quaternion)O)abcdxyzt)_MinimalMatrix_CastableMatrix)MatrixSubspacesMatrixReductionsc                       ] tR t^GtRtR# )SubspaceOnlyMatrix N__name__
__module____qualname____firstlineno____static_attributes__rm       `/var/www/html/photoedit/myenv/lib/python3.14/site-packages/sympy/matrices/tests/test_matrices.pyrl   rl   G       rt   rl   c                       ] tR t^LtRtR# )ReductionsOnlyMatrixrm   Nrn   rm   rt   ru   rx   rx   L   rv   rt   rx   c                     \        W R  4      # )c                     \        W8H  4      # N)intijs   &&ru   <lambda> eye_Reductions.<locals>.<lambda>Q   s
    3qv;rt   rx   ns   &ru   eye_Reductionsr   P   s    &>??rt   c                     \        W R  4      # )c                     ^ # r   rm   r}   s   &&ru   r   "zeros_Reductions.<locals>.<lambda>U   s    1rt   r   r   s   &ru   zeros_Reductionsr   T   s    n55rt   c                     \        \        4       EF=  w  rVP                  ^^4      pVP                  R8X  dO   \        ;QJ d&    R VP                   4       F  '       d   K   RM	  RM! R VP                   4       4      '       g   Q hVP
                  ^8X  d   \        VP
                  4      \        J g   Q hVP                  ^8X  d   \        VP                  4      \        J g   Q hV ^,          '       g3   \        VP                  4       4      \        \        \        39   g   Q hEK  \        VP                  4       4      \        J d   EK>  Q h	  R# )   c              3   D   "   T F  p\        V4      \        J x  K  	  R # 5ir{   )typer|   ).0r~   s   & ru   	<genexpr>test_args.<locals>.<genexpr>\   s     (IAaCs    FTN)r      )	enumerateclassesrB   shapeallrowsr   r|   colsflatlisttuplerK   todokdict)r   clsms      ru   	test_argsr   X   s    G$IIaOww& SS(I(ISSS(I(I%I%IIIvv{tAFF|s222vv{tAFF|s2221uu>dE5%9999	?d*** %rt   c                     \         \        3 FS  p V P                  ^^4      p\        4       ;_uu_ 4        VP                  pRRR4       XVP                  4       8X  d   KS  Q h	  \        \        3 FS  p V P                  ^^4      p\        4       ;_uu_ 4        VP                  pRRR4       XVP                  4       8X  d   KS  Q h	  R#   + '       g   i     L; i  + '       g   i     L@; i)r   N)
r5   r3   rB   rW   _matr   r6   r4   _smatr   )r   r   matsmats       ru   test_deprecated_mat_smatr   e   s    &IIaO#%%&&C &affh	 '
 22IIaO#%%77D &qwwy   	 3 &%
 &%s   C	C	CC,c                      \        ^^\        \        .4      p V \        ,          \        ^^\        \        ,          \        \        ,          .4      8X  g   Q hR#    N)r5   rc   rd   re   )vs    ru   test_divisionr   r   s;    q!aVAQ3&A!QqSz****rt   c            	        a a \        . RO\        \        \        .^\        ,          R\        \        ,          ..4      o S S ,           \        . RO^\        ,          ^\        ,          ^\        ,          .^\        ,          R^\        ,          \        ,          ..4      8X  g   Q h\        ^^^^.4      o\	        \
        V V3R l4       R# )r   c                     < S S,           # r{   rm   r   r   s   ru   r   test_sum.<locals>.<lambda>{   s	    q1urt   Nr   r   r   i)r         i)r5   rc   rd   re   rR   r+   r   s   @@ru   test_sumr   w   s    	Aq!9qsC1o67Aq5FI!QqS!A#1dAaCE8JKLLLLq!aVA
:}%rt   c                      \        ^^R\        .4      p \        ^^^\        \        4      .4      p\        V 4      V8X  g   Q hR# )r   N)r5   rc   r   absr   s     ru   test_absr   }   s:    q!b!WAq!aQ[!Aq6Q;;rt   c                      \        R4      p \        R4      pW,           V P                  ! V4      u;8X  d   \        ^^.^^..4      8X  g   Q h Q hR# )r   N)r   r   r   r   r   )r   r   )r5   addr_   r`   s     ru   test_additionr      sY      	A
 	  	A
 5AEE!H8AA'7 888888rt   c            
      ,   \         \        3 EF  p V ! ^^\        ^	4      4      pWR,          8X  g   Q hVR,          \        ^^. RO4      8X  g   Q hVR,          \        . RO4      8X  g   Q hV^ ^.R3,          \        . RO. RO.4      8X  g   Q hV^ ^.^3,          V^ ^.^.3,          8X  g   Q hV^^ ^.3,          V^.^ ^.3,          8X  g   Q hVR^ ^.3,          \        ^ ^.^^.^^..4      8X  g   Q hVR	,          ^ 8X  g   Q hVR
,          \        . RO. RO.4      8X  g   Q hVR,          \        ^ ^.^^.^^..4      8X  g   Q hVR,          V^ ^.^3,          8X  g   Q hVR,          V^^ ^.3,          8X  g   Q hV ! ^^\        ^	4      4      pV. ROR3,          \        . RO. RO. RO. RO. RO.4      8X  g   Q hVR. RO3,          \        . RO. RO. RO.4      8X  d   EK  Q h	  \        P                  ! ^4      p^VR&   ^VR&   ^VR&   VP                  ! ^^.^.4      \        ^.^..4      8X  g   Q hVP                  ! ^^ .. RO4      \        . RO. RO.4      8X  g   Q hVP                  ! . RO. RO4      \        . RO. RO. RO. RO.4      8X  g   Q hR# )r   NNNNr   r   r   r   r   r      r   r   )r   r      r   r   r   r   r   )r   r   Nr   )r   r   )NNr   r   )r   r   )r   r   r   r   r   )r   r      )r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   )r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   r   r   )r   r   r   r   )r5   r6   rangerB   extract)Mr_   s     ru   test_fancy_index_matrixr      s   l#aE!HdG||w&Ay1111w&++++!Q|vy)&<====!Q|q!Q!~---QF|q!q!f~---QF|v1v1v1v&>????w!||yFIy#9::::yFQFQFQF#;<<<<yAq!faiL(((yAa!QiL(((aE!H!#$0 )  	  K F, %  	 + $4 	1AAdGAdGAdG99aVaS!VCC-	 &
 
 
 
 99aVY'63 ,    99\<0F	< 5   rt   c                    aa \        R4      o\        R4      oSS,          p V R,          ^8X  g   Q hV R	,          ^8X  g   Q hV R
,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^ 8X  g   Q h \        R4       V R,          ^8X  g   Q hV R	,          ^8X  g   Q hV R
,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^ 8X  g   Q h \        SV 4      pVSP                  ! V 4      8X  g   Q hVR,          ^8X  g   Q hVR	,          ^8X  g   Q hVR
,          ^8X  g   Q hVR,          ^8X  g   Q hVR,          ^ 8X  g   Q hVR,          ^ 8X  g   Q h\        \        VV3R l4       S\        R4      ,          p \        V \         4      '       g   Q hV R,          \        8X  g   Q hV R	,          ^\        ,          8X  g   Q hV R
,          ^\        ,          8X  g   Q hV R,          ^ 8X  g   Q h\        S,          pW8X  g   Q h^S,          p \        V \         4      '       g   Q hV R,          ^8X  g   Q hV R	,          ^
8X  g   Q hV R
,          ^8X  g   Q hV R,          ^ 8X  g   Q h \        R4       \        V \         4      '       g   Q hV R,          ^8X  g   Q hV R	,          ^
8X  g   Q hV R
,          ^8X  g   Q hV R,          ^ 8X  g   Q hR#   \         d     ELi ; i  \         d     R# i ; i)r   z	c = a @ bc                     < \        S S4      # r{   )r;   r   s   ru   r   %test_multiplication.<locals>.<lambda>   s    :1a@rt   rc   z	c = 5 @ bN)r   r   r   r   r   r   r   r   r   r   r   r   r   r   )
r5   evalSyntaxErrorr;   multiply_elementwiserR   r+   r   
isinstancerc   )ra   hc2r_   r`   s      @@ru   test_multiplicationr      s     	A 	  	A
 	
!AT7a<<T7a<<T7a<<T7a<<T7b==T7a<<
[ w!||w!||w!||w!||w"}}w!||#Aq)A&&q))))T7a<<T7a<<T7b==T7a<<T7a<<T7a<<
:@A	F3KAa    T7a<<T7ac>>T7ac>>T7a<<	
QB7N7	AAa    T7a<<T7c>>T7c>>T7a<<	[ !V$$$$w!||w#~~w#~~w!||]  N  s$   K) K; )K87K8;L
	L
c                  
  aa \        \        ^ .^ \        ..4      p V ^,          V 8X  g   Q h\        \        \        .^ ^ ..4      p V ^,          \        \        \        .\        \        ..4      8X  g   Q h\        ^ ^ ^ \        ^4      ) ^,          .. RO. RO\        ^4      ) ^,          ^ ^ ^ ..4      o\        \        \
        ^,          ,          ^ \        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           \        \
        ^,          ,          ^,          \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \        ,          \        ,          \
        ^,          ,          ^0,          ,           \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \
        ^,          4      ,           .^ \        \
        ^,          ,          ^,          \        \
        ^,          4      \        \
        ^,          4      .\        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           \        \
        ^,          4      \        \        ^,          ,          \
        ^,          ,          ^ ,          \        \        ,          \        ,          \
        ^,          ,          ^0,          ,           \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \
        ^,          4      ,           \        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           .\        \
        ^,          ,          ^,          \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \        ,          \        ,          \
        ^,          ,          ^0,          ,           \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \
        ^,          4      ,           \        \
        ^,          4      \        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           \        \
        ^,          ,          ^,          ^\        ,          \        ^,          ,          \
        ^,          ,          ^@,          ,           \        \        ,          \        ,          \
        ^,          ,          ^ ,          ,           ^\        ,          \        ^,          ,          \
        ^,          ,          ^@,          ,           \        \
        ^	,          4      ,           ..4      o\        \        ) \
        ^,          ,          ^,          \        \        ^,          ,          \
        ^,          ,          ^@,          ,
          \        \        ,          \        ,          \
        ^,          ,          ^`,          ,
          \        \        ^,          ,          \
        ^,          ,          ^@,          ,
          \        \
        ^,          4      ,           \        \
        ^,          4      \        ) \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,
          \        \
        ^,          4      ,           \        ) \
        ^,          ,          ^,          ^\        ,          \        ^,          ,          \
        ^,          ,          ^,          ,
          \        \        ,          \        ,          \
        ^,          ,          ^@,          ,
          ^\        ,          \        ^,          ,          \
        ^,          ,          ^,          ,
          \        \
        ^	,          4      ,           .^ \        \
        ^,          ,          ^,          \        \
        ^,          4      \        \
        ^,          4      .\        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           \        \
        ^,          4      \        \        ^,          ,          \
        ^,          ,          ^ ,          \        \        ,          \        ,          \
        ^,          ,          ^0,          ,           \        \        ^,          ,          \
        ^,          ,          ^ ,          ,           \        \
        ^,          4      ,           \        \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,           \        \
        ^,          4      ,           .\        ) \
        ^,          ,          ^,          ^ \        ) \        ,          \
        ^,          ,          ^,          \        \        ,          \
        ^,          ,          ^,          ,
          \        \
        ^,          4      ,           \        ) \
        ^,          ,          ^,          \        \        ^,          ,          \
        ^,          ,          ^@,          ,
          \        \        ,          \        ,          \
        ^,          ,          ^`,          ,
          \        \        ^,          ,          \
        ^,          ,          ^@,          ,
          \        \
        ^,          4      ,           ..4      p\        ^^VV3R l4      pSS,          Tu;8X  d	   V8X  g   Q h Q hR# )r   c                 F   <a a \        VVV V3R  l\        ^4       4       !  # )c              3   Z   <"   T F   pSSV3,          SVS3,          ,          x  K"  	  R # 5ir{   rm   )r   kABr~   r   s   & ru   r   Atest_multiplication_inf_zero.<locals>.<lambda>.<locals>.<genexpr>  s'     (I11Q3!A#s   (+)r   r   )r~   r   r   r   s   ffru   r   .test_multiplication_inf_zero.<locals>.<lambda>  s    3(Ia(I#Jrt   Nr   r   r   r   r   r   r   r   )	r5   r   r   r   r   rc   r`   r_   r^   )r   CC2r   r   s      @@ru   test_multiplication_inf_zeror     s.   Q!R!"A6Q;;R1a&!"A6Vc3Z#s45555	
Aq1Q4%'
A$q!Q	 	A 		AqD!RT!Q$Yq[2a419Q;.1a48"QT'!)bAgaQRdlSUo:UXZ[\X\]^X^_`bc_cXcdfXf:fiklmoplpipqrtuquiuvxix:x{|}~  AB  ~B  |C  ;C  	D	
Bq!tGAIqAw!Q$(	Aad1r!tAqDy{	"Qq!tW	,a1gr!Q$wq!t|BAaPQSTPTUW7WZ\]^`a]aZabcefbfZfgiZi7ilmnoqrnrls7suwxyuyz{}~z~u~  @B  vB  EG  HI  EI  JK  MN  JN  EN  OQ  EQ  vQ  TU  VW  YZ  VZ  T[  v[  	\	AqDR1WQT\"_	$r!tAvad{2~	51a41R	G!AqD'	QSTUVXYUYSZ\^_`\`abdeae\efh\hkmnokopqstptktuwkw\wz{|}  @A  }A  {B  ]B  DF  GH  JK  GK  DK  LM  DM  PQ  RT  PT  UV  XY  UY  PY  Z[  ]^  Z^  P^  _a  Pa  Da  df  gh  dh  ij  dj  kl  no  ko  do  pr  dr  Dr  uv  wy  uy  z{  }~  z~  u~  @  BC  C  uC  DF  uF  DF  IJ  KL  NO  KO  IP  DP  	Q	R 	SA 	SAXaZ"QT'!Q$,r/!BqDF1a4KN2R1WQT\"_DqAwNRSTUWXTXRYrtqtuvqvwxz{w{q{|~q~  BD  EF  BF  GH  JK  GK  BK  LN  BN  rN  QR  ST  VW  SW  QX  rX  []  Z]  ^_  ab  ^b  Zb  cd  Zd  gh  ik  gk  lm  op  lp  gp  qr  tu  qu  gu  vy  gy  Zy  |~  @  |@  AB  |B  CD  FG  CG  |G  HJ  |J  ZJ  MN  OQ  MQ  RS  UV  RV  MV  WX  Z[  W[  M[  \_  M_  Z_  bc  de  gh  dh  bi  Zi  jMNPRSTVWSWPWXYPY  QR  ST  VW  SW  QX  bc  de  gh  dh  bi  j+-a419Q;Aad1+DqAw+NRSTUWXTXRY[]^_ab^b[bcdfgcg[ghj[jmopqmqrsmstuwxtxmxy{m{[{  A  BC  EF  BF  F  GH  JK  GK  K  LN  N  \N  QR  ST  VW  SW  QX  \X  DF  GH  DH  IJ  LM  IM  DM  NP  DP  SU  VW  SW  XY  [\  X\  S\  ]_  S_  D_  bc  de  gh  dh  bi  Di  jEGC1HQJXYrtqtuvqvwxz{w{q{|~q~  BD  EF  BF  GH  JK  GK  BK  LN  BN  rN  QR  ST  VW  SW  QX  rX  ac  `c  de  gh  dh  `h  ij  `j  mo  pq  st  pt  mt  uv  xy  uy  my  z|  m|  `|  A  BC  C  DE  E  FG  IJ  FJ  J  KM  M  `M  PR  ST  VW  SW  PW  XY  [\  X\  P\  ]_  P_  `_  bc  de  gh  dh  bi  `i  j	k 	lA 
1J	KBQ3!>r>>>rt   c                    aaaa \        \        R  4       \        p \        ^^.^^..4      oSR,          R,          V ! R4      ^,          V ! ^4      ^,          V ! ^34      ^,          V ! R4      ^,          .8X  g   Q hS^,          R,          . RO8X  g   Q h\        . R O. R!O. R"O.4      oS^,          R,          . R#O8X  g   Q hS^ ,          \	        ^4      8X  g   Q hS^,          S8X  g   Q h\        ^..4      ^d,          R$,          ^^d,          8X  g   Q h\	        ^4      R,          \	        ^4      8X  g   Q h\        ^^.^^..4      \        ^4      ,          \        ^^
.^^..4      8X  g   Q h\        ^!^.^0^9..4      oS\        P                  ,          R,          . R%O8X  g   Q h\        ^ ^.R&^..4      oS\        P                  ,          ^,          S8X  g   Q h\        ^^ .^^..4      \        P                  ,          \        ^^ .\        P                  ^..4      8X  g   Q h\        ^^ .^^..4      R,          \        ^^ .R^..4      8X  g   Q h^ RIH	o \        ^\        .^ ^..4      S,          \        ^\        S,          .^ ^..4      8X  g   Q h\        \        \        .^ \        ..4      S,          \        \        S,          \        \        S^,
          ,          ,          S,          .^ \        S,          ..4      8X  g   Q h\        \        S,          \        S^,
          ,          S,          \        S,          S^,          ,          \        S,          S,          ,
          ^\        ^,          ,          ,          .^ \        S,          \        S^,
          ,          S,          .^ ^ \        S,          ..4      '       g   Q h\        \        ^^ .^ \        ^ .^ ^ \        ..4      S,          \        \        S,          \        S^,
          ,          S,          ^ .^ \        S,          ^ .^ ^ \        S,          ..4      8X  g   Q h\        ^^ .^^..4      oSP                  \        ^4      4      SP                  ^4      8X  g   Q h\        ^..4      oS^
,          \        R'..4      u;8X  d8   SP                  \        ^
4      4      u;8X  d   SP                  ^
4      8X  g   Q h Q h\        . R(O. R)O. R*O. R+O. R,O.4      o\        \        V3R l4       \        \        R 4       \        . R-O. R.O. R/O.4      oSR,          \        . R/O. R/O. R/O.4      8X  g   Q h\        \        V3R l4       \        \        V3R	 l4       \        ^^.^^..4      oSR,          \        R
R.RR..4      8X  g   Q h\        . R.O. R.O. R.O.4      oSR,          \        . R.O. R.O. R.O.4      8X  g   Q h\        . R-O. R.O. R.O.4      oSR,          \        . R.O. R.O. R.O.4      8X  g   Q h\!        RRR7      o\#        SS,          \$        4      '       g   Q h\!        RRRR7      o\        \        VV3R l4       \!        RRRR7      oSS,          \        \'        ^ S4      \'        ^S4      \'        ^ S4      ) \'        ^S4      ,
          ^,           .^ \'        ^ S4      ^\'        ^ S4      ,
          .. R.O.4      8X  g   Q hSS^,           ,          \        . R.O. R.O. R.O.4      8X  g   Q h\        \        V3R l4       \        . R.O. R0O. R1O.4      oSR,          \        . R2O. R3O. R4O.4      8X  g   Q hSR,          S^,          8X  g   Q h\        . R-O. R5O. R/O.4      o\!        R4      oSS,          oSP)                  S^4      P+                  4       S^,          8X  g   Q h\        \        VV3R l4       SS,          S^S,          ,          8X  g   Q h\        . R/O. R/O. R/O.4      o\!        RRRR7      oSS,          S8X  g   Q h\!        RRRR7      oSS,          \-        ^ S,          ^ S,          ^ S,          4      8X  g   Q hSS,          P)                  S^ 4      \	        ^4      8X  g   Q hSS,          P)                  S^4      \/        ^4      8X  g   Q h\        . R6O. R7O. R8O.4      oSR,          \-        R9R9R94      8X  g   Q hS\0        ,          \-        ^\0        ,          ^\0        ,          ^\0        ,          4      8X  g   Q h\        . R-O. R.O. R.O.4      o\        \        V3R l4       \        \        V3R l4       \        \        P                  \        P                  .\        P                  \        P                  ..4      oS\        P                  ,          S8X  g   Q h\        ^^.^^..4      oS\        P                  ,          \        \        P                  \        P                  .^\        P                  ,          ^\        P                  ,          ..4      8X  g   Q hR# ):c                  &    \        R4      ^,          # r   r   r5   rm   rt   ru   r   test_power.<locals>.<lambda>#  s    ):rt   r   i g      ?r   c                  8   < S P                  \        ^
4      4      # 
   )_matrix_pow_by_jordan_blocksr   r   s   ru   r   r   L  s     > >qu Ert   c                  T    \        ^^.^^..4      P                  \        R4      4      # )r   )r5   r   r   rm   rt   ru   r   r   O  s(    AA'7 8 U UVWX[V\ ]rt         $@c                     < S R ,          #  @rm   r   s   ru   r   r   R  	    q#vrt   c                  *   < S \        ^^4      ,          # r   r   r   s   ru   r   r   S      q(1a.0rt   ihiO<iv0ir   T)integer)r  negativec                     < S S,          # r{   rm   )r   r   s   ru   r   r   ]  s	    q!trt   )r  nonnegativec                  *   < S \        ^^4      ,          # r  r  r   s   ru   r   r   d  r  rt   g      @c                  D   < S P                  SR4      P                  4       # )r   )subsdoit)Anr   s   ru   r   r   l  s    rwwq"~224rt   )r  positiver   c                     < S R ,          # r   rm   r   s   ru   r   r   {  r   rt   c                     < S \         ,          # r{   r
   r   s   ru   r   r   |  s	    q!trt   Nr   ii)i  i  i,*  i7  )r   r   r   )r   r   r   )r      r   )	i"  i     i  i  ip  i  i  i  r   )r   r   r   r   i   )r   r   r   r   r   )r   r   r   r   r   )r   r   r   r   r   )r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   )r   r   r   )   H   Y   )i#        )i<  i  iI  r  r   r   r   r   r   )r   r   r   r   r   r   gd%%@)rR   r-   r   r5   r9   r   r   Half	sympy.abcr   r_   r`   r   _eval_pow_by_recursionr,   
ValueErrorr   r   r\   rM   r  r  r8   rB   r
   )Rr   r  r   r   s    @@@@ru   
test_powerr%  "  s   
!:;AAA ArE1:!D'!)QsVAXquQw#qAAAAqD!92222	9j12AqD!9EEEEa43q6>>a4199QC5MS $'1c6111q68s1v%%%Aq6Aq6"#WQZ/6Ar7RH:M3NNNNR2r(#$AqvvIq>\)))AQ !AqvvI>QAq6Aq6"#QVV+v1v{6K/LLLLAq6Aq6"#S(FQFS!H3E,FFFFAq6Aq6"#Q&&1ac(QF1C*DDDDAq6Aq6"#Q&&1a41qs8A2FAqD	1R*SSSS	
Aq1q5z!|ad1a4i!Q$q&01QT6:;q!ta!a%jl;!ad;= > > > > Aq!9q!Qi!Q34a76	
Aq1Q3xz1	
AqD!	
Aq!tC <   
 	AA A))!A$/13K3KA3NNNNuAb5FUG9% %)G)G")N %	  $% % % % % 	 "3_oWfghA
;EF ;]^	9i01Ad7fiI>????
:~&
:01AA Ad7fz95	97MNOOOO	9i01Ad7fiI>????	9i01Ad7fiI>????sD!AadF####sD40A
:|$sDd3Aa46	1	~a3nQ6J5J^\]_`Ma5ade5ef~a3QQ_`acdQeMeffh i i i i q1u:Iy ABBBB
:01	9i01Ac6V_oOPPPPc6QT>>	*Y/0AsA	
AB771a=1a4'''
:457a!A#h 	()AsD40Aa4199sDd3Aa441adAqD))))qD;;q!A&&&qD;;q!a((()*Ac6T6662222a44Aq!tQT****	9i01A
:~&
:|$ 166166"234Aaff9>>A1vAaff9!&&!&&!1AaffHah3G HIIIIrt   c                     \        ^\        ,           ^\        ,
          .^\        ,
          ^\        ,           ..4      p \        R4      ;_uu_ 4        V P                  4       P	                  4       \        \        ^\        ,          4      \        ^4      ,           ^,          \        ^\        ,          4      ) \        ^4      ,           ^,          .\        ^\        ,          4      ) \        ^4      ,           ^,          \        ^\        ,          4      \        ^4      ,           ^,          ..4      8X  g   Q h RRR4       R#   + '       g   i     R# ; ir   TN)r5   rc   rF   r   r   r   s    ru   $test_issue_17247_expression_blowup_1r)    s    1ac
QqS!A#J'(A	T		uuw~~61Q3x#a& !#s1Q3xi#a&&8!%;<1Q3xi#a& !#s1Q3x#a&'8!&;<+> $? ? 	? ? 
			s   C#D??E	c                  P   \        ^\        ,           ^\        ,
          .^\        ,
          ^\        ,           ..4      p \        R4      ;_uu_ 4        V P                  4       w  rW,          VP	                  4       ,          '       g   Q h RRR4       R#   + '       g   i     R# ; ir'  )r5   rc   rF   jordan_forminvr   PJs      ru   $test_issue_17247_expression_blowup_2r0    sg    1ac
QqS!A#J'(A	T		}}s1557{{{ 
			s   8BB%	c            
         \        ^\        ,           ^\        ,
          .^\        ,
          ^\        ,           ..4      p \        R4      ;_uu_ 4        V ^d,          \        R\        ^d,          ,          R,           RR\        ^d,          ,          ,
          .RR\        ^d,          ,          ,
          R\        ^d,          ,          R,           ..4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   Tl                N)r5   rc   rF   r(  s    ru   $test_issue_17247_expression_blowup_3r2    s    1ac
QqS!A#J'(A	T		#v+AsF25SSUs  wU  VW  Y\  V\  w\  V\  ]+.LQPSV.SSUstuwztzUz  ~\  V\  ]!^ _ _ 	_ _ 
			s   BC  C1	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V ^
,          \        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)a1  [
        [             -3/4,       45/32 - 37*I/16,         1/4 + I/2,      -129/64 - 9*I/64,      1/4 - 5*I/16,      65/128 + 87*I/64],
        [-149/64 + 49*I/32, -177/128 - 1369*I/128,  125/64 + 87*I/64, -2063/256 + 541*I/128,  85/256 - 33*I/16,  805/128 + 2415*I/512],
        [          1/2 - I,         9/4 + 55*I/16,              -3/4,       45/32 - 37*I/16,         1/4 + I/2,      -129/64 - 9*I/64],
        [   -5/8 - 39*I/16,   2473/256 + 137*I/64, -149/64 + 49*I/32, -177/128 - 1369*I/128,  125/64 + 87*I/64, -2063/256 + 541*I/128],
        [            1 + I,         -19/4 + 5*I/4,           1/2 - I,         9/4 + 55*I/16,              -3/4,       45/32 - 37*I/16],
        [         21/8 + I,    -537/64 + 143*I/16,    -5/8 - 39*I/16,   2473/256 + 137*I/64, -149/64 + 49*I/32, -177/128 - 1369*I/128]]Ta  [
            [     7369525394972778926719607798014571861/604462909807314587353088 - 229284202061790301477392339912557559*I/151115727451828646838272,   -19704281515163975949388435612632058035/1208925819614629174706176 + 14319858347987648723768698170712102887*I/302231454903657293676544,      -3623281909451783042932142262164941211/604462909807314587353088 - 6039240602494288615094338643452320495*I/604462909807314587353088,    109260497799140408739847239685705357695/2417851639229258349412352 - 7427566006564572463236368211555511431*I/2417851639229258349412352, -16095803767674394244695716092817006641/2417851639229258349412352 + 10336681897356760057393429626719177583*I/1208925819614629174706176,    -42207883340488041844332828574359769743/2417851639229258349412352 - 182332262671671273188016400290188468499*I/4835703278458516698824704],
            [50566491050825573392726324995779608259/1208925819614629174706176 - 90047007594468146222002432884052362145*I/2417851639229258349412352,  74273703462900000967697427843983822011/1208925819614629174706176 + 265947522682943571171988741842776095421*I/1208925819614629174706176, -116900341394390200556829767923360888429/2417851639229258349412352 - 53153263356679268823910621474478756845*I/2417851639229258349412352, 195407378023867871243426523048612490249/1208925819614629174706176 - 1242417915995360200584837585002906728929*I/9671406556917033397649408,   -863597594389821970177319682495878193/302231454903657293676544 + 476936100741548328800725360758734300481*I/9671406556917033397649408, -3154451590535653853562472176601754835575/19342813113834066795298816 - 232909875490506237386836489998407329215*I/2417851639229258349412352],
            [   -1715444997702484578716037230949868543/302231454903657293676544 + 5009695651321306866158517287924120777*I/302231454903657293676544,     -30551582497996879620371947949342101301/604462909807314587353088 - 7632518367986526187139161303331519629*I/151115727451828646838272,           312680739924495153190604170938220575/18889465931478580854784 - 108664334509328818765959789219208459*I/75557863725914323419136,    -14693696966703036206178521686918865509/604462909807314587353088 + 72345386220900843930147151999899692401*I/1208925819614629174706176,  -8218872496728882299722894680635296519/1208925819614629174706176 - 16776782833358893712645864791807664983*I/1208925819614629174706176,      143237839169380078671242929143670635137/2417851639229258349412352 + 2883817094806115974748882735218469447*I/2417851639229258349412352],
            [   3087979417831061365023111800749855987/151115727451828646838272 + 34441942370802869368851419102423997089*I/604462909807314587353088, -148309181940158040917731426845476175667/604462909807314587353088 - 263987151804109387844966835369350904919*I/9671406556917033397649408,   50259518594816377378747711930008883165/1208925819614629174706176 - 95713974916869240305450001443767979653*I/2417851639229258349412352,  153466447023875527996457943521467271119/2417851639229258349412352 + 517285524891117105834922278517084871349*I/2417851639229258349412352,  -29184653615412989036678939366291205575/604462909807314587353088 - 27551322282526322041080173287022121083*I/1208925819614629174706176,   196404220110085511863671393922447671649/1208925819614629174706176 - 1204712019400186021982272049902206202145*I/9671406556917033397649408],
            [     -2632581805949645784625606590600098779/151115727451828646838272 - 589957435912868015140272627522612771*I/37778931862957161709568,     26727850893953715274702844733506310247/302231454903657293676544 - 10825791956782128799168209600694020481*I/302231454903657293676544,      -1036348763702366164044671908440791295/151115727451828646838272 + 3188624571414467767868303105288107375*I/151115727451828646838272,     -36814959939970644875593411585393242449/604462909807314587353088 - 18457555789119782404850043842902832647*I/302231454903657293676544,      12454491297984637815063964572803058647/604462909807314587353088 - 340489532842249733975074349495329171*I/302231454903657293676544,      -19547211751145597258386735573258916681/604462909807314587353088 + 87299583775782199663414539883938008933*I/1208925819614629174706176],
            [  -40281994229560039213253423262678393183/604462909807314587353088 - 2939986850065527327299273003299736641*I/604462909807314587353088, 331940684638052085845743020267462794181/2417851639229258349412352 - 284574901963624403933361315517248458969*I/1208925819614629174706176,      6453843623051745485064693628073010961/302231454903657293676544 + 36062454107479732681350914931391590957*I/604462909807314587353088,  -147665869053634695632880753646441962067/604462909807314587353088 - 305987938660447291246597544085345123927*I/9671406556917033397649408,  107821369195275772166593879711259469423/2417851639229258349412352 - 11645185518211204108659001435013326687*I/302231454903657293676544,     64121228424717666402009446088588091619/1208925819614629174706176 + 265557133337095047883844369272389762133*I/1208925819614629174706176]]N)r5   r   rF   r(  s    ru   $test_issue_17247_expression_blowup_4r4    sh    8 	q K L 	MA 
T		"uq "H  I J J 	J J 
			s   #AA'	c            	         \        ^^R 4      p \        R4      ;_uu_ 4        V P                  R4      \        \        ^,          R^\
        ,          ,
          \        ^,          ,          ,           ^$\
        ,          \        ^,          ,          ,           \        RR7      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   c                 D    ^RW,           ,          \         ,          ,           # r   r  r  r}   s   &&ru   r   6test_issue_17247_expression_blowup_5.<locals>.<lambda>      !rQSk!m"3rt   Trc   EXdomainN)r5   rF   charpolyr$   rc   r
   r(  s    ru   $test_issue_17247_expression_blowup_5r?    su    q!34A	T		zz#(1a42!8QT/+ABqDAI+MqY]"^^^^ 
			s   BB00C	c            	         \        ^^\        ^@4       U u. uF  p \        V ,           NK  	  up 4      p\        R4      ;_uu_ 4        VP	                  R4      ^ 8X  g   Q h RRR4       R# u up i   + '       g   i     R# ; i)r   TbareissNr5   r   rc   rF   detr~   r   s     ru   $test_issue_17247_expression_blowup_6rE    sb    q!5":.:aacc:./A	T		uuY1$$$ 
	 /			   A)A..A?	c                      \        ^^R 4      p \        R4      ;_uu_ 4        V P                  R4      ^ 8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   c                 D    ^RW,           ,          \         ,          ,           # r7  r  r}   s   &&ru   r   6test_issue_17247_expression_blowup_7.<locals>.<lambda>  r9  rt   T	berkowitzN)r5   rF   rC  r(  s    ru   $test_issue_17247_expression_blowup_7rK    sB    q!34A	T		uu[!Q&&& 
			s   AA	c            	         \        ^^\        ^@4       U u. uF  p \        V ,           NK  	  up 4      p\        R4      ;_uu_ 4        VP	                  R4      ^ 8X  g   Q h RRR4       R# u up i   + '       g   i     R# ; i)r   TluNrB  rD  s     ru   $test_issue_17247_expression_blowup_8rN    sa    q!5":.:aacc:./A	T		uuT{a 
	 /			rF  c                  H   \        ^^\        ^@4       U u. uF  p \        V ,           NK  	  up 4      p\        R4      ;_uu_ 4        VP	                  4       \        . RO. RO. RO. RO. RO. RO. RO. RO.4      R38X  g   Q h RRR4       R# u up i   + '       g   i     R# ; i)r   TN)r   r   r  r  r   r=  r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r5   r   rc   rF   rrefrD  s     ru   $test_issue_17247_expression_blowup_9rU    s    q!5":.:aacc:./A	T		vvxF********$, - /56 6 	6 6 
	 /			s   B;BB!	c                      \        ^^R 4      p \        R4      ;_uu_ 4        V P                  ^ ^ 4      ^ 8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   c                 D    ^RW,           ,          \         ,          ,           # r7  r  r}   s   &&ru   r   7test_issue_17247_expression_blowup_10.<locals>.<lambda>  r9  rt   TN)r5   rF   cofactorr(  s    ru   %test_issue_17247_expression_blowup_10rZ    sD    q!34A	T		zz!Q1$$$ 
			s   AA	c            	          \        ^^R 4      p \        R4      ;_uu_ 4        V P                  4       \        ^^^ .^$,          4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   c                 D    ^RW,           ,          \         ,          ,           # r7  r  r}   s   &&ru   r   7test_issue_17247_expression_blowup_11.<locals>.<lambda>  r9  rt   TN)r5   rF   cofactor_matrixr(  s    ru   %test_issue_17247_expression_blowup_11r_    sQ    q!34A	T		  "fQA3r6&:::: 
			s   +AA)	c            	          \        ^^R 4      p \        R4      ;_uu_ 4        V P                  4       ^^^\        ,          ^^ ^/8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   c                 D    ^RW,           ,          \         ,          ,           # r7  r  r}   s   &&ru   r   7test_issue_17247_expression_blowup_12.<locals>.<lambda>  r9  rt   TN)r5   rF   	eigenvalsr
   r(  s    ru   %test_issue_17247_expression_blowup_12rd    sP    q!34A	T		{{}AqsAq! 4444 
			s   )AA'	c            
         \        ^ ^\        ,
          \        ^,           ^\        ,
          .^\        ,
          \        ^,           ^ \        ^,           .^ ^\        ,
          \        ^,           ^\        ,
          .^ ^ ^\        ,
          ^ ..4      p V P                  4       pV^ ,          ^ ^\        . RO4      .38X  g   Q hV^,          ^ ,          \        \        ^4      \        ^,
          ,          ,
          ^,           8X  g   Q hV^,          ^,          ^8X  g   Q hV^,          ^,          ^ ,          P	                  RRR7      \        \        ) \        ^4      \        ^,
          ,          ,           ^,
          \        ^,
          ,          .R\        ,          \        ^,          ^\        ,          ,
          ^,           ,          \        ^,           \        \        ^4      \        ^,
          ,          ,
          ^,           ,          \        ^,          ^\        ,          ,
          ^,           ,          ,           .\        ) \        ^4      \        ^,
          ,          ,           ^,
          \        ^,
          ,          .^..4      8X  g   Q hV^,          ^ ,          \        \        ^4      \        ^,
          ,          ,           ^,           8X  g   Q hV^,          ^,          ^8X  g   Q hV^,          ^,          ^ ,          P	                  RRR7      \        \        ) \        ^4      \        ^,
          ,          ,
          ^,
          \        ^,
          ,          .R\        ,          \        ^,          ^\        ,          ,
          ^,           ,          \        ^,           \        \        ^4      \        ^,
          ,          ,           ^,           ,          \        ^,          ^\        ,          ,
          ^,           ,          ,           .\        ) \        ^4      \        ^,
          ,          ,
          ^,
          \        ^,
          ,          .^..4      8X  g   Q hR# )r   FT)deepnumerN)r   r  r   r   rP  )r5   rc   
eigenvectsr   r   )r   evs     ru   %test_issue_17247_expression_blowup_13rj    s   Aq1ua!e$	
QA1a!e$Aq1ua!e$A1q5!$	& 	'A 
Ba5QF=123333a58q47AE?*Q....a58q==a58A;55"tAwA

"QU	+,	Aq!taczA~	!a%!d1gq1uo*=*A!BAqD1Q3JQRN!S	ST"tAwA

"QU	+,	
	A :    a58q47AE?*Q....a58q==a58A;55"tAwA

"QU	+,	Aq!taczA~	!a%!d1gq1uo*=*A!BAqD1Q3JQRN!S	ST"tAwA

"QU	+,	
	A :   rt   c                     \        ^^^\        ,           ^\        ,
          .^,          ^\        ,
          ^\        ,           .^,          ,           ^,          4      p \        R4      ;_uu_ 4        V P                  4       \        \        ^,           ^\        ,
          \        ^,           ^\        ,
          \        ^,           ^\        ,
          \        ^,           ^\        ,
          .^ ^\        ,          ^ ^\        ,          ^ ^\        ,          ^ ^\        ,          .. RO. RO. RO. RO. RO. RO.4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   TNrS  )r5   rc   rF   echelon_formr(  s    ru   %test_issue_17247_expression_blowup_14rm    s    q!qsAaCjlac1Q3Z\1145A	T		~~6UAE1q5!a%QAq1ua!eDacqAaCQ!A#1!DDDDDDD+F $G G 	G G 
			s   -C	EE	c                     \        ^^^\        ,           ^\        ,
          .^,          ^\        ,
          ^\        ,           .^,          ,           ^,          4      p \        R4      ;_uu_ 4        V P                  4       \        \        ^,           ^\        ,
          \        ^,           ^\        ,
          \        ^,           ^\        ,
          \        ^,           ^\        ,
          ..4      \        ^ ^\        ,          ^ ^\        ,          ^ ^\        ,          ^ ^\        ,          ..4      .8X  g   Q h RRR4       R#   + '       g   i     R# ; ir   TN)r5   rc   rF   rowspacer(  s    ru   %test_issue_17247_expression_blowup_15rq    s   q!qsAaCjlac1Q3Z\1145A	T		zz|QAq1ua!eQUAPQESTWXSXZ[^_Z_(`'a bdjmnpqrspsuvxyz{x{}~  AB  CD  AD  FG  IJ  KL  IL  mM  lN  eO   P  P  	P  P 
			s   -CD::E	c                     \        ^^^\        ,           ^\        ,
          .^,          ^\        ,
          ^\        ,           .^,          ,           ^,          4      p \        R4      ;_uu_ 4        V P                  4       \        \        ^,           .^\        ,
          .\        ^,           .^\        ,
          .\        ^,           .^\        ,
          .\        ^,           .^\        ,
          ..4      \        ^\        ,
          .\        ^,           .^\        ,
          .\        ^,           .^\        ,
          .\        ^,           .^\        ,
          .\        ^,           ..4      .8X  g   Q h RRR4       R#   + '       g   i     R# ; iro  )r5   rc   rF   columnspacer(  s    ru   %test_issue_17247_expression_blowup_16rt  !  s   q!qsAaCjlac1Q3Z\1145A	T		}}6AE7AE7AE7AE7APQE7TUXYTYSZ\]`a\a[bdehidicj*k#lntwx{|w|v}  @A  DE  @E  F  HI  LM  HM  GN  PQ  TU  PU  OV  XY  \]  X]  W^  `a  de  `e  _f  hi  lm  hm  gn  pq  tu  pu  ov  vw  ox  #y  y  	y  y 
			s   -C<E44F	c                  :   \        ^^\        ^@4       U u. uF  p \        V ,           NK  	  up 4      p\        R4      ;_uu_ 4        VP	                  4       \        ^.R.^.^ .^ .^ .^ .^ ..4      \        ^.R.^ .^.^ .^ .^ .^ ..4      \        ^.R.^ .^ .^.^ .^ .^ ..4      \        ^.R.^ .^ .^ .^.^ .^ ..4      \        ^.R.^ .^ .^ .^ .^.^ ..4      \        ^.R.^ .^ .^ .^ .^ .^..4      .8X  g   Q h RRR4       R# u up i   + '       g   i     R# ; i)	r   TNr  r   rP  rQ  r=  )r5   r   rc   rF   	nullspacerD  s     ru   %test_issue_17247_expression_blowup_17rx  &  s^   q!5":.:aacc:./A	T		{{}QCaS!aS!aS!56QCaS!aS!aS!56QCaS!aS!aS!56QCaS!aS!aS!56QCaS!aS!aS!56QCaS!aS!aS!56!8 8 	8 8 
	 /			s   DB4D		D	c                  J   \        ^^^\        ,           ^\        ,
          .^,          ^\        ,
          ^\        ,           .^,          ,           ^,          4      p \        R4      ;_uu_ 4        V P                  4       '       d   Q h RRR4       R#   + '       g   i     R# ; i)r   TN)r5   rc   rF   is_nilpotentr(  s    ru   %test_issue_17247_expression_blowup_18r{  1  sj    q!qsAaCjlac1Q3Z\1145A	T		>>##### 
			s   -BB"	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  4       '       d   Q h RRR4       R#   + '       g   i     R# ; i)ay  [
        [             -3/4,                     0,         1/4 + I/2,                     0],
        [                0, -177/128 - 1369*I/128,                 0, -2063/256 + 541*I/128],
        [          1/2 - I,                     0,                 0,                     0],
        [                0,                     0,                 0, -177/128 - 1369*I/128]]TN)r5   r   rF   is_diagonalizabler(  s    ru   %test_issue_17247_expression_blowup_19r~  6  sO    q a b 	cA
 
T		&&((((( 
			s   AA	c                  h   \        \        ^,           ^\        ,
          ^ ^ .^\        ,
          \        ^,           ^ \        ^,           .^ ^\        ,
          \        ^,           ^ .^ ^ ^ \        ^,           ..4      p \        R4      ;_uu_ 4        V P                  4       \        ^^^ \        ^,           \        ^,
          ,          .. RO. RO. RO.4      \        . RO^ ^\        ,          ^ ^ .^ ^ \        ^,           ^ .^ ^ ^ \        ^,           ..4      38X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   TN)r   r  r   r   r   r   r   r   r   r   r   r   r   )r5   rc   rF   diagonalizer(  s    ru   %test_issue_17247_expression_blowup_20r  ?  s	   UQU#UQUQU#	
QUQU#	
QU#	% 	&A
 
T		}}6AAA''''	+) $*
 "!"!QU"!AE"	$ %	#& 	& 		& 	& 
			s   B
D  D1	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  RR7      \        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)  [
        [             -3/4,       45/32 - 37*I/16,                   0,                     0],
        [-149/64 + 49*I/32, -177/128 - 1369*I/128,                   0, -2063/256 + 541*I/128],
        [                0,         9/4 + 55*I/16, 2473/256 + 137*I/64,                     0],
        [                0,                     0,                   0, -177/128 - 1369*I/128]]TGEmethodv  [
            [-26194832/3470993 - 31733264*I/3470993, 156352/3470993 + 10325632*I/3470993, 0, -7741283181072/3306971225785 + 2999007604624*I/3306971225785],
            [4408224/3470993 - 9675328*I/3470993, -2422272/3470993 + 1523712*I/3470993, 0, -1824666489984/3306971225785 - 1401091949952*I/3306971225785],
            [-26406945676288/22270005630769 + 10245925485056*I/22270005630769, 7453523312640/22270005630769 + 1601616519168*I/22270005630769, 633088/6416033 - 140288*I/6416033, 872209227109521408/21217636514687010905 + 6066405081802389504*I/21217636514687010905],
            [0, 0, 0, -11328/952745 + 87616*I/952745]]Nr5   r   rF   r,  r(  s    ru   %test_issue_17247_expression_blowup_21r  Q  g    q c d 	eA
 
T		uuDu!VA /: -; &< < 	< < 
			   ,AA0	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  RR7      \        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r  TLUr  r  Nr  r(  s    ru   %test_issue_17247_expression_blowup_22r  ^  r  r  c                     \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  RR7      P	                  4       \        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r  TADJr  r  N)r5   r   rF   r,  r   r(  s    ru   %test_issue_17247_expression_blowup_23r  k  sp    q c d 	eA
 
T		uuEu"))+va 9: 7; 0< < 	< < 
			s   :A--A>	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  RR7      \	        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r  TCHr  r  Nr6   r   rF   r,  r5   r(  s    ru   %test_issue_17247_expression_blowup_24r  x  sg    Q c d 	eA
 
T		uuDu!VA /: -; &< < 	< < 
			r  c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  RR7      \	        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r  TLDLr  r  Nr  r(  s    ru   %test_issue_17247_expression_blowup_25r    sg    Q c d 	eA
 
T		uuEu"fQ 0: .; '< < 	< < 
			r  c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  4       ^8X  g   Q h RRR4       R#   + '       g   i     R# ; i)a  [
        [             -3/4,       45/32 - 37*I/16,         1/4 + I/2,      -129/64 - 9*I/64,      1/4 - 5*I/16,      65/128 + 87*I/64,         -9/32 - I/16,      183/256 - 97*I/128],
        [-149/64 + 49*I/32, -177/128 - 1369*I/128,  125/64 + 87*I/64, -2063/256 + 541*I/128,  85/256 - 33*I/16,  805/128 + 2415*I/512, -219/128 + 115*I/256, 6301/4096 - 6609*I/1024],
        [          1/2 - I,         9/4 + 55*I/16,              -3/4,       45/32 - 37*I/16,         1/4 + I/2,      -129/64 - 9*I/64,         1/4 - 5*I/16,        65/128 + 87*I/64],
        [   -5/8 - 39*I/16,   2473/256 + 137*I/64, -149/64 + 49*I/32, -177/128 - 1369*I/128,  125/64 + 87*I/64, -2063/256 + 541*I/128,     85/256 - 33*I/16,    805/128 + 2415*I/512],
        [            1 + I,         -19/4 + 5*I/4,           1/2 - I,         9/4 + 55*I/16,              -3/4,       45/32 - 37*I/16,            1/4 + I/2,        -129/64 - 9*I/64],
        [         21/8 + I,    -537/64 + 143*I/16,    -5/8 - 39*I/16,   2473/256 + 137*I/64, -149/64 + 49*I/32, -177/128 - 1369*I/128,     125/64 + 87*I/64,   -2063/256 + 541*I/128],
        [               -2,         17/4 - 13*I/2,             1 + I,         -19/4 + 5*I/4,           1/2 - I,         9/4 + 55*I/16,                 -3/4,         45/32 - 37*I/16],
        [     1/4 + 13*I/4,    -825/64 - 147*I/32,          21/8 + I,    -537/64 + 143*I/16,    -5/8 - 39*I/16,   2473/256 + 137*I/64,    -149/64 + 49*I/32,   -177/128 - 1369*I/128]]TNr5   r   rF   rankr(  s    ru   %test_issue_17247_expression_blowup_26r    sI    q z { 	|A 
T		vvx1}} 
			   AA	c            	      J   \        ^ ^\        ,
          \        ^,           ^\        ,
          .^\        ,
          \        ^,           ^ \        ^,           .^ ^\        ,
          \        ^,           ^\        ,
          .^ ^ ^\        ,
          ^ ..4      p \        R4      ;_uu_ 4        V P                  4       w  rVP	                  4       \        \        R4      4      P	                  4       8X  g   Q hV\        \        R4      4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   Ta  [
            [    0,  4*x/(x**2 - 2*x + 1), -(-17*x**4 + 12*sqrt(2)*x**4 - 4*sqrt(2)*x**3 + 6*x**3 - 6*x - 4*sqrt(2)*x + 12*sqrt(2) + 17)/(-7*x**4 + 5*sqrt(2)*x**4 - 6*sqrt(2)*x**3 + 8*x**3 - 2*x**2 + 8*x + 6*sqrt(2)*x - 5*sqrt(2) - 7), -(12*sqrt(2)*x**4 + 17*x**4 - 6*x**3 - 4*sqrt(2)*x**3 - 4*sqrt(2)*x + 6*x - 17 + 12*sqrt(2))/(7*x**4 + 5*sqrt(2)*x**4 - 6*sqrt(2)*x**3 - 8*x**3 + 2*x**2 - 8*x + 6*sqrt(2)*x - 5*sqrt(2) + 7)],
            [x - 1, x/(x - 1) + 1/(x - 1),                       (-7*x**3 + 5*sqrt(2)*x**3 - x**2 + sqrt(2)*x**2 - sqrt(2)*x - x - 5*sqrt(2) - 7)/(-3*x**3 + 2*sqrt(2)*x**3 - 2*sqrt(2)*x**2 + 3*x**2 + 2*sqrt(2)*x + 3*x - 3 - 2*sqrt(2)),                       (7*x**3 + 5*sqrt(2)*x**3 + x**2 + sqrt(2)*x**2 - sqrt(2)*x + x - 5*sqrt(2) + 7)/(2*sqrt(2)*x**3 + 3*x**3 - 3*x**2 - 2*sqrt(2)*x**2 - 3*x + 2*sqrt(2)*x - 2*sqrt(2) + 3)],
            [    0,                     1,                                                                                            -(-3*x**2 + 2*sqrt(2)*x**2 + 2*x - 3 - 2*sqrt(2))/(-x**2 + sqrt(2)*x**2 - 2*sqrt(2)*x + 1 + sqrt(2)),                                                                                            -(2*sqrt(2)*x**2 + 3*x**2 - 2*x - 2*sqrt(2) + 3)/(x**2 + sqrt(2)*x**2 - 2*sqrt(2)*x - 1 + sqrt(2))],
            [1 - x,                     0,                                                                                                                                                                                               1,                                                                                                                                                                                             1]]a  [
            [0, 1,                       0,                       0],
            [0, 0,                       0,                       0],
            [0, 0, x - sqrt(2)*(x - 1) + 1,                       0],
            [0, 0,                       0, x + sqrt(2)*(x - 1) + 1]]N)r5   rc   rF   r+  r   r   r-  s      ru   %test_issue_17247_expression_blowup_27r    s   Aq1ua!e$	
QA1a!e$Aq1ua!e$A1q5!$	& 	'A
 
T		}}xxzVA 'o %p q rx  rx  rz	z 	z z
 F1 I J K K 	K K 
			s   !A%DD"	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  4       \        R4      8X  g   Q h RRR4       R#   + '       g   i     R# ; i)r  Ta  [
            sqrt(14609315/131072 + sqrt(64789115132571/2147483648 - 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3) + 76627253330829751075/(35184372088832*sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))) - 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)))/2 + sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))/2),
            sqrt(14609315/131072 - sqrt(64789115132571/2147483648 - 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3) + 76627253330829751075/(35184372088832*sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))) - 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)))/2 + sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))/2),
            sqrt(14609315/131072 - sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))/2 + sqrt(64789115132571/2147483648 - 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3) - 76627253330829751075/(35184372088832*sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))) - 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)))/2),
            sqrt(14609315/131072 - sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))/2 - sqrt(64789115132571/2147483648 - 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3) - 76627253330829751075/(35184372088832*sqrt(64789115132571/4294967296 + 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)) + 2*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3))) - 3546944054712886603889144627/(110680464442257309696*(25895222463957462655758224991455280215303/633825300114114700748351602688 + sqrt(1213909058710955930446995195883114969038524625997915131236390724543989220134670)*I/22282920707136844948184236032)**(1/3)))/2)]N)r5   r   rF   singular_valuesr(  s    ru   %test_issue_17247_expression_blowup_28r    sc    q c d 	eA
 
T		  "a )c 'd d 	d d 
			s   !AA%	c                      \        \        R 4      4      p \        R4      ;_uu_ 4        V P                  4       ^8X  g   Q h RRR4       R#   + '       g   i     R# ; i)a	  [
        [1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I,1/4+1/2*I,-129/64-9/64*I,1/4-5/16*I,65/128+87/64*I,-9/32-1/16*I,183/256-97/128*I,3/64+13/64*I,-23/32-59/256*I,15/128-3/32*I,19/256+551/1024*I],
        [21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I,125/64+87/64*I,-2063/256+541/128*I,85/256-33/16*I,805/128+2415/512*I,-219/128+115/256*I,6301/4096-6609/1024*I,119/128+143/128*I,-10879/2048+4343/4096*I,129/256-549/512*I,42533/16384+29103/8192*I],
        [-2,17/4-13/2*I,1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I,1/4+1/2*I,-129/64-9/64*I,1/4-5/16*I,65/128+87/64*I,-9/32-1/16*I,183/256-97/128*I,3/64+13/64*I,-23/32-59/256*I],
        [1/4+13/4*I,-825/64-147/32*I,21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I,125/64+87/64*I,-2063/256+541/128*I,85/256-33/16*I,805/128+2415/512*I,-219/128+115/256*I,6301/4096-6609/1024*I,119/128+143/128*I,-10879/2048+4343/4096*I],
        [-4*I,27/2+6*I,-2,17/4-13/2*I,1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I,1/4+1/2*I,-129/64-9/64*I,1/4-5/16*I,65/128+87/64*I,-9/32-1/16*I,183/256-97/128*I],
        [1/4+5/2*I,-23/8-57/16*I,1/4+13/4*I,-825/64-147/32*I,21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I,125/64+87/64*I,-2063/256+541/128*I,85/256-33/16*I,805/128+2415/512*I,-219/128+115/256*I,6301/4096-6609/1024*I],
        [-4,9-5*I,-4*I,27/2+6*I,-2,17/4-13/2*I,1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I,1/4+1/2*I,-129/64-9/64*I,1/4-5/16*I,65/128+87/64*I],
        [-2*I,119/8+29/4*I,1/4+5/2*I,-23/8-57/16*I,1/4+13/4*I,-825/64-147/32*I,21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I,125/64+87/64*I,-2063/256+541/128*I,85/256-33/16*I,805/128+2415/512*I],
        [0,-6,-4,9-5*I,-4*I,27/2+6*I,-2,17/4-13/2*I,1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I,1/4+1/2*I,-129/64-9/64*I],
        [1,-9/4+3*I,-2*I,119/8+29/4*I,1/4+5/2*I,-23/8-57/16*I,1/4+13/4*I,-825/64-147/32*I,21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I,125/64+87/64*I,-2063/256+541/128*I],
        [0,-4*I,0,-6,-4,9-5*I,-4*I,27/2+6*I,-2,17/4-13/2*I,1+I,-19/4+5/4*I,1/2-I,9/4+55/16*I,-3/4,45/32-37/16*I],
        [0,1/4+1/2*I,1,-9/4+3*I,-2*I,119/8+29/4*I,1/4+5/2*I,-23/8-57/16*I,1/4+13/4*I,-825/64-147/32*I,21/8+I,-537/64+143/16*I,-5/8-39/16*I,2473/256+137/64*I,-149/64+49/32*I,-177/128-1369/128*I]]TNr  r(  s    ru   test_issue_16823r    sI    q F G 	HA 
T		vvx1}} 
			r  c                  <   \        . RO^\        ^4      ,           R\        ^4      ,           ^\        ^4      ,
          \        ^4      ) ^,
          ^^R^^.	R^\        ^4      ,          ,           R^\        ^4      ,          ,
          R^\        ^4      ,          ,
          R^\        ^4      ,          ,           R^RR^ .	R\        ^4      ,          ^,
          ^^\        ^4      ,          ,
          R^\        ^4      ,          ,           ^^\        ^4      ,          ,           RR^R^.	^^\        ^4      ,          ,
          ^\        ^4      ,          ^,           ^\        ^4      ,          ^,           ^^\        ^4      ,          ,
          ^R	^^^ .	R^\        ^4      ,          ,           ^^\        ^4      ,          ,           R\        ^4      ,          ^,
          ^^\        ^4      ,          ,
          ^RRR^.	R^\        ^4      ,          ,           R^\        ^4      ,          ,
          R^\        ^4      ,          ,
          R^\        ^4      ,          ,           R^RR^ .	^\        ^4      ,
          \        ^4      ) ^,
          ^\        ^4      ,           R\        ^4      ,           R^^^^.	.4      p \        R4      ;_uu_ 4        V P                  4       \        ^^ ^ ^ ^ ^ ^ ^ \	        ^4      ^,          .	^ ^^ ^ ^ ^ ^ ^ \	        ^4      ) ^,          .	^ ^ ^^ ^ ^ ^ ^ \	        ^4      ^,          .	^ ^ ^ ^^ ^ ^ ^ \	        ^4      ) ^,          .	. R
O^ ^ ^ ^ ^ ^^ ^ \	        ^4      ) ^,          .	. RO^ ^ ^ ^ ^ ^ ^ ^\	        ^4      ) ^,          .	.4      R38X  g   Q h RRR4       R#   + '       g   i     R# ; i)r   TN)	r   r   r   r   r   r   r   r   r   r  rQ  rv  r  r   )	r   r   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   r   r   rR  )r5   r   rF   rT  r   r(  s    ru   test_issue_18531r    s   #	
T!Wb47lAQK$q'Aq!RAN	aQiaQiaQiaQiQPRTVXYZ	DGaQtAwYQtAwYAd1gIr2qRTVWX	
QtAwY$q'	Aqay1}a!DG)mQQPRTUV	aQiQtAwY47
QAd1gIq"bRTVWX	aQiaQiaQiaQiQPRTVXYZ	
T!WtAwhlAQKd1gr1aAN	
 		A 
T		vvxF1aAq!ad1f-1aAq!adU1W-1aAq!ad1f-1aAq!adU1W-*1aAq!adU1W-*1aAq!adU1W-$/ 0 2JK K 	K K 
			s   6C	N

N	c            	         \        \        R  4       \        \        R 4       \        \        R 4       \        \        4      ;_uu_ 4        ^\        R
4      ^&   RRR4       \        4       \        . 4      u;8X  d   \        ^ ^ . 4      8X  g   Q h Q h\        . .4      \        ^^ . 4      8X  g   Q h\        . . .4      \        ^^ . 4      8X  g   Q h\	        4       ;_uu_ 4        \        ^.R..4      P                  4       ^.R..8X  g   Q h RRR4       \	        4       ;_uu_ 4        \        ^.R..4      P                  P                  4       ^..R..8X  g   Q h RRR4       \        ^ ..4      p \	        4       ;_uu_ 4        \        P                  V R&   RRR4       \        \        ^ .^ ^ ..4      pTpVP                  VP                  8X  g   Q hVP                  ^8X  g   Q hVR,          \        ^ ^ ^ .8X  g   Q h\        ^^\        ^ ^ ^ .4      pTpVP                  VP                  8X  g   Q hVP                  ^8X  g   Q hVR,          \        ^ ^ ^ .8X  g   Q hW8X  g   Q h\        V4      V8X  g   Q h\        ^^\        ^^4      4      p\        ^^\        ^^
4      4      p\        WE.4      pVP                  ^8X  g   Q hVP                  ^8X  g   Q hVR,          . RO8X  g   Q h\        \        ^4      4      \        ^4      8X  g   Q h\        \        \        ^4      4      4      \        \        ^4      4      8X  g   Q h\        V4      VP                  ! 4       8X  g   Q h\        \        V4      4      \        V4      P!                  4       8X  g   Q hV\        V4      Jg   Q h\#        ^^4      \#        ^^4      ^,          .\#        ^^4      ^,          \#        ^^4      ^,          ..p\        V4      p V \        . RO. RO. RO. RO. RO.4      8X  g   Q hV P                  4       V8w  g   Q h\        VRR7      P                  4       V8X  g   Q h\%        R^^4      p\#        ^4      V.p\        V4      \        ^^.^^.VR,          VR,          .VR,          VR,          ..4      8X  g   Q h\	        4       ;_uu_ 4        \        VRR7      P                  4       V U	u. uF  q.NK  	  up	8X  g   Q h RRR4       \        \#        ^4      \#        ^ 4      .4      \        \#        ^4      .4      8X  g   Q h\        \        R 4       \        \        R	 4       \        ^^.^^..4      p \        V R.4      p
V
\        ^^.^^.^^..4      8X  g   Q hR#   + '       g   i     EL?; i  + '       g   i     EL; i  + '       g   i     ELg; i  + '       g   i     EL@; iu up	i   + '       g   i     L; i)c                  .    \        ^^\        ^4      4      # r   )r5   r   rm   rt   ru   r   test_creation.<locals>.<lambda>  s    vaE"I6rt   c                      \        ^R. 4      # )r   r  r   rm   rt   ru   r   r    s    vaR0rt   c                  &    \        R4      ^,          # r   r   rm   rt   ru   r   r    s    vf~a0rt   Nr   Fevaluater   c                  B    \        \        ^4      \        ^ ^4      .4      # r   r5   r<   rm   rt   ru   r   r  4      vtAwQ
&;<rt   c                  B    \        \        ^4      \        ^^ 4      .4      # r  r  rm   rt   ru   r   r  5  r  rt   r   r  r   	r   r   r   r   r   r   r   r   	   )r   r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   )rR   r#  
IndexErrorr5   rW   tolistTr   EmptySetrc   r   r   r   r9   r3   as_immutable
as_mutabler<   rE   )r   r_   r   r`   c23c13ra   datr   r~   M2s              ru   test_creationr    s2   
:67
:01
:01	
		vq 
 8vbz5VAq"%5555552$<6!Q++++2r(vaB//// 
 	!	!T{m$++-A3+>>> 
"		!	!T{m$&&--/aSED6?BBB 
"uA		!	!**$ 
" 	AA A	A66QVV66Q;;Q4Aq!Q<q!aAq\"A	A66QVV66Q;;Q4Aq!Q<6M6!9>>
AuQ{
#C
AuQ|
$CzA66Q;;66Q;;Q4....#a&>SV###?3q623s1v7NNNN1!1111/!$%);)F)F)HHHHF1I1ItAay{#d1Qik4!9Q;%?
@CsA     88:#&--/3666S!QA7A,C#;&AAtWagtWag	"    
 
 	!	!cE*113S7ISS7IIII 
" 47DG$%a	)::::
:<=
:<= 	AA A	F	B!Q!Q!Q01111W 
		 
"	!	!	!	!	! 
"	!	!h 8J 
"	!sN   V!)V$%4V8W>W%W 
(	W%V!	$V5	8W		W	 W%%W5	c                  X   \         P                  ! ^\        ^^4      \        ^^4      ^,          \        ^^4      ^,          \        ^^4      ^,          \        ^^4      ^,          \        ^^4      ^,          \        ^^4      ^,          4      \        . RO. RO. RO. RO.4      8X  g   Q hR# )r   N)r   r   r   r   r   r   )r   r   r   r   r   r   )r   r   r   r   r   r   )r5   	irregularr<   rm   rt   ru   test_irregular_blockr  =  s    AtAay$q)A+tAay{Qq	!T!AYq[$q)A+tAay{<?E	G @  rt   c                  P   \         P                  \         P                  \        \        ,          \         P
                  .\        \        \        \        ^,          .\        \         P                  ) \        \        ,          ^..p \        V 4      pVP                  4       V 8X  g   Q hR# r   N)	r   Oner   rc   rd   Zerore   r5   r  )lstr   s     ru   test_tolistr  F  sc    EE1661Q3'!Q1a41quufac1:M
NCsA88:rt   c                  ,   \        ^ ^4      P                  4       \        ^ ^4      8X  g   Q h\        ^ ^4      P                  4       \        \        ^ ^4      4      8X  g   Q h\        ^^ 4      P                  4       \        \        ^^ 4      4      8X  g   Q hR# r   N)rB   r  r  r3   rm   rt   ru   test_as_mutabler  L  su    A;!!#uQ{222A;##%q!)EEEEA;##%q!)EEEErt   c                     \        ^4      p V R,          \        ^4      8X  g   Q hV R,          \        ^4      8X  g   Q h\        ^^R 4      pVR,          \        ^^R4      8X  g   Q hVR,          \        ^^R	4      8X  g   Q h\        . R
O. RO. RO. RO.4      pVR,          \        ^^. RO4      8X  g   Q hVRR1R3,          \        . RO. RO.4      8X  g   Q hR# )r   c                     W,           # r{   rm   r}   s   &&ru   r   test_slicing.<locals>.<lambda>W      15rt   r   N)Nr   Nr  )r   r   Nr   r   r   r   ):r   r   Nr   r   r   )r   r   r   r   r   r   r   r   r   r  r      r           r   r  )r   r   r  r  r  )r9   rB   r5   )m0m1m2s      ru   test_slicingr  R  s    	QBf:Qh<58###	1(	)Bd8vaI....f:1f----	|^=MN	OBe9q!^4444bc1f:1A BCCCCrt   c                     \        ^4      p \        ^4      V R&   V \        R4      8X  g   Q h\        ^4      V R&   V \        ^4      8X  g   Q h\        ^^R4      V R&   V \        R4      8X  g   Q h\        ^4      V R&   V \        ^4      8X  g   Q h. R	OV R&   V \        R	4      8X  g   Q h^ ^ .V R
&   V \        R4      8X  g   Q hR# )r   N)r  r  )r   r   r   r   r  r   r  )Nr   Nr  r   r   r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   r   r  r   r   r   r   r  r   r  r  r  r  r     )r  r   ))r   r   r   r   )r   r   r   r   r  r  )rB   r9   r5   r   s    ru   test_submatrix_assignmentr  `  s    aAa&AhK & ' ' ' ' AAfIA;;Q<(AdG & ' ' ' ' AhAdGa==MAdG * + + + + 1vAeH * + + + +rt   c                  r  a  \        ^^R 4      o S P                  . RO^ ^.4      \        ^^. RO4      8X  g   Q hS P                  ^ ^.. RO4      \        ^^. RO4      8X  g   Q hS P                  \        ^4      \        ^4      4      S 8X  g   Q h\        \        V 3R l4       \        \        V 3R l4       R# )	r   c                 "    V ^,          V,           # r  rm   r}   s   &&ru   r   test_extract.<locals>.<lambda>}  s    !A#'rt   c                  ,   < S P                  ^.^ .4      # )r   r   r   s   ru   r   r        qyy!qc2rt   c                  ,   < S P                  ^ .^.4      # r   r  r   s   ru   r   r    r  rt   N)r   r   r   )r   r   r   r   r  r   r  )r   r   r   r  r  r  )r5   r   r   rR   r  r   s   @ru   test_extractr  |  s    q!)*A99YA'6!Q8K+LLLL99aVY'6!Q8K+LLLL99U1XuQx(A---
:23
:23rt   c                     \        ^4      p V P                  ^^	4      \        ^^	R4      8X  g   Q h\        ^^R 4      pVP                  ^^4      \        R4      8X  g   Q hVP                  ^^4      \        R4      8X  g   Q hR# )r   c                     W,           # r{   rm   r}   s   &&ru   r   test_reshape.<locals>.<lambda>  r  rt   N)	r   r   r   r   r   r   r   r   r   )r   r   r   r   )r   r   r   r   ))r   r   r   r   r   r   )r   r   r   r   r   r   )r9   reshaper5   r  r  s     ru   test_reshaper    s    	QB::ava,GHHHH	1(	)B::	1DEF F F::av&NOOOOrt   c                      \        ^4      p V P                  R 4      \        ^4      ^,          8X  g   Q hV P                  R 4      \        ^4      8X  g   Q hR# )r   c                     ^V ,          # r  rm   rc   s   &ru   r    test_applyfunc.<locals>.<lambda>  s    !A#rt   c                     ^ # r   rm   r  s   &ru   r   r    s    !rt   N)r9   	applyfuncrB   )r  s    ru   test_applyfuncr    sE    	QB<<&#a&(222<<$a000rt   c            	         \        \        \        \        ,           ,          ^.\        \        ,           \        ,          \        ,          \        \        \        \        \        ,           ,          ,           ,          ..4      p V P                  4       pV\        \        \        ,          \        ^,          ,           ^.\        \        ^,          ,          \        \        ^,          ,          ,           \        \        ,          \        \        ^,          ,          ,           \        ^,          ,           ..4      8X  g   Q h\	        RRR7      p\        \        \        V,          4      .4      P                  RR7      \        \        V4      \        \        V4      ,          ,           .4      8X  g   Q h\        . RO. RO. RO.4      P                  4       \        ^^\        ^^4      .. R	O. R
O.4      8X  g   Q hR# )r   r_   Treal)complexNr   )r   r   r  r   )r   r   r  r  )
r5   rc   rd   r   r   r   r
   r   r   r   )r  r  r_   s      ru   test_expandr
    sW   	!QU)QAE19a-A1q5	M1B!CD	EB	B
A#1*a1QT6AadF?AaC!AqD&L1a4,?@AC C C C 	sA3qs8*$$T$2A3q6!"#$ $ $ 9j)4599;v	
Ax1~G @   rt   c            	         \        \        \        4      ^,          \        \        ^,          4      .\        \        ^,          4      \        \        4      ^,          ,          \        \        ^,          4      \        \        4      ^,          ,          ..4      p V P                  \        P                  ! \        4      \        P                  ! \        4      ,          4      pV\        \        ^,          \        \        4      .\        ^,          \        \        4      ,          \        ^,          \        \        4      ,          ..4      8X  g   Q hV P                  \        P                  ! \        4      \        P                  ! \        4      ,          4      pV\        \        ^,          \        .\        \        ^,          ,          \        ^,          \        ,          ..4      8X  g   Q hV P                  \        P                  ! \        4      \        P                  ! \        4      ,          4      pV\        \        ^,          \        ) .\        ) \        ^,          ,          \        ^,          ) \        ,          ..4      8X  g   Q hR# r  )
r5   r   rc   r   rd   refinerY   r  r  r  r  s     ru   test_refiner    st   	#a&!)T!Q$Z(adCFAI%tAqDz#a&!)';<> 
?B	166!9qvvay(	)B!Q$A!Q$s1v+q!tCF{)CDEEEE	1::a=1::a=0	1B!Q$Qq!tVQT!V$456666	1::a=1::a=0	1B!Q$qbAg1uQw%789999rt   c            	         \        ^^4      p \        ^^^R7      p V \        ^^^R7      8X  g   Q h\        ^^^ ^4      p \        ^^RR7      p V \        ^^RR7      8X  g   Q hV P                  4       pVP                  ! 4        W8X  g   Q h\        P                  ! ^4      pV \        ^RVR7      8X  g   Q hR F,  pR	 F#  p\        VRWBR7      p W P
                  8X  d   K#  Q h	  K.  	  \        ^
^^FR7      p ^ p\        V P                  ^ ,          4       F@  p\        V P                  ^,          4       F  pWV3,          ^ 8X  g   K  V^,          pK  	  KB  	  V^8X  g   Q hR# )
r   )seedT)r  	symmetric)r  prng)r  percentr  )minr  N)r   r  )d   F      )r=   copyr)   randomRandomr  r   r   )r   r   rngsizer  
zero_countr~   r   s           ru   test_randomr    sD   1aA1aa A
1aa((((1aC A11-A
115555	AJJL6M6
--
C
137777 $G44KA8O8 % 
 	21b)AJ1771:qwwqz"AAw!|a
 #  rt   c                  "  aa \        ^4      p V P                  4       \        ^4      8X  g   Q hV P                  RR7      \        ^4      8X  g   Q hV P                  RR7      \        ^4      8X  g   Q hV P                  RR7      \        ^4      8X  g   Q hV P                  RR7      \        ^4      8X  g   Q hV P                  RR7      \        ^4      8X  g   Q h\        . RO. RO. RO.4      p V P                  4       pW,          \        ^4      8X  g   Q hV P                  RR7      V8X  g   Q hV P                  RR7      V8X  g   Q hV P                  RR7      V8X  g   Q hV P                  RR7      V8X  g   Q hV P                  RR7      V8X  g   Q h\        . RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. RO. R O. R!O. R"O. R#O. R$O. R%O. R&O. R'O. R(O. R)O.4      pVP                  RR7      V,          \        VP                  ^ ,          4      8X  g   Q h\        oS! . R*O. R+O. R,O.4      o\
        ;QJ d.    VV3R lR	P                  4        4       F  '       d   K   R
M#	  RM! VV3R lR	P                  4        4       4      '       g   Q h\        oS! . R*O. R+O. R,O.4      o\
        ;QJ d.    VV3R lR	P                  4        4       F  '       d   K   R
M#	  RM! VV3R lR	P                  4        4       4      '       g   Q hR# )-r   r  r  r  r  r  QRBLOCKc              3   \   <"   T F!  p\        SP                  V4      4      SJ x  K#  	  R # 5ir{   r   r,  r   sr   r   s   & ru   r   test_inverse.<locals>.<genexpr>  $     L.KtAEE!H~$.K   ),zGE ADJ LU CH LDL QRFTc              3   \   <"   T F!  p\        SP                  V4      4      SJ x  K#  	  R # 5ir{   r"  r#  s   & ru   r   r%    r&  r'  N)r   r   r   )r   r   r   )r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )0   1      )r  G   ^   );      A   )r9   r,  r5   r   r3   r   splitr4   )r   AinvAAr   r   s      @@ru   test_inverser4    s   AA557c!f555Q'''555#a&(((555Q'''555#a&(((555Q'''	 	A 557D6SV555%%%555$&&&555%%%555$&&&555%%%	\WWWWWWWWWWWWWWWWWWWWWWWW1Y 
ZB2 666!B&#bhhqk*::::
C\ 	A 3L.C.I.I.KL333L.C.I.I.KLLLLL
C\ 	A 3L.C.I.I.KL333L.C.I.I.KLLLLLrt   c                     \        ^^\        ^,          \        ,          ^\        ^,          ,          \        \        ,          ,           .4      p \        \        .pV P                  V4      \        ^\        ,          \        ,          \        ^,          .\        ^\        ,          \        ,           ..4      8X  g   Q h\        ^^\        \        ^,          \        ^,          ,          .4      p V P                  V4      \        ^^ .^\        ,          \        ^,          ,          \        ^,          ^,          \        ^,          ,          ..4      8X  g   Q h\        ^,          \        ,          p\        \        .p\	        W!4      \        ^\        ,          ^\        ,          .^\        ,          ^ ..4      8X  g   Q h\        ^,          \        ^,          ,          p\	        W!4      \        ^\        ^,          ,          ^\        ,          \        ^,          ,          .^\        ,          \        ^,          ,          ^\        ^,          ,          \        ,          ..4      8X  g   Q h\
        \        \        ^,          ,          ,           p\        ^,          ^\        ^,          ,          ,           p\        ^ ^\        ,          .^\        ,          ^\        ,          ..4      pV\	        V\        \        \        .4      4      8X  g   Q hV\	        V\        \        \        .4      P                  4      8X  g   Q h\	        V\        \        3V.4      \        ^ ^\        ^,          ,          ^\        ,          .^\        ^,          ,          ^\        ,          ^\        ,          .^\        ,          ^\        ,          ^ ..4      8X  g   Q hR# r   )r5   rc   rd   jacobianr:   re   r  )Lsymsfganss        ru   test_jacobian_hessianr<  
  s   q!ad1fa1fqsl+,Aq6D::dv!Aq!t}q!A#'l&CDDDDq!aAad^$A::dv1v!AqD!Q$q&A+/F&GHHHH	1QAq6D1v!QqSzAaC8&<====	1QT	A11a41QT"QqSAXqAvax$89:; ; ; 	
AadF
A	1qAvA
1!*Q3!* C'!VQF^,,,,'!VQF^--....1q!fqc"f1a41	
1a4AaC1aCAaC1. '       rt   c                     \        \        \        4      \        \        4      .\        4      \        \        4      ^,          \        \        4      ^,          ,           8X  g   Q h\        \	        \        4      \	        ^\        ,          4      .\        4      \	        ^\        ,          4      8X  g   Q h\        \	        \        4      \        .\        4      \	        \        4      \        \	        \        4      ,          ,
          8X  g   Q h\        ^\        \        ^,          .\        4      ^8X  g   Q hR\	        \        4      ,          \        \        4      ,          \        ,          ^\        \        4      ,          \	        \        4      ,          \        ^,          ,          ,           ^\	        \        4      ,          \        \        4      ,          \        ,          ,
          \	        \        4      \        \        4      ,          \        ^,          ,          ,
          \	        \        4      \        \        4      ,          \        ^,          ,          ,           p \        \	        \        4      \        \        4      \        ^,          .\        4      P                  4       V 8X  g   Q h\        \	        \        4      \        \        4      \        ^,          .\        RR7      P                  4       V 8X  g   Q h\        ^,          ) \        \        4      ^,          ,          \        ^,          \        \        4      ^,          ,          ,
          ^\        ,          \        \        4      ^,          ,          ,
          ^\        ,          \        \        4      ^,          ,          ,
          p\        \        \        4      \        \        4      \        ^,          .\        4      P                  4       V8X  g   Q h\        \        \        4      \        \        4      \        ^,          .\        RR7      P                  4       V8X  g   Q h\        . \        4      ^8X  g   Q hR# )r   rJ  r  Nr=  )rA   r   rc   r   r   r   )w1w2s     ru   test_wronskianr@  &  s   c!fc!f%q)SVQYQ-BBBBc!fc!A#h'+s1Q3x777c!fa[!$A3q6(9999aAqD\1%***	CF3q6	!	aAhs1voad2	2Qs1vXc!f_Q5F	FAs1vad
 VCF]1a4/
0Bc!fc!fad+Q/668B>>>c!fc!fad+Q{CJJL  
Q$s1vqy1a4A	>	)AaCA	M	9AaCA	M	IBc!fc!fad+Q/668B>>>c!fc!fad+Q{CJJL  Rq   rt   c                     \        ^\        .\        ^..4      P                  \        ^4      \        ^^.^^..4      8X  g   Q h\        \        ^.\        \        ,           ^..4      P                  \        R.\        R..4      \        R^.R^..4      8X  g   Q h\        \        ^.\        \        ,           ^..4      P                  \        R3\        R3.4      \        R^.R^..4      8X  g   Q h\        \        ^.\        \        ,           ^..4      P                  \        R\        R/4      \        R^.R^..4      8X  g   Q h\        \        \        ,          .4      P                  \        \        ^,
          \        \        ^,
          /RR7      \        \        ^,
          \        ^,
          ,          .4      8X  g   Q h\         F:  p \        ^^ .^ ^..4      V P                  ^4      P                  ^^4      8X  d   K:  Q h	  R# )r   T)simultaneousNr  r  r   )r5   rc   r  rd   r   r9   r   s    ru   	test_subsrD  7  s   Aq6Aq6"#((A.&1a&1a&9I2JJJJAq6AE1:&',,q"g2w-?@Q"a!"# # #Aq6AE1:&',,q"g2w-?@Q"a!"# # #Aq6AE1:&',,aQ^<Q"a!"# # #1Q3%=q!a%AE2FQQ !" " " 1v1v&'3771:??1a+@@@@ rt   c                     \        ^\        .\        ^..4      P                  \        ^/4      \        ^^.^^..4      8X  g   Q h\        \        ^.\        \        ,           ^..4      P                  \        R\        R/4      \        R^.R^..4      8X  g   Q h\         F;  p \        ^^ .^ ^..4      V P                  ^4      P                  ^^/4      8X  d   K;  Q h	  R# )r   Nr  r  r   )r5   rc   xreplacerd   r   r9   rC  s    ru   test_xreplacerG  E  s    Aq6Aq6"#,,aV4AA ! ! !Aq6AE1:&'00!RB@Q"a!"# # #1v1v&'3771:+>+>1v+FFFF rt   c                     \        R 4      p \        R4      p\        ^\        ,          ^\        ,          ,           \        \        \        ,          ,           \        ,          .V! \        4      \        V! \        4      ,          ,           V! \        4      ,          ^^V ,          \        V \        ,          4      V ,          ,
          ,          \        ,          ..4      pVP                  4        V\        \        \        ,           \        \        ,          ,          ^\        ,           .^\        ,           ^^^\        \        V ,          4      ,          ,
          \        V ,          ,          ,          ..4      8X  g   Q h^\        ,           ^,          p\        V..4      pVP                  4        V\        V..4      8X  g   Q hVP                  \        R7       V\        VP                  \        R7      ..4      8X  g   Q hR# )r   r9  )ratioN)	r   r   r5   rc   rd   r   r   r)   r   )r   r9  r   eqs       ru   test_simplifyrK  M  s\   sAAQqS1Q3YQqSAHQ4!AaD&=!A$&QqS3q2v;q=-@(AB(FHJ 	KAJJL1q51q5/!a%I$%E1q1SAY;A.F+GIK L L L L
a%!BvAJJLvJJRJ2./01111rt   c                     \        . RO. RO.4      p V P                  \        ^^.^^.^^.^^.^^.^^.^^.^^.^	^	.^ ^ ..
4      8X  g   Q hV P                  P                  V 8X  g   Q hV P                  V P                  4       8X  g   Q hR# )r   N)
r   r   r   r   r   r   r   r   r  r   )r5   r  r   r(  s    ru   test_transposerM  ^  s    ..0 	1A33&QFFFFFFFFFF	% 	& 	& 	& 	& 3355A::33!++-rt   c                     \        ^ \        ^.. RO.4      p V P                  \        ^ ^.\        ^.^^ ..4      8X  g   Q hV P                  \        ^ \        ) ^.. RO.4      8X  g   Q hV P                  V P	                  4       8X  g   Q hV P
                  V P                  P                  8X  g   Q hV P
                  \        ^ ^.\        ) ^.^^ ..4      8X  g   Q hR# )r   N)r   r   r   )r5   r
   r  r   	conjugateHr(  s    ru   test_conjugaterQ  o  s    Aq	 	A 33&1a&a&a&" # # # # 33&1qb!*$& ' ' ' '33!++-33!##%%<<33&Aq'2q'q'# $ $ $ $rt   c            
      (   \        \        R  4       \        ^\        \        \        .^ ^\        \        .^ ^ ^\        .. RO.4      p V P                  \        . RO\        ) ^^ ^ .\        ) \        ) R^ .\        ) \        ) \        R..4      8X  g   Q hR# )c                  ,    \        ^4      P                  # r  )r9   Drm   rt   ru   r   !test_conj_dirac.<locals>.<lambda>  s    3q688rt   Nr  r  r  )rR   AttributeErrorr5   r
   rT  r(  s    ru   test_conj_diracrW    s    
>+,Aq!Aq!Aq! 	A
 33&*2A*2r2*2rAr*, - - - -rt   c                  \    \        . RO. RO. RO.4      p V P                  4       ^8X  g   Q hR# )r   Nr   r   r   r   r   r   )r   r   r   r5   tracer(  s    ru   
test_tracer]    s-    	 	A 779??rt   c                  b    \        \        ^ ^ .^ \        ^ ..4      p V P                  R8X  g   Q hR# )r   Nr  )r5   rc   rd   r   r(  s    ru   
test_shaper_    s4    Aq	Aq	 	A77frt   c                  t   \        \        ^ ^ .^ \        ^ ..4      p V P                  ^R 4       V \        \        ^ ^ .^\        ^,           ^..4      8X  g   Q hV P	                  ^ R 4       V \        \        ^,           ^ ^ .^\        ,           \        ^,           ^..4      8X  g   Q hV P                  ^ 4      \        \        ^,           ^ ^ ..4      8X  g   Q hV P                  ^ 4      p^*V^ &   V R,          \        ^,           8X  g   Q hV ^ RR13,          p^*V^ &   V R,          \        ^,           8X  g   Q hV P                  ^ 4      pV\        \        ^,           ^\        ,           .4      8X  g   Q h^ V^ &   V R,          \        ^,           8X  g   Q hV R,          p^*V^ &   V R,          \        ^,           8X  g   Q hR# )r   c                      W,           ^,           # r   rm   )rr   s   &&ru   r   !test_col_row_op.<locals>.<lambda>  s
    QUQYrt   c                 *    V \         V,          ,           # r{   )rd   )ra   r   s   &&ru   r   rd    s    QAXrt   Nr   r  r  )r5   rc   rd   row_opcol_oprowcol)r   r1c1s      ru   test_col_row_oprl    s   Aq	Aq	 	AHHQ&'1AE1' ( ( ( ( HHQ%&QAq)QAq)+ , , , ,
 558vAq!}o....	
qBBqET7a!e	
1crc6BBqET7a!e	
qBQA''''BqET7a!e	
4BBqET7a!ert   c                      \        . RO. RO.4      p V P                  ^^4       V R,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^8X  g   Q hR# )r   Nr   r   r   r   r   r   r   )r5   row_multr(  s    ru   test_row_multrp    sS     	AJJqOS6R<<S6Q;;S6R<<rt   c                      \        . RO. RO. RO.4      p V P                  ^^ ^4       V R,          ^8X  g   Q hV R,          ^8X  g   Q hV R,          ^8X  g   Q hR# )r   Nr   rn  r   r   r   r   r   r   r   )r5   row_addr(  s    ru   test_row_addru    s[     	A IIa!S6Q;;S6Q;;S6Q;;rt   c                  :   \         R ,           F  p V P                  ^4      pVP                  ^^ R 4       W! . RO. RO. RO.4      8X  g   Q hV P                  ^4      ^,          pRVR&   VP                  ^^ R 4       V W! . R	O. R
O. RO.4      8X  d   K  Q h	  R# )r  c                 "    V ^V,          ,           # r  rm   r   us   &&ru   r   !test_zip_row_op.<locals>.<lambda>      AaCrt   c                 "    V ^V,          ,           # r  rm   rx  s   &&ru   r   rz    r{  rt   NrY  r   r   r   r  r  r   )r   r  r   )r   r   r   r  )r   r9   
zip_row_op)r   r   s     ru   test_zip_row_opr    s    r{{GGAJ	Q/0C""$ % % 	% % GGAJqL$	Q/0!C##% & & 	& & rt   c                      \        . RO4      p \        . RO4      p\        . RO4      pV R9   d   Q hV R9   d   Q hV ^8w  g   Q hW8X  g   Q hW8w  g   Q hR# )r   Nr   )r   r   r   rm   rb  r   )r   r_   r`   s      ru   test_issue_3950r    sU    yAyAyAR=S>6M66M66M6rt   c                      ! R  R4      p  ! R R4      pV ! 4       pV! 4       p\        . RO4      pWC,          ^8X  g   Q h^WC&   V^,          ^8X  g   Q h\        . RO. RO.4      pWBV3,          ^8X  g   Q hV^V3,          ^8X  g   Q hWB^3,          ^8X  g   Q h^WBV3&   VR,          ^8X  g   Q h^V^V3&   VR,          ^8X  g   Q h^WB^3&   VR,          ^8X  g   Q hR# )c                   &   a  ] tR tRt o R tRtV tR# )test_issue_3981.<locals>.Index1i  c                    ^# rb  rm   selfs   &ru   	__index__)test_issue_3981.<locals>.Index1.__index__      rt   rm   Nro   rp   rq   rr   r  rs   __classdictcell____classdict__s   @ru   Index1r         	 	rt   r  c                   &   a  ] tR tRt o R tRtV tR# )test_issue_3981.<locals>.Index2i  c                    ^# r  rm   r  s   &ru   r  )test_issue_3981.<locals>.Index2.__index__  r  rt   rm   Nr  r  s   @ru   Index2r    r  rt   r  Nr   rn  r   r   )r  r  index1index2r   s        ru   test_issue_3981r    s       XFXFyA9>>AIQ4199	9%&AV^!!!QY<1QY<1AfnT7a<<AaiLT7a<<AaiLT7a<<rt   c                  ,  a  \        \        ^4      ^.4      o \        ;QJ d(    V 3R l\        ^4       4       F  '       d   K   RM	  RM! V 3R l\        ^4       4       4      '       g   Q h\        ;QJ d(    V 3R l\        ^4       4       F  '       d   K   RM	  RM! V 3R l\        ^4       4       4      '       g   Q h\        ;QJ d(    V 3R l\        ^4       4       F  '       d   K   RM	  RM! V 3R l\        ^4       4       4      '       g   Q hR# )r   c              3      <"   T F5  pSP                   ! 4       V,          SV,          P                  4       8H  x  K7  	  R # 5ir{   evalfr   r~   r_   s   & ru   r   test_evalf.<locals>.<genexpr>  s*     >Xqwwy|qtzz|+Xs   =A FTc              3      <"   T F7  pSP                   ! ^4      V,          SV,          P                  ^4      8H  x  K9  	  R# 5ir  r  r  s   & ru   r   r    s.     @x!qwwqz!}!

1-x   ?Ac              3      <"   T F7  pSP                   ! ^4      V,          SV,          P                  ^4      8H  x  K9  	  R# 5ir  r   r  s   & ru   r   r  	  s.     8x!qss1vayAaDFF1I%xr  N)r5   r   r   r   r_   s   @ru   
test_evalfr    s    Q|A3>U1X>333>U1X>>>>>3@uQx@333@uQx@@@@@38uQx83338uQx88888rt   c                  h   \        \        \        .\        \        ..4      p V P                  ! 4       R J g   Q h\        . RO. RO.4      p V P                  ! 4       RJ g   Q h\        . RO^^\        ^..4      p V P                  ! 4       R J g   Q h\        ^\        ^..4      p V P                  ! 4       R J g   Q h\        . RO.4      p V P                  ! 4       RJ g   Q h\        ^.\        .^..4      p V P                  ! 4       R J g   Q h\        ^.^.^..4      p V P                  ! 4       RJ g   Q hR# )TFNr  r  r   )r5   rc   is_symbolicr  s    ru   test_is_symbolicr    s   AA A==?d"""l+,A==?e###q!Ql+,A==?d"""Aq	{A==?d"""	{A==?e###aS1#A==?d"""aS1#A==?e###rt   c                      \        . RO.4      p V P                  RJ g   Q h\        ^.^.^..4      p V P                  RJ g   Q h\        ^^4      p V P                  RJ g   Q hR# )r   TFNr   )r5   is_upperrB   r  s    ru   test_is_upperr    se    	{A::aS1#A::aA::rt   c                      \        . RO.4      p V P                  RJ g   Q h\        ^.^.^..4      p V P                  RJ g   Q hR# )r   FTNr   )r5   is_lowerr  s    ru   test_is_lowerr  &  sH    	{A::aS1#A::rt   c                      \        ^^. RO4      p V P                  ! 4       '       g   Q h\        ^^ .^ ^..4      p V P                  ! 4       '       d   Q h\        . 4      p V P                  ! 4       '       g   Q hR# )r   N)r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r5   rz  r  s    ru   test_is_nilpotentr  -  sj    q!EFA>>AA A~~r
A>>rt   c                     ^^r\        W4      pVP                  ! ^4       ^\        W4      ,          pW#8X  g   Q hVP                  VP                  u;8X  d	   ^8X  g   Q h Q hVP                  VP                  u;8X  d	   ^8X  g   Q h Q hVP
                  VP
                  u;8X  d	   R8X  g   Q h Q h\        ^4      \        ^^4      8X  g   Q h\        ^4      \        ^^4      8X  g   Q h\        ^^4      \        ^^^ .^,          4      8X  g   Q h\        ^^4      \        ^^^.^,          4      8X  g   Q hVP                  ! ^ 4       V\        W4      8X  g   Q hR# )r   N)r   r   )rB   fillr<   r   r   r   r5   )r   r   r_   r`   s       ru   test_zeros_ones_fillr  6  s)   aqaAFFAK	DJA6M666QVV q     66QVV q     77agg''''''8uQ{"""7d1aj   A;&As1u----1:1qc!e,,,,FF1Iart   c                     \        ^ 4      p V \        4       8X  g   Q h\        ^ ^4      p V P                  ^ 8X  g   Q hV P                  ^8X  g   Q h\        ^^ 4      p V P                  ^8X  g   Q hV P                  ^ 8X  g   Q hR# r  )rB   r5   r   r   r  s    ru   test_empty_zerosr  K  so    aA==aA66Q;;66Q;;aA66Q;;66Q;;rt   c                     \        \        ^,          \        \        ,          .\        \        \        4      ,          \        \	        \        4      ,          ..4      p V P
                  ! \        4      \        ^\        ,          \        .\        \        4      \	        \        4      ..4      8X  g   Q h\        \        \        ) \        ^,          .\        \        4      ^\        ,          \        \        ) 4      ,
          \        ^\        ,          ,           ..4      P                  \        \        4      \        \        \        ) \        .\        ^ \        ..4      8X  g   Q h\        \        \        4      ^,
          \        ,          ^\        ,          \        \        ,          ,           \        \        ,          .^\        ,          \        \        4      \        \        \        ^,           4      4      ..4      P                  \        ^ 4      \        . RO\        ^ \        ^4      ..4      8X  g   Q hV P                  ! \        4      \        \        ^^4      \        ^,          ,          \        \        ^,          ,          ^,          .\        ^,          \        \        4      ,          ^,          \        ^,          \	        \        4      ,          ^,          ..4      8X  g   Q hR# )r   N)r   r   r   )r5   rc   rd   r   r   r   r   limitr   r   r   r   r  s    ru   test_issue_3749r  V  s   AqsaAh#a&123A66!91a3q63q6*:;<<<<	
QB1	Q1sA2wAaC(* ++05B<bS"Ar{+,- - - 
a&1*a1qsAqD*	
1c!fc#a!e*o&( )).q!	B3q6?+,- - - ;;q>V	!Q1	a1fQh'	
Ac!fQ1SVA&%( ) ) ) )rt   c                      \        ^4      p V P                  ^ ^4       R F.  pV P                  VR R7      V P                  RR7      8X  d   K.  Q h	  R# )r   c                     V ^ 8H  # r   rm   r  s   &ru   r   %test_inv_iszerofunc.<locals>.<lambda>j  s    art   )r  
iszerofuncr  r  N)r  r  )r9   col_swapr,  )r   r  s     ru   test_inv_iszerofuncr  f  sQ    AAJJq!uuF/?u@EEE  	    rt   c                  ^   \        R 4      w  r\        V \        V4      ,          V \        V4      ,          .4      p\        W.4      pVP	                  V4      pWBP	                  VP
                  4      8X  g   Q hWBP
                  P	                  V4      8X  g   Q hWBP
                  P	                  VP
                  4      8X  g   Q hVP
                  \        VP                  ^ ,          4      ,          V,          pVP                  \        4      pV\        ^^ .^ V ^,          ..4      8X  g   Q hR# zrho,phiN)
r   r5   r   r   r6  r  r9   r   r  r*   )rhophiXYr/  r:  s         ru   test_jacobian_metricsr  n  s    y!HCCHc#c(l+,AzA	

1A

133q!!!!qss####	C
OAA	HAACF,----rt   c                  n   \        R 4      w  r\        V \        V4      ,          V \        V4      ,          V ^,          .4      p\        W.4      p\        \        V4      V ) \        V4      ,          .\        V4      V \        V4      ,          .^V ,          ^ ..4      pVP	                  V4      V8X  g   Q hR# r  )r   r5   r   r   r6  )r  r  r  r  r/  s        ru   test_jacobian2r  {  s    y!HCCHc#c(lCF34AzA	SC4C=!	SCCL!cEq! 	A
 ::a=Art   c                  ^   \        \        \        \        ,           \        ,           4      \        \        \        ,           \        ,           4      \        \        \        ,           \        ,           4      .4      p \        \        \        \        .4      p\        ^^4       F  p\        ^^4       Fx  pV RV1R3,          pVRV1R3,          pVP                  V4      pVP                  V8X  g   Q hVP                  V8X  g   Q h\        V4       F  pVRV3,          V8X  d   K  Q h	  Kz  	  K  	  R# )r   Nr   )	r5   r   rc   rd   re   r   r6  r   r   )r  r  r~   r   X_sliceY_slicer/  r   s           ru   test_issue_4564r    s    AEAIAEAIAEAI?@A1ayA1a[q!AAhGAhG  )A66Q;;66Q;;1XAw')))   rt   c                    a a \        \        \        \        ,           \        ,           4      \        \        \        ,           \        ,           4      .\        \        \        ,           \        ,           4      \        \        \        ,           \        ,           4      ..4      o \        \        V 3R  l4       S R,          o \        \        \        .\        \        ..4      o\        \        V V3R l4       \        \        V 3R l4       R# )c                  V   < S P                  \        \        \        \        .4      4      # r{   r6  r5   rc   rd   re   r  s   ru   r   (test_nonvectorJacobian.<locals>.<lambda>  s    ajjAq	):;rt   c                  &   < S P                  S4      # r{   )r6  r  r  s   ru   r   r    s    ajjmrt   c                  d   < S P                  \        \        \        .\        \        ..4      4      # r{   r  r  s   ru   r   r    s     ajj1a&1a&0B)CDrt   Nr  )r5   r   rc   rd   re   rR   	TypeErrorr  s   @@ru   test_nonvectorJacobianr    s    QUQYQUQY0QUQYQUQY02 	3A
9;<	$AAA A
9+,
9DErt   c                      \        ^^.^^..4      p V P                  4       pVP                  ^8X  g   Q h\        ^4       F  pW,          V^,           8X  d   K  Q h	  R# r   )r5   vecr   r   )r   m_vecr~   s      ru   test_vecr    sV    AA AEEGE::??1Xx1q5    rt   c                     \        ^^.^^..4      p V P                  4       pVP                  ^8X  g   Q h\        ^4       F  pW,          V^,           8X  d   K  Q h	  V P                  RR7      pV^ ,          ^8X  g   Q h\        ^\        \        \
        ,           ,          .\
        \        ,          \        ^,          ,           ^..4      p V P                  RR7      pV^ ,          \
        \        ,          \        ^,          ,           8X  g   Q h\        ^\        \        \
        ,           ,          .\
        \        ,          ^..4      p V P                  RRR7      pV^ ,          \
        \        ,          8X  g   Q h\        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       R# )	r   F)diagonal)r  check_symmetryc                  :    \        ^^..4      P                  4       # rb  r5   vechrm   rt   ru   r   test_vech.<locals>.<lambda>      v1vh/446rt   c                  @    \        ^^.^^..4      P                  4       # rb  r  rm   rt   ru   r   r        v1v1v&67<<>rt   c                  :    \        ^^..4      P                  4       # rb  r  rm   rt   ru   r   r    r  rt   c                  @    \        ^^.^^..4      P                  4       # rb  r  rm   rt   ru   r   r    r  rt   N)	r5   r  r   r   rc   rd   rR   r+   r#  )r   m_vechr~   s      ru   	test_vechr    sH   AA AVVXF;;!1XyAE!!! VVUV#F!9>>Aq1uI1q!tQ01AVVUV#F!9!ad
"""Aq1uI1a)*AVVU5V9F!9!
:67
:>?
:67
:>?rt   c                     a  \        . RO4      \        . RO4      8X  g   Q h^^^..o \        \        V 3R l4       \        S RR7      \        . RO4      8X  g   Q hR# )r   c                     < \        S 4      # r{   )r8   r   s   ru   r   test_diag.<locals>.<lambda>  s	    tAwrt   F)strictNr   )r8   r5   rR   r#  r   s   @ru   	test_diagr    sK    	?fY////	
AsA
:'% F9$5555rt   c                  N   \        ^^.^^..4      p \        ^\        .\        ^..4      p\        ^\        ^.\        ^\        .\        \        \        ..4      pV P                  ! 4       V .8X  g   Q hVP                  ! 4       V.8X  g   Q hVP                  ! 4       V.8X  g   Q hR# r   )r5   rc   rd   re   get_diag_blocksr_   r`   ra   s      ru   test_get_diag_blocks1r    s    AA AAA AAq	Aq!9q!Qi01A1#%%%1#%%%1#%%%rt   c                     \        ^^.^^..4      p \        ^\        .\        ^..4      p\        ^\        ^.\        ^\        .\        \        \        ..4      p\	        WV4      P                  4       WV.8X  g   Q h\	        WV4      P                  4       WV.8X  g   Q h\	        WV4      P                  4       WV.8X  g   Q h\	        W"V4      P                  4       W"V.8X  g   Q hR# r   )r5   rc   rd   re   r8   r  r  s      ru   test_get_diag_blocks2r    s    AA AAA AAq	Aq!9q!Qi01Aa=((*qQi777a=((*qQi777a=((*qQi777a=((*qQi777rt   c                     \        ^^.^^..4      p \        ^\        .\        ^..4      p\        ^\        ^.\        ^\        .\        \        \        ..4      p\	        WV4      pVP                  RR7      \	        V P
                  ! 4       VP
                  ! 4       VP
                  ! 4       4      8X  g   Q h\	        WV4      pVP                  RR7      \	        V P
                  ! 4       VP
                  ! 4       VP
                  ! 4       4      8X  g   Q h\	        WV4      pVP                  RR7      \	        V P
                  ! 4       VP
                  ! 4       VP
                  ! 4       4      8X  g   Q h\	        W WW 4      pVP                  RR7      \	        V P
                  ! 4       V P
                  ! 4       VP
                  ! 4       V P
                  ! 4       VP
                  ! 4       V P
                  ! 4       4      8X  g   Q hVP                  RRR7      \	        V P
                  ! RR7      V P
                  ! RR7      VP
                  ! RR7      V P
                  ! RR7      VP
                  ! RR7      V P
                  ! RR7      4      8X  g   Q hR# )r   T)try_block_diagr  )r  r  r  N)r5   rc   rd   re   r8   r,  )r_   r`   ra   r   s       ru   test_inv_blockr    s   AA AAA AAq	Aq!9q!Qi01AQ1A555%aeegquuw)HHHHQ1A555%aeegquuw)HHHHQ1A555%aeegquuw)HHHHQ1A555%	!%%'1557AEEGQUUW*> > > >55U53t	UQUU%0!%%u2E	UQUU%0!%%u2E8G G G Grt   c                 r   \        \        R 4       \        \        R 4       \        \	        ^4      4      \        ^4      8X  g   Q h\        \        ^4      4      \        ^4      8X  g   Q h\        \        R 4       \        \	        ^4      4      \        ^4      8X  g   Q h\        \        ^4      4      \        ^4      8X  g   Q h\        \        R 4       \        \	        ^4      \        ^4      4      \        ^^4      8X  g   Q h\        \        R 4       \        \        R 4       \        \        R 4       R# )	z[
Check that matrix dimensions can be specified using any reasonable type
(see issue 4614).
c                      \        ^R4      # )r   r  rB   rm   rt   ru   r   $test_creation_args.<locals>.<lambda>  s
    uQ|rt   c                      \        ^^^^4      # rb  r  rm   rt   ru   r   r    s    eAq!Q/rt   c                      \        R 4      #       @r  rm   rt   ru   r   r    s    uRyrt   c                      \        R 4      # r  r9   rm   rt   ru   r   r    s    s2wrt   c                      \        ^4      # r  r   rm   rt   ru   r   r    s    fQirt   c                      \        ^^4      # rb  r   rm   rt   ru   r   r    s
    fQlrt   c                      \        ^^..4      # rb  r   rm   rt   ru   r   r    s    vq1#h/rt   N)rR   r#  r  rB   r|   r   r9   r<   rm   rt   ru   test_creation_argsr    s    
 :+,
9/0Q=E!H$$$a(((
:()s1v;#a&   wqz?c!f$$$
:'A
#tAqz111
9'(
9*+
:/0rt   c                  <   \        ^^. RO4      p V P                  4       '       d   Q hV P                  4       '       g   Q hV P                  RR7      '       g   Q h\        ^^. RO4      p V P                  4       '       g   Q h\        ^^^4      p V P                  4       '       g   Q hV P                  4       '       g   Q h\        ^^. RO4      p V \        ^^^4      8X  g   Q h\        ^^\	        ^^4      4      p V P                  4       '       d   Q hV P                  4       '       g   Q h\        R4      p V P                  4       '       g   Q h\        R4      p V P                  4       '       g   Q h\        ^^^\
        ^,          ^\
        ,          ,           ^,           \        \
        ^,           ^,          ^^ \        ^ ^.	4      p V P                  4       '       g   Q hV P                  RR7      '       d   Q hV P                  4       P                  RR7      '       g   Q hR# )	r   Fr)   Nr   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   r   r   ))r   r   r   r   )r   r   r   )r   r   r   )r5   is_diagonalis_symmetricr8   rB   rc   rd   r   r   s    ru   test_diagonal_symmetricalr
    s   q!\"A}}>>>>5>))))q!\"A==???Q1A==???>>q!01AQ1q!U1a[!A~~==???'(A==???%&A==???q!aA!aQUQJ1aAFGA>>~~u~----88:""E"2222rt   c                  ,	  a \        ^^\        ,           .^\        ,
          ^..4      oSP                  4       '       g   Q h\        ^^. RO4      oSP                  4       '       d   Q hSP                  4       '       d   Q h\	        \
        V3R l4       \        ^^^4      oSP                  4       w  rV \        ^4      8X  g   Q hVS8X  g   Q h\        ^^. R	O4      oSP                  4       '       g   Q hSP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q h\        ^^. R
O4      oSP                  4       '       g   Q hSP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q hV \        ^4      8X  g   Q hVS8X  g   Q h\        ^^. RO4      oSP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q h\        ^^. RO4      oSP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q hV  F!  pVP                  4       ^,          ^8X  d   K!  Q h	  \        ^^. RO4      oSP                  4       '       g   Q hSP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q hV \        ^ ^.^^ ..4      8X  g   Q h\        ^^. RO4      oSP                  R4      '       d   Q h\	        \        V3R l4       SP                  4       '       g   Q hSP                  4       w  rV P                  4       S,          V ,          V8X  g   Q h\        ^^. RO4      oSP                  4       '       d   Q h\	        \        V3R l4       \        ^^. RO4      oSP                  4       '       d   Q h\	        \        V3R l4       \        R4      w  r4rV\        ^^W5WT.4      oSP                  4       '       g   Q hSP                  4       '       g   Q hR# )r   c                  $   < S P                  4       # r{   r  r   s   ru   r   &test_diagonalization.<locals>.<lambda>%      rt   Tc                  &   < S P                  R 4      # )Tr  r   s   ru   r   r  Q  s    d 3rt   c                  $   < S P                  4       # r{   r  r   s   ru   r   r  Y      rt   c                  $   < S P                  4       # r{   r  r   s   ru   r   r  ]  r  rt   a b c dNr   r   r      r   r   r  )r   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   rP  r   r  )r   r   r  r   r   )	r   r   r   r  r   r   r   r  r   )r5   r
   r}  r	  rR   r-   r8   r  r9   r,  as_numer_denomr  r,   r   )r.  rT  r~   r_   r`   ra   rb   r   s          @ru   test_diagonalizationr    s   AaC1Q3(#$A    q!+,A""$$$$~~
!89 	Q1A]]_FQA;;6M6q!\"A>>    ]]_FQ557Q;?aq!\"A>>    ]]_FQ557Q;?aA;;6M6q!\"A    ]]_FQ557Q;?aq!12A    ]]_FQ557Q;?a!!$)))  	q!\"A==???    ]]_FQ557Q;?aAA'(((( 	q!]#A""4((((
;34    ]]_FQ557Q;?a 	q!\"A""$$$$
;/0q!56A""$$$$
;/0 #JA!q!aA\"A>>    rt   c                    a  \        ^ ^.^^ ..4      o S P                  ! 4       RJ g   Q h^ S R&   S P                  ! 4       RJ g   Q h\        ^ ^.^^ ..4      o S P                  ! 4        ^ S R&   \        \        V 3R l4       R# )r   TFc                  &   < S P                   ! 4       # r{   r  r  s   ru   r   "test_issue_15887.<locals>.<lambda>p  s    rt   Nr   )rC   r}  r  rR   r,   r  s   @ru   test_issue_15887r  f  s    QFQF+,A D(((AdG E)))QFQF+,AMMOAdG
;/0rt   c            	        a \        ^^. RO4      o\        \        V3R l4       \        ^^. R	O4      o\        ^^. R
O4      p SP                  4       w  rW8X  g   Q hV SP	                  4       ^,          8X  g   Q h\        ^^. RO4      o\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        ^^. RO4      o\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        ^^. RO4      o\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        ^^. RO4      o\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        ^^. RO4      o\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        ^^. RO4      oSP                  4       '       d   Q h\        ^^. RO4      p SP                  4       w  rW8X  g   Q h\        \        R^nR7      \        R^nR7      .\        R^nR7      \        R^nR7      ..4      oSP                  4       w  rVP                  4        F.  p\        V\        4      '       g   K  VP                  ^n8X  d   K.  Q h	  R# )r   c                  $   < S P                  4       # r{   )r+  r   s   ru   r   "test_jordan_form.<locals>.<lambda>v  r  rt   z1.0)	precisionz2.0z3.14159265358979323846264338327z4.0Nr  )	r   r  r   r   ir   r  ir  )	r  r   r   r   r   r   r   r   r   	r   r   r   rP  r   r   r  r   r   )	r   r   r   r   r   r   r   r   r   )	r   r   r   r   rQ  r   r   r=  )	r  r   r   r   r  r   r   r   r  )	r   rQ  r   r   rv  r   r   r   )	r   r   r   r   r   r   r   r   r   )r   r   r  r   r   r  r   r   r   r   r  r   r  r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r=  r   r   r  r   r   r  r%  r=  r  r   r   r   )r  r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r  r  r  r  r   r   r   r   r  r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )
r5   rR   r-   r+  r  r}  r	   valuesr   _prec)Jmustr.  r/  termr   s       @ru   test_jordan_formr*  s  sV   q!+,A
!89 	q!:;A1a56E==?DA::AMMOA&&&& 	q!23A 1a45E==?DA:: 	q!45A 1a78E==?DA:: 	q!34A1a45E==?DA::q!LMA1a  E
 ==?DA::q!LMA 1a ' (E ==?DA::q!JKA""$$$$1aIJE==?DA:: 	u,eES.IJ8CH%PUadJefh 	iA==?DA
dE""::$$$ rt   c            	         \        . RO. RO. RO. RO.4      p ^^\        ,          ,
          p^^\        ,          ,           p\        V^^ ^ .^ V^ ^ .^ ^ V^.^ ^ ^ V..4      p\        V^^ ^ .^ V^ ^ .^ ^ V^.^ ^ ^ V..4      pV P                  4       w  rVWc8X  g	   Wd8X  g   Q h\        WV,          VP	                  4       ,          4      V 8X  g   Q hR# )r   N)r   r   r   r   )rP  r   r   r   )r   r   r   r   )r   r   rP  r   )r5   r
   r+  r)   r,  )r   pqJmust1Jmust2r.  r/  s          ru   #test_jordan_form_complex_issue_9274r0    s        " 	#A 	
AaCA	AaCAaAq\Aq\Aq\Aq\# $F aAq\Aq\Aq\Aq\# $F ==?DA;!+%%ACK A%%%rt   c                      \        . RO. RO. RO. RO.4      p V P                  4       w  rV\        . RO. RO. RO. R	O.4      8X  g   Q hV\        . R
O. RO. R	O. RO.4      8X  g   Q hR# )r   Nr  r  )r   r   r   r   r  r   r   r   r   r  r   r   r  )r5   r+  r-  s      ru   test_issue_10220r3    s     	A ==?DA$$$& ' ' ' ' $$$$	& ' ' ' 'rt   c                  $   \        . RO. RO. RO. RO.4      p V P                  4       w  rVP                  4       \        \        ) \        ) ^,          \        \        ^,          .R\        ,           ^ R\        ,
          ^ .^ \	        ^4      ) ^,          \        ^,          ,
          ^ \	        ^4      ) ^,          \        ^,          ,           .. RO.4      8X  g   Q hV\        \        ) ^^ ^ .^ \        ) ^ ^ .^ ^ \        ^.^ ^ ^ \        ..4      8X  g   Q hR# )r   Nr  r  r  r  r   r  r   r   r  r  r   r   r   r   r2  )r5   r+  r   r
   r   )r   r.  r/  s      ru   test_jordan_form_issue_15858r8    s    	 	A
 ]]_FQ88:R1"Q$Q!A#6	aQQA6!A$q1Q3Q1a!A#66	!8 9 9 9 9
 
Q1	
QB1	
Aq!	
Aq!	    rt   c                  x   \        R 4      w  rp\        V) V ,
          V^,          W,           ,          ,           W,          W,           ,          .W,          W,           ,          V) V^,          W,           ,          ,           ..4      pVP                  \        4      pV\	        \        ^,          W,          W ,          ,           ^V,          V,          ,           W,           ,          \        ,          ,           W,          V ,          W,           ,          ,           \        RR7      8X  g   Q h\        V4      \        J g   Q h\        ^^.^^ ..4      pVP                  4       VP                  \        4      u;8X  d/   \        \        ^,          \        ,
          ^,
          4      8X  g   Q h Q h\        ^^.\        ^ ..4      pVP                  \        4      pVP                  \        8w  g   Q hVP                  4       P                  VP                  \        4      \        ^,          ^\        ,          ,
          8X  g   Q hR# )z
UA K_i K_wzZZ(K_i,K_w,UA)r;  N)
r   r5   r>  rc   r#   r   r$   genas_exprr  )UAK_iK_wr   r>  r,  s         ru   test_Matrix_berkowitz_charpolyr?    s   <(LBS#S!VSY//sw	7JKGSY/#Q	8J1JKM 	NA zz!}HQTSVcf_quSy039=a??
CI !*:	<< < < >X%%%AA A::<1::a=BHQTAX\,BBBBBBAA A	

1A55A::99;AEE1%A!333rt   c            
         \        R 4      p \        P                  ! ^V 4      pVP                  4       \        \	        V 4      ..4      8X  g   Q h\        P                  ! ^V 4      pVP                  4       \        \	        V 4      \	        V 4      \	        V 4      ^,          .^ \	        V 4      \	        V 4      .^ ^ \	        V 4      ..4      8X  g   Q hR# )lamdaN)r   r5   jordan_block_eval_matrix_exp_jblockr   lr   s     ru   test_exp_jordan_blockrF    s    wAAq!A$$&&3q6(*<<<<Aq!A$$&VSVSVAX&AA3q6N 	  rt   c            	      
   \        ^^.^ R..4      p \        \        ^4      R\        R4      ,          ^,          ^\        ^4      ,          ^,          ,           .^ \        R4      ..4      pV P                  4       V8X  g   Q h\        V 4      V8X  g   Q h\        ^^ .^ ^..4      p V P                  4       \        \        ^ .^ \        ..4      8X  g   Q h\        V 4      \        \        ^ .^ \        ..4      8X  g   Q h\        ^R.^^..4      p V P                  4       \        \        \        ^4      ,          \        ) \	        ^4      ,          .\        \	        ^4      ,          \        \        ^4      ,          ..4      8X  g   Q hR# )r   Nr  rP  r  )r5   r   r   r   r   )r   m_exps     ru   test_exprI    s<   AB !ASVRBZ\Ac!fHQJ67!SWFGE557eq6U??AA A557fq!fq!f-....q6VaVaV,----B!Q !A557fqQx!CF3aAh#a&5IJKKKKrt   c                  x  a \        R 4      p \        P                  ! ^V 4      oSP                  4       \        \	        V 4      ..4      8X  g   Q h\        P                  ! ^V 4      oSP                  4       \        \	        V 4      ^V ,          R^V ^,          ,          ,          ^^V ^,          ,          ,          .^ \	        V 4      ^V ,          R^V ^,          ,          ,          .^ ^ \	        V 4      ^V ,          .^ ^ ^ \	        V 4      ..4      8X  g   Q h\        . RO. RO. RO.4      o\        \        V3R l4       R# )rA  c                  $   < S P                  4       # r{   )r   r   s   ru   r   test_log.<locals>.<lambda>=  s    rt   Nr  r  r   r  )r   r5   rB  _eval_matrix_log_jblockr   rR   r,   rD  s    @ru   test_logrN  '  s   wAAq!A$$&&3q6(*<<<<Aq!A$$&Q1b!AqD&k1a1f:6CFAaCQq!tV-As1vqs#Aq#a&!		

 
 
 				A
 ;(rt   c                  H   \        \        \        3R34      p V P                  \        4      '       g   Q hV P                  \        4      '       d   Q hV P                  \
        4      '       g   Q hV P                  \        ^4      p V P                  \        4      '       d   Q hR# )r   Nr  )r5   rc   rd   hasre   r   r  r   s    ru   test_hasrQ  @  sm    A A5588O8uuQxx<55===	q!AuuQxx<xrt   c            	          \        R 4      p \        ^^V \        V 4      ^,          \        V 4      ^,          ,           \        P
                  .4      p\        V4      w  r#rEV\        P
                  8X  g   Q hR# )rc   N)r   r5   r   r   r   r   r0   rc   columnpivot_offset	pivot_valpivot_assumed_nonzero
simplifieds         ru   :test_find_reasonable_pivot_naive_finds_guaranteed_nonzero1rY  J  s_     	sAAq1c!fai#a&!)3QVV<=F$V, ?L2rt   c            
      .   \        R 4      p \        ^^V \        V 4      ^,          \        V 4      ^,          ,           V ^,          ,           \        V 4      ^,          \        V 4      ^,          ,           .4      p\	        V\
        R7      w  r#rEV^8X  g   Q hR# rc   )simpfuncN)r   r5   r   r   r0   r1   rS  s         ru   :test_find_reasonable_pivot_naive_finds_guaranteed_nonzero2r]  V  s~     	sAAqVQYs1vqy(A-VQYs1vqy(*+F
 	%Vi@ ?L2>>rt   c            
         \        R 4      p \        ^^V \        V 4      ^,          \        V 4      ^,          ,           V ,           \        V 4      ^,          \        V 4      ^,          ,           .4      p\	        V\
        R7      w  r#rE\        V4      ^8X  g   Q hV^ ,          ^ ,          ^8X  g   Q hV^ ,          ^,          ^V ,           8X  g   Q hV^,          ^ ,          ^8X  g   Q hV^,          ^,          ^8X  g   Q hR# r[  )r   r5   r   r   r0   r1   lenrS  s         ru   +test_find_reasonable_pivot_naive_simplifiesr`  e  s     	sAAqVQYs1vqy(*VQYs1vqy(*+F
 	%Vi@ ?L2 z?aa=q   a=qs"""a=q   a=q   rt   c                  X  a a \        \        R  4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \
        R 4       \        \        R 4       \        \        R	 4       \        \        R
 4       \        \        R 4       \        \
        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        \
        R 4       \        \        R 4       \        \        R 4       \        \
        R 4       \        \        R 4       \        \
        R 4       \        \        R 4       \        \
        R 4       \        \
        R 4       \        \
        R 4       \        \        R  4       \        \        R! 4       \        \        R" 4       \        \        R# 4       \        \        R$ 4       \        \        R% 4       \        \        R& 4       \        \        R' 4       \        R,4      o \        \        V 3R( l4       \        . R-O.4      o\        . R.O. R/O. R0O.4      o \        \        V V3R) l4       \        . R.O. R/O. R0O.4      o \        \        V V3R* l4       R+# )1c                  "    \        ^^.^..4      # rb  r   rm   rt   ru   r   test_errors.<locals>.<lambda>x  s    v1vsm4rt   c                  ,    \        ^^..4      R,          # )r   )g333333?r   r   rm   rt   ru   r   rc  y      v1vh/7rt   c                  ,    \        ^^..4      R,          # )r   )r   g@r   rm   rt   ru   r   rc  z  re  rt   c                      \        ^^RR7      # )r   T)ra   r  )r=   rm   rt   ru   r   rc  {  s    z!qDArt   c                  <    \        ^^.4      P                  ^^4      # rb  )r5   r  rm   rt   ru   r   rc  |  s    vq!f~55a;rt   c                  ^    \        ^^.^^..4      P                  ^^ .\        ^^.4      4      # rb  )r5   copyin_matrixrm   rt   ru   r   rc  ~  s/    AA'(661vvq!f~Nrt   c                  X    \        ^^.^^..4      P                  ^ ^.\        4       4      # rb  )r5   copyin_listsetrm   rt   ru   r   rc    s1    fq!fq!f%56BBADurt   c                  @    \        . RO. RO.4      P                  4       # )r   r   )r   r   r   r5   r,  rm   rt   ru   r   rc    s    I0F)G)K)K)Mrt   c                  \    \        ^^^^.4      P                  \        ^^.^^..4      4      # rb  )r5   row_joinrm   rt   ru   r   rc    s.    q!aV$--fq!fq!f5E.FGrt   c                  X    \        ^^.4      P                  \        ^^.^^..4      4      # rb  )r5   col_joinrm   rt   ru   r   rc    s)    FAq6N33FQFQF;K4LMrt   c                  X    \        ^.4      P                  ^\        ^^.^^..4      4      # rb  )r5   
row_insertrm   rt   ru   r   rc    6    vqc{55a!B1vA :  rt   c                  X    \        ^.4      P                  ^\        ^^.^^..4      4      # rb  )r5   
col_insertrm   rt   ru   r   rc    rv  rt   c                  8    \        ^^.4      P                  4       # rb  r[  rm   rt   ru   r   rc    s    A)=)=)?rt   c                  8    \        ^.4      P                  ^4      # rb  )r5   r  rm   rt   ru   r   rc    s    faSk33A6rt   c                  D    \        ^^.^^..4      P                  ^^4      # rb  )r5   minorrm   rt   ru   r   rc    s#    v1v1v&67==aCrt   c                  D    \        ^^.^^..4      P                  ^^4      # rb  )r5   minor_submatrixrm   rt   ru   r   rc    s#    v1v1v&67GG1Mrt   c                  :    \        . RO4      P                  ^4      # r   r   r5   crossrm   rt   ru   r   rc    s    fY/55a8rt   c                  :    \        . RO4      P                  ^4      # r  r5   dotrm   rt   ru   r   rc    s    fY/33A6rt   c                  P    \        . RO4      P                  \        ^^.4      4      # r  r  rm   rt   ru   r   rc    s    vi044VQF^Drt   c                  :    \        ^^.4      P                  . 4      # rb  r  rm   rt   ru   r   rc    s    vq!f~11"5rt   c                  :    \        ^^.4      P                  R4      # )r   r_   r  rm   rt   ru   r   rc    s    faVn005rt   c                  >    \        ^^.4      P                  . RO4      # r  r  rm   rt   ru   r   rc    s    vq!f~11)<rt   c                  8    \        . RO4      P                  4       # r  )r5   r   rm   rt   ru   r   rc    s    	):)>)>)@rt   c                  @    \        ^^.^^..4      P                  4       # rb  )r5   
normalizedrm   rt   ru   r   rc        v1v1v&67BBDrt   c                  <    \        ^^.4      P                  RR7      # )r   znot a methodr  ro  rm   rt   ru   r   rc    s    vq!f~111Hrt   c                  8    \        ^^.4      P                  4       # rb  r5   
inverse_GErm   rt   ru   r   rc        A)B)B)Drt   c                  @    \        ^^.^^..4      P                  4       # rb  r  rm   rt   ru   r   rc    r  rt   c                  8    \        ^^.4      P                  4       # rb  r5   inverse_ADJrm   rt   ru   r   rc    s    A)C)C)Ert   c                  @    \        ^^.^^..4      P                  4       # rb  r  rm   rt   ru   r   rc    s    v1v1v&67CCErt   c                  8    \        ^^.4      P                  4       # rb  )r5   
inverse_LUrm   rt   ru   r   rc    r  rt   c                  8    \        ^^.4      P                  4       # rb  r  rm   rt   ru   r   rc    s    A)D)D)Frt   c                  8    \        ^^.4      P                  4       # rb  r5   rC  rm   rt   ru   r   rc    s    A););)=rt   c                  D    \        ^^.^^..4      P                  RR7      # )r   zNot a real methodr  r  rm   rt   ru   r   rc    s%    AA'(,,4G,Hrt   c                  P    \        . RO. RO. RO. RO.4      P                  RR7      # )r   zNot functionr  r  r  r  r  r  rm   rt   ru   r   rc    s*    l)+ ,,/C>C,JKrt   c                  P    \        . RO. RO. RO. RO.4      P                  RR7      # )r   Fr  r  r  r  r  r  rm   rt   ru   r   rc    s*    l)+ ,,/C5C,ABrt   c                  V    \        \        ^^.^^..4      \        ^^.^^..4      4      # rb  r:   r5   rm   rt   ru   r   rc    s/    AA/0&1a&1a&9I2JKrt   c                  8    \        \        ^^.^^..4      . 4      # rb  r  rm   rt   ru   r   rc    s    wv1v1v.>'?Drt   c                  :    \        \        R 4      ^,          R4      # )rc   r_   )r:   r   rm   rt   ru   r   rc    s    wvc{A~s;rt   c                  &    \        ^4      R,          # )r   )r   r   r  rm   rt   ru   r   rc        s1vd|rt   c                  &    \        ^4      R,          # )r   )r   r   r  rm   rt   ru   r   rc    r  rt   c                  &   < S P                  R 4      # )zmethod=LU_decomposition())rC  r(  s   ru   r   rc    s    quu%@Art   c                  (   < S P                  R S4      # )g@)ru  r   Vs   ru   r   rc    s    q||C3rt   c                  (   < S P                  RS4      # )g@g)rx  r  s   ru   r   rc    s    q||D!4rt   Nr  )r   r   r   r   r   r   r   r   )rR   r#  r  r+   r  r-   r5   r  s   @@ru   test_errorsr  w  sI   
:45
:78
:78
:AB
:;<
:NP
9  
!MN
:GI
MO
:  
:  
!?@
967
:CD
:MN
989
967
:DE
:56
956
:<=
!@A
:DE
:HI
!DE
:DE
!EF
:EF
!DE
!FG
!=>
:HJ
:	KL :	BC :KM
:DE
:;<
:+,
:+,NOA
:AB~A	9i01A
:34	9i01A
:45rt   c                  
   \        \        4       4      ^ 8X  g   Q h\        \        ^^..4      4      \        \        ^.^..4      4      u;8X  d	   ^8X  g   Q h Q h\        \        ^ ^R 4      4      \        \        ^^ R 4      4      u;8X  d	   ^ 8X  g   Q h Q h\        \        . RO. RO.4      4      ^8X  g   Q h\        ^.4      \        ^..4      8X  g   Q h\        4       '       d   Q h\        4       \        . 4      8X  g   Q hR# )r   c                     ^ # r   rm   r}   s   &&ru   r   test_len.<locals>.<lambda>  s    rt   c                     ^ # r   rm   r}   s   &&ru   r   r    s    art   Nr   r   )r_  r5   rm   rt   ru   test_lenr    s    vx=Av1vh CaSz(:$;@q@@@@@vaN+,F1a()/-./ / / / /vy),-.!3331#;&1#-'''xx<8vbz!!!rt   c            
         \        ^^\        3\        ^^3^
^\        ^,          334      p V P                  \        4      \        \        ^\        ,          \        ^,          ^,          3\        \        ,          ^\        ,          ^\        ,          3^
\        ,          ^\        ,          \        ^,          ^,          334      8X  g   Q hV P                  \        4      \        \        ^\        ,          \        \        ,          3\        ^,          ^,          ^\        ,          ^\        ,          3^
\        ,          ^\        ,          \        \        ^,          ,          334      8X  g   Q hR# r   )r5   rc   rd   r   r   s    ru   test_integrater    s    Aq	Aq!9r1adm45A;;q>AaCAa 1Q3!QqS/BqD!A#q!tAv3FGHI I I;;q>AaC11Q!QqS1BqD!A#qAv3FGHI I Irt   c                      \        ^^\        \        4      \        ,          3\        ^^3^
^\        ^,          ^,           334      p V P	                  \        ^ 4      \        R\        ^^3R34      8X  g   Q hR# )r   N)r   r   r   )r   r   r   )r5   r   rc   rd   r  r   s    ru   
test_limitr    sa    As1vax 1a)b!QTAX->?@A771a=FI1ay*#EFFFFrt   c            	         \        ^^\        3\        ^^3^
^\        ^,          ^,           334      p \        V P	                  \        4      \        V 4      4      '       g   Q hV P	                  \        4      \        RR^ ^ ^\        ,          334      8X  g   Q hV P	                  \        4      \        R4      8X  g   Q h\	        V \        4      \        RR^ ^ ^\        ,          334      8X  g   Q h\	        V \        4      \        R4      8X  g   Q hV P                  4       p\        VP	                  \        4      \        V4      4      '       g   Q hVP	                  \        4      \        RR^ ^ ^\        ,          334      8X  g   Q hVP	                  \        4      \        R4      8X  g   Q h\	        V\        4      \        RR^ ^ ^\        ,          334      8X  g   Q h\	        V\        4      \        R4      8X  g   Q hV P	                  \        RR7      \        V \        RR7      8X  g   Q h\	        V \        RR7      \        V \        RR7      8X  g   Q hR# )r   Fr  Nr  r   )r   rY  r   )	rC   rc   rd   r   r   r   r  rD   r[   )r   A_imms     ru   	test_diffr    s   Q1I1ay2q!Q$(2CDEAaffQia))))66!9*Iy1a1++NOOOO66!9*+LMMMM1:+Y	Aq!A#;,OPPPP1:+,MNNNNNNEejjmT%[1111::a=0)YAqQRs1TUUUU::a=01RSSSSq>19i!QPQRSPS2UVVVVq>12STTTT66!e6$1u(MMMM1u%A)NNNNrt   c                  p   \        \        \        .\        \        ..4      p V P                  V 4      \        ^^ .^ ^ ..^ ^.^ ^ ...^ ^ .^^ ..^ ^ .^ ^....4      8X  g   Q h\        W 4      \        ^^ .^ ^ ..^ ^.^ ^ ...^ ^ .^^ ..^ ^ .^ ^....4      8X  g   Q hV P                  4       pVP                  V4      \        ^^ .^ ^ ..^ ^.^ ^ ...^ ^ .^^ ..^ ^ .^ ^....4      8X  g   Q h\        W4      \        ^^ .^ ^ ..^ ^.^ ^ ...^ ^ .^^ ..^ ^ .^ ^....4      8X  g   Q hV P                  \        4      \        ^ ^ .^ ^ ..4      8X  g   Q h\        \        \        .4      pV P                  V4      \        P                  ! ^^^^4      8X  g   Q hV P                  V 4      \        ^^ .^ ^ ..^ ^.^ ^ ...^ ^ .^^ ..^ ^ .^ ^....4      8X  g   Q hVP                  \        \        ..4      pVP                  R8X  g   Q hV\        ^.^ ..^ .^...4      8X  g   Q h\        R4      pV! \        \        \        4      pVP                  \        \        \        ..4      \        VP                  \        4      VP                  \        4      VP                  \        4      .4      8X  g   Q hVP                  \        \        \        .^34      \        VP                  \        ^4      VP                  \        \        4      VP                  \        \        4      .VP                  \        \        4      VP                  \        ^4      VP                  \        \        4      .VP                  \        \        4      VP                  \        \        4      VP                  \        ^4      ..4      8X  g   Q h\        \        4      \        \        4      ,          pVP                  \        \        ..4      \        \!        \        4      \        \        4      ,          \        \        4      \        \        4      ,          .4      8X  g   Q hVP                  \        \        \        334      \        \!        \        4      \        \        4      ,          \        \        4      \        \        4      ,          .4      8X  g   Q hVP                  \        \        \        334      \        \        \        4      ) \        \        4      ,          \!        \        4      \        \        4      ,          .4      8X  g   Q hVP                  \        \        33\        \        ..4      \        \!        \        4      \        \        4      ,          \        \        4      ) \        \        4      ,          .\        \        4      \        \        4      ,          \!        \        4      \        \        4      ,          ..4      8X  g   Q hVP                  \        4      P                  \        4      P                  \        4      VP                  \        \        \        33^4      R,          8X  g   Q hVP                  \        4      P                  \        4      P                  \        4      VP                  \        \        \        33^4      R,          8X  g   Q hVP                  \        \        \        ..\        \        \        334      \        \        \        \        3 UUu. uFA  p\        \        \        3 Uu. uF"  qP                  V4      P                  V4      NK$  	  upNKC  	  upp4      8X  g   Q h\        P
                  ! \#        \        \        ..4      4      p	\%        V	\        4      '       g   Q hV	\#        ^^ ..4      8X  g   Q h\        ^,          P                  \#        \        \        ..4      4      p	\%        V	\        4      '       g   Q hV	\#        ^\        ^,          ,          ^ ..4      8X  g   Q hR# u upi u uppi )r   r9  N)r   r   r   r  r}  )rC   rc   rd   re   rf   r   rZ   r  r_   rD   r`   rB   r   r   r   r   r   r5   r   )
r   r  r   dBr9  fxyzexprr   r~   ress
             ru   test_diff_by_matrixr    s    	QFQF+,A66!9!Q!Q 0Aq6Aq62BC1vPQSTvFVZ[]^Y_bcefagXhEijkkkk:1a&1a&!1QFQF3C DAQRTUPVGW[\^_Z`cdfgbhYiFjkllllNNE::e!Q!Q(8Aq6Aq6:J'KPQSTvXY[\W]N^bcefagjkmnio`pMq&r ssss1a&1a&)9QFQF;K(LQRTUPVYZ\]X^O_cdfgbhklnojpaqNr's!tttt 66!9*QFQF+;<<<<aV$A66!9Aq!Q////66!9!Q!Q 0Aq6Aq62BC1vPQSTvFVZ[]^Y_bcefagXhEijkkkk 
!Q	B88y   !qc
aS1#J/0000AQ1:D99q!Qi[!UDIIaL$))A,		RS+U%VVVV99q!Qi^$	1a$))Aq/499Q?;	1a$))Aq/499Q?;	1a$))Aq/499Q?;/ )    q6#a&=D99q!fX%QAAs1v(F"GGGG99Q!Q	"eSVCF]CF3q6M,J&KKKK99Q!Q	"ec!fWSV^SVCF],K&LLLL99q!fY!Q)USVCF]SVGCPQFN4SVYZ[V\]`ab]cVcehijeklopqlrerUs3t-uuuu 99Q<Q$$Q'499q!Qi\1+Eg+NNNN99Q<Q$$Q'499q!Qi\1+Eg+NNNN99q!Qi[Aq!9,/5noqrtumv9wmvhiZ[]^`aYb:cYbTU99Q<;L;LQ;OYb:cmv9w3xxxx &&!Q!
"Cc/0000&1a&""""a4++fq!fX&
'Cc/0000&1QT61+'''' ;d9ws   `2((`-`2-`2c            	         a  \        ^^\        3\        ^^3^
^\        ^,          ^,           334      o \        \        V 3R l4       \        S R4      ! \        4      \        RR^ ^ ^\        ,          334      8X  g   Q hR# )r   c                     < S P                   # r{   )nonexistantattributer   s   ru   r   test_getattr.<locals>.<lambda>  s    1#9#9rt   r   Nr  r   )r5   rc   rd   rR   rV  getattrr   s   @ru   test_getattrr    sj    Aq	Aq!9r1adQh&789A
>9:1fa FIy1a1++N$OOOOrt   c                  T   \        . RO. RO. RO.4      p V P                  '       g   Q hV P                  p V P                  '       g   Q h^V R&   V P                  RJ g   Q h\        . RO. RO. RO.4      p V P                  '       d   Q h\	        ^^4      p V P                  '       g   Q hR# )r   FN)r   r   r   )r   r   r   r   )r   r  r  )r5   is_upper_hessenbergr  is_lower_hessenbergrB   r   s    ru   test_hessenbergr    s    	9i01A     	A     AeH  E)))	9i01A$$$$$aA     rt   c            	      r   \        \        R  4       \        \        R 4       \        \        R 4       \        \        R 4       \        \        R 4       \        ^\        ,           ^ 3R34      P                  RR7      \        \        ^\        ,           4      ^ .^ ^..4      8X  g   Q h\        R4      p V P                  RR7      pV\        ^^ .^^\        ^4      ,          \        ,          ..4      8X  g   Q hWP                  ,          V 8X  g   Q h\        R4      p V P                  4       pWP                  ,          V 8X  g   Q hVP                  '       g   Q hV\        . RO. RO. RO.4      8X  g   Q h\        ^R\        ,          ^^\        ,          ,           3^\        ,          ^R\        ,           3^^\        ,          ,
          R\        ,
          ^334      p V P                  4       P                  4       \        R\        ^^ 3^\        ,
          ^ ^334      8X  g   Q h\        \        R 4       \        \        R 4       \        \        R	 4       \        \        R
 4       \        \        R 4       \        ^\        ,           ^ 3R34      P                  RR7      \        \        ^\        ,           4      ^ .^ ^..4      8X  g   Q h\        R4      p V P                  RR7      pV\        ^^ .^^\        ^4      ,          \        ,          ..4      8X  g   Q hWP                  ,          V 8X  g   Q h\        R4      p V P                  4       pWP                  ,          V 8X  g   Q hVP                  '       g   Q hV\        . RO. RO. RO.4      8X  g   Q h\        ^R\        ,          ^^\        ,          ,           3^\        ,          ^R\        ,           3^^\        ,          ,
          R\        ,
          ^334      p V P                  4       \        R\        ^^ 3^\        ,
          ^ ^334      8X  g   Q hR# )c                  4    \        R4      P                  4       # r   r5   choleskyrm   rt   ru   r   test_cholesky.<locals>.<lambda>.  s    )@)@)Brt   c                  4    \        R4      P                  4       # r   r   )r   r   r  rm   rt   ru   r   r  /      v&67@@Brt   c                  R    \        ^\        ,           ^ 3R34      P                  4       # r   r   )r5   r
   r  rm   rt   ru   r   r  0  s    vAqz6&:;DDFrt   c                  4    \        R4      P                  4       # r   )r   r   )r   r   r  rm   rt   ru   r   r  1  r  rt   c                  8    \        R4      P                  RR7      # r   F	hermitianr  r  rm   rt   ru   r   r  2  s    v&67@@5@Qrt   Fr  c                  4    \        R4      P                  4       # r   r6   r  rm   rt   ru   r   r  A  s    f)=)F)F)Hrt   c                  4    \        R4      P                  4       # r  r  rm   rt   ru   r   r  B      |,<=FFHrt   c                  R    \        ^\        ,           ^ 3R34      P                  4       # r  )r6   r
   r  rm   rt   ru   r   r  C  s    |a!eQZ,@AJJLrt   c                  4    \        R4      P                  4       # r  r  rm   rt   ru   r   r  D  r  rt   c                  8    \        R4      P                  RR7      # r  r  rm   rt   ru   r   r  E  s    |,<=FFQVFWrt   Nr   r  ))   r  rQ  )r     r   )rQ  r   r  r  )r   r   r   r  )r  r   r   r  r  )rR   r-   r#  r5   r
   r  r   r  r  r   r6   )r   r7  s     ru   test_choleskyr  -  sg   
!BC
:BC
:FG
:BC
:QRAE1:v&'0050AV	a!ea1a&M" F# # # # A	

U
#AAAd1gIaK 012222SS5A::78A	

Ass7a<<::::	9j9::::BqD!ac'"QqS!R!V$4q1Q3wQ6KLMA::< FI1ay1q5!Q-+P$QQQQ
!HI
:HI
:LM
:HI
:WX!a%V,-666G6	a!ea1a&S" L# # # #%&A	

U
#AAAd1gIaK 012222SS5A::=>A	

Ass7a<<::::	9j9::::q"Q$AaC(1Q326*:Q1Wb1fb<QRSA::<69q!Qi!a%A"GHHHHrt   c            	      z   \        R RR7      p \        \        V 4      \        V 4      .4      p\	        VP                  ^4      4      ^8X  g   Q hVP                  ^
4      \        \        V 4      ^
,          \        V 4      ^
,          ,           \        ^^
4      4      8X  g   Q h\        ^\        ^^4      ..4      pVP                  4       \        ^\        ^	^4      ,           \        P                  4      8X  g   Q hVP                  \        4      \        V4      8X  g   Q hVP                  \        ) 4      \        V4      8X  g   Q h\        ^^.^^..4      pVP                  ^4      ^8X  g   Q hVP                  R
4      ^ 8X  g   Q hVP                  R4      ^8X  g   Q h\        ^
4      P                  ^4      \        ^
4      P                  R
4      u;8X  d	   ^8X  g   Q h Q hVP                  \        4      ^8X  g   Q h\        ^\        \        .V \        P                  \        ) ..4      pVP                  R4      \!        \        ^%^4      ^\#        \        4      ^,          ,          ,           \        ^,          ,           V ^,          ,           4      8X  g   Q h\        . RO^^\        ^^4      ..4      pVP                  ^4      \!        \        R^4      \!        R4      ^,          ,           4      8X  g   Q hVP                  R
4      \        P$                  J g   Q hVP                  R4      \!        R4      ^,          8X  g   Q h\        ^^.^^..4      p\        ^^.R
^..4      p\        ^ \&        ) .\&        ^ ..4      p\        ^^ .^ R..4      pW#WE.p\        RRR7      pR F  p\)        ^4      P                  V4      \        P$                  J g   Q hV FX  p	V FO  p
V	P                  V4      V
P                  V4      ,           W,           P                  V4      ,
          pV^ 8  d   KO  Q h	  KZ  	  W#WE3 FR  p\+        W|,          P                  V4      \#        V4      VP                  V4      ,          ,
          4      pV^ 8X  d   KR  Q h	  K  	  \        ^^^\&        ,          ,
          R.4      p\        \        P                  ^\&        ,          ^.4      p\        . RO4      p\        ^^\&        .4      p\        \-        R4      \        ^R4      ^.4      pWVVV.p\        RRR7      p^^RR
\        P.                  \        P0                  \        3 EF<  pV^ 8  d1   \        . RO4      P                  V4      \        P$                  J g   Q hV^8  dx   V Fq  p	V Fh  p
V	P                  V4      V
P                  V4      ,           W,           P                  V4      ,
          p\+        V^ 8  4      \        P2                  J d   Kh  Q h	  Ks  	  V^^RR
\        P.                  \        P0                  39   g   K  V FR  p	\+        Wy,          P                  V4      \#        V4      V	P                  V4      ,          ,
          4      pV^ 8X  d   KR  Q h	  EK?  	  \        ^^. RO4      pVP                  ^4      ^8X  g   Q hR	# )rc   Tr  	frobeniusfroi  iI3 alphag      Y@Nr  r   )r   r   r   r  )r  r   r  )r  r  r  r   )	r   r   r   r  r  r   r   r  r   )r   r5   r   r   r*   normr   r   r   r   r   maxr  r9   rd   r   r   r   r  r
   rB   r)   r   InfinityNegativeInfinitytrue)rc   r   r   r   r   rT  r7  r  orderr  r  difr   r_   r`   ra   rb   es                     ru   test_matrix_normr  U  s~    	sAAA AAFF1I!###66":SVRZ#a&"*4hq"oFFFF 	HQN#$%A668s2A.777766":Q662#;#a&    	AA A66!9>>66":??66+!###r7<<?c"gll2.3!3333366":?? 	Aq	Aqvvs+,-AFF5M8B?Qs1vqy[02q581a4?@A B A 	
Q8B?345A66!9Xc1-UA=>>>>66":66+$s)A+---
 	AA AAQ !AQB!Q !AAB !A	
qA7&EQx}}U#qvv---Avve}qvve}4ULL'(q!   AEG>>%0J./ 0C!8O8   $ 	1qs7B A!Q A|A1ayAhq#.23A	
q!QA7&EQB

A,>,>C19)$))%0AFF:::A:A66%=166%=8U+,C#C1H-777   Q2r1::q/A/ABB~~e4Z!&&-/ 1 2axx  D& 	q!23A66!9??rt   c                    aa \        R RR7      o\        ^4      p ^
V R&   \        ^^
4      V R&   V P                  4       ^d8X  g   Q hSV R&   V P                  4       \	        ^
\        S4      4      \        \        ^^
4      \        S4      4      ,          8X  g   Q h\        \        S4      \        S4      .\        S4      ) \        S4      ..4      pVP                  4       o\        ;QJ dl    VV3R l\        ^^4      \        P                  \        ^^
4      \        ^,          \        \        \        ^^4      ,          3 4       F  '       d   K   RMa	  RM]! VV3R l\        ^^4      \        P                  \        ^^
4      \        ^,          \        \        \        ^^4      ,          3 4       4      '       g   Q h\        . 4      P                  4       ^ 8X  g   Q hR# )	rc   Tr  c              3      <"   T F=  p\        R 4      P                  SP                  SV4      P                  4       4      x  K?  	  R# 5i)      ?N)r	   
epsilon_eqr  r  )r   valMcrc   s   & ru   r   (test_condition_number.<locals>.<genexpr>  s@      QO BEuRy##BGGAsO$9$9$;<<Os   AAFNr   r   r   r   )r   r9   r   condition_numberr   r   r   r5   r   r   r   r   r   r   )r   r   r  rc   s     @@ru   test_condition_numberr    sv   sAAAAdGq"oAdG3&&&AdG3r3q6?S!R#a&5Q#QQQQQQ CF7CF"345A	
			B3 Q	!Q!R"Q$Bx1~<MOQ333 Q	!Q!R"Q$Bx1~<MOQ Q Q Q Q ":&&(A---rt   c                      \        R4      p \        R4      pW R,          8X  g   Q hW R,          8w  d   Q hW8X  d   Q hW8w  g   Q hV ^
8w  g   Q hV ^
8X  d   Q h\        R4      p\        R4      pW#8X  g   Q hW#8w  d   Q hR# )r   N)r   rn  r   r   r  ))r  r   r   )r   r   r   )r   r   r   r   )rY  r  r  )r5   r6   )r   r   r   rT  s       ru   test_equalityr     s    01A01A$<<dG|v:6M67N7Bw; 	67A01A6M6v:vrt   c                      \        ^4      P                  \        . RO.4      4      \        . RO. RO. RO. RO.4      8X  g   Q hR# )r   N)r   r   r   rY  r  r  )r9   rs  r5   rm   rt   ru   test_col_joinr    s@    q6??69+./	 	  rt   c                     \        . RO.4      p \        R^4       F`  p. ROpVP                  V^4       \        \	        ^4      P                  W4      P                  ^ 4      P                  4       4      V8X  d   K`  Q h	  R# )r   Nr   r   r   rP  rY  )r5   r   insertrN   r9   ru  ri  r  )r4r~   rE  s      ru   test_row_insertr    sf    		B2q\	As1v((/33A6==?@AEEE rt   c                      \        . RO4      p \        R^4       F`  p. ROpVP                  V^4       \        \	        ^4      P                  W4      P                  ^ 4      P                  4       4      V8X  d   K`  Q h	  R# )r   Nr  rP  r   )r5   r   r  rN   rB   rx  rh  r  )c4r~   rE  s      ru   test_col_insertr
    sd    			B2q\	AuQx**1155a8??ABaGGG rt   c                  F   \        ^^.4      P                  4       \        \        ^^4      \        ^^4      .4      8X  g   Q h\        . RO4      P                  4       \        . RO4      8X  g   Q h\        . RO4      p V P                  R R7      \        . RO4      8X  g   Q hR# )r   c                 >    V P                   ! ^
RR7      P                  # )r   T)r   chop)r  is_zeror  s   &ru   r   !test_normalized.<locals>.<lambda>  s    2D199rt   r  Nr   )r   r   g0.++)r5   r  r   r   s    ru   test_normalizedr    s    1a&>$$&AA/01 1 1 )'')VI->>>> 	}A<<9  			  rt   c                  R    \        R  4      R8X  g   Q h\        R 4      R8X  g   Q hR# )c                  4    \        ^4      P                  4       # r  r9   print_nonzerorm   rt   ru   r   $test_print_nonzero.<locals>.<lambda>  s    3q6//1rt   z[X  ]
[ X ]
[  X]
c                  6    \        ^4      P                  R4      # )r   .r  rm   rt   ru   r   r    s    3q6//4rt   z[.  ]
[ . ]
[  .]
N)rO   rm   rt   ru   test_print_nonzeror    s8    12     45     rt   c                  D   \         P                  ! ^4      \        ^4      8X  g   Q h\         P                  ! ^4      \        ^4      8X  g   Q h\        ^^4      \        ^^^.^,          4      8X  g   Q h\        ^^ .^ ^..4      p \        ^ ^ .^ ^ ..4      p\         F  pVP                  ^4      pW8X  g   Q hV \        ^VR7      8X  g   Q h\        V4      V8X  g   Q hVP                  ^4      pW8X  g   Q hV\        ^VR7      8X  g   Q h\        V4      V8X  d   K  Q h	  R# )r   rC  N)r5   r9   rB   r<   r   r   )r~   re   r   r   s       ru   test_zeros_eyer  	  s   ::a=CF"""<<?eAh&&&1:1qc"f----AA AAA AGGAJvvCsO###Aw#~~IIaLvvE!%%%%Aw#~~ rt   c                     \        4       P                  '       g   Q h\        ^ ^ .^ ^ ..4      P                  '       g   Q h\        ^^4      P                  '       g   Q h\        ^4      P                  '       d   Q h\        \        ^ .^ ^ ..4      P                  R8X  g   Q h\        \        ^ .^ ^ ..4      P                  R8X  g   Q h\        \        ^ .^ ^ ..4      P                  R8X  g   Q h\        \        ^ .^ ^ ..4      P                  R8X  g   Q h\        \        ^.^ ^ ..4      P                  R8X  g   Q h\        RRR7      p \        V ^ .^ ^ ..4      P                  R8X  g   Q hR# )r   NFr_   T)nonzero)	r5   is_zero_matrixrB   r9   rc   r6   r3   r4   r   r  s    ru   test_is_zeror  	  sX   8"""""Aq6Aq6"#22222A;%%%%%1v$$$$$Aq6Aq6"#22d:::!Q!Q()88D@@@QFQF+,;;tCCC 1a&1a&!12AATIIIAq6Aq6"#22e;;;sD!AAq6Aq6"#22e;;;rt   c                  8   \         ^,          p \        V 4      p\        V ) 4      p\        V 4      p\        V ) 4      p\        V 4      p\        V ) 4      pW!,          \	        ^4      ,          \	        ^4      8X  g   Q hWC,          \	        ^4      ,          \	        ^4      8X  g   Q hWe,          \	        ^4      ,          \	        ^4      8X  g   Q hVP                  4       ^^\        V 4      ,          ,           8X  g   Q hVP                  4       ^^\        V 4      ,          ,           8X  g   Q hVP                  4       ^^\        V 4      ,          ,           8X  g   Q h\        ^ 4      \	        ^4      8X  g   Q h\        ^ 4      \	        ^4      8X  g   Q h\        ^ 4      \	        ^4      8X  g   Q h\        P                  ! . RO\         ^,          4      p\        P                  ! . RO\         ^,          4      p\        P                  ! . RO\         ^,          4      p	\        \         ) ^,          4      VP                  4       8X  g   Q h\        \         ) ^,          4      VP                  4       8X  g   Q h\        \         ) ^,          4      V	P                  4       8X  g   Q h\        \         5^,          4      VP                  4       8X  g   Q h\        \         5^,          4      VP                  4       8X  g   Q h\        \         5^,          4      V	P                  4       8X  g   Q hR# )r   NrY  r  r  )r   r@   r?   r>   r9   r\  r   r]   from_axis_angleto_rotation_matrixrG   rH   rI   )
thetar3_plusr3_minusr2_plusr2_minusr1_plusr1_minusq1q2q3s
             ru   test_rotation_matricesr,  #	  s`   qDEG% HG% HG% HCF"c!f,,,CF"c!f,,,CF"c!f,,, ==?a!CJ,....==?a!CJ,....==?a!CJ,.... Q<3q6!!!Q<3q6!!!Q<3q6!!! 
	#	#IrAv	6B		#	#IrAv	6B		#	#IrAv	6BrTAX""7"7"9999rTAX""7"7"9999rTAX""7"7"99992"b&;&;&====2"b&;&;&====2"b&;&;&====rt   c                      \        \        R 4      ^,          4      R8X  g   Q h\        \        R4      4      \        R4      8X  g   Q h\        \        R 4       \        \        R4      4      R8X  g   Q h\        \        R4      4      R8X  g   Q hR# )vectorz	vector[4]rb   c                  &    \        R 4      R,          # )rb   r  r.   rm   rt   ru   r   %test_DeferredVector.<locals>.<lambda>K	  s    ~c226rt   testzDeferredVector('test')N)strr/   r   rR   r  reprrm   rt   ru   test_DeferredVectorr4  H	  sx    ~h'*+{:::>#&'>#+>>>>
:67~c"#s***v&'+CCCCrt   c                  >    \        \        R 4      4      '       d   Q hR# )r  N)rP   r/   rm   rt   ru    test_DeferredVector_not_iterabler6  O	  s    s+,,,,,rt   c                  (    \        \        R  4       R# )c                  *    \        \        R 4      4      # )r  )r5   r/   rm   rt   ru   r   ,test_DeferredVector_Matrix.<locals>.<lambda>S	  s    f^C%89rt   N)rR   r  rm   rt   ru   test_DeferredVector_Matrixr:  R	  s    
99:rt   c            
      &   \         p \        ^^^^.4      p\        ^^^^.4      p\        W.4      \        ^^^^.4      \        ^^V ! ^4      ^,          V ! R4      ^,          .4      .8X  g   Q h\        VP                  VP                  .4      \        ^^^^.4      \        ^^V ! ^4      ^,          V ! R4      ^,          .4      .8X  g   Q h\        \        ^^.4      \        ^^.4      .R4      \        ^\	        ^
4      ,          ^
,          \	        ^
4      ^
,          .4      \        \	        ^
4      ) ^
,          ^\	        ^
4      ,          ^
,          .4      .8X  g   Q h\        \        ^.4      4      p\        V4      \        ^..4      .8X  g   Q hR# )r   TNr  )r   r5   r2   r  r   r'   )r$  r  r  r7  s       ru   test_GramSchmidtr<  U	  s{   A	1q!f	B	1q!f	Bx 	1q!f	vaQqT!VQrU1W,=>?@ @ @bdd|$	1q!f	vaQqT!VQrU1W,=>?@ @ @ 1v1v7>$r(
2tBx{+,b	"aRjm,-C/ / / / 	&!+Aq>fqcUm_,,,rt   c                  ^    \        . RO^4      ^ 8X  g   Q h\        . RO^RR7      ^ 8X  g   Q hR# )r   F)zeroNr  )r7   rm   rt   ru   test_casoratianr?  f	  s.    lA&!+++lAE2a777rt   c                     \        4       \        ^ ^4      ,          P                  R8X  g   Q h\        ^^ 4      \        ^ ^4      ,          \        ^^4      8X  g   Q h\        ^ ^4      \        ^^ 4      ,          \        4       8X  g   Q hR# )r   N)r   r   )r5   rB   r   rm   rt   ru   test_zero_dimension_multiplyrA  k	  si    HU1a[ ''6111A;uQ{"eAqk111A;uQ{"fh...rt   c                    a  \        ^^\        ^4      4      o S R,          \        ^^..4      8X  g   Q hS R,          \        ^^..4      8X  g   Q hS R,          \        ^^..4      P                  8X  g   Q hS R,          \        ^^..4      P                  8X  g   Q h\        \        V 3R l4       \        \        V 3R l4       R# )r   c                     < S R,          # )r   )r   r   rm   r   s   ru   r   'test_slice_issue_2884.<locals>.<lambda>w	  	    qwrt   c                     < S R,          # )r   r  rm   r   s   ru   r   rD  x	  rE  rt   Nr   )r  r   r   r  )r5   r   r  rR   r  r   s   @ru   test_slice_issue_2884rG  q	  s    q!U1XAT7fq!fX&&&&U8v1vh''''T7fq!fX&(((((U8v1vh')))))
:'
:'rt   c                      \        ^ ^4      R,          P                  R8X  g   Q h\        ^^ 4      R,          \        ^^ . 4      8X  g   Q hR# )r   Nr  r   r  )rB   r   r5   rm   rt   ru   test_slice_issue_3401rI  {	  sC    A;u##v---A;tq!R 0000rt   c                     \        ^^4      p ^V ^&   V R,          \        . RO4      8X  g   Q hV ^,          ^8X  g   Q hV R,          ^.8X  g   Q h^*V R&   V R,          ^*8X  g   Q hV R,          \        ^*^ ..4      8X  g   Q h\        ^^..4      V R&   V R,          \        . RO.4      8X  g   Q h^*^+..V R&   V R,          \        . R	O.4      8X  g   Q h^V R
&   V R,          \        . RO4      8X  g   Q h. ROV R
&   V R,          \        . RO4      8X  g   Q h\        . RO4      V R
&   V R,          \        . RO4      8X  g   Q h\        . RO4      V R
&   V R,          \        . RO4      8X  g   Q hR# )r   :r   r   NNr  r  r   )r   :r   NNr   )r   r   r   )r   *   +   r   )r   :Nr   N)   r   r   rr  )rB   r5   r6   )r$  s    ru   test_copyinrN  	  sd   aAAaDT7fY''''Q4199T7qc>>AdGT7b==U8vAwi((((1vhAeHT7fi[))))RzAeHT7fk]++++AdGU8vj))))AdGT7fY''''YAdGT7fY''''9%AdGT7fY''''rt   c                  V  a  \        ^^.^^..4      P                  4       \        ^^.^ ^ ..4      R38X  g   Q h\        \        R 4       \        . R	O\        ^^.^\        R..4      o \        S P                  4       ^,          4      S P                  8w  g   Q hS P                  4       ^ ,          \        ^4      8w  g   Q hS P                  \        R7      ^ ,          \        ^4      8w  g   Q h\        \        V 3R l4       \        \        V 3R l4       \        \        V 3R l4       R# )
r   c                  @    \        ^^.^^..4      P                  4       # rb  ro  rm   rt   ru   r   'test_invertible_check.<locals>.<lambda>	  s    v1v1v&67;;=rt   r  c                  (   < S P                  R R7      # )r  r  r,  r   s   ru   r   rQ  	  s    quuEu2rt   c                  (   < S P                  R R7      # )r  r  rS  r   s   ru   r   rQ  	      quuDu1rt   c                  (   < S P                  R R7      # )r  r  rS  r   s   ru   r   rQ  	  rU  rt   Nr   r  )r  r  r   )	r5   rT  rR   r#  rc   r_  r   r9   r(   r   s   @ru   test_invertible_checkrW  	  s     Aq6Aq6"#((*v1v1v6F/G.NNNN
:=>
a!
a 	A
 qvvx{qvv%%% 668A;#a&   6686$Q'3q6111
:23
:12
:12rt   c                      \        R 4      w  rW,          pVP                  V \        . RO4      4      \        . RO4      V,          8X  g   Q hR# )zx, yN)r   r   r   )r   r  r5   )rc   rd   r  s      ru   test_issue_3959rY  	  s<    6?DA	A66!VI&'6)+<Q+>>>>rt   c                  H    \        \        ^^.^^..4      4      R8X  g   Q hR# )r   zMatrix([[1, 2], [3, 4]])N)r2  r5   rm   rt   ru   test_issue_5964r[  	  s(    v1v1v&'(,FFFFrt   c                      \        R 4      w  r\        \        V ^V,          .V^,          V ^,           ..4      4      R8X  g   Q hR# )x yz'Matrix([
[   x,   2*y],
[y**2, x + 3]])N)r   r&   r5   )rc   rd   s     ru   test_issue_7604r^  	  sE    5>DAAaC1a4Q-01234 4 4rt   c                  |   \        ^4      P                  '       g   Q h\        ^4      P                  4       P                  '       g   Q h\        ^4      P                  '       d   Q h\	        ^4      P                  '       d   Q h\        . RO.4      P                  '       d   Q h\        ^^R^R^R^/4      P                  '       g   Q h\        ^^\        ^4      4      P                  '       d   Q h\        ^^R^R^/4      P                  '       d   Q h\        ^^R^R^R^R^R^/4      P                  '       d   Q hR# )r   NrY  r   r   r  r   rs  )r9   is_Identityr  rB   r<   r5   r6   r   rm   rt   ru   test_is_Identityra  	  s   q6q6 ,,,,,Qx#####Aw"""""yk".....!eAuQa89EEEEEAaq*66666Aa%5!34@@@@@Aa%5E!U1eA!NO[[[[[[rt   c                     \        ^^4      P                  \        ^^4      4      ^8X  g   Q h\        ^^4      P                  . RO4      ^8X  g   Q h\        . RO4      P                  \        . RO4      4      ^8X  g   Q h\        ^^^\        ,          .4      P                  \        \        ^^\        ,          .4      4      R\        ,           8X  g   Q h\        ^^^\        ,          .4      P                  \        \        ^^\        ,          .4      RR7      R\        ,           8X  g   Q h\        ^^^\        ,          .4      P                  \        \        ^^\        ,          .4      RR7      ^\        ,           8X  g   Q h\        ^^^\        ,          .4      P                  \        \        ^^\        ,          .4      RRR7      ^\        ,
          8X  g   Q h\        ^^^\        ,          .4      P                  \        ^^\        ,          ^.4      RRR7      ^^\        ,          ,           8X  g   Q h\        ^^^\        ,          .4      P                  \        ^^\        ,          ^.4      RRR7      ^^\        ,          ,
          8X  g   Q h\        \        ^\        ,          .4      P                  \        \        ^\        ,          .4      RRR7      R8X  g   Q h\        \        ^\        ,          .4      P                  \        \        ^\        ,          .4      RR7      ^8X  g   Q h\	        \
        R	 4       R
# )r   Fr  Tphysicsr  conjugate_conventionrightleft)re  c                  V    \        ^^.4      P                  \        ^^.4      RRR7      # )r   Tr1  rd  r  rm   rt   ru   r   test_dot.<locals>.<lambda>	  s'    vq!f~11&!Q.Dgm1nrt   Nrr  r   rQ  )r<   r  r5   r
   rR   r#  rm   rt   ru   test_dotrj  	  s   1:>>$q!*%***1:>>)$))))  	!23r9991a1+""61a1+#6726AAA1a1+""61a1+#6%"HBQRFRRR1a1+""61a1+#6$"G2PQ6QQQ1a1+""61a1+#6$]f"gkmpqkqqqq1a1+""61ac1+#6$]d"eijmnopmpipppp1a1+""61ac1+#6$]c"dhilmnolohoooo1ac(1Q3x 0EX^_ceeee1ac(1Q3x 0vNRSSSS
:nort   c            	         \        R RR7      w  rr#rE\        ^ W4V3V) ^ W!) 3V) V) ^ V 3V) W) ^ 334      p\        ^ V ) V) V) 3V ^ WT) 3W) ^ V3W$V) ^ 334      pVP                  4       P                  V4      '       g   Q h\	        ^4      P                  4       P                  \        ^4      4      '       g   Q hVP                  4       P                  4       P                  V) 4      '       g   Q hR# )zB_x B_y B_z E_x E_y E_zTr  N)r   r5   dualequalsr9   rB   )B_xB_yB_zE_xE_yE_zFFds           ru   	test_dualrv  	  s   #*!$. CcC 
!c4 
tC 
Ta 	 	A 
sdSD3$	S$	dq3	SDQ	 
B 668??2q6;;=a))))668==?!!1"%%%%rt   c                     \        ^^.4      P                  4       RJ g   Q h\        ^^^ \        ^,          ^\        ,          ,           ^,           \        \        ^,           ^,          ) ^ \        \        ,          \        ) \        ) \        ,          ^ .	4      p V P                  4       RJ g   Q hV P                  RR7      e   Q hV P                  R R7      e   Q hV R,          ) V R&   V P                  4       e   Q hV R,          ) V R&   V P	                  4       p V P                  RR7      RJ g   Q h^V R&   V P                  4       RJ g   Q hR# )r   FTr  Nc                     V # r{   rm   r  s   &ru   r   %test_anti_symmetric.<locals>.<lambda>	  s    !rt   r   r   )r5   is_anti_symmetricrc   rd   r   r   s    ru   test_anti_symmetricr{  	  s/   1a&>++-666q!aA!aa!eaZKAaC!aRT1MNA D(((.6664<<< whAdG (((whAdG	
A.$666AdG E)))rt   c                     \        R4      p V P                  RR7      w  rWP                  ,          VP                  V,          u;8X  d   \        VP                  4      8X  g   Q h Q hV P                  RRR7      w  rWP                  ,          VP                  V,          u;8X  d   \        VP                  4      8X  g   Q h Q hW,          VP                  4       ,          V 8X  g   Q hR# )r   T)	normalize)r}  sortN)r   r   )r5   r  r  r9   r   r,  )r   r.  rY   s      ru   #test_normalize_sort_diogonalizationr  	  s     A==4=(DASS5ACCE(S[(((((==4d=3DASS5ACCE(S[(((((3quuw;!rt   c                  (    \        \        R  4       R# )c                  6    \        . RO\        ^ ^. 4      .4      # r  r   rm   rt   ru   r   !test_issue_5321.<locals>.<lambda>	
  s    vy&Ar2B&CDrt   N)rR   r#  rm   rt   ru   test_issue_5321r  
  s    
:DErt   c            	         \         P                  ! \        ^4      ^\        ^4      ,          4      \        . RO. RO.4      8X  g   Q h\         P                  ! \        ^4      ^\        ^4      ,          4      \        ^^ .^ ^.^^ .^ ^..4      8X  g   Q h\        p V P                  V ! \        ^4      4      V ! ^\        ^4      ,          4      4      \        . RO. RO.4      8X  g   Q hR# )r   N)r   r   r   r   )r   r   r   r   )r5   hstackr9   vstackr6   rC  s    ru   test_issue_5320r  
  s    ==Q3q6*f6 /    ==Q3q6*f	
A	
A	
A	
A	6 /    C::c#a&k3qQx=1V= 6   rt   c                      \        ^..4      p \        V 4      p\         P                  ! W4      \        ^^..4      8X  g   Q h\         P                  ! W4      \        ^.^..4      8X  g   Q hR# r   )r5   r   r  r  )r   AIms     ru   test_issue_11944r  
  s^    uA
!*C== FQF8$4444== FQC!:$6666rt   c                    a . ROp . ROp\        . RO4      pVP                  pV3R lp\         F  oS! V 4      pS! V4      pV! VP                  V4      V4       V! VP                  VP                  4      V4       V! VP                  P                  VP                  4      V4       V! VP                  P                  V4      V4       K  	  \	        \
        R 4       R# )r   c                 <   < W8X  g   Q h\        V 4      S8X  g   Q hR # r{   )r   )r   r;  r   s   &&ru   r1  test_cross.<locals>.test)
  s    xxAw#~~rt   c                  X    \        ^^^^.4      P                  \        ^^^^.4      4      # rb  r  rm   rt   ru   r   test_cross.<locals>.<lambda>3
  s*    q!aV""6!QA#78rt   Nr   r   )r  r   r  )r5   r  r   r  rR   r+   )r_   r`   ri  rh  r1  r   r   r   s          @ru   
test_crossr  #
  s    AA

C
%%C FFQWWQZQWWQSS\3QSSYYqss^S!QSSYYq\3  : 9 :rt   c                  ^   \        \        \        \        .4      p \        \        \
        \        .4      pV P                  4       V,          V P                  V4      8X  g   Q hV P                  4       P                  4       '       g   Q hV P                  4       P                  4       V 8X  g   Q hR # r{   )r5   rc   rd   re   r_   r`   ra   hatr  rz  vee)v1v2s     ru   test_hat_veer  6
  sw    	Aq		B	Aq		B668b=BHHRL(((668%%''''668<<>Rrt   c                  L   \         RR  F[  p V P                  ^4      V P                  ^4      0p\        V4      ^8X  d'   VP                  4       V P                  ^4      8X  d   K[  Q h	  \         R,           F*  p \	        V P                  ^4      \
        4      '       g   K*  Q h	  R# )r   Nr  r  )r   r9   r_  popr   r   )r   r$  s     ru   	test_hashr  =
  s    rs|WWQZ$1v{quuw#''!*444  r{{cggaj(3333 rt   c                  J   a  \         ^ ,          o \        \        V 3R l4       R# )r   c                  8   < \        S P                  ^4      4      # rb  )hashr9   rC  s   ru   r   !test_issue_3979.<locals>.<lambda>K
  s    4
#3rt   N)r   rR   rV  rC  s   @ru   test_issue_3979r  F
  s     !*C
>34rt   c                      ^ \         .^^ ..p \        ^ ^.\         ) ^ ..4      p\         F  pW! V 4      P                  4       8X  d   K  Q h	  R# r  )r
   r5   r   adjoint)r  r;  r   s      ru   test_adjointr  N
  sR    q6Aq6
C
1a&A2q'"
#Cc#h&&(((( rt   c                  h   ^ RI p V P                  P                  P                  P	                  R4      pV P                  P                  P                  V4      p^ \        V,          .^ V..p\        ^ ^ .\        ) V,          V..4      p\         F  pWE! V4      P                  4       8X  d   K  Q h	  R# )r   Nr_   )
sympy.physics.quantumrc  quantumoperatorr!   r"   r
   r5   r   r  )sympyr_   a_dagr  r;  r   s         ru   test_adjoint_with_operatorr  U
  s     &&//4AMM!!((+Eq1u:5z
"C
1a&A2:q/*
+Cc#h&&(((( rt   c                      \        \        \        \        4      ^,          \	        \        4      ^,          ,           ..4      4      \        ^..4      8X  g   Q hR# r  )r)   r3   r   rc   r   rm   rt   ru   test_simplify_immutabler  `
  sE    Oc!fai#a&!)&;%<$=>?#aSE*+ + +rt   c                     aa \        R \        R7      w  oo\        ^^V3R l4      p \        ^^V3R l4      pVP                  SS4      pW 8X  g   Q hR# )zF, GrC  c                     < S! W,           4      # r{   rm   )r~   r   Gs   &&ru   r   test_replace.<locals>.<lambda>f
      !AC&rt   c                     < S! W,           4      # r{   rm   )r~   r   rt  s   &&ru   r   r  g
  r  rt   N)r   r   r5   replace)Kr   Nrt  r  s      @@ru   test_replacer  d
  sI    6x(DAqq!()Aq!()A			!QA6M6rt   c                     \        ^^.\        ^^\        ,          ,
          ..4      p V P                  4       \        P                  \        ^4      \        P
                  \        08X  g   Q hV P                  \        4      \        08X  g   Q hR# r   )r5   rc   atomsr   r  NegativeOner   r   s    ru   
test_atomsr  l
  sb    AA!G%&A779qtAMM15555776?qc!!!rt   c            	        a	 \        \        \        .\        \        ..4      p \        V P                  R R7      4      \        V P                  4       4      8X  g   Q h\        ^^h.R^.R^..4      \        . RO. RO.4      \        \        \        .4      .pV F{  pVP                  R R7      pW#,          pW2,          p\        WB,          4      V8X  g   Q h\        WS,          4      V8X  g   Q hVP                  V8X  g   Q hVP                  V8X  d   K{  Q h	  V F  p\        VP                  RR7      4      pW#,          pW2,          p\        WB,          4      V8X  g   Q h\        WS,          4      V8X  g   Q hVP                  V8X  g   Q hVP                  V8X  d   K  Q h	  ^ RI	H
o	 V P                  RR7      pV P                  4       p\        R\        R\        R\        R/p\        ;QJ dG    V	3R l\        VP                  V4      VP                  V4      4       4       F  '       d   K   R	M<	  R
M8! V	3R l\        VP                  V4      VP                  V4      4       4       4      '       g   Q hR# )RDr  i  ED)compib,  i  ij  c              3   r   <"   T F,  w  rS! VP                  4       VP                  4       4      x  K.  	  R # 5ir{   r   )r   r~   r   r  s   &  ru   r   test_pinv.<locals>.<genexpr>
  s0      
3DA 	QSSUACCE3s   47FTNr   )r   r   r  )r  rM     i_)r5   r_   r`   ra   rb   r)   pinvr,  rP  sympy.core.numbersr  r   zipr  )
A1Asr   A_pinvAApApAr-  wrepsr  s
            @ru   	test_pinvr  r
  s	   	!Q!Q 	!BBGG4G()Xbffh-???? 2s)dAYQ0
1
)\*
+
!Q.
B t$jj A%%%%///uu||uu||  !&&&-.jj A%%%%///uu||uu||  (
tA
Avq%E1e4D3 
taffTl3
333 
taffTl3
 
 
 
 
rt   c            	          \        . RO. RO. RO. RO. R	O. R
O.4      .p V  FV  pVP                  RR7      pW,          pW!,          pVP                  V8X  g   Q hR pV! VP                  V4      '       d   KV  Q h	  R# )=   r  r  c                 &   V P                  \        4      pV Uu/ uF  q3VP                  4       bK  	  ppW,
          P                  V4      P                  4       p\        ;QJ d    R  V 4       F  '       d   K   R# 	  R# ! R  V 4       4      # u upi )c              3   >   "   T F  p\        V4      R 8  x  K  	  R# 5i)g|=Nr   r   r  s   & ru   r   Xtest_pinv_rank_deficient_when_diagonalization_fails.<locals>.allclose.<locals>.<genexpr>
  s     ;{!s1v~{s   FT)r  r%   r  rF  r   )M1r  rootofsrc  rootofs_approxdiff_approxs   &&    ru   allcloseEtest_pinv_rank_deficient_when_diagonalization_fails.<locals>.allclose
  su    hhv&G4;<GqlGN<7,,^<BBDK3;{;33;3;3;{;;; =s   BN)r  r  7   r  r,  r   )>   `   U   r  r  r   )E   8   rM  r   6   r   )r   r  [   )   r,  r   )r   r  r   r)  Z   r   )r   r   r   r   r   r   )r5   r  rP  )r  r   r  r  r  r  s         ru   3test_pinv_rank_deficient_when_diagonalization_failsr  
  s    
 	#####  	!
B t$jjuu||	< s####- rt   c                      \        ^ ^4      \        ^ ^4      ,           \        ^ ^. 4      8X  g   Q h\        ^^ 4      \        ^^ 4      ,           \        ^^ . 4      8X  g   Q hR# r  )r<   r5   rm   rt   ru   test_issue_7201r  
  sR    1:Q
"fQ2&66661:Q
"fQ2&6666rt   c                      \         \        \        \        3 F)  p V ! \        .^ ..4      P
                  \        08X  d   K)  Q h	  R# r  )r3   r4   r5   r6   rc   free_symbolsr(  s    ru   test_free_symbolsr  
  s9    3V\I1#s}))aS000 Jrt   c                 *  a   ^ RI Ho  \	        S ! . RO4      4      \	        . RO4      8X  g   Q h\	        S ! . RO.4      4      \	        . RO.4      8X  g   Q h\	        S ! . RO. RO.4      4      \	        . RO. RO.4      8X  g   Q h\	        S ! \
        \        \        .4      4      \	        \
        \        \        .4      8X  g   Q h\        \        V 3R l4       \	        S ! ^^.4      S ! ^^.4      .4      \	        ^^.^^..4      8X  g   Q h\	        S ! ^^.4      ^^..4      \	        ^^.^^..4      8X  g   Q h\	        S ! . 4      S ! . 4      .4      \	        ^^ . 4      u;8X  d   \	        ^ ^ . 4      8w  g   Q h Q hR#   \         d    \        R4        ELi ; i)zSee issue 7465.array?NumPy must be available to test creating matrices from ndarraysc                  B   < \        S ! ^^.^^..^^.^^...4      4      # rb  r   r  s   ru   r   #test_from_ndarray.<locals>.<lambda>
  s.    u1v1v.!Q!Q0@ABCrt   Nr   rn  )
numpyr  ImportErrorrU   r5   rc   rd   re   rR   NotImplementedErrorr  s   @ru   test_from_ndarrayr  
  s{   P %	"#vi'8888%$%)<<<<%I./0	9%&' ' '%Aq	"#vq!Qi'8888
CE5!Q=%A-01VaVaV<L5MMMM5!Q=1a&)*fq!fq!f5E.FFFF59eBi()VAq"-=Q1bAQQQQQQ  PNOPs   E8 8FFc            	         ^ RI Hp   ^ RIHpHp V ! X! ^^.^^..4      4      pV^,          ^8X  g   Q h\        V4      . RO8X  g   Q h\        \        4      ;_uu_ 4        V ! X! ^^.^^..4      4      pRRR4       V^,          ^8X  g   Q h\        V4      . RO8X  g   Q hR#   \
         d    \        R4        Li ; i  + '       g   i     LR; i)r   
_matrixify)r  matrixzPNumPy must be available to test indexing matrixified NumPy ndarrays and matricesNr  )
sympy.matrices.commonr  r  r  r  r  rU   r   rQ   PendingDeprecationWarning)r  r  r  r   s       ru   test_17522_numpyr  
  s    0a' 	51a&1a&)*+AQ41997l"""	2	3	3v1v1v./0 
4Q41997l"""  a_`a 
4	3s   B" B>"B;:B;>C	c                      ^ RI Hp   ^ RIHp V ! X! ^^.^^..4      4      pV^,          R8X  g   Q h\        V4      . RO8X  g   Q hR#   \         d    \        R4        LRi ; i)r   r  )r  zEmpmath must be available to test indexing matrixified mpmath matrices      @N)r  g       @r  r  )r  r  mpmathr  r  rU   r   )r  r  r   s      ru   test_17522_mpmathr  
  sk    0V! 	6Aq6Aq6*+,AQ43;;7****  VTUV   A A! A!c                      ^ RI Hp   ^ RIHp V ! X! ^^.^^..4      4      pV^,          ^8X  g   Q h\        V4      . RO8X  g   Q hR#   \         d    \        R4        LRi ; i)r   r  )
csr_matrixzJSciPy must be available to test indexing matrixified SciPy sparse matricesNr  )r  r  scipy.sparser  r  rU   r   )r  r  r   s      ru   test_17522_scipyr  
  sj    0[+ 	:1v1v./0AQ41997l"""  [YZ[r  c                  r   \        ^\        .\        ) ^..4      p V P                  '       g   Q h^\        ,          V R&   V P                  RJ g   Q h\        V R&   V P                  e   Q hV R,          \        ,          V R&   V P                  RJ g   Q h\	        R4      p\        R4      p\        V..4      P                  RJ g   Q h\        W.\        V4      V..4      P                  RJ g   Q h\        W.W!..4      P                  RJ g   Q h\        W.\        V4      V..4      P                  RJ g   Q hR# )	r   FNr`   ra   Tr   r   r   )r5   r
   is_hermitianrc   r    r!   r"   r   r  s      ru   test_hermitianr    s&   A!Q !A>>>>cAdG>>U"""AdG>>!!!gaiAdG>>U"""#AAA3%=%%---A6F1Iq>*+88D@@@A6A6"#00E999A6IaL!,-.;;uDDDrt   c                      \        \        \        \        R R7      ..4      p V ^ ,          ^\        ,          8w  g   Q hV P                  ! 4       \        ^\        ,          ..4      8X  g   Q hR# )Fr  N)r5   r   rc   r  r  s    ru   	test_doitr    sO    Qq5)*+,AQ41Q3;;668v!ug&&&rt   c                    aa \        . RO. RO. RO.4      p V P                  ^4       V \        . RO. RO.4      8X  g   Q h\        . RO. RO. RO.4      pVP                  R	4       V\        . RO. RO.4      8X  g   Q h\        . RO. R
O. RO.4      pVP                  ! R4       V\        . RO. R
O.4      8X  g   Q h\        . RO. RO. RO.4      o\        \        V3R l4       \        . RO. RO. RO.4      o\        \        V3R l4       \        . RO. RO. RO.4      p V P	                  ^4       V \        ^^.^^.^^..4      8X  g   Q h\        . RO. RO. RO.4      pVP	                  R	4       V\        ^^.^^.^^..4      8X  g   Q h\        . RO. RO. RO.4      o\        \        V3R l4       \        . RO. RO. RO.4      o\        \        V3R l4       R# )r   c                  &   < S P                  ^
4      # r   row_delr.  s   ru   r   +test_issue_9457_9467_9876.<locals>.<lambda>&      qyy}rt   c                  (   < S P                   ! R4      # r   r   r	  rX   s   ru   r   r  (      qyy~rt   c                  &   < S P                  ^
4      # r   col_delr  s   ru   r   r  2  r  rt   c                  (   < S P                   ! R4      # r  r  rX   s   ru   r   r  4  r  rt   Nr   r  r   r  )r   r   r   )r  r   r  r  )r5   r
  rR   r  r  )r   r  r^   r.  rY   s      @@ru   test_issue_9457_9467_9876r    s   	9i01AIIaL	9-....	9i01AIIbM	9-....	9k23AIIbM	9-....	9i01A
:,-	9i01A
:-. 		9i01AIIaLAAA/0000	9i01AIIbMAAA/0000	9i01A
:,-	9i01A
:-.rt   c                     \        R RR7      w  r\        R4      w  r#\        ^4      p\        ^^WV\        .4      pW,          V,          W,          V,          8w  g   Q hW2,          V,          W#,          V,          8X  g   Q hW,          WP,          8w  g   Q hW%,          WR,          8X  g   Q hW,          V,          \        W,          ^ .^ W,          ..4      8X  g   Q hR# )r]  F)commutativeza bN)r   r9   r5   re   )rc   rd   r_   r`   r   r  s         ru   test_issue_9422r  6  s    5e,DA5>DAAA	1qQl	#B3q5ACE>>3q5ACE>>424<<424<<3q5FQS!Hq!#h/0000rt   c                  J   \        . 4      p R R.\        . RO4      3pRR.V^,          P                  3pRR .\        ^^.^^..4      3pWV3 FR  w  rEV FG  p\        W4      pRV9   d	   V! V4      MV! ^*V4      pW8X  d   \        V4      \        V4      8w  d   KG  Q h	  KT  	  R# )rx  rq  ru  rs  joinN)r  r   r   )r5   r  r  id)	r   r_   r`   ra   opsr   opr9  news	            ru   test_issue_10770r  B  s    r
A	z"F9$55A	z"AaDFF*A	|$fq!fq!f-=&>>A)BA B,!A$Ab!HC832a5 000  rt   c                     \        . RO. RO. RO.4      p V P                  . RO. RO4      \        ^^.^^.^^..4      8X  g   Q hV P                  . RO. RO4      \        ^.^.^..4      8X  g   Q hV P                  . RO. RO4      \        . RO.4      8X  g   Q hV P                  . RO. RO4      \        . RO. RO.4      8X  g   Q hV P                  . RO. R	O4      \        ^^ . 4      8X  g   Q hV P                  . R	O. RO4      \        ^ ^. 4      8X  g   Q hV P                  . RO. R
O4      \        ^.^..4      8X  g   Q hR# )r   Nr   rn  r  r   )TTF)TFF)TFT)FFF)FTF)r5   r   r   s    ru   test_issue_10658r!  N  s8   	9i01A99Y 34AAA'() ) )99Y 45!qcA39PPPP99)959LLLL99()4	9%&' ' '99Y 56&Ar:JJJJ99*I6&Ar:JJJJ99(*>?aSz  rt   c                     \        R^\        ^4      ,          ,           R.R\        ^4      ,          ^,          ^,           R\        ^4      ,          ^,          ..4      p V P                  4       ^8X  g   Q h\        ^^\        ^4      ,          \        ,          ,           R.^R^\        ^4      ,          \        ,          ,           ..4      p \	        V P                  4       ^ ,          \        ^R^^\        ^4      ,          \        ,          ,           ,          .^ ^ ..4      ,
          4      \        ^^4      8X  g   Q h\        R4      w  rr4rVrxrr\        WW,          W+,          ^ .W4W;,          WK,          ^ .WVW[,          Wk,          ^.WxW{,          W,          ^.W^V	,          V,          W,          ,
          ^V
,          V,          W,          ,
          ^ ..4      p V P                  4       ^8X  g   Q hR# )r   /a_x a_y b_x b_y c_x c_y d_x d_y e_x e_y t_0 t_1NrQ  r$  r   )r5   r   r  r
   r)   rT  rB   r   )r   axaybxbycxcydxdyexeyt0t1s                ru   !test_opportunistic_simplificationr0  [  s    	aQi$r$q'z!|a'7DGA&FGHA668q== 	1T!W9Q;#Qr!DG)A+~$678AAFFHQK&1b!aQik/.B*CaV)L"MMNRWXY[\R]]]] +22c*d'B"bB"ruRU1%rRU25&;R25q<QSUY[Y^_a_defRgikopqsostvovwyw|o|}~  @B  ~B  CE  ~E  FH  FK  ~K  LM  iN  O  	PA668q==rt   c                     \        . RO. RO.4      p V P                  4       ^ ,          \        . RO. RO.4      ,
          P                  4       R8  g   Q h\        . RO. RO. R	O.4      p\        . R
O. RO. RO.4      p\        . RO. RO. R	O.4      pVP                  4       V,
          P                  4       R8  g   Q hVP                  4       V,
          P                  4       R8  g   Q hR# )~jth?V瞯<N)r2  gRM@g(\M@)gDl)@gQgp=
cG@)r  r   r   )r   r  r  )q7K<r  r   )      r   r   r  )r4  r  r  )r5                 @)r6  r6  r  )r   r5  r7  )r  r4  g      )r5   rT  r  r,  )mmm_mixedm_floatm_invs       ru   test_partial_pivotingr<  k  s    
 $&;<=BGGIaL6#3#3#5 6 67;tvF F F o%%' (G ''') *G (((* +E
 KKME!'')E111KKME!'')E111rt   c                     \        . RO. RO. RO.4      p V P                  R R7      ^ ,          p\        . RO. R	O. R
O.4      pW,
          pVP                  4       R8  g   Q hVR,          ^ 8X  g   Q hR# )zWhen doing numerical computations, all elements that pass
the iszerofunc test should be set to numerically zero if they
aren't already. c                     \        V 4      R 8  # )g [n<r  r  s   &ru   r   *test_iszero_substitution.<locals>.<lambda>  s    Qrt   r  r3  N)?皙gɿr   )皙r@  gٿr   )rA  rB  g333333?r   )r  r   g4iҤIӿr   )r   r  glٲer   r  r  )r5   rT  r  )r   m_rref	m_correctm_diffs       ru   test_iszero_substitutionrF    sp     	$%9:NOPAVV5V6q9F9:[\}~IF;;=5   #;!rt   c            	         ^ RI Hp  ^\        \        \	        ^^-4      ,          4      ,          \        \        \	        ^^-4      ,          4      \        ^4      ,           ,          pR\        ^4      ,          \        \        \	        ^^-4      ,          4      ^,          ,          ^\        \        \	        ^^-4      ,          4      ,          ,           R\        \        \	        ^^-4      ,          4      ^,          ,           ,          pV ! ^ ^ 4      pV ! ^\        ^4      ) 4      pV ! W4      p\        V\        V4      ,
          V\        V4      ,
          .4      p\        W5,
          WE,
          .4      p\        \        W5,
          4      \        WE,
          4      .4      pR p	VP                  RV	R7      ^8X  g   Q hVP                  RV	R7      ^8X  g   Q hVP                  RV	R7      ^8X  g   Q hR# )r   )Pointc                 <    \        V P                  ! 4       4      R 8  # )g#B;)r   r   r  s   &ru   r   "test_issue_11238.<locals>.<lambda>  s    #acce*u$rt   T)r)   r  Nr%  r   )	sympy.geometry.pointrH  r   r   r   r   r5   r)   r  )
rH  xxyyp1p2p0r  r  m3Zs
             ru   test_issue_11238rS    sx   *	
3r(2r""#	#SHR,<)<%=Q%G	HB
T!W*SHR,,-q0
02c"Xb"=M:M6N3N
NQSVYZ\]efhjl]mZmVnpqVqQq	rB	q!B	q47(	B	rB	hrl"B"$56	7B	"'"	#B	"'"HRW$56	7B
 	%A77DQ7/144477DQ7/144477DQ7/1444rt   c                     \        ^^. RO4      p V \        P                  ,          pW,           p\         FQ  pV! V4      P	                  4       w  rE\        V4      \        V 4      8X  g   Q h\        V4      \        V 4      8X  d   KQ  Q h	  R# )r   Nr  )r5   r   ImaginaryUnitr   as_real_imagr   )r  r  rQ  klsr_   r`   s         ru   test_as_real_imagrX    so    	!I	B	AOO	B	B"g""$Aw$r("""Aw$r(""" rt   c                      \        ^^. RO4      p V P                  4       w  rV^,          \        ^^^.4      8X  g   Q hV^ ,          \        ^^. RO4      8X  g   Q hR# )r   Nr"  )r   r   r   r   )r5   jordan_cells)r   r.  Jcellss      ru   test_deprecatedr\    sY    
 	q!23A IA!9q!aS))))!9q!\2222rt   c                      ^ RI Hp  \        . RO4      p\        . RO4      pV ! V^4      \        . RO4      8X  g   Q hV ! V^4      \        . RO4      8X  g   Q hR# )r   )ModN)r  r   r   )r   r  r#  )r   r   r   )r   r   r   )sympy.core.modr^  r5   )r^  r   r   s      ru   test_issue_14489r`    sK    "zA}Aq!9y))))q!9y))))rt   c                       ^ RI Hp  \	        ^^.^^..4      pX ! V\
        R7      P                  P                  R8X  g   Q hR#   \         d    \        R4        LTi ; i)r   r  r  )dtypefloat64N)r  r  r  rU   r5   floatrb  name)r  r   s     ru   test_issue_14943rf    s`    P 	1!u~A% &&++y888	  PNOPs   A AAc                  r    \        R ^^4      p \        R4      pV R,          ^ 8  p\        V4      R8X  g   Q hR# )r   r_   zm[0, 0] > 0Nr   )rE   r   r2  )r   r_   s     ru   test_case_6913rh    s8    S!QAsA	$	Aq6]"""rt   c                  h    \        R ^^4      p \        R4      pV P                  V4      W/8X  g   Q hR# )r   r_   N)rE   rL   match)r   r_   s     ru   test_issue_11948rk    s1    S!QAS	A771:!rt   c                     \        ^\        .4      \        ^\        ) .4      .p \         P                  ! V !  \        ^.\        ..4      \        ^.\        ) ..4      .8X  g   Q h\        ^\        ^ .4      \        \        ^ \        ) .4      .p \         P                  ! V !  \        ^.\        .^ ..4      \        \        ^,          .\        ^4      ^,          .\        ) ..4      .8X  g   Q h\        ^\        .^\        ) ..4      pVP	                  4       w  r#W"P
                  ,          \         P                  ! ^4      8X  g   Q hR# r   )r5   r
   orthogonalizer   QRdecompositionrP  r9   )vecsr   rY   r$  s       ru   test_gramschmidt_conjugate_dotrp    s#   Aq6NFAr7O,D&	!qc
	VaSA2$K012 2 2 Aq!9vq!aRj12D&	!qcA3	 &1Q3%!A$q&QB4)@"ABC C C 1a&1qb'"
#C DAss7fjjm###rt   c                      \        \        R ^^4      4      p \        \        R^^4      4      pV P                  ! V4      p\        W R,          4      p\        W0R,          4      pW1R,          8X  g   Q hV^ 8X  g   Q hR# )r_   r`   Nr   )r5   rE   r  r   )r_   r`   ra   rb   r  s        ru   test_issue_8207rr    sj    |CA&'A|CA&'A	aAQ$AQ$A$<<6M6rt   c            	      V  a ^ RI Hp  \        ^^.^ ^..4      oSP                  \	        \
        \        ,          4      \
        4      \        \	        \        4      \	        ^\        ,          4      \	        \        4      ,
          .^ \	        ^\        ,          4      ..4      8X  g   Q h\        ^^.^^..4      o\        S,          ^,          P                  \        \
        4      \
        4      \        \        ^4      ^,          \        ^4      ) ^,          .\        ^4      ) ^,          \        ^4      ^,          ..4      8X  g   Q h\        \        R 4       \        \        V3R l4       \        . RO. RO. RO. RO.4      oSP                  \        \
        4      \
        4      SP                  4       8X  g   Q h\        \        V3R l4       \        ^)^.^^"..4      o\        SP                  \        \
        4      \
        4      ^,          4      S8X  g   Q h\        . R	O. R
O. RO.4      o\        SP                  V ! \
        ^4      \
        4      ^,          4      S8X  g   Q h\        . RO. RO. RO. RO.4      oSP                  \        \
        4      \
        4      SP                  4       8X  g   Q h\        . RO. RO. RO. RO.4      oSP                  \        \
        \        ,          4      \
        4      \        \        S\        ,          P                  4       4      4      8X  g   Q hR# )r   )nthrootc                  Z    \        ^4      P                  \        \        4      \        4      # r  )rB   analytic_funcr   rc   rm   rt   ru   r   test_func.<locals>.<lambda>  s    a 6 6s1vq Art   c                  `   < S \         ,          P                  \        \         4      \         4      # r{   )rc   rv  r   r   s   ru   r   rw    s    1 3 3CFA >rt   c                  J   < S P                  \        \        4      \        4      # r{   )rv  r   rc   r   s   ru   r   rw    s    Q ;rt   N)r   r  r  r   )r   r   r   r  )r   r   r  r   )r   r  r   )r  r   r  )r  r   r  r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   r  )sympy.simplify.simplifyrt  r5   rv  r   rc   rf   r   r   r   rR   r#  r   r)   r   )rt  r   s    @ru   	test_funcr{    sD   /A1vA??3qs8Q'6CFC!Hs1v<M3NQRTWXYZ[X[T\P]2^+____A1vAFQJ%%c!fa0FT!WQYaQR
<SW[\]W^V^_`V`bfghbijkbkUl;m4nnnn :AB
:>?}MNA??3q61%000
:;<R"b"#AAOODGQ/23q888[+67AAOOGAqM15q89Q>>>lL,GHA??3q61%000lL,GHA??3qs8Q'6(AaC99;2G+HHHHrt   z$Cannot create threads under pyodide.c            	      D   R  p \        R4      ;_uu_ 4        \        P                  P                  4       ;_uu_ 4       pVP	                  V 4      pVP                  4       '       g   Q h RRR4       RRR4       R#   + '       g   i     L; i  + '       g   i     R# ; i)c                  ^    \         P                  R 8X  g   Q h\        ^..4      p W ,          p R# )NT)rJ   stater5   r   s    ru   r9  test_issue_19809.<locals>.f  s-    !''4///QC5MErt   TN)rF   
concurrentfuturesThreadPoolExecutorsubmitresult)r9  executorfutures      ru   test_issue_19809r  
  si     
T		2244__Q'F==??"? 5 
	44 
		s"   'B*A;	)B;BBB	c                      \        \        \        .4      p \        V \        ^ ^3\        ^ ^34      \        \        P
                  .\        P
                  ..4      8X  g   Q hR# r  )r5   rc   rd   r   r   r   r(  s    ru   test_issue_23276r    sT    1vAQAq	Aq!9-	
	
9 2   rt   c                  t   \        . RO. RO. RO. RO.4      p V P                  4       pV^ ,          \        . RO4      8X  g   Q hV^,          \        . RO4      8X  g   Q hV^,          \        . RO4      8X  g   Q h\        V4      ^8X  g   Q h\        P                  ! V .VO5!  P                  4       V8X  g   Q hR# )	r   Nr   r   r   r   r   r  rQ  r   r  r%  r   r   r   r   r   r   r   r   rv  r   )r   r  r   r   )r   rQ  r   r   )r   r  r   rv  )rl   rs  r5   r_  r  r   basiss     ru   test_columnspace_oner  !  s    00002 	3A
 MMOE8vm,,,,8vn----8vn----u:??==#U#//1U:::rt   c                      \        . RO. RO. RO. RO.4      p V P                  4       pV^ ,          \        . RO.4      8X  g   Q hV^,          \        . RO.4      8X  g   Q hV^,          \        . RO.4      8X  g   Q h\        V4      ^8X  g   Q hR# )r   Nr  r  r  r  )r   r  r   r   r   )r   r   r   r   r   )rl   rp  r5   r_  r  s     ru   test_rowspacer  0  s    00002 	3A
 JJLE8v/00008v/011118v/0000u:??rt   c                     \        . RO. RO. RO. R	O.4      p V P                  4       pV^ ,          \        . R
O4      8X  g   Q hV^,          \        . RO4      8X  g   Q h\        ;QJ d)    R W^ ,          ,           4       F  '       d   K   RM	  RM! R W^ ,          ,           4       4      '       g   Q h\        ;QJ d)    R W^,          ,           4       F  '       d   K   RM	  RM! R W^,          ,           4       4      '       g   Q hR# )r   c              3   8   "   T F  qP                   x  K  	  R # 5ir{   r  r  s   & ru   r   %test_nullspace_one.<locals>.<genexpr>H       -*Qyy*   FTc              3   8   "   T F  qP                   x  K  	  R # 5ir{   r  r  s   & ru   r   r  I  r  r  Nr  r  r  r  )r  r   r   r   r   )r  r  r   r  r   )rl   rw  r5   r   r  s     ru   test_nullspace_oner  >  s    00002 	3A
 KKME8v.////8v011113-!!H*-333-!!H*-----3-!!H*-333-!!H*-----rt   c                    a \        ^4      o\        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R	 l4       \        \        V3R
 l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       SP                  R^ ^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^R7      \	        . RO. RO. RO.4      8X  g   Q h\        ^^^ .^,          4      p V P                  ! R^^4      \	        ^^^ .^,          4      8X  g   Q hV P                  ! R^ ^4      \	        ^^^ .^,          4      8X  g   Q hV P                  ! R^ ^^4      \	        ^^^ .^,          4      8X  g   Q hR# )r   c                  &   < S P                  R 4      # abcelementary_row_opr  s   ru   r   test_row_op.<locals>.<lambda>P      q2259rt   c                  $   < S P                  4       # r{   r  r  s   ru   r   r  Q      q224rt   c                  ,   < S P                  R ^^R7      # n->knrh  r   r  r  s   ru   r   r  R      q227Q2Grt   c                  ,   < S P                  R R^R7      # )r  r  rQ  r  r  s   ru   r   r  S      q227a2Hrt   c                  ,   < S P                  R ^^R7      # n<->mrow1row2r  r  s   ru   r   r  T      q2272Krt   c                  ,   < S P                  R ^^R7      # r  r  r  s   ru   r   r  U  r  rt   c                  ,   < S P                  R R^R7      # r  r  rQ  r  r  s   ru   r   r  V      q227!2Lrt   c                  ,   < S P                  R ^RR7      # r  r  r  s   ru   r   r  W      q2272Lrt   c                  .   < S P                  R ^^^R7      # n->n+kmr  r  r   r  r  s   ru   r   r  X      q22911PQ2Rrt   c                  .   < S P                  R ^^^R7      # r  r  r  s   ru   r   r  Y  r  rt   c                  .   < S P                  R R^^R7      # r  r  rQ  r  r  s   ru   r   r  Z      q2292AQR2Srt   c                  .   < S P                  R ^R^R7      # r  r  r  s   ru   r   r  [      q22912QR2Srt   c                  .   < S P                  R ^^^R7      # r  r  r  s   ru   r   r  \  r  rt   r  r  )r  r   r  r  )rh  r  r  )rh  r   r  )r  r   r  Nr  r  r  rY  rZ  )r   r   r   )r   rR   r#  r  r5   rx   r_   r  s    @ru   test_row_opr  M     qA
:9:
:45
:GH
:HI
:KL
:KL
:LM
:LM
:RS
:RS
:ST
:ST
:RS w1-Iy8Y1ZZZZw1-Iy8Y1ZZZZwA3vy)U^>_7````wQ!4	9V_?`8aaaaw1-Iy8Y1ZZZZwQQ769iYbBc;ddddwAA6&)YXaAb:ccccy!Q2fiT]=^6____ya11=T]_hHiAjjjjyqAA>&)U^`iIjBkkkk 	QA3q5)Aw1-1qc!e1DDDDw1-1qc!e1DDDDy!Q2fQA3q56IIIIrt   c                    a \        ^4      o\        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R	 l4       \        \        V3R
 l4       \        \        V3R l4       \        \        V3R l4       \        \        V3R l4       SP                  R^ ^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^4      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^R7      \	        . RO. RO. RO.4      8X  g   Q hSP                  R^ ^^R7      \	        . RO. RO. RO.4      8X  g   Q h\        ^^^ .^,          4      p V P                  ! R^^4      \	        ^^^ .^,          4      8X  g   Q hV P                  ! R^ ^4      \	        ^^^ .^,          4      8X  g   Q hV P                  ! R^ ^^4      \	        ^^^ .^,          4      8X  g   Q hR# )r   c                  &   < S P                  R 4      # r  elementary_col_opr  s   ru   r   test_col_op.<locals>.<lambda>t  r  rt   c                  $   < S P                  4       # r{   r  r  s   ru   r   r  u  r  rt   c                  ,   < S P                  R ^^R7      # )r  ri  r   r  r  s   ru   r   r  v  r  rt   c                  ,   < S P                  R R^R7      # )r  r  rQ  r  r  s   ru   r   r  w  r  rt   c                  ,   < S P                  R ^^R7      # r  col1col2r  r  s   ru   r   r  x  r  rt   c                  ,   < S P                  R ^^R7      # r  r  r  s   ru   r   r  y  r  rt   c                  ,   < S P                  R R^R7      # r  r  rQ  r  r  s   ru   r   r  z  r  rt   c                  ,   < S P                  R ^RR7      # r  r  r  s   ru   r   r  {  r  rt   c                  .   < S P                  R ^^^R7      # r  r  r  r   r  r  s   ru   r   r  |  r  rt   c                  .   < S P                  R ^^^R7      # r  r  r  s   ru   r   r  }  r  rt   c                  .   < S P                  R R^^R7      # r  r  rQ  r  r  s   ru   r   r  ~  r  rt   c                  .   < S P                  R ^R^R7      # r  r  r  s   ru   r   r    r  rt   c                  .   < S P                  R ^^^R7      # r  r  r  s   ru   r   r    r  rt   r  r  )r  r   r  r  )ri  r  r  )ri  r   r  )r  r   r  Nr  r  r  rY  rZ  )r   r   r   )r   rR   r#  r  r5   rx   r  s    @ru   test_col_opr  q  r  rt   c                  :   \        ^4      p \        ^4      pV P                  '       g   Q hVP                  '       g   Q h\        ^ ^ . 4      pVP                  '       g   Q h\        ^^. RO4      pVP                  '       g   Q h\        ^^. RO4      pVP                  '       d   Q h\	        R4      p\        ^^V^ ^ .4      pVP                  '       g   Q h\        ^^W3^ .4      pVP                  '       d   Q h\        ^^. RO4      pVP                  '       d   Q hR# )r   rc   N)r   r   r   r   r   r   )r   r   r   r   r   r   	r   r   r   r   r   r   r   r   r   )r   r   
is_echelonrx   r   )zroidentr_   rc   s       ru   test_is_echelonr    s    
1
C1E>>>>Q2&A<<<<Q#56A<<<<Q#56A|||sAQAq!9-A<<<<QA!9-A|||Q#>?A||||rt   c                  |  a \        ^4      p \        ^4      pV P                  ! 4       V 8X  g   Q hVP                  4       V8X  g   Q h\        ^ ^ . 4      p V P                  ! 4       V 8X  g   Q h\        ^^^.4      p V P                  ! 4       V 8X  g   Q hV3R lp\        ^^. RO4      p \	        ^.R.^..4      .p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       \        ^^. RO4      p . p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       \        ^^. RO4      p \	        \        R^4      .^.^ ..4      \	        \        R	^4      .^ .^..4      .p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       \        ^^. R
O4      p \	        ^ .R	.^..4      .p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       \        ^^. RO4      p \	        ^.^ .^ ..4      \	        ^ .R.^..4      .p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       \        ^^. RO4      p \	        R.^.^ ..4      .p\        V P                  4       Uu. uF  q@VR3,          NK  	  ppV P                  ! 4       oSP                  '       g   Q hV! WV4       R# u upi u upi u upi u upi u upi u upi )r   c                   < V FM  p\         ;QJ d#    R  SV,           4       F  '       d   K   RM	  RM! R  SV,           4       4      '       d   KM  Q h	  V F  p\         ;QJ d    R V 4       F  '       d   K   RM	  RM! R V 4       4      '       d   K?  \         ;QJ d1    R SVP                  4       ,           4       F  '       d   K   RM&	  RM"! R SVP                  4       ,           4       4      '       g   K  Q h	  R# )c              3   8   "   T F  qP                   x  K  	  R # 5ir{   r  r   rf   s   & ru   r   Ctest_echelon_form.<locals>.verify_row_null_space.<locals>.<genexpr>  s     6+Qyy+r  FTc              3   8   "   T F  qP                   x  K  	  R # 5ir{   r  r  s   & ru   r   r    s     ,!Qyy!r  c              3   8   "   T F  qP                   x  K  	  R # 5ir{   r  r  s   & ru   r   r    s     J2IQyy2Ir  N)r   r   )r   r   nullsr   	a_echelons   &&& ru   verify_row_null_space0test_echelon_form.<locals>.verify_row_null_space  s    A36)A+63336)A+66666 A3,!,333,!,,,3J)AKKM2IJ333J)AKKM2IJJJJJ rt   r   Nr  r  )	r   r   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   r   r   r  r   )	r   r   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   )	r   r   rl  rx   r5   r   r   r  r   )r_   r  r  r  r~   r   r  s         @ru   test_echelon_formr    s   
 	AqA >>q   >>q   Q2&A>>q   QA3'A>>q   K 	Q#>?ATTT  E #166]+]adGG]D+ I!5) 	Q#>?AE"166]+]adGG]D+ I!5)Q#>?Ar1oVV  r1oVV 	E #166]+]adGG]D+ I!5) 	Q#>?AVVV  E #166]+]adGG]D+ I!5)Q#>?ASSS  TTT 	E #166]+]adGG]D+ I!5)Q#56ATSS  E #166]+]adGG]D+ I!5)y , , , , , ,s$   N ;N%N*N/ N4N9c                  X   \        ^ ^ . 4      p V P                  RR7      V 8X  g   Q h\        ^^^.4      p \        ^^^.4      pV P                  RR7      VP                  ! RR7      u;8X  d	   V 8X  g   Q h Q h\        ^^. RO4      pVP                  ! RR7      \        ^.^ .^ ..4      8X  g   Q h\        ^^. RO4      pVP                  ! RR7      \        . RO.4      8X  g   Q h\        ^^. RO4      pVP                  ! RR7      \        . RO. RO. R	O.4      8X  g   Q h\        ^^. R
O4      p\        ^^. RO4      p\        ^^. RO4      p\        ^^. RO4      pVP                  ! RR7      VP                  ! RR7      u;8X  d;   VP                  ! RR7      u;8X  d"   VP                  ! RR7      u;8X  d	   V8X  g   Q h Q h\        ^4      p \	        ^4      pV P                  RR7      V 8X  g   Q hVP                  ! RR7      V8X  g   Q h\        . RO. RO. RO. RO.4      pVP                  ! 4       w  rgV\        . RO. RO. RO. RO.4      8X  g   Q hVR8X  g   Q h\        \        ^^4      \        ^^4      ^^.. RO. RO. RO.4      pVP                  ! RR7      \        ^^ ^ \        R^4      .^ ^^ \        R^4      .^ ^ ^\        ^^4      .. RO.4      8X  g   Q h\        R4      p\        ^^V^^\        V4      V^.4      p\        VP                  ! RR7      ^^ \        V4      V) ^,           ,          V\        ^^4      ,          ) V,           ,          ^ ^^\        V4      V,           ^,           ,          .4       F(  w  r\        W,
          4      P                  '       d   K(  Q h	  R# )r   F)pivotsrc   Nr   r  )r   r   r  r   r   )	r   r   r   r   r   r   r   r   r   )	r   r   r   r   r   r   r   r   r   r  )	r   r   r   r   r   r   r   r   r   rQ  )r   r   r   r   r   rQ  r   )r  r   r   r   r   rv  r   )r   r   r  r   r   r   rQ  )r  r   r   r  r  r   r   )r   rQ  r   r   r   r   r  )r   r   r   r   r   r  r   )r   r   r   r   r   r  r   )r   r   r   r   r   r   r   )r   r   r   r  r  r  ir  )rx   rT  r5   r   r   r   r   r   r  r)   r  )r  r_   r`   ra   rb   re   r   pivot_offsetsrc   r~   r   s              ru   	test_rrefr    s   Q2&A6661$$$QA3'AQA3'A666166#7<1<<<<<Q9-A666A3aS/#::::Q9-A6669+#6666Q#>?A666%/%/%/+1 $2 2 2 2
 	Q#>?AQ#>?AQ#>?AQ#>?A66FF% &FF% & FF% & %&& & & & &
 	qAA6661$$$661$$$''''	) 	*A
 C&----	/ 0 0 0 0
 I%%%x2!QA!L0002 	3A 666*+Q8C3E)F*+QHR4E)F*+Q8C3E)F);	+= $> > > > 	sAQAq!T!Wa#;<AAFF%(47QBF#a!Q&7%7!%;<1aa1q)+, &&&&&,rt   c                  >   \        R 4      w  rr#\        ^ ^ .^ ^ .^^.^^..4      p\        WW#.4      pVP                  V4      \        ^^ .^ ^.^ ^ .^ ^ ..4      \        RV,          V,           .^V,          ^,          V^,          ,
          .V .V..4      38X  g   Q hR# )r  Nr  )r   r5   rref_rhs)r_   r`   ra   rb   r   r   s         ru   test_rref_rhsr  P  s    #JA!AAAA/0Aa|A::a=VFFFF	%  
1qMqSUQqS[MMM	 	   rt   c            	        a \        . RO. RO. RO. RO. RO. RO.4      oSP                  R^^R7      p SP                  R^^RR7      pSP                  R^^R7      pV R,          \        . RO.4      8X  g   Q hVR,          \        . RO.4      8X  g   Q hVR,          \        . RO.4      8X  g   Q h\        \        V3R l4       \        \        V3R l4       \        \        V3R	 l4       R
# )r   r  r  r  r  r  r  c                  ,   < S P                  R ^^R7      # r  r  r   s   ru   r   "test_issue_17827.<locals>.<lambda>p  r  rt   c                  ,   < S P                  R ^^R7      # r  r  r   s   ru   r   r  q  r  rt   c                  .   < S P                  R R^^R7      # )r  r  r  r  r   s   ru   r   r  r  r  rt   N)r   r   r  r   )r  r  r   r   )r   r   r   r   )r   r   r   r  )r   r   r   rQ  )r   r  r   r   rP  )r   r   )r   r   r   r  )r  r  r   r  )r5   r  rR   r#  )rT  r   rt  r   s      @ru   test_issue_17827r  _  s     	A 	
G!!4A	IAA<A	Ga0AT7fl^,,-,T7fm_--.-T7fo.//0/
:KL
:GH
:STrt   c                  "   \        ^^.\        ^^\        ,          ,
          ..4      p V P                  4       ^8X  g   Q h\        ^^\        ^^
4      4      pVP                  4       ^8X  g   Q h\	        ^4      pVP                  4       ^ 8X  g   Q hR# r   )r5   rc   r  r   rB   )r   r   r,  s      ru   	test_rankr  t  st    AA!G%&A668q==q!U1b\"A668q==aA668q==rt   c                  l   \        R 4      w  rr#rErgrr\        WW
,          W,          ^ .W#W*,          W:,          ^ .WEWJ,          WZ,          ^.WgWj,          Wz,          ^.W^V,          V,          W,          ,
          ^V	,          V,          W,          ,
          ^ ..4      pVP                  4       ^8X  g   Q hR# )r#  N)r   r5   r  )r$  r%  r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r   s                ru   test_issue_11434r  |  s    AB 3BBBBBBq)q)q)q)2b25!B$r'BE/1=	? 	@A
 668q==rt   c            
          \        R 4      w  r\        RV ,          ^^ ^ .^V ,          RV ,          ^,
          ^^ .^ ^V ,          RV ,          V,
          ^,
          ^.^ ^ W,           R..4      p\        ^^ ^ ^V ^,          V) V ,
          ,          ,          .^ ^^ ^W) V ,
          ,          ,          .^ ^ ^^V) V ,
          ,          .. RO.4      pRpVP                  4       w  rV\        W5,
          4      \	        VP
                  !  8X  g   Q hWd8X  g   Q hR# )z
nu, lambdaNr   r  r  r  r   )r   r5   rT  r)   rB   r   )nulambr   expected_reducedexpected_pivotsreducedr  s          ru   test_rank_regression_from_sor    s    |$HBBA:B$2	A:QrTBrET>A#5:BIr:< 	=A 1aBED52:,>)?@ !1aAr52:,?@ !1aTEBJ@@B C  OffhOG$./5!''?BBB$$$rt   c                  *   \        . RO. RO. RO. RO.4      p V \         P                  ! ^4      \        ,          ,
          pVP                  4       ^8X  g   Q hV^,          P                  4       ^8X  g   Q hV^,          P                  4       ^8X  g   Q hR# )r   Nr  r5  r6  r7  )r5   r9   r
   r  )r   r   s     ru   test_issue_15872r    sn    o~|LMA	FJJqMAA668q==qD;;=AqD;;=Art   ([  r  concurrent.futuresr  collections.abcr   sympy.core.addr   sympy.core.functionr   r   r   r  r   r	   r
   r   r   r   r   r   sympy.core.powerr   sympy.core.singletonr   sympy.core.symbolr   r   sympy.core.sympifyr   $sympy.functions.elementary.complexesr   &sympy.functions.elementary.exponentialr   r   (sympy.functions.elementary.miscellaneousr   r   r   (sympy.functions.elementary.trigonometricr   r   r   sympy.integrals.integralsr   $sympy.matrices.expressions.transposer   sympy.physics.quantum.operatorr    r!   r"   sympy.polys.polytoolsr#   r$   sympy.polys.rootoftoolsr%   sympy.printing.strr&   sympy.sets.setsr'   rz  r(   r)   sympy.simplify.trigsimpr*   sympy.matrices.exceptionsr+   r,   r-   sympy.matrices.matrixbaser/   sympy.matrices.determinantr0   sympy.matrices.utilitiesr1   sympy.matricesr2   r3   r4   r5   r6   r7   r8   r9   r:   r;   r<   r=   r>   r?   r@   rA   rB   rC   rD   rE   rF   rG   rH   rI   rJ   
sympy.corerK   rL   (sympy.functions.special.tensor_functionsrM   sympy.utilities.iterablesrN   rO   rP   sympy.utilities.exceptionsrQ   sympy.testing.pytestrR   rS   rT   rU   rV   rW   sympy.assumptionsrY   sympy.tensor.arrayrZ   $sympy.tensor.array.array_derivativesr[   sympy.matrices.expressionsr\   sympy.algebrasr]   r  r^   r!  r_   r`   ra   rb   rc   rd   re   rf   r   r  rg   rh   sympy.matrices.matricesri   rj   rl   rx   r   r   r   r   r   r   r   r   r   r   r   r%  r)  r0  r2  r4  r?  rE  rK  rN  rU  rZ  r_  rd  rj  rm  rq  rt  rx  r{  r~  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r
  r  r  r4  r<  r@  rD  rG  rK  rM  rQ  rW  r]  r_  rl  rp  ru  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r
  r  r  r*  r0  r3  r8  r?  rF  rI  rN  rQ  rY  r]  r`  r  r  r  r  r  r  r  r  r  r  r  r   r  r  r
  r  r  r  r  r,  r4  r6  r:  r<  r?  rA  rG  rI  rN  rW  rY  r[  r^  ra  rj  rv  r{  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r!  r0  r<  rF  rS  rX  r\  r`  rf  rh  rk  rp  rr  r{  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  rm   rt   ru   <module>r6     s     $  6 6 L L L   " / & 4 = E E D D / : N N 2 * # % 8 ,  4 C .L L L L L L L 8 " C @ @ 6: :  $ @ - %  , , , <2G
H B E ^_o  
 ~@P  
@6
+
!+
&
9)XDN>^JB?_*JX_
%
'
 
6%
;
5
8GP
y
	8$
)&$<<<<<K&d&K0O2dFD+84P1$
:<7Mt 8!"AG2" "$$-:&>9$"*)  
.	*F!@.6&8G$1&3@E!P
1H%V&&'"$4.L)2
!$;6z"IG
O*4(nP!%IPZz.&"FH $<">JD-;-"8
/(1
(.3(?G4\p&(*&F"7:& 4 5 5))+"*
Z "$ "$J71R$# 	+	#E '
/8	1	1 205&#3*9# 
$I< :;# <#;.!JH!JH6\*~>'BU*%(ij 
 s   R3R$R!	$R5	