| 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/linalg/tests/__pycache__/ |
Upload File : |
�
d�c�� � � � d Z dZddlZddlZddlZddlZddlmZm Z m
Z
mZmZm
Z
ddlZddlmZ ddlmZmZmZmZmZmZmZmZmZmZmZmZmZmZm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/ ddl0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z: ddl;m<Z< dd l=m>Z> dd
l?m@Z@ ddlmAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZHmIZImJZJmKZKmLZLmMZMmNZNmOZOmPZPmQZQmRZRmSZSmTZTmUZUmVZVmWZW ddlXmYZYmZZZ dd
l[m\Z\ ddl]m^Z^ ddl_m`Z` ddlambZb decfd�ZdejI eje gZfejJ ejg gZhefehz ZidVd�Zjd� Zkd� Zld� Zm G d� d� � Zn G d� d� � Zo G d� d� � Zp G d� d� � Zq G d � d!� � Zr G d"� d#� � Zs G d$� d%es� � Zt G d&� d'� � Zu G d(� d)� � Zv G d*� d+ev� � Zw G d,� d-� � Zx G d.� d/� � Zy G d0� d1� � Zz G d2� d3� � Z{ G d4� d5� � Z| G d6� d7� � Z} G d8� d9� � Z~d:� Z G d;� d<� � Z� G d=� d>� � Z� G d?� d@� � Z�dA� Z�ej� �� ej� � � dBk dC�D� � dE� � � Z�dF� Z�dG� Z�ej� �� ddH�I� � dJ� � � Z� G dK� dL� � Z�dWdM�Z�ej� j� ej� �� ejU ej� � � j� dNk dO�D� � dP� � � � � Z�dQ� Z�dR� Z�dS� Z� G dT� dU� � Z�dS )Xz* Test functions for linalg.decomp module
z~
Build linalg:
python setup_linalg.py build
Run tests if scipy is installed:
python -c 'import scipy;scipy.linalg.test()'
� N)�assert_equal�assert_almost_equal�assert_array_almost_equal�assert_array_equal�assert_�assert_allclose)�raises)�eig�eigvals�lu�svd�svdvals�cholesky�qr�schur�rsf2csf�lu_solve� lu_factor�solve�diagsvd�
hessenberg�rq�
eig_banded�eigvals_banded�eigh�eigvalsh�qr_multiply�qz�orth�ordqz�subspace_angles�hadamard�eigvalsh_tridiagonal�eigh_tridiagonal�
null_space�cdf2rdf�LinAlgError)
�dgbtrf�dgbtrs�zgbtrf�zgbtrs�dsbev�dsbevd�dsbevx�zhbevd�zhbevx�get_lapack_funcs)�norm)�_select_function)�ortho_group)�array�diag�ones�full�linalg�argsort�zeros�arange�float32� complex64�ravel�sqrt� iscomplex�shape�sort�sign�asarray�isfinite�ndarray�eye�dtype�triu�tril)�seed�random)�assert_no_overwrite)�matrix)�check_free_memory)� HAS_ILP64Fc � � |t v rdt j � | | � � t j � | | � � dz z }||� � � j z dz }n-t j � | | � � }||j z dz }|r*|t
d| z � � t j | � � z z
}|� |� � S )z7Generate random sym/hermitian array of the given size n� �?� ) �COMPLEX_DTYPES�nprM �rand�conj�Tr@ rH �astype)�n�posdefrI �As �@/usr/lib/python3/dist-packages/scipy/linalg/tests/test_decomp.py�_random_hermitian_matrixr_ 3 s� � ������I�N�N�1�a� � �2�9�>�>�!�Q�#7�#7��#<�<��
�������^�Q�����I�N�N�1�a� � ��
���W�a�K��
� !� �T�!�A�#�Y�Y�r�v�a�y�y�
� ���8�8�E�?�?�� � c � � t j | � � }t j |� � r,t |j � � dt |j � � z z S t j |j � � j }t j |� � \ }}d||z
z }||z }t j
||�� � ||z }t j ||� � S )zCClears trailing `fuss_binary_bits` of mantissa of a floating numberrS � @)�out)rV �
asanyarray�iscomplexobj�
clear_fuss�real�imag�finforI �nmant�frexp�rint�ldexp)�ar�fuss_binary_bits�x�significant_binary_bits�x_mant�x_exp�fs r^ rg rg G s� � �
�
�b���A� ��q��� <��!�&�!�!�B��A�F�);�);�$;�;�;� �h�q�w�/�/�5���H�Q�K�K�M�F�E��%�(8�8�9�A�
�a�K�F��G�F������
�a�K�F�
�8�F�E�"�"�"r` c � � t | t � � r| j } nt | � � } t |t � � r|j }nt |� � }t | |k d� | |� � � � d S )Nz%dtype mismatch: "{}" (should be "{}"))�
isinstancerG rI r �format)�act�dess r^ �assert_dtype_equalr{ X s� � ��#�w��� ��i����C�j�j���#�w��� ��i����C�j�j���C�3�J�3�:�:�3��D�D�F� F� F� F� Fr` c � � t | t � � r| }t |� � dz dz
}nLt | t � � r(t | j � � dk r| j d }| }nt
d� � �t j |� � }|j �
� � t |� � z |z }d|j |z z }|S )a Return a random symmetric (Hermitian) matrix.
If 'dim_or_eigv' is an integer N, return a NxN matrix, with eigenvalues
uniformly distributed on (-1,1).
If 'dim_or_eigv' is 1-D real array 'a', return a matrix whose
eigenvalues are 'a'.
rT � r zinput type not supported.� �?)rw �intrM rG �lenrB � TypeErrorr4 �rvsrY rX r6 )�dim_or_eigv�dim�d�v�hs r^ �symrandr� i s� � � �+�s�#�#� 5����3�K�K��M�A����
�[�'�
*�
*� 5�
�k��
�
�A�
%�
%����"������3�4�4�4������A� ����
�
�T�!�W�W��q� �A��Q�S��U��A��Hr` c � � t | � � t | � � }}|dt |� � t |� � z
z z }|� |� � S )NrS )r� rJ rK rZ )r� rI �a1�a2�as r^ �_complex_symrandr� � sG � �
�S�\�\�7�3�<�<��B�
�R��b���$�r�(�(�"�
#�#�A��8�8�E�?�?�r` c �&