o
    GZh                     @   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 )   )
MatrixExpr)Identity    )S)ExprBuilder)cacheit)Pow)_sympify)
MatrixBase)NonSquareMatrixErrorc                       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dd Zdd Z  ZS )MatPowFc                    sX   t |}|jstd|jdu rtd| t |}t | ||}|r*|jdd}|S )NzMatPow base should be a matrixFzPower of non-square matrix %s)deep)r	   Z	is_Matrix	TypeErrorZ	is_squarer   super__new__doit)clsbaseexpevaluateoptionsobj	__class__ P/var/www/auris/lib/python3.10/site-packages/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   Z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   Z
is_IntegerZis_positiverangeZ_is_shape_symbolicr$   Z"sympy.matrices.expressions.matexprr/   )r    ijkwargsr&   r/   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                 3   s     | ]}|j di  V  qd S )Nr   )r   )r)   arghintsr   r   	<genexpr>>   s    zMatPow.doit.<locals>.<genexpr>r   r   InverseZ_eval_powerr   )getr   r1   r   r
   r   ZOneZZeror   rowsZNegativeOner0   r<   r   getattr)r    r9   r   r   r<   Z
eval_powerr   r8   r   r   <   s(   









zMatPow.doitc                 C      | j \}}t| |S r"   )r   r   Z	transposer    r   r   r   r   r   _eval_transposeZ      
zMatPow._eval_transposec                 C   r@   r"   )r   r   ZadjointrA   r   r   r   _eval_adjoint^   rC   zMatPow._eval_adjointc                 C   r@   r"   )r   r   	conjugaterA   r   r   r   _eval_conjugateb   rC   zMatPow._eval_conjugatec                 C   s   t | |S r"   )r   _eval_derivative)r    xr   r   r   rG   f   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%   r;   )r   r   )r   rJ      )         )	validatorrL   Tc                    r'   r   r(   r)   r3   r   r   r   r-      r.   z8MatPow._eval_derivative_matrix_lines.<locals>.<listcomp>c                    s   g | ]} j qS r   r(   rQ   r<   r    r   r   r-      s    z cannot evaluate %s derived by %s)Z0sympy.tensor.array.expressions.array_expressionsrI   Z*tensor.array.expressions.array_expressionsrK   matmulr&   Zinverser<   r   r   r#   has_eval_derivative_matrix_linesr   r   _lines	_validater   Z_first_pointer_parentZ_first_pointer_indexZ_second_pointer_parentZ_second_pointer_indexZfromiterr2   r   NotImplementedError)
r    rH   rI   rK   r&   r   lrr3   ZsubexprZnewexprr   rS   r   rV   i   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$   r6   r   rB   rD   rF   rG   rV   r[   __classcell__r   r   r   r   r      s$    



,r   N)Zmatexprr   Zspecialr   Z
sympy.corer   Zsympy.core.exprr   Zsympy.core.cacher   Zsympy.core.powerr   Zsympy.core.sympifyr	   Zsympy.matricesr
   Zsympy.matrices.exceptionsr   r   r   r   r   r   <module>   s    