o 5ήc;@sNddlmmZddlmZddlZddlZdZ ej j dZ e e dZeeZee e dZeeeZeejeZdgZddZejeejejdZd Ze d e e d ejdd ejd dej!ded<ej"eeedd\Z#Z$edZ%ejee#eeZ&e dede%de&e$dkre dded<ej'eeedddd\Z(Z$edZ)ejee(eeZ*e dedde)de*e$dkre dded<ej'eeedddd\Z+Z$edZ,ejee+eeZ-e dedde,de-e$dkr%e ddSdS)NzSPARSKIT/drivcav/e05r0200z(https://math.nist.gov/pub/MatrixMarket2/z.mtx.gzz _rhs1.mtx.gzcCs.tdd7<tjtddt|S)Nr )countsysstderrwriteAm)vr T/tmp/pip-target-vg8gfxp4/lib/python/scipy/sparse/linalg/_isolve/tests/demo_lgmres.pymatvecsr )r shapedtypedzMatrixMarket problem %szInvert z x rz matrix; nnz = g+=)restrttolzGMRES(z): z matvecs, relative residual: zDidn't converge )inner_mouter_krzLGMRES(z, 6) [same memory req.]: z, 6) [same subspace size]: ).scipy.sparse.linalgsparselinalglascipy.ioionumpynprproblemlib _datasource Repositorymmopenfmmreadtocsrrclosearrayravelbnormbnormrr LinearOperatorr rAMprintnnzgmresx0infocount_0err0lgmresx1count_1err1x2count_2err2r r r r sZ   (