| Server IP : 217.160.0.135 / Your IP : 216.73.217.85 Web Server : Apache System : Linux www 6.18.52-i1-ampere #1203 SMP Mon Sep 14 18:29:59 CEST 2026 aarch64 User : sws1074145052 ( 1074145052) PHP Version : 8.3.32 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : ON | Perl : ON | Python : OFF | Sudo : OFF | Pkexec : OFF Directory : /lib/python3/dist-packages/scipy/sparse/linalg/_isolve/__pycache__/ |
Upload File : |
�
d�c � � � d Z g ZddlmZmZmZmZ ddlmZm Z m
Z
i dd�dd�dd �d
d �d d�d
d�dd �dd �dd �dd �dd �dd �dd �dd �dd �dd �Z d� Z d� Z
d� ZdS ) zrestructuredtext en� )�
asanyarray�asarray�array�zeros)�aslinearoperator�LinearOperator�IdentityOperator)�fr
r
)r
�dr )r
�Fr )r
�Dr
)r r
)r r )r r )r r
)r r
)r r )r r )r r
)r
r
)r
r )r
r )r
r
c �8 � | dvrd} |dvrd}t | |f S )N�fdFDr )�
_coerce_rules)�x�ys �C/usr/lib/python3/dist-packages/scipy/sparse/linalg/_isolve/utils.py�coercer s- � ����� ������ �� ��1�� � � c � � | S �N� �r s r �idr s � � �Hr c �h � | }t | � � } | j d | j d k rt d| j f� �� � �| j d }t |� � }|j |dfk s,|j |fk s t d| j � d|j � d�� � �|j j dvr|� d� � }d � }t | d
� � r
| j j }n | � |� � j j }t ||j j � � }t ||� � � }|� � � }|��t |d
� � r|j }nt }t |d� � r|j } nt } |t u r%| t u r t | j | j �� � }nLt! | j || | j �� � }n.t |� � }| j |j k rt d� � �|�t# ||� � � }
n�t% |t& � � r0|dk r)|� � � } |� | � � }
n^t+ ||� � � }
|
j |dfk s,|
j |fk s t d| j � d|
j � d�� � �|
� � � }
| ||
||fS )a� Make a linear system Ax=b
Parameters
----------
A : LinearOperator
sparse or dense matrix (or any valid input to aslinearoperator)
M : {LinearOperator, Nones}
preconditioner
sparse or dense matrix (or any valid input to aslinearoperator)
x0 : {array_like, str, None}
initial guess to iterative method.
``x0 = 'Mb'`` means using the nonzero initial guess ``M @ b``.
Default is `None`, which means using the zero initial guess.
b : array_like
right hand side
Returns
-------
(A, M, x, b, postprocess)
A : LinearOperator
matrix of the linear system
M : LinearOperator
preconditioner
x : rank 1 ndarray
initial guess
b : rank 1 ndarray
right hand side
postprocess : function
converts the solution vector to the appropriate
type and dimensions (e.g. (N,1) matrix)
r � z&expected square matrix, but got shape=z shapes of A z and b z are incompatibler r c � � | S r r r s r � postprocessz make_system.<locals>.postprocessQ s � ��r �dtype)r N�psolve�rpsolve)�shaper )�matvec�rmatvecr z/matrix and preconditioner have different shapes�Mbz and x0 )r r" �
ValueErrorr r �char�astype�hasattrr# r r �ravelr r r! r r r �
isinstance�str�copyr ) �A�M�x0�b�A_�Nr �xtyper r! r �bCopys r � make_systemr6 s� � �B
�B�����A��w�q�z�Q�W�Q�Z� � ��N�1�7�*�N�N�O�O�O� ���
�A��1�
�
�A�
�G��!�u� � ���A�4���� (��� (� (��� (� (� (� )� )� )� �w�|�6�!�!�
�H�H�S�M�M��� � � �q���� '��� ������� � �!�&�� �5�!�'�,� '� '�E���� � � �A� ��� � �A� �y� �2�h� � � ��Y�F�F��F� �2�i� � � ��j�G�G��G� �R�<�<�G�r�M�M� �q�w�a�g�>�>�>�A�A� �q�w�v�w�%&�W�.� .� .�A�A�
�Q� � �� �7�a�g� � � � N�O�O� O�
�z� �!�5� !� !� !��� �B�� � � �
��:�:��F�F�H�H�E�������A�� �"�E� "� "� "����A�q�6�!�!�Q�W���_�_� � >�A�G� >� >�#$�7� >� >� >� ?� ?�
?�
�G�G�I�I�� �a��A�{� "�"r N)�
__docformat__�__all__�numpyr r r r � scipy.sparse.linalg._interfacer r r r r r r6 r r r �<module>r; sU ��%�
�
�� 4� 3� 3� 3� 3� 3� 3� 3� 3� 3� 3� 3�� � � � � � � � � � ��3� � �#� �y�� ��3� � )�#� �/8�� ��3� � )�#� �/8�� � �3� � !*�#� � 09�� � �3� � !*�#� � 09�� �
�3� �
� � � �
�
�
�`#� `#� `#� `#� `#r