ó î&]\c@`sddlmZmZmZdZgZddlmZmZm Z m Z m Z m Z ddl mZmZmZidd 6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6dd6Zd „Zd „Zd „Zd S(i(tdivisiontprint_functiontabsolute_importsrestructuredtext en(t asanyarraytasarraytasmatrixtarraytmatrixtzeros(taslinearoperatortLinearOperatortIdentityOperatortftdtFtDcC`s8|dkrd}n|dkr*d}nt||fS(NtfdFDR (t _coerce_rules(txty((s?lib/python2.7/site-packages/scipy/sparse/linalg/isolve/utils.pytcoerces     cC`s|S(N((R((s?lib/python2.7/site-packages/scipy/sparse/linalg/isolve/utils.pytidsc `s|}t|ƒ}|jd|jdkrEtd|jfƒ‚n|jd}tˆƒ‰ˆj|dfkp‚ˆj|fks”tdƒ‚nˆjjdkr¸ˆjdƒ‰n‡fd†}t|dƒrå|jj}n|jˆƒjj}t |ˆjjƒ}t ˆd|ƒ‰ˆj ƒ‰|dkrNt |d|ƒ}nTt|d|ƒ}|j|dfkp„|j|fks–td ƒ‚n|j ƒ}|dkrMt|d ƒrÉ|j} nt} t|d ƒrê|j} nt} | tkr&| tkr&td |jd|jƒ}qzt|jd | d| d|jƒ}n-t|ƒ}|j|jkrztdƒ‚n|||ˆ|fS(s8Make a linear system Ax=b Parameters ---------- A : LinearOperator sparse or dense matrix (or any valid input to aslinearoperator) M : {LinearOperator, Nones} preconditioner sparse or dense matrix (or any valid input to aslinearoperator) x0 : {array_like, None} initial guess to iterative method b : array_like right hand side Returns ------- (A, M, x, b, postprocess) A : LinearOperator matrix of the linear system M : LinearOperator preconditioner x : rank 1 ndarray initial guess b : rank 1 ndarray right hand side postprocess : function converts the solution vector to the appropriate type and dimensions (e.g. (N,1) matrix) iis(expected square matrix, but got shape=%ss$A and b have incompatible dimensionsRR c`s.tˆtƒrt|ƒ}n|jˆjƒS(N(t isinstanceRRtreshapetshape(R(tb(s?lib/python2.7/site-packages/scipy/sparse/linalg/isolve/utils.pyt postprocessPstdtypes$A and x have incompatible dimensionstpsolvetrpsolveRtmatvectrmatvecs/matrix and preconditioner have different shapesN(R Rt ValueErrorRRtchartastypethasattrRRRtraveltNoneRRRRRR R ( tAtMtx0RtA_tNRtxtypeRRR((Rs?lib/python2.7/site-packages/scipy/sparse/linalg/isolve/utils.pyt make_system!sL   '  '     N(R R (R R (R R(R R(R R (R R (R R(R R(RR (RR (RR(RR(RR (RR (RR(RR(t __future__RRRt __docformat__t__all__tnumpyRRRRRRtscipy.sparse.linalg.interfaceR R R RRRR,(((s?lib/python2.7/site-packages/scipy/sparse/linalg/isolve/utils.pyts.