+ 0i R^RIt^RIHt^RIHtR.tRRRRRR RR RR RR R /RlltR#)N)_get_atol_rtol) make_systemtfqmrrtolgh㈵>atolmaxiterMcallbackshowFcpVPp \P!V \P4'd\p VP V 4p\P!VP\P4'dVP V 4p\ WW!4wrr\PPV4R8XdVP4p V ^3#VP^,p Vf\RV ^ ,4pVfVP4p MWPV 4, p T p V P4pT pVPVPV 44pTpR;p;pp\P!VP4V 4PpTp\P !V4pTpV^8XdV ^3#\#RVWC4wpp\%V4EF6pV^,^8HpV'dL\P!VP4V4pVR8XdV R3u#VV, pV VV,, pVXV,,pV V^,V, V,V,,p\PPV4V, p\P !R^V^,,, 4pVVV,,pV^,V,pVPV4p V VV ,, p Ve V!V 4V\P !V^,4,V8d$V'd\'RV^, 24V ^3u#V'g\P!VP4V4pVV, p!VV!V ,,p V!V,V!^,V,,pVPVPV 44pVV, pEKVPVPX44pTp TpEK9 V'd\'RX^, 24W3#)a) Solve ``Ax = b`` with the Transpose-Free Quasi-Minimal Residual method. Parameters ---------- A : {sparse array, ndarray, LinearOperator} The real or complex N-by-N matrix of the linear system. Alternatively, `A` can be a linear operator which can produce ``Ax`` using, e.g., `scipy.sparse.linalg.LinearOperator`. b : {ndarray} Right hand side of the linear system. Has shape (N,) or (N,1). x0 : {ndarray} Starting guess for the solution. rtol, atol : float, optional Parameters for the convergence test. For convergence, ``norm(b - A @ x) <= max(rtol*norm(b), atol)`` should be satisfied. The default is ``rtol=1e-5``, the default for ``atol`` is ``0.0``. maxiter : int, optional Maximum number of iterations. Iteration will stop after maxiter steps even if the specified tolerance has not been achieved. Default is ``min(10000, ndofs * 10)``, where ``ndofs = A.shape[0]``. M : {sparse array, ndarray, LinearOperator} Inverse of the preconditioner of A. M should approximate the inverse of A and be easy to solve for (see Notes). Effective preconditioning dramatically improves the rate of convergence, which implies that fewer iterations are needed to reach a given error tolerance. By default, no preconditioner is used. callback : function, optional User-supplied function to call after each iteration. It is called as ``callback(xk)``, where ``xk`` is the current solution vector. show : bool, optional Specify ``show = True`` to show the convergence, ``show = False`` is to close the output of the convergence. Default is `False`. Returns ------- x : ndarray The converged solution. info : int Provides convergence information: - 0 : successful exit - >0 : convergence to tolerance not achieved, number of iterations - <0 : illegal input or breakdown Notes ----- The Transpose-Free QMR algorithm is derived from the CGS algorithm. However, unlike CGS, the convergence curves for the TFQMR method is smoothed by computing a quasi minimization of the residual norm. The implementation supports left preconditioner, and the "residual norm" to compute in convergence criterion is actually an upper bound on the actual residual norm ``||b - Axk||``. References ---------- .. [1] R. W. Freund, A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems, SIAM J. Sci. Comput., 14(2), 470-482, 1993. .. [2] Y. Saad, Iterative Methods for Sparse Linear Systems, 2nd edition, SIAM, Philadelphia, 2003. .. [3] C. T. Kelley, Iterative Methods for Linear and Nonlinear Equations, number 16 in Frontiers in Applied Mathematics, SIAM, Philadelphia, 1995. Examples -------- >>> import numpy as np >>> from scipy.sparse import csc_array >>> from scipy.sparse.linalg import tfqmr >>> A = csc_array([[3, 2, 0], [1, -1, 0], [0, 5, 1]], dtype=float) >>> b = np.array([2, 4, -1], dtype=float) >>> x, exitCode = tfqmr(A, b, atol=0.0) >>> print(exitCode) # 0 indicates successful convergence 0 >>> np.allclose(A.dot(x), b) True ri'rg?z:TFQMR: Linear solve converged due to reach TOL iterations z@TFQMR: Linear solve not converged due to reach MAXIT iterations )dtypenp issubdtypeint64floatastyperlinalgnormcopyshapeminmatvecinner conjugaterealsqrtrrangeprint)"Abx0rrr r r r rxndofsruwrstarvuhatdthetaetarhorhoLastr0normtau_iterevenvtrstaralphauNextczbetas"&&&$$$$$$ _/var/www/html/photoedit/myenv/lib/python3.14/site-packages/scipy/sparse/linalg/_isolve/tfqmr.pyrr sYh GGE }}UBHH%% HHUO }}QWWbhh'' HHUOQ2)JA! yy~~aB FFH1v GGAJEeURZ( z FFH  O A A E !A DA ((5??$a ( - -CG WWS\F C {1v Wfd9GD!gax1} hhu0!4G"}2w'ME ME UT\ E!S(1, , q!C' GGB!eQh,' ( uqy!tun HHQK S1W    QK a 4 '$$(F8-.q6M((5??,a0C=DD1H At tQw!m+A88AHHQK(D IA88AHHUO,DAGSV   F8% & <)N)numpyr iterativerutilsr__all__rr=r<rCsO% )jjBjjjj#jr=