(phQItSrSSKrSSKJr SSKJr SSKrSSKJ r Sr Sr Sr S r SS jrS rS rSS jrg)zn Routines for removing redundant (linearly dependent) equations from linear programming equality constraints. N)svd) interp_decompdtrsmcSn[R"[U5U:RSS95R 5$)a^ Counts the number of nonzeros in each row of input array A. Nonzeros are defined as any element with absolute value greater than tol = 1e-13. This value should probably be an input to the function. Parameters ---------- A : 2-D array An array representing a matrix Returns ------- rowcount : 1-D array Number of nonzeros in each row of A gvIh%<=axis)nparrayabssumflatten)Atols T/var/www/html/venv/lib/python3.13/site-packages/scipy/optimize/_remove_redundancy.py _row_countrs7" C 88SVc\&&A&. / 7 7 99cH[U5n[R"X!-5$)a Returns the index of the densest row of A. Ignores rows that are not eligible for consideration. Parameters ---------- A : 2-D array An array representing a matrix eligibleRows : 1-D logical array Values indicate whether the corresponding row of A is eligible to be considered Returns ------- i_densest : int Index of the densest row in A eligible for consideration )rr argmax)r eligibleRows rowCountss r _get_densestr#s&1 I 99Y- ..rcSnSn[U5S:HnU[R"U5SS24n[R"XS5(dSnSnU[R"U5nXX#4$)a Eliminates trivial equations from system of equations defined by Ax = b and identifies trivial infeasibilities Parameters ---------- A : 2-D array An array representing the left-hand side of a system of equations b : 1-D array An array representing the right-hand side of a system of equations Returns ------- A : 2-D array An array representing the left-hand side of a system of equations b : 1-D array An array representing the right-hand side of a system of equations status: int An integer indicating the status of the removal operation 0: No infeasibility identified 2: Trivially infeasible message : str A string descriptor of the exit status of the optimization. rNzbThere is a zero row in A_eq with a nonzero corresponding entry in b_eq. The problem is infeasible.)rr logical_notallclose)rbstatusmessagei_zeros r_remove_zero_rowsr#:sq4FG ]a F ".. ! #$A ;;qy! $ $> ".. !A   rcUupEX!n[SXFSSS9nUSUS-USUS-2U4'XsS-SnXCU4n XCS-S2U4==X- - ss'XE4$)Nr)lowerdiagr) pluperm_rvjLUpvpermulpivs rbg_update_denser1`sv EB IE a!!,ADQqS'Bt!tQwK A#$A T(CstQwKAEK 5Lrc SnSnSnSnSn[X5uppVUS:waXXV4$URup[[U55n [U 5n /n Sn Un[R "XU -4SS9n[R "US5 XSS2US24'[R "U5n[R"XU -[S 9n[R"UR[S 9n[R"USS9[R"U54nUSn U GHnSUU'US:aSUUS- 'U US- n[UXSS2U4US- 5n[RR%UUSS 9nUUnS n[S['U5U5HnUU[)UU-['U55n[+USS2U4R-5R/U55nUU:R15(dMgUU[R2"U5-nUU U'S UUU- ' GM UR4R/UR7S S55n[RR9U5n [+U5SU -- U:a SnUnXXV4s $U R;U5 UcGM['U 5X- :XdGM O [=[U55n![U![=U 5- 5n!UU!SS24UU!XV4$![ai [RR!USS2U 45nUunn[[U55n [#U5HunnU UU UsU U'U U'M GNf=f) Eliminates redundant equations from system of equations defined by Ax = b and identifies infeasibilities. Parameters ---------- A : 2-D sparse matrix An matrix representing the left-hand side of a system of equations rhs : 1-D array An array representing the right-hand side of a system of equations Returns ------- A : 2-D sparse matrix A matrix representing the left-hand side of a system of equations rhs : 1-D array An array representing the right-hand side of a system of equations status: int An integer indicating the status of the system 0: No infeasibility identified 2: Trivially infeasible message : str A string descriptor of the exit status of the optimization. References ---------- .. [2] Andersen, Erling D. "Finding all linearly dependent rows in large-scale linear programming." Optimization Methods and Software 6.3 (1995): 219-227. :0yE>rrThere is a linear combination of rows of A_eq that results in zero, suggesting a redundant constraint. However the same linear combination of b_eq is nonzero, suggesting that the constraints conflict and the problem is infeasible.NF)orderr)dtype)trans2Fr)r#shapelistranger zeros fill_diagonalarangeintonesbooleyer1 Exceptionscipylinalg lu_factor enumeratelu_solvelenminr transposedotanyrTreshapenormappendset)"rrhs true_ranktolapiv tolprimalr r! inconsistentmnr)rdr(A_orige js_candidatesjs_maskluir*r+r,i1i2pijsbatchj_index j_indicescbibarbnormkeeps" r_remove_redundancy_pivot_denserolsP@GI FG5L 07AF {v&& 77DA U1XA QA A F F !U3'AQaeH  AIIa1C0Mggm))6G  ryy| +B UF ! q5AacF @!A#A Vq!tWac:B\\ " "2q " 2 7 #QB/G7C s2w$?@IAalO--/33B78AG   ""w1-.!$! 0DDHHS[[Q/0EIINN3'E5z1U7#i/&F33 (SVq}-DUX uQx=D s1v D $'?CIv 66K @''!Q$0BEB%(^F#A,B)/VBZ&r F2J'  @s/K22A/M%$M%cSnSnSnSnSn[X5uppEUS:waXXE4$URupx[[U55n [U 5n [ [XwU-55n /n Un [ R R[ R RU5U45R5n[R"U5nU GHbnUSS2U 4nSX'US:aSXS- '[ R RRUR5U5RSS5n[U [ U 5- 5n[R "USS2U4R5R#U55U:R%5Sn['U5S:aUUSnUX'MUR(R#URSS55n[RR+U5n[!U5SU-- U:a SnUnXXE4s $U R-U5 GMe [ [U55n[U[ U 5- 5nU USS24UUXE4$) r3r4rrr5Nrr;r)r#r<r=r>rUrGsparsehstackrEtocscr r?rHspsolverNrRr rOnonzerorLrQrSrT)rrVrXrYr r!rZr[r\r)rkr]r^r_rcBrfrgrkr*rlrmrns r_remove_redundancy_pivot_sparserxs2BGI FG5L 07AF {v&& 77DA U1XA QA E!qSMA A F U\\--a0!45;;=A  A, adG q5AcF \\ ( ( : B B2q I !CF(^VVAaeH&&(,,R0 1G ; D D Fq I q6A:1Q4AADDDHHS[[Q/0EIINN3'E5z1u9% 1&F33 QT uQx=D s1v D $'?CIv 66rc[X5upp#US:waXX#4$[U5upEn[R"[5R nUR 5[ UR5-U-nURupX::aUSOSn [U 5U:Ga?USS2S4n [R"U 5US-:n [R"U 5(aA[R"[R"U RU55U:5(aSnSnO[R"[R"U RU55US-:5(aSnS nOp[X 5n[R"XSS 9n[R"X5n[U5upEnURupX::aUSOSn [U 5U:aGM?XX#4$) ar Eliminates redundant equations from system of equations defined by Ax = b and identifies infeasibilities. Parameters ---------- A : 2-D array An array representing the left-hand side of a system of equations b : 1-D array An array representing the right-hand side of a system of equations Returns ------- A : 2-D array An array representing the left-hand side of a system of equations b : 1-D array An array representing the right-hand side of a system of equations status: int An integer indicating the status of the system 0: No infeasibility identified 2: Trivially infeasible message : str A string descriptor of the exit status of the optimization. References ---------- .. [2] Andersen, Erling D. "Finding all linearly dependent rows in large-scale linear programming." Optimization Methods and Software 6.3 (1995): 219-227. rr;NgcAzDue to numerical issues, redundant equality constraints could not be removed automatically. Try providing your constraint matrices as sparse matrices to activate sparse presolve, try turning off redundancy removal, or try turning off presolve altogether.drr5r ) r#rr finfofloatepsmaxr<r rPrOrdelete)rrr r!UsVhr~rr[r\s_minr)ri_removes r_remove_redundancy_svdrjsB.a3A& {V$$1vHA" ((5/  C %%'CL 3 &C 77DAVAbEE" e*s  aeHvvay3:- vvl##rvvbffQUU1X.>.D'E'EF%G  66"&&q"S3Y. / /F8G 0 IIa * IIa "q6bww"Q7 e*s :   rcdSnSnSn[X5uppEUS:waXXE4$URupxUn Uc[RR U5n [ UR XS9up[R"XSU U -XU S5(dSnUn[U SU 5n X SS24n Xn XXE4$)agEliminates redundant equations from a system of equations. Eliminates redundant equations from system of equations defined by Ax = b and identifies infeasibilities. Parameters ---------- A : 2-D array An array representing the left-hand side of a system of equations rhs : 1-D array An array representing the right-hand side of a system of equations rank : int, optional The rank of A randomized: bool, optional True for randomized interpolative decomposition Returns ------- A : 2-D array An array representing the left-hand side of a system of equations rhs : 1-D array An array representing the right-hand side of a system of equations status: int An integer indicating the status of the system 0: No infeasibility identified 2: Trivially infeasible message : str A string descriptor of the exit status of the optimization. rrr5N)randr) r#r<r rH matrix_rankrrQrsorted)rrVrank randomizedr r!rZr[r\rvidxprojA2rhs2s r_remove_redundancy_idrs@FG5L 07AF {v&& 77DA A | II ! !! $acc16IC ;;sr7|d*CABL 9 9 Ra/C 6B 8D V $$r)N)NT)__doc__numpyr scipy.linalgrscipy.linalg.interpolativerrGscipy.linalg.blasrrrr#r1rorxrrrrrsK 4 #:*/.#!L |7~|7~Y!xD%r