o
    5ήc	                     @   s   d dl Zd dlmZmZmZ d dlmZ d dl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)assert_equalassert_allcloseassert_)minres)norm)raises   )assert_normclosec                  C   s8   t jd t jdd} | | j } t jd}| |fS )N  
   )nprandomseedrandT)matrixvector r   T/tmp/pip-target-vg8gfxp4/lib/python/scipy/sparse/linalg/_isolve/tests/test_minres.pyget_sample_problem
   s
   
r   c                  C   sJ   t  \} }d| d< d|d< t| |\}}t|d t| ||dd d S )Nr   )r   h㈵>tol)r   r   r   r	   dotAbxpinfor   r   r   test_singular   s   

r   c                     sv   t  \} }tjd tj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%   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   

r3   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=Truer   )r      r4   )r4   r   T)checkN)r   assert_raises
ValueErrorr   r   r   r   r   test_asymmetric_fail8   s   

"r8   c                  C   sp   t jd d} t jdd}t ||j}t jd}t jd}t|||| dd }t|||| d d S )Nr
   ư>   )r&   r   r   r   r   r   r   randnr   r   r   r	   )r   ar   cxr   r   r   test_minres_non_default_x0A   s   r@   c                  C   s   t jd d} t jdd}t ||j}t jd}t jd}t jdd}t ||j}t||||| dd }t|||| d d S )Nr    r9   r:   Mr&   r   r   r   r;   r   r=   r   r>   mr?   r   r   r   "test_minres_precond_non_default_x0L   s   rE   c                  C   sx   t jd d} t d}t d}t d}t jdd}t ||j}t||||| dd }t	|||| d d S )Nr
   r9   r   rA   r   r   )
r   r   r   r-   onesr<   r   r   r   r	   rC   r   r   r   test_minres_precond_exact_x0Y   s   


rG   )numpyr   numpy.testingr   r   r   scipy.sparse.linalg._isolver   scipy.linalgr   pytestr   r6   test_iterativer	   r   r   r*   r3   r8   r@   rE   rG   r   r   r   r   <module>   s    
	
	