o
    Yi	                     @   s   d dl Zd dlmZ d dlmZmZmZ d dlm	Z	 d dl
m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dS )    N)norm)assert_equalassert_allcloseassert_)minres)raisesc                  C   s4   t jd} | dd}||j }| d}||fS )N  
   )nprandomRandomStaterandT)rngmatrixvector r   T/tmp/pip-target-1s0edx8b/lib/python/scipy/sparse/linalg/_isolve/tests/test_minres.pyget_sample_problem	   s
   

r   c                  C   sV   t  \} }d| d< d|d< t| |\}}t|d t| | | dt| ks)J d S )Nr   )r   h㈵>)r   r   r   r   Abxpinfor   r   r   test_singular   s   

$r   c                     st   t  \} }tjd}|d}g   fdd}t| |||d  }g  t| ||d tt|d  d   d S )N90  r	   c                    s     |  d S )N)append)xktracer   r   trace_iterates#   s   z*test_x0_is_used_by.<locals>.trace_iterates)x0callback)r#   r   )r   r
   r   r   r   r   r   array_equal)r   r   r   r"   r!   trace_with_x0r   r   r   test_x0_is_used_by   s   

r&   c                  C   s\   t  \} }d}| |td  }t| ||d\}}t||\}}t|d t||dd d S )Ng      ?r	   )shiftr   r   )rtol)r   r
   eyer   r   r   )r   r   r'   	shifted_Ax1info1x2info2r   r   r   
test_shift-   s   

r/   c                  C   sZ   t  \} }d| d< d| d< tt t| |dd\}}W d   dS 1 s&w   Y  dS )z;Asymmetric matrix should raise `ValueError` when check=True   )r0      r1   )r1   r0   T)checkN)r   assert_raises
ValueErrorr   r   r   r   r   test_asymmetric_fail7   s   

"r5   c                  C   sv   t jd} d}| dd}t ||j}| d}| d}t||||dd }t|| | |t| ks9J d S )Nr   ư>   )r"   r(   r   r
   r   r   randndotr   r   r   )r   r(   ar   cxr   r   r   test_minres_non_default_x0@   s   

$r>   c                  C   s   t jd} d}| dd}t ||j}| d}| d}| dd}t ||j}t|||||dd }t|| | |t| ksGJ d S )Nr   r6   r7   Mr"   r(   r   r8   r   r(   r;   r   r<   mr=   r   r   r   "test_minres_precond_non_default_x0K   s   

$rC   c                  C   s   t jd} d}t d}t d}t d}| dd}t ||j}t|||||dd }t	|| | |t	| ks?J d S )Nr   r6   r	   r?   r   )
r
   r   r   r)   onesr9   r:   r   r   r   rA   r   r   r   test_minres_precond_exact_x0X   s   


$rE   )numpyr
   numpy.linalgr   numpy.testingr   r   r   scipy.sparse.linalg._isolver   pytestr   r3   r   r   r&   r/   r5   r>   rC   rE   r   r   r   r   <module>   s    
	
	