kh-{rddlmZmZmZmZddlmZddlZddlZdZ GddeZ e Z da Gd d eZ y) )LpSolver LpSolver_CMD subprocessPulpSolverError) constantsNc~|jxs,t|jxst|jS)zyCheck whether lp is a MIP. From an XPRESS point of view, a problem is also a MIP if it contains SOS constraints.)isMIPlensos1sos2)lps z/var/www/html/SchoolMeal/SchoolMeal/pds_admin_SchoolMeal/Backend/venv/lib/python3.12/site-packages/pulp/apis/xpress_api.py_ismipr!s* 88: 5RWW 5RWW5cleZdZdZdZ d dZdZdZdZe dZ dZ e d Z y) XPRESSzOThe XPRESS LP solver that uses the XPRESS command line tool in a subprocessNc Htj|||||||||| | | |  y)a Initializes the Xpress solver. :param bool mip: if False, assume LP even if integer variables :param bool msg: if False, no log is shown :param float timeLimit: maximum time for solver (in seconds) :param float gapRel: relative gap tolerance for the solver to stop (in fraction) :param heurFreq: the frequency at which heuristics are used in the tree search :param heurStra: heuristic strategy :param coverCuts: the number of rounds of lifted cover inequalities at the top node :param preSolve: whether presolving should be performed before the main algorithm :param options: Adding more options, e.g. options = ["NODESELECTION=1", "HEURDEPTH=5"] More about Xpress options and control parameters please see https://www.fico.com/fico-xpress-optimization/docs/latest/solver/optimizer/HTML/chapter7.html :param bool warmStart: if True, then use current variable values as start ) gapRelmipmsg timeLimitoptionspath keepFilesheurFreqheurStra coverCutspreSolve warmStartN)r__init__) selfrrrrrrrrrrr r!s rr"zXPRESS.__init__/s9>   rc$|jdS)N optimizer)executableExtensionr#s r defaultPathzXPRESS.defaultPath^s'' 44rc8|j|jS)True if the solver is available) executablerr's r availablezXPRESS.availableastyy))rc  |j|jstd|jz|j|jddddd\}}}}}|j |d|j }|jjd d rL|Dcgc]/}|j|j|jf1} }|j|| g} t|r|j r| jgd tjtjtjtjtjtjtj tjd } d} n| jddgtj"tj tjtjtjtj$tjtj"tjd } d} t'|d5} |j(s| j+d| j+d|j-|zdz|j.| j+d|j.z|jjd}|| j+d|dd|jjd}|| j+d|z|jjd}|| j+d|z|jjd}|| j+d|z|jjd }|| j+d!|z|jjd d r&| j+d"|j-|zdz|j0D]}| j+|dzt|r|j r| j+d#n| j+d$| j+d%|j-|zdz| j+d&|j-|d'| D]}| j+d(|d)|d*| j+d+| j+d,d d d t'|5} d }d }d }|j(sVt2j4d-k\r!t6j8}t6j8}n"t6j:}t6j<}d.}t7j>|j|jgd.| ||d./}|r|j@D]}|jCd0k7rtd1|jz d d d |jE||\}}}}}|jG||||| j|j| d2tj}|jI||jK||jM||jO||jQ||Scc}w#1swYxYw#1swYxYw)3"Solve a well formulated lp problemzPuLP: cannot execute rprtcmdattrslxr)writeSOSrr!FN) mipobjval bestbound mipstatusr rrr6lpobjvallpstatus r rrr8r9r:r;r<wz OUTPUTLOG=0 z readprob  z MAXTIME=%d rz MIPRELSTOP=frz HEURFREQ=%d rzHEURSTRATEGY=%d rz COVERCUTS=%d r z PRESOLVE=%d z readslxsol z mipoptimize z lpoptimize z writeprtsol z set fh [open z w]; list z puts $fh "z=$z" z close $fh zQUIT iT)shellstdinstdoutstderruniversal_newlinesr zPuLP: Error while executing ))r+rrcreate_tmp_filesnamewriteLPr optionsDictgetvalue writeslxsolrextendrLpStatusUndefinedLpStatusInfeasibleLpStatusOptimalLpStatusNotSolvedLpStatusUnboundedopenrwrite quote_pathrrsys hexversionrDEVNULLPIPESTDOUTPopenrFwaitreadsoldelete_tmp_filesassignVarsVals assignVarsDjassignConsSlack assignConsPi assignStatus) r#rtmpLptmpSoltmpCmdtmpAttrtmpStart variablesvstart attrNames statusmap statuskeyr0 targetGaprrrr optionr1consumesuboutsuberrxpress_valuesredcostslacksdualsattrsstatuss r actualSolvezXPRESS.actualSolveestyy)!"9DII"EF F373H3H GGT5%4 0vvwJJuqdhhJ?      U 32;UQqwwy?Taffaggi(UEU   Xu - ":$((   D E..........//,,.. I$I   j*5 6..,,//............ I#I &# ' #88 /* IIkDOOE$::TA B~~) .4>>9:((,,X6I$ K !}B78''++J7H# /H45''++J7H# -89((,,[9I$ *Y67''++J7H# /H45##K7 -$//(*CCdJK,, ) &4-( )bzdhh /* .) IIntv'>>E F II 89E " : JtfBtfC89 : IIm $ IIh O' P&\ RSGFF88>>Z/'//F'//F(__F'..F"G%%BGG$#' F A{{}!%&Dtyy&PQQ"= R@15 VW0M- eVVW=uyyB79T9TU &!   6"   aV@' ' P R Rs&Y-.Y-5JY2CY?2Y<?Zc(i}i}i}i}t|5}t|D]\}}|dkr |dk(r-|j}t|d} t|d} >|dkrD|dk(rJ|j} t | dkDsi| ddk(r(| d} t | d|| <t | d|| <t | ddk(s| dd vs| d} t | d|| <t | d|| < d d d t } t|5}|D]} | jjd }t |dk(s1|dj|dk(sK|dj} t|dj}|| |dj< d d d ||||| fS#1swYxYw#t$r0 t |dj}n#t$rYnwxYwYlwxYw#1swYXxYw) zRead an XPRESS solution filer;rr: rr Cr9LGREN=) rW enumeratesplitintr floatdictstriplower ValueError)filenameattrfilerzr{r|r}rClineno_linerowscolslinerKr~fieldsrOs rrazXPRESS.readsols1 (^" 9q!*1! 9 A:q[!KKMEuQx=DuQx=Db[r\ !;;=D4y1}7c>$(7D+0a>F4L,1$q'NGDM a\Q.47f3D$(7D+0a>F4L*/Q.E$KC! 9" 9H (^ 5q 5++C0v;!#q (9VAY(F"1IOO-E! #F1IOO$5 6 05E&)//+, 5 5wu44e" 9" 9X&!!$)&)//*;$G? ?HHHHHc t|d5}t|D]<\}}|jd|z|D]\}}|jd|d|dd >|jddddy#1swYyxYw) a Write a solution file in SLX format. The function can write multiple solutions to the same file, each solution must be passed as a list of (name,value) pairs. Solutions are written in the order specified and are given names "solutionN" where N is the index of the solution in the list. :param string name: file name :param list values: list of lists of (name,value) pairs rAzNAME solution%d z C  z.16frBzENDATA N)rWrrX)r#rKrzr2isolrOs rrPzXPRESS.writeslxsols$_ "#F+ @3 -12#&@KD%IIad|2>?@ @ IIj !  " " "s AA22A;c<dtjdd|zdzS)z Quotes a path for the Xpress optimizer console, by wrapping it in double quotes and escaping the following characters, which would otherwise be interpreted by the Tcl shell: \ $ " [ "z ([\\$"[])z\\\1)resub)rs rrYzXPRESS.quote_path-s!RVVL'4883>>r) TTNNNFNNNNNF) __name__ __module__ __qualname____doc__rKr"r(r,r staticmethodrarPrYrrrr)sy D   - ^5*yv8585t"$??rrcxeZdZdZdZ d dZdZddZdZddZ dZ dd Z e d Z d Zd Zy) XPRESS_PYz0The XPRESS LP solver that uses XPRESS Python APINc || t| }tj|||||| |||||  d|_| |_y)a Initializes the Xpress solver. :param bool mip: if False, assume LP even if integer variables :param bool msg: if False, no log is shown :param float timeLimit: maximum time for solver (in seconds) :param float gapRel: relative gap tolerance for the solver to stop (in fraction) :param heurFreq: the frequency at which heuristics are used in the tree search :param heurStra: heuristic strategy :param coverCuts: the number of rounds of lifted cover inequalities at the top node :param preSolve: whether presolving should be performed before the main algorithm :param bool warmStart: if set then use current variable values as warm start :param string export: if set then the model will be exported to this file before solving :param options: Adding more options. This is a list the elements of which are either (name,value) pairs or strings "name=value". More about Xpress options and control parameters please see https://www.fico.com/fico-xpress-optimization/docs/latest/solver/optimizer/HTML/chapter7.html N) rrrrrrrrr r!)absrr" _available_export) r#rrrrrrrr r!exportrs rr"zXPRESS_PY.__init__AsW@  YI   rc|j- ddlatjdd|_|jS|jS#d|_Y|jSxYw)r*Nr FT)rrxsetOutputEnabledr's rr,zXPRESS_PY.availablezsX ?? " (''."&t ("'s A Ac |j}|jD] }d|_ |jj D] }d|_ |j a|j j jdr|j|j dn|j|j |||t|r|js|jy|jy#tjtjtj f$r}t#t%|d}~wwxYw)aPerform the actual solve from actualSolve() or actualResolve(). :param prepare: a function that is called with `lp` as argument and allows final tweaks to `lp.solverModel` before the low level solve is started. FNz.lpl) solverModelrmmodified constraintsrzrrendswithrXrr lpoptimizesolverx ModelErrorInterfaceError SolverErrorrstr)r#rpreparemodelrncerrs r callSolverzXPRESS_PY.callSolvers ,NNE\\^ #"  #^^**, #"  #||'<<%%'007KK c2KK -" bz$((  " !!6#8#8&:L:LM ,!#c(+ + ,sC,D/D2E 2EE c  |j}t|r|jrggddf\}}}} |j||tj tj tj tj tj tj tjtj d}d}nggggf\}}}} |j||||tjtjtj tj tj tjtj tjtj d }d}|H|j|jD cic]} | j|| jd!c} |H|j|jD cic]} | j|| jd!c} |M|j!|j"j%D cic]\} } | || jdc} } |M|j'|j"j%D cic]\} } | || jdc} } |j)|j+|tj } |j-| | S#d\}}YxYw#d\}}}}YxYwcc} wcc} wcc} } wcc} } w#t.j0t.j2t.j4f$r} t7t9| d} ~ wwxYw)N)NNr7r6)NNNNr?r>r )rrr getmipsolrrRrSrTgetlpsolrUrVrcrmrK_xprsrdrfritemsrerN getAttribrgrxrrrrr)r#rrxr|r}djsrqrrrnnrrrs rfindSolutionValueszXPRESS_PY.findSolutionValuess> ,NNEbzdhh(*Bd(:%65#+OOAv.!22 22 22 22 22 33 00 22  ( )+BB%65#CNN1feS9 !22 00 33 22 22 22 22 22 22  ' }!!",,."QQ1661QWWQZ=#8"QR",,. QQQWWQZ!8 QR 8:8L8L8NOfq!Qaggaj))O!""9;9M9M9OPv1Qqwwqz**P]] *I,G,GF OOF #Mg+ *IAv"C,B)Avuc#R QPQ!!6#8#8&:L:LM ,!#c(+ + ,sy.LK!B L K.!B+L $K=0&L$L:1L+L 2L:L A L!K+(L. K:7L2MMMc|jsd} ddl}t ||j ||j |||j|S#t$r}t|}Yd}~t |d}~wwxYw)r.zXPRESS Python API not availabler N)r,rx ImportErrorrrbuildSolverModelrr)r#rrmessagerxrs rrzXPRESS_PY.actualSolvesy~~7G #"'* * b! G$&&r**  #c(!'* * #sA A? A::A?c|j| |j}dD]6\}}|jj|d}|%|j ||8|j D]}t |tr |d}|d}n]|jdd}t|dk7rtdt|z|dj}|dj} |j |t||jjdd rt!}t!} t#|j%d D]P} | j'|j)| j'| j)| j*dRt-|rT|j.rHt||j0j2k(r|j5|n(|j7|| d n|j9|ddd|j:rd } |j=| yy#t$rYnwxYw |j |t|#t$rYnwxYw|j ||#t>j@t>jBt>jDf$r} tt| d} ~ wwxYw)zP Takes the pulp lp model and translates it into an xpress model ))r MIPRELSTOP)rMAXTIME)rHEURFREQ)r HEURSTRATEGY)r COVERCUTS)r PRESOLVENr rrrzInvalid option r!Fc |jdSNr rrs rz,XPRESS_PY.buildSolverModel..2saggajrkey warmstartc&|dkDr t|yyr)print)probdatarmsgtypes rrz+XPRESS_PY.buildSolverModel..messageCs{c #r)#_extractrrMrN setControlr isinstancetuplerr rrrrrrlistsortedrmrOappendrrr attributesr loadmipsol addmipsol loadlpsolr addcbmessagerxrrr) r#rrrrKrOrtrsolvalcolindrnrrs rrzXPRESS_PY.buildSolverModels bK ,NNE 2 T((,,S$7$$$T51 2,, .fe,!!9D"1IE#\\#q1F6{a'-.?#f+.MNN!!9??,D"1IOO-E$$T3u:6 .,##K7 4HI2Awwy, aggi0 aggaj12":$((6{e&6&6&;&;;((0 DOOFD$=xx#""7+ 7"$$T5<8!  u-6!!6#8#8&:L:LM ,!#c(+ + ,sn2JBJ"H3=A!JCJ3 H?<J>H??JI J I,)J+I,,J2K5K  Kc t}t}t|jD]o}|j|}|jst |ds t d|j |jd|j |j qt|dkDr|jj||t}t}t} |jD]} | jst | ds t d|j | jd|j d| j | jtj n | j|j | jd|j d| j | j tjn | j t|dkDr|jj#||| |j%|||j'|S#tj(tj*tj,f$r} t t/| d} ~ wwxYw)z>Resolve a problem that was previously solved by actualSolve().rzCannot add new constraintsr zCannot add new variablesLNG)rrrrhasattrrrrconstantr rchgrhsrmlowBoundrxinfinityupBound chgboundsrrrrrr) r#rrrhsindrhsvalrKconbndindbndtypebndvalrnrs r actualResolvezXPRESS_PY.actualResolveKs& ,VFVFr~~. -nnT*||sG,)*FGG ciil+ s||m, -6{Q%%ff5VFfGVF\\^ Szzq'*)*DEE aggaj)s# !**2Dv.!**U aggaj)s# 1Bfoo R S7|a((&A OOB (**2. .!!6#8#8&:L:LM ,!#c(+ + ,sH6H992J+I??Jct|dr t|d|jD]}t|dst|d|jj D]}t|dst|dy)z,Reset any XPRESS specific information in lp.rrN)rdelattrrmrrz)rrnrs r_resetzXPRESS_PY._resetuss 2} % B & $Aq'"7# $&&( $Aq'"7# $rc  |j| tj}|jtj k(r|j tjt}t}t}t}t}|jD]'}|j|jtj n |j|j|jtjn |j|j|jj|d|j tj"k(r|jdn@|j tj$k(r|jdn|jd|j|j&*|j)|dgt+|dzzgg||||t-t/|j|j1D]\} \}} | | f|_t} t-t5|j6D]\} } |j6| }tj9dt5|j;d D}|j< }|jtj>k(r'tjA|tjB| }n|jtjDk(r'tjA|tjF| }ne|jtjHk(r'tjA|tjJ| }n!tMd tO|jz| j| ||ft+| d kDsj|jQ| Dcgc]\}}}| c}}| D]\} }}| |f|_t} t+| dkDr9|jQ| Dcgc]\}}}| c}}| D]\} }}| |f|_d}|||jRd|||jTd||_+ycc}}wcc}}w#tjXtjZtj\f$r*}|j|tMtO|d}~wwxYw)aExtract a given model to an XPRESS Python API instance. The function stores XPRESS specific information in the `solverModel` property of `lp` and each variable and constraint. These information can be removed by calling `_reset`. NgIBrr rc3FK|]\}}||jdzyw)rNr).0ras r z%XPRESS_PY._extract..s(!1 N!s!c&|djdSrrrs rrz$XPRESS_PY._extract..s!A$**Q-rr)bodysenserhszUnsupprted constraint type dc ^g}t|D]}g}g}||jD]4\}}|j|jd|j|6|jtj |||t |t|r|j|yy)z"Extract sos constraints from PuLP.r N) rrrrrxsosrr addSOS) msosdictsostypesoslistrKindicesweightsrnvals raddsosz"XPRESS_PY._extract..addsoss#7OUD G G")$-"5"5"7,3qwwqz2s+,NN6::gwT#ST Uw<HHW% rr)/rrxproblemr r LpMaximize chgobjsensemaximizerrmrrrr objectiverNcat LpIntegerLpBinaryrKaddcolsr rzip getVariablerrrSumrrLpConstraintLE constraintleqLpConstraintGEgeqLpConstraintEQeqrr addConstraintr rrrrr)r#rrobjlbubctypenamesrnjrconsrrKrlhsr rryrrs rrzXPRESS_PY._extracts BW ,NN$Exx9///!!&//2 &CBBFEFE\\^ % ajj.@6??*ajjQ QYY->&//AIIN 2<<++As3455I///LL%UUi000LL%LL% QVV$ % MM#sc#hl3RRUE R&s2<<>5;L;L;N'OP ! 6Aqa& !6D$VBNN%;< "4nnT*jj! &syy{8O P!||m99 8 88))s&**#)NAYY)":"::))s&**#)NAYY)":"::))s&)))MA)5CIIF Q3K(t9s?''$(?wq!Q(?@%)+ 1c%&F +6D1 "24y1}##d$;71aQ$;N:RR R ,AR3 R AR R2S7 %S22S7c8|jj|S)z`Get an arbitrary attribute for the model that was previously solved using actualSolve().)rr)r#rwhichs r getAttributezXPRESS_PY.getAttributes~~''..r) TTNNNNNNNNN)N)rrrrrKr"r,rrrrrrrrr7rrrrr<sv: D  7r ,@?,B +P,d(,T $ $_,B/rr)corerrrrrrZrrr XPRESS_CMDrxrrrrr;sE6FE 6K?\K?\ i/i/r