_MhwHddlmZmZmZddlZddlmZddlmZm Z m Z m Z m Z m Z ddlmZmZmZmZddlmZddlZddlZddlmZd Z ddlmZn #e$rd ZYnwxYwejagd ZGd d e Z!dZ"dZ#ddZ$ ddZ% ddZ&dZ' ddZ(dZ)dZ*dS))warncatch_warnings simplefilterN)asarray)issparseSparseEfficiencyWarning csr_array csc_array eye_array diags_array)is_pydata_spmatrixconvert_pydata_sparse_to_scipyget_index_dtypesafely_cast_index_arrays) LinAlgError)_superluFT) use_solverspsolvespluspilu factorizedMatrixRankWarningspsolve_triangularis_sptriangular spbandwidthceZdZdS)rN)__name__ __module__ __qualname__d/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/scipy/sparse/linalg/_dsolve/linsolve.pyrrsDr"rc d|vr|dt_tjr!d|vrtj|ddSdSdS)aq Select default sparse direct solver to be used. Parameters ---------- useUmfpack : bool, optional Use UMFPACK [1]_, [2]_, [3]_, [4]_. over SuperLU. Has effect only if ``scikits.umfpack`` is installed. Default: True assumeSortedIndices : bool, optional Allow UMFPACK to skip the step of sorting indices for a CSR/CSC matrix. Has effect only if useUmfpack is True and ``scikits.umfpack`` is installed. Default: False Notes ----- The default sparse solver is UMFPACK when available (``scikits.umfpack`` is installed). This can be changed by passing useUmfpack = False, which then causes the always present SuperLU based solver to be used. UMFPACK requires a CSR/CSC matrix to have sorted column/row indices. If sure that the matrix fulfills this, pass ``assumeSortedIndices=True`` to gain some speed. References ---------- .. [1] T. A. Davis, Algorithm 832: UMFPACK - an unsymmetric-pattern multifrontal method with a column pre-ordering strategy, ACM Trans. on Mathematical Software, 30(2), 2004, pp. 196--199. https://dl.acm.org/doi/abs/10.1145/992200.992206 .. [2] T. A. Davis, A column pre-ordering strategy for the unsymmetric-pattern multifrontal method, ACM Trans. on Mathematical Software, 30(2), 2004, pp. 165--195. https://dl.acm.org/doi/abs/10.1145/992200.992205 .. [3] T. A. Davis and I. S. Duff, A combined unifrontal/multifrontal method for unsymmetric sparse matrices, ACM Trans. on Mathematical Software, 25(1), 1999, pp. 1--19. https://doi.org/10.1145/305658.287640 .. [4] T. A. Davis and I. S. Duff, An unsymmetric-pattern multifrontal method for sparse LU factorization, SIAM J. Matrix Analysis and Computations, 18(1), 1997, pp. 140--158. https://doi.org/10.1137/S0895479894246905T. Examples -------- >>> import numpy as np >>> from scipy.sparse.linalg import use_solver, spsolve >>> from scipy.sparse import csc_array >>> R = np.random.randn(5, 5) >>> A = csc_array(R) >>> b = np.random.randn(5) >>> use_solver(useUmfpack=False) # enforce superLU over UMFPACK >>> x = spsolve(A, b) >>> np.allclose(A.dot(x), b) True >>> use_solver(useUmfpack=True) # reset umfPack usage to default useUmfpackassumeSortedIndices)r&N)r%uumfpack configure)kwargss r#rr sd|vl+ |M-77f5J.KLLLLLLMM77r"ctjtjfdtjtjfdtjtjfdtjtjfdi}t t|jj}t t|jjj} |||f}n+#t$r}d|d|d}t||d}~wwxYw|d d z}tj |}tj |j tj |_ tj |jtj |_||fS) z8Get umfpack family string given the sparse matrix dtype.dizidlzlz]only float64 or complex128 matrices with int32 or int64 indices are supported! (got: matrix: z , indices: )Nrldtype)npfloat64int32 complex128int64getattrr3nameindicesKeyError ValueErrorcopyrindptr)A _familiesf_typei_typefamilyemsgA_news r#_get_umf_familyrHcs9 RX !4 RX !4 IR & &FR- . .F%FF+, %%%T7=TTJPTTToo1$%AY_F IaLLE:ahbh777ELJqy999EM 5=s' B22 C<CCc  t|}|r|jnd}t|}t|}t|r |jdvs&t |}t dtdt|}|st|}|j dkp|j dko|j ddk}| | }tj|j|j}|j|kr||}|j|kr||}|j \} } | | krt#d| | fd| |j d kr&t#d |j d |j d dt%t&d st( t&_|o t&j}|r|r|r|} n|} t| |j } t(rt1d|jjdvrt#dt5|\} }t7j| } | t6j|| d}n|r|r|}d}|s|jdkrd}nd }|jtj d}|j!tj d}tE|}tGj$| |j%|j&|||||\}}|d kr6t dtNd|(tj)|r|}ntU|}|jdks5t|s&t dtdt |}g}g}g}tW|j dD]}|dd|f}||}tj,|}|j d }|-||-tj.||t^ |-tj|||j tj0|}tcte|j }tj0|| }tj0|| } |||| ff|j |j}|r|3|}|S)aSolve the sparse linear system Ax=b, where b may be a vector or a matrix. Parameters ---------- A : ndarray or sparse array or matrix The square matrix A will be converted into CSC or CSR form b : ndarray or sparse array or matrix The matrix or vector representing the right hand side of the equation. If a vector, b.shape must be (n,) or (n, 1). permc_spec : str, optional How to permute the columns of the matrix for sparsity preservation. (default: 'COLAMD') - ``NATURAL``: natural ordering. - ``MMD_ATA``: minimum degree ordering on the structure of A^T A. - ``MMD_AT_PLUS_A``: minimum degree ordering on the structure of A^T+A. - ``COLAMD``: approximate minimum degree column ordering [1]_, [2]_. use_umfpack : bool, optional if True (default) then use UMFPACK for the solution [3]_, [4]_, [5]_, [6]_ . This is only referenced if b is a vector and ``scikits.umfpack`` is installed. Returns ------- x : ndarray or sparse array or matrix the solution of the sparse linear equation. If b is a vector, then x is a vector of size A.shape[1] If b is a matrix, then x is a matrix of size (A.shape[1], b.shape[1]) Notes ----- For solving the matrix expression AX = B, this solver assumes the resulting matrix X is sparse, as is often the case for very sparse inputs. If the resulting X is dense, the construction of this sparse result will be relatively expensive. In that case, consider converting A to a dense matrix and using scipy.linalg.solve or its variants. References ---------- .. [1] T. A. Davis, J. R. Gilbert, S. Larimore, E. Ng, Algorithm 836: COLAMD, an approximate column minimum degree ordering algorithm, ACM Trans. on Mathematical Software, 30(3), 2004, pp. 377--380. :doi:`10.1145/1024074.1024080` .. [2] T. A. Davis, J. R. Gilbert, S. Larimore, E. Ng, A column approximate minimum degree ordering algorithm, ACM Trans. on Mathematical Software, 30(3), 2004, pp. 353--376. :doi:`10.1145/1024074.1024079` .. [3] T. A. Davis, Algorithm 832: UMFPACK - an unsymmetric-pattern multifrontal method with a column pre-ordering strategy, ACM Trans. on Mathematical Software, 30(2), 2004, pp. 196--199. https://dl.acm.org/doi/abs/10.1145/992200.992206 .. [4] T. A. Davis, A column pre-ordering strategy for the unsymmetric-pattern multifrontal method, ACM Trans. on Mathematical Software, 30(2), 2004, pp. 165--195. https://dl.acm.org/doi/abs/10.1145/992200.992205 .. [5] T. A. Davis and I. S. Duff, A combined unifrontal/multifrontal method for unsymmetric sparse matrices, ACM Trans. on Mathematical Software, 25(1), 1999, pp. 1--19. https://doi.org/10.1145/305658.287640 .. [6] T. A. Davis and I. S. Duff, An unsymmetric-pattern multifrontal method for sparse LU factorization, SIAM J. Matrix Analysis and Computations, 18(1), 1997, pp. 140--158. https://doi.org/10.1137/S0895479894246905T. Examples -------- >>> import numpy as np >>> from scipy.sparse import csc_array >>> from scipy.sparse.linalg import spsolve >>> A = csc_array([[3, 2, 0], [1, -1, 0], [0, 5, 1]], dtype=float) >>> B = csc_array([[2, 0], [-1, 0], [2, 0]], dtype=float) >>> x = spsolve(A, B) >>> np.allclose(A.dot(x).toarray(), B.toarray()) True Ncsccsrz.spsolve requires A be CSC or CSR matrix format stacklevelrz!matrix must be square (has shape r0rz!matrix - rhs dimension mismatch (z - r'r2Scikits.umfpack not installed.dD\convert matrix data to double, please, using .astype(), or set linsolve.useUmfpack.u = FalseT autoTransposeFrK)r>)ColPerm)optionszMatrix is exactly singularzCspsolve is more efficient when sparse b is in the CSC matrix format)maxval)shaper3)4r __class__rrformatr rrrndimrXsum_duplicates _asfptyper4 promote_typesr3astyper=hasattrr%noScikitr'toarrayravel RuntimeErrorcharrHr(UmfpackContextlinsolve UMFPACK_Ar;intcr?dictrgssvnnzdatarfillnanrrange flatnonzeroappendfullint concatenatermaxfrom_scipy_sparse)!r@b permc_spec use_umfpackis_pydata_sparsepydata_sparse_cls b_is_sparse b_is_vector result_dtypeMNb_vec umf_familyumfxflagr;r?rVinfo Afactsolve data_segsrow_segscol_segsjbjxjwsegment_length sparse_data idx_dtype sparse_row sparse_cols! r#rrszd*!,,'7A T&q))A&q))A QKK4AH66 aLL = $ 4 4 4 41++K  AJJFaKEQVq[%DQWQZ1_K A#AGQW55Lw, HH\ " "w, HH\ " " 7DAq QFaVFFFGGGAGAJVQWVVQRVVVWWW :s # #$#| .*,KF;{F;  IIKKEEEQW---3355  A?@@ @ 7>>GAtqyy13DQRSSSSrv GGII$AJH%%);A)>)>%3,<<<<aLLIHH171:&& C Cqqq!tW__&&,,..Z^^N2&&!"""" E E EFFF  BqE!A!A!ABBBB.33K's17||<<>> import numpy as np >>> from scipy.sparse import csc_array >>> from scipy.sparse.linalg import splu >>> A = csc_array([[1., 0., 0.], [5., 0., 2.], [0., -1., 0.]], dtype=float) >>> B = splu(A) >>> x = np.array([1., 2., 3.], dtype=float) >>> B.solve(x) array([ 1. , -3. , -1.5]) >>> A.dot(B.solve(x)) array([ 1., 2., 3.]) >>> B.solve(A.dot(x)) array([ 1., 2., 3.]) clsc:|t|SNrwr ras r#csc_construct_funcz splu..csc_construct_func((A77 7r"rK&splu converted its input to CSC formatrMrNcan only factor square matricesSuperLU)DiagPivotThreshrU PanelSizeRelaxNrUNATURALT SymmetricModeFrilurVr typeto_scipy_sparsetocscr rrZrrr\r]rXr=rr4rirjupdatergstrfrlrm) r@rydiag_pivot_threshrelax panel_sizerVA_clsrrrr;r?_optionss r#rrHsH!'Q', 8 8 8 8 8     % % ' '& QKK4AH-- aLL 5 $ 4 4 4 4 A 7DAq Q:;;;.q"'9EEOGV$5z(777H    y(($(! >!QUAFGV-?#X 7 7 77r"c t|r= 1.0) for ILU. (default: 10) drop_rule : str, optional Comma-separated string of drop rules to use. Available rules: ``basic``, ``prows``, ``column``, ``area``, ``secondary``, ``dynamic``, ``interp``. (Default: ``basic,area``) See SuperLU documentation for details. Remaining other options Same as for `splu` Returns ------- invA_approx : scipy.sparse.linalg.SuperLU Object, which has a ``solve`` method. See also -------- splu : complete LU decomposition Notes ----- To improve the better approximation to the inverse, you may need to increase `fill_factor` AND decrease `drop_tol`. This function uses the SuperLU library. Examples -------- >>> import numpy as np >>> from scipy.sparse import csc_array >>> from scipy.sparse.linalg import spilu >>> A = csc_array([[1., 0., 0.], [5., 0., 2.], [0., -1., 0.]], dtype=float) >>> B = spilu(A) >>> x = np.array([1., 2., 3.], dtype=float) >>> B.solve(x) array([ 1. , -3. , -1.5]) >>> A.dot(B.solve(x)) array([ 1., 2., 3.]) >>> B.solve(A.dot(x)) array([ 1., 2., 3.]) rc:|t|Srrrs r#rz!spilu..csc_construct_funcrr"rKz'spilu converted its input to CSC formatrMrNrr) ILU_DropRule ILU_DropTolILU_FillFactorrrUrrNrUrTrrr)r@drop_tol fill_factor drop_ruleryrrrrVrrrrr;r?rs r#rrsv!'Q', 8 8 8 8 8     % % ' '& QKK4AH-- aLL 6 $ 4 4 4 4 A 7DAq Q:;;;.q"'9EEOGV#.$5z(777H    y(($(! >!QUAFGV-?"H 6 6 66r"ctr&ttdst t_tjrt rtdtr j dks&ttdtd jjdvrt!dt#\}t%j|fd }|St+jS) a Return a function for solving a sparse linear system, with A pre-factorized. Parameters ---------- A : (N, N) array_like Input. A in CSC format is most efficient. A CSR format matrix will be converted to CSC before factorization. Returns ------- solve : callable To solve the linear system of equations given in `A`, the `solve` callable should be passed an ndarray of shape (N,). Examples -------- >>> import numpy as np >>> from scipy.sparse.linalg import factorized >>> from scipy.sparse import csc_array >>> A = np.array([[ 3. , 2. , -1. ], ... [ 2. , -2. , 4. ], ... [-1. , 0.5, -1. ]]) >>> solve = factorized(csc_array(A)) # Makes LU decomposition. >>> rhs1 = np.array([1, -2, 0]) >>> solve(rhs1) # Uses the LU factors. array([ 1., -2., -2.]) r'rPrKrrMrNrQrRctjdd5tj|d}dddn #1swxYwY|S)Nignore)divideinvalidTrS)r4errstatesolver(rh)rxresultr@rs r#rzfactorized..solveLsHh??? P P7#4a$OO P P P P P P P P P P P P P P PMs$AA A )r rrr`r%rar'rdrrZr rrr]r3rer=rHr(rfnumericrr)r@rrrs` @r#rrsT<!(     % % ' ' :s # #$#| |  A?@@ @  8E 1 1! A 9(Q 8 8 8 8 KKMM 7>> import numpy as np >>> from scipy.sparse import csc_array >>> from scipy.sparse.linalg import spsolve_triangular >>> A = csc_array([[3, 0, 0], [1, -1, 0], [2, 0, 1]], dtype=float) >>> B = np.array([[2, 0], [-1, 0], [2, 0]], dtype=float) >>> x = spsolve_triangular(A, B) >>> np.allclose(A.dot(x), B) True rrLTrKz?CSC or CSR matrix format is required. Converting to CSC matrix.rMrNz+A must be a square matrix but its shape is .rrNrz&A is singular: zero entry on diagonal.)rrMz)b must have 1 or 2 dims but its shape is zlThe size of the dimensions of A must be equal to the size of the first dimension of b but the shape of A is z and the shape of b is r2)r3rZzA is singular.)$r rrrrZrrrr r>rXr=rrsetdiagdiagonalr4anyrr r\ asanyarrayr[r^r3float32r_r rgstrsrlrmr;r?reshapelen)r@rxlower overwrite_A overwrite_b unit_diagonaltransrrdiaginvdiagrLUrrs r#rrXsp!(     % % ' ' E{{qx5(( C  QKKAH-- N $ 4 4 4 4 aLL  FFHH 7DAqAvv D!' D D DFF F /      #: ; ; ; IIaLLL               zz|| 6$!)   :8:: :D& C<<K(((AA{7+++.A aAvV B B B BDD DAGAJ :w : :/0w : : :   #B$4QWbj$I$I17SSLw, HH\ " "w, HH\ " "  FFHH   q!fL 1 1 1 a|E : : :  ! nUqvqy!(qvqy!(  GAt ,*+++ !'/"Bs17||a/?(@BBB K Hs-+D$$D(+D(c*t|r |jdvs&tdtdt |}|jdkr:|jdk|jdkfS|jdkr:|j\}}||k ||k fS|jdkrVtd | Dtd | DfS|jd krZtd t|j D}td t|j D}||fS|j |j}}t|dz }d\}}|dzddfD]V}|||||dz}|o||k }|o||k }|s|sdSWt!jt!j|t!j|}|}|o||k }|o||k }|jdkr||fS||fS)aReturns 2-tuple indicating lower/upper triangular structure for sparse ``A`` Checks for triangular structure in ``A``. The result is summarized in two boolean values ``lower`` and ``upper`` to designate whether ``A`` is lower triangular or upper triangular respectively. Diagonal ``A`` will result in both being True. Non-triangular structure results in False for both. Only the sparse structure is used here. Values are not checked for zeros. This function will convert a copy of ``A`` to CSC format if it is not already CSR or CSC format. So it may be more efficient to convert it yourself if you have other uses for the CSR/CSC version. If ``A`` is not square, the portions outside the upper left square of the matrix do not affect its triangular structure. You probably want to work with the square portion of the matrix, though it is not requred here. Parameters ---------- A : SciPy sparse array or matrix A sparse matrix preferrably in CSR or CSC format. Returns ------- lower, upper : 2-tuple of bool .. versionadded:: 1.15.0 Examples -------- >>> import numpy as np >>> from scipy.sparse import csc_array, eye_array >>> from scipy.sparse.linalg import is_sptriangular >>> A = csc_array([[3, 0, 0], [1, -1, 0], [2, 0, 1]], dtype=float) >>> is_sptriangular(A) (True, False) >>> D = eye_array(3, format='csr') >>> is_sptriangular(D) (True, True) )rKrLcoodiadoklilzCis_sptriangular needs sparse and not BSR format. Converting to CSR.rMrNrrrrc3(K|] \}}||kVdSrr!.0rcs r# z"is_sptriangular..s*//da16//////r"c3(K|] \}}||kVdSrr!rs r#rz"is_sptriangular..s*4Q4Q1Q!V4Q4Q4Q4Q4Q4Qr"rc32K|]\}}|D] }||kV dSrr!rrowcolscols r#rz"is_sptriangular..7RR93TRRcC3JRRRRRRRr"c32K|]\}}|D] }||kV dSrr!rs r#rz"is_sptriangular..rr"r)TTr)FFrL)rrZrrr offsetsrvmincoordsallkeys enumeraterowsr?r;rr4repeatarangediff) r@rrrupperr?r;rrs r#rrsR QKKAH(RRR R $ 4 4 4 4 aLL x5y}}!#QY]]__%999 U  X d !!##ddl%7%7%9%999 U  //affhh/////4Q4Q4Q4Q4Q1Q1QQQ U  RR !&0A0ARRRRRRR !&0A0ARRRRRe|h GF F aALE5Q2  vc{6#'?23-3$;++---3$;++-- U << 9RYq\\276?? 3 3D D  *tt|((**E  *tt|((**Ex5e| %<r"ct|r |jdvs&tdtdt |}|jdkrst d|j t d|jfS|jdvrk|j |j }}t|dz }tj tj|tj||z }|jd kr| }ns|jd kr|jd|jdz }nL|jd krAd |Ddgz}t| t |fSt tj| dt tj|dfS) abReturn the lower and upper bandwidth of a 2D numeric array. Computes the lower and upper limits on the bandwidth of the sparse 2D array ``A``. The result is summarized as a 2-tuple of positive integers ``(lo, hi)``. A zero denotes no sub/super diagonal entries on that side (tringular). The maximum value for ``lo``(``hi``) is one less than the number of rows(cols). Only the sparse structure is used here. Values are not checked for zeros. Parameters ---------- A : SciPy sparse array or matrix A sparse matrix preferrably in CSR or CSC format. Returns ------- below, above : 2-tuple of int The distance to the farthest non-zero diagonal below/above the main diagonal. .. versionadded:: 1.15.0 Examples -------- >>> import numpy as np >>> from scipy.sparse.linalg import spbandwidth >>> from scipy.sparse import csc_array, eye_array >>> A = csc_array([[3, 0, 0], [1, -1, 0], [2, 0, 1]], dtype=float) >>> spbandwidth(A) (2, 0) >>> D = eye_array(3, format='csr') >>> spbandwidth(D) (0, 0) )rKrLrrrzCspbandwidth needs sparse format not LIL and BSR. Converting to CSR.rMrNrrrJrrLrrcg|] \}}||z  Sr!r!rs r# zspbandwidth..gs ,,,41aA,,,r")rrZrrr rvrritemr?r;rr4rrrrr)r@r?r;rgaps r#rr2sH QKKAH(KKK R $ 4 4 4 4 aLL x51qy}}++---..Aqy}}7K7K7M7M0N0NNNx>!!(AI KK!Oi ! bgfoo66@ 8u  $C U  hqkAHQK' U  ,,16688,,,s2Cy#c(("" s   """A & &BF3KK,<,<,>,>(B(B BBr")NT)NNNNN)NNNNNNNN)TFFF)+warningsrrrnumpyr4r scipy.sparserrr r r r scipy.sparse._sputilsr rrr scipy.linalgrr> threadingrrascikits.umfpackr( ImportErrorlocalr%__all__ UserWarningrrrHrrrrrrrr!r"r#r s7777777777HHHHHHHHHHHHHHHHNNNNNNNNNNNN$$$$$$  %%%%%%%HHHY_    X X X        AMAMAMF   D@ @ @ @ F04.2f7f7f7f7RJNGK_6_6_6_6DBBBJIN%*G G G G TMMM`7C7C7C7C7CsAAA