o
    *ήc/                     @   s   d dl mZ ddlmZ ddlmZ d dlmZ d dlm	Z	 d dl
mZ d dlmZ d d	lmZ d d
lmZ G dd deZdS )    )NonSquareMatrixError   )
MatrixExpr)Identity)S)ExprBuilder)cacheit)Pow)_sympify)
MatrixBasec                       s~   e Zd Zd fdd	Zedd Zedd Zedd	 Zed
d Z	dd Z
dd Zdd Zdd Zdd Zdd Z  ZS )MatPowFc                    sT   t |}|jstd|jstd| t |}t | ||}|r(|jdd}|S )NzMatPow base should be a matrixzPower of non-square matrix %sF)deep)r
   	is_Matrix	TypeError	is_squarer   super__new__doit)clsbaseexpevaluateoptionsobj	__class__ H/tmp/pip-target-vg8gfxp4/lib/python/sympy/matrices/expressions/matpow.pyr      s   zMatPow.__new__c                 C   
   | j d S )Nr   argsselfr   r   r   r         
zMatPow.basec                 C   r   )Nr   r   r!   r   r   r   r   !   r#   z
MatPow.expc                 C   s   | j jS N)r   shaper!   r   r   r   r%   %   s   zMatPow.shapec                 C   s   | j  | j S r$   )r   as_explicitr   r!   r   r   r   _get_explicit_matrix)   s   zMatPow._get_explicit_matrixc                    s   ddl m} |   t tr> jjr& jjr&| fddt jD   n| 	 s2 
 ||f S ddlm} || ||S  ||f S )Nr   MatMulc                       g | ]} j qS r   r   ).0kAr   r   
<listcomp>3       z!MatPow._entry.<locals>.<listcomp>)MatrixElement)sympy.matrices.expressionsr)   r   
isinstancer   r   
is_Integeris_positiverange_is_shape_symbolicr'   "sympy.matrices.expressions.matexprr2   )r"   ijkwargsr)   r2   r   r.   r   _entry-   s   
zMatPow._entryc                    s     ddr fdd| jD \}}n| j\}}t|tr.||jd 9 }|jd }t|tst|tr7|| S |tjkr>|S |tjkrHt|j	S |tj
kr]ddlm} ||jd	i  S t|dd }|d urk||S t||S )
Nr   Tc                    s   g | ]
}|j d i  qS )r   )r   )r,   arghintsr   r   r0   >   s    zMatPow.doit.<locals>.<listcomp>r   r   Inverse_eval_powerr   )getr    r4   r   r   r   OneZeror   rowsNegativeOner3   rB   r   getattr)r"   r@   r   r   rB   
eval_powerr   r?   r   r   <   s(   









zMatPow.doitc                 C   s   | j \}}t|j|S r$   )r    r   T)r"   r   r   r   r   r   _eval_transposeZ   s   
zMatPow._eval_transposec                 C   s   t | |S r$   )r	   _eval_derivative)r"   xr   r   r   rM   ^   s   zMatPow._eval_derivativec           
         s  ddl m} ddlm} ddlm} ddlm  j}j	j
dkro||soj	|}|D ]>}t|t|td|jd |j	|d   |jd tdgd	d
g|jd}|jd j|_d|_|jd j|_d|_|g|_q.|S |dkdkr|fddt|D }	n:|dkdkr j	|S |dk dkr| fddt| D }	n|dkdkr |S td|f |	|S )Nr   )ArrayContraction   )ArrayTensorProductr   r(   rA   )r   r   )r   rP      )         )	validatorrR   Tc                    r*   r   r+   r,   r:   r!   r   r   r0      r1   z8MatPow._eval_derivative_matrix_lines.<locals>.<listcomp>c                    s   g | ]} j qS r   r+   rW   rB   r"   r   r   r0      s    z cannot evaluate %s derived by %s)0sympy.tensor.array.expressions.array_expressionsrO   *tensor.array.expressions.array_expressionsrQ   matmulr)   inverserB   r   r   r%   has_eval_derivative_matrix_linesr   r   _lines	_validater    _first_pointer_parent_first_pointer_index_second_pointer_parent_second_pointer_indexfromiterr7   r   NotImplementedError)
r"   rN   rO   rQ   r)   r   lrr:   subexprnewexprr   rY   r   r_   a   sN   

"
z$MatPow._eval_derivative_matrix_linesc                 C   s   t | j| j S r$   )r   r   r   r!   r   r   r   _eval_inverse   s   zMatPow._eval_inverse)F)__name__
__module____qualname__r   propertyr   r   r%   r   r'   r=   r   rL   rM   r_   rk   __classcell__r   r   r   r   r      s     



,r   N)sympy.matrices.commonr   matexprr   specialr   
sympy.corer   sympy.core.exprr   sympy.core.cacher   sympy.core.powerr	   sympy.core.sympifyr
   sympy.matricesr   r   r   r   r   r   <module>   s    