| Server IP : 93.86.61.54 / Your IP : 216.73.216.60 Web Server : Apache/2.4.62 (Ubuntu) System : Linux rasin.ddns.net 6.8.0-124-generic #124~22.04.1-Ubuntu SMP PREEMPT_DYNAMIC Tue May 26 21:05:19 UTC x86_64 User : www-data ( 33) PHP Version : 8.4.22 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : ON | Perl : ON | Python : OFF | Sudo : ON | Pkexec : ON Directory : /usr/lib/python3/dist-packages/sympy/polys/matrices/__pycache__/ |
Upload File : |
o
�8Va] � @ s� d dl mZ d dlmZ d dlmZ d dlmZ d dlm Z d dl
m Z ddl m
Z
mZmZmZ d d
� Zd d � Zd
d� Zdd� Zdd� Zdd� ZdS )� )� defaultdict��Add)�Mul)�S)�construct_domain)�PolyNonlinearError� )�SDM� sdm_irref�sdm_particular_from_rref�sdm_nullspace_from_rrefc C sF t |�}t| |�\}}t|||�}|j}|js|jr$|�� �� d �� }t |�\}}} |r5|d |kr5dS t
||d |�}
t ||j ||| �\} } t
t�}
|
�� D ]\}}|
|| �|�|�� qPt| | �D ] \}}|| }|�� D ]\}}|
|| �||�|� � qrqfdd� |
�� D �}
tj}t|�t|
� D ]}||
|< q�|
S )a Solve a linear system of equations.
Examples
========
Solve a linear system with a unique solution:
>>> from sympy import symbols, Eq
>>> from sympy.polys.matrices.linsolve import _linsolve
>>> x, y = symbols('x, y')
>>> eqs = [Eq(x + y, 1), Eq(x - y, 2)]
>>> _linsolve(eqs, [x, y])
{x: 3/2, y: -1/2}
In the case of underdetermined systems the solution will be expressed in
terms of the unknown symbols that are unconstrained:
>>> _linsolve([Eq(x + y, 0)], [x, y])
{x: -y, y: y}
r ���Nr c S � i | ] \}}|t |� �qS � r )�.0�s�termsr r �?/usr/lib/python3/dist-packages/sympy/polys/matrices/linsolve.py�
<dictcomp>j � z _linsolve.<locals>.<dictcomp>)�len�_linear_eq_to_dict�sympy_dict_to_dm�domain� is_RealField�is_ComplexField�to_ddm�rref�to_sdmr r r
�oner �list�items�append�to_sympy�zipr �Zero�set)�eqs�syms�nsyms�eqsdict�rhs�Aaug�K�Arref�pivots�nzcols�P�V� nonpivots�sol�i�v�npi�Vi�sym�zeror r r r � _linsolve. s0 �
r< c
s� t |�jdd� | D �� }t|ddd�\}}tt||��� t| �}t|�}tt|t|����g }t| |�D ] \} }
� �fdd�| �� D �} |
rM� |
| |< | rT|�| � q4t t
|�||d f|�} | S )z?Convert a system of dict equations to a sparse augmented matrixc s s � | ]}|� � V qd S �N)�values)r �er r r � <genexpr>w s � z#sympy_dict_to_dm.<locals>.<genexpr>T)�field� extensionc s i | ] \}}�| � | �qS r r )r r �c��elem_map� sym2indexr r r s z$sympy_dict_to_dm.<locals>.<dictcomp>r ) r' �unionr �dictr% r �ranger" r# r
� enumerate)
�
eqs_coeffs�eqs_rhsr) �elemsr. �elems_K�neqsr* r+ �eqr, �eqdict�sdm_augr rD r r u s
�r c s dd� � � fdd�| D �S )z�Use expand to cancel nonlinear terms.
This approach matches previous behaviour of linsolve but should be
deprecated.
c S s | j r | j| j } | �� S r= )� is_Equality�lhsr, �expand)rP r r r � expand_eq� s z)_expand_eqs_deprecated.<locals>.expand_eqc s g | ]}� |��qS r r )r rP �rV r r �
<listcomp>� s z*_expand_eqs_deprecated.<locals>.<listcomp>r )r( r rW r �_expand_eqs_deprecated� s rY c C s2 zt | |�W S ty t| �} t | |� Y S w �z1Convert a system Expr/Eq equations into dict form)�_linear_eq_to_dict_innerr rY )r( r) r r r r � s �r c C sF t |�}g g }}| D ]}t||�\}}|�|� |�|� q ||fS rZ )r' � _lin_eq2dictr# )r( r) r+ rL rP r, rQ r r r r[ � s
r[ c sH | |v r t j| t jifS | jrHtt�}g }| jD ] }t||�\}}|�|� |� � D ] \}}|| �|� q*qt
|� � dd� |� � D �} � | fS | j r�d } }
g }| jD ]}t||�\}}|se|�|� qT| du rn|} |}
qTt �t
|� � | du r|� i fS � fdd�| � � D �} � |
| fS | jr�t| j| j |�S | j|@ s�| i fS t �)z?Efficiently convert a linear equation to a dict of coefficientsc S r r r )r r: �coeffsr r r r � r z _lin_eq2dict.<locals>.<dictcomp>Nc s i | ] \}}|� | �qS r r )r r: rC ��coeffr r r � r )r r& �One�is_Addr r! �argsr\ r# r" r �is_Mulr r rS rT r, � free_symbols) �a�symset�
terms_list�
coeff_list�ai�ci�ti�mij�cijr � terms_coeffr r^ r r\ � sF
�
r\ N)� collectionsr �sympy.core.addr �sympy.core.mulr �sympy.core.singletonr �sympy.polys.constructorr �sympy.polys.solversr �sdmr
r r r
r< r rY r r[ r\ r r r r �<module> s G