3
Od4                 @   s.  d dl jjZd dljZd dlZd dlZdZ	ej
jjdZejde	 Zejej Zej  ejde	 Zejejej Zej  d gZdd ZejeejejdZd	Zed
e	  edejd  ejd ejf  d ed < ej eeedd\Z!Z"ed  Z#ejj$ee! e ejj$e Z%ede e#de% e"d krBed d ed < ej&eeed ddd\Z'Z"ed  Z(ejj$ee' e ejj$e Z)eded  e(de) e"d kred d ed < ej&eeed ddd\Z*Z"ed  Z+ejj$ee* e ejj$e Z,eded  e+de, e"d kr*ed dS )    NzSPARSKIT/drivcav/e05r0200z&ftp://math.nist.gov/pub/MatrixMarket2/z	%s.mtx.gzz%s_rhs1.mtx.gzc             C   s,   t d  d7  < tjjdt d   t|  S )Nr      z%d)countsysstderrwriteAm)v r	   `/var/www/html/virt/lib64/python3.6/site-packages/scipy/sparse/linalg/isolve/tests/demo_lgmres.pymatvec   s    r   )r   shapedtyped   zMatrixMarket problem %szInvert %d x %d matrix; nnz = %dr   g+=)Zrestrttolz
GMRES(%d):zmatvecs, residualzDidn't converge      )Zinner_mZouter_kr   z LGMRES(%d,6) [same memory req.]:zmatvecs, residual:z"LGMRES(%d,6) [same subspace size]:   r   )-Zscipy.sparse.linalgsparseZlinalgZlaZscipy.ioioZnumpynpr   problemlibZ_datasourceZ
RepositorymmopenfZmmreadZtocsrr   closearrayZravelbr   r   ZLinearOperatorr   r   AMprintZnnzZgmresZx0infoZcount_0ZnormZerr0Zlgmresx1Zcount_1Zerr1Zx2Zcount_2Zerr2r	   r	   r	   r
   <module>   sP   
  
 
 
