o
    à8Va‹
  ã                   @   st   d dl mZ d dlmZmZ d dlmZ d dlmZ G dd„ deƒZ	d dl
mZmZ d dlmZ d	d
„ Zeed< dS )é    )Ú_sympify)ÚSÚBasic)ÚNonSquareMatrixError)ÚMatPowc                   @   s`   e Zd ZdZdZejZejfdd„Ze	dd„ ƒZ
e	dd„ ƒZd	d
„ Zdd„ Zdd„ Zdd„ ZdS )ÚInversea  
    The multiplicative inverse of a matrix expression

    This is a symbolic object that simply stores its argument without
    evaluating it. To actually compute the inverse, use the ``.inverse()``
    method of matrices.

    Examples
    ========

    >>> from sympy import MatrixSymbol, Inverse
    >>> A = MatrixSymbol('A', 3, 3)
    >>> B = MatrixSymbol('B', 3, 3)
    >>> Inverse(A)
    A**(-1)
    >>> A.inverse() == Inverse(A)
    True
    >>> (A*B).inverse()
    B**(-1)*A**(-1)
    >>> Inverse(A*B)
    (A*B)**(-1)

    Tc                 C   s6   t |ƒ}|jstdƒ‚|jstd| ƒ‚t | ||¡S )Nzmat should be a matrixzInverse of non-square matrix %s)r   Z	is_MatrixÚ	TypeErrorZ	is_squarer   r   Ú__new__)ÚclsZmatÚexp© r   úD/usr/lib/python3/dist-packages/sympy/matrices/expressions/inverse.pyr	   #   s   zInverse.__new__c                 C   s
   | j d S ©Nr   )Úargs©Úselfr   r   r   Úarg-   s   
zInverse.argc                 C   s   | j jS ©N)r   Úshaper   r   r   r   r   1   s   zInverse.shapec                 C   s   | j S r   )r   r   r   r   r   Ú_eval_inverse5   s   zInverse._eval_inversec                 C   s   ddl m} d|| jƒ S )Nr   )Údeté   )Z&sympy.matrices.expressions.determinantr   r   )r   r   r   r   r   Ú_eval_determinant8   s   zInverse._eval_determinantc                 K   sB   d|v r|d dkr| S | j }| dd¡r|jdi |¤Ž}| ¡ S )NZ
inv_expandFZdeepTr   )r   ÚgetÚdoitZinverse)r   Zhintsr   r   r   r   r   <   s   zInverse.doitc                 C   sB   | j d }| |¡}|D ]}| j| j 9  _| j| 9  _q|S r   )r   Ú_eval_derivative_matrix_linesZfirst_pointerÚTZsecond_pointer)r   Úxr   ÚlinesÚliner   r   r   r   F   s   

z%Inverse._eval_derivative_matrix_linesN)Ú__name__Ú
__module__Ú__qualname__Ú__doc__Z
is_Inverser   ZNegativeOner   r	   Úpropertyr   r   r   r   r   r   r   r   r   r   r      s    



r   )ÚaskÚQ)Úhandlers_dictc                 C   sT   t t | ¡|ƒr| jjS t t | ¡|ƒr| j ¡ S t t | ¡|ƒr(td| j ƒ‚| S )zÌ
    >>> from sympy import MatrixSymbol, Q, assuming, refine
    >>> X = MatrixSymbol('X', 2, 2)
    >>> X.I
    X**(-1)
    >>> with assuming(Q.orthogonal(X)):
    ...     print(refine(X.I))
    X.T
    zInverse of singular matrix %s)	r%   r&   Z
orthogonalr   r   ZunitaryÚ	conjugateZsingularÚ
ValueError)ÚexprZassumptionsr   r   r   Úrefine_InverseS   s   

r+   N)Zsympy.core.sympifyr   Z
sympy.corer   r   Zsympy.matrices.commonr   Z!sympy.matrices.expressions.matpowr   r   Zsympy.assumptions.askr%   r&   Zsympy.assumptions.refiner'   r+   r   r   r   r   Ú<module>   s    G