ó ʽ÷Xc@`sddlmZmZmZdZgZddlmZddlm Z m Z m Z m Z m Z mZddlmZmZmZidd6dd6dd6d d6dd6dd6d d6d d6dd6d d6dd6d d6d d6d d6d d6d d6Zd „Zd „Zd d „Zd S(i(tdivisiontprint_functiontabsolute_importsrestructuredtext en(twarn(t asanyarraytasarraytasmatrixtarraytmatrixtzeros(taslinearoperatortLinearOperatortIdentityOperatortftdtFtDcC`s8|dkrd}n|dkr*d}nt||fS(NtfdFDR(t _coerce_rules(txty((s?/tmp/pip-build-7oUkmx/scipy/scipy/sparse/linalg/isolve/utils.pytcoerces     cC`s|S(N((R((s?/tmp/pip-build-7oUkmx/scipy/scipy/sparse/linalg/isolve/utils.pytidsc `sß|}t|ƒ}|jd|jdkrEtd|jfƒ‚n|jd}tˆƒ‰ˆj|dfkp‚ˆj|fks”tdƒ‚nˆjjdkr¸ˆjdƒ‰n‡fd†}|dkrt|dƒrñ|jj}n|j ˆƒjj}t |ˆjjƒ}nCt d t ƒ|dkrFˆjj}n|dkratd ƒ‚nt ˆd|ƒ‰ˆjƒ‰|dkr t|d|ƒ}nTt|d|ƒ}|j|dfkpÖ|j|fksètd ƒ‚n|jƒ}|dkrŸt|d ƒr|j} nt} t|d ƒr<|j} nt} | tkrx| tkrxtd|jd|jƒ}qÌt|jd| d| d|jƒ}n-t|ƒ}|j|jkrÌtdƒ‚n|||ˆ|fS(s‡Make 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 xtype : {'f', 'd', 'F', 'D', None}, optional dtype of the x vector 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 dimensionsRRc`s.tˆtƒrt|ƒ}n|jˆjƒS(N(t isinstanceRRtreshapetshape(R(tb(s?/tmp/pip-build-7oUkmx/scipy/scipy/sparse/linalg/isolve/utils.pyt postprocessSstdtypesTUse of xtype argument is deprecated. Use LinearOperator( ... , dtype=xtype) instead.s#xtype must be 'f', 'd', 'F', or 'D's$A and x have incompatible dimensionstpsolvetrpsolveRtmatvectrmatvecs/matrix and preconditioner have different shapesN(R Rt ValueErrorRRtchartastypetNonethasattrRRRtDeprecationWarningRtravelR RRRRR R ( tAtMtx0RtxtypetA_tNRRRR((Rs?/tmp/pip-build-7oUkmx/scipy/scipy/sparse/linalg/isolve/utils.pyt make_system"sZ!   '     '     N(R R (R R(R R(R R(RR (RR(RR(RR(RR (RR(RR(RR(RR (RR(RR(RR(t __future__RRRt __docformat__t__all__twarningsRtnumpyRRRRRR tscipy.sparse.linalg.interfaceR R R RRRR$R.(((s?/tmp/pip-build-7oUkmx/scipy/scipy/sparse/linalg/isolve/utils.pyts.