khJmddlZddlZddlZddlZddlmZmZmZmZddl m Z ddl m Z m Z mZmZdZdZGd d eZeZGd d eZeZGd deZy)N) LpSolver_CMDLpSolver subprocessPulpSolverError) constants)DictListOptionalTuplescipfscipceZdZdZdZ d dZejejejejejejejejejejejejejejejdZ ejejejhZ dZdZdZedZy) SCIP_CMDzThe SCIP optimization solverNc Ftj|||||||||| | |  y)a+ :param bool mip: if False, assume LP even if integer variables :param bool msg: if False, no log is shown :param list options: list of additional options to pass to solver :param bool keepFiles: if True, files are saved in the current directory and not deleted after solving :param str path: path to the solver binary :param float timeLimit: maximum time for solver (in seconds) :param float gapRel: relative gap tolerance for the solver to stop (in fraction) :param float gapAbs: absolute gap tolerance for the solver to stop :param int maxNodes: max number of nodes during branching. Stops the solving when reached. :param int threads: sets the maximum number of threads :param str logPath: path to the log file mipmsgoptionspath keepFiles timeLimitgapRelgapAbsmaxNodesthreadslogPathNr__init__) selfrrrrrrrrrrrs x/var/www/html/SchoolMeal/SchoolMeal/pds_admin_SchoolMeal/Backend/venv/lib/python3.12/site-packages/pulp/apis/scip_api.pyr zSCIP_CMD.__init__/66   )unknownzuser interruptznode limit reachedztotal node limit reachedzstall node limit reachedztime limit reachedzmemory limit reachedzgap limit reachedzsolution limit reachedz"solution improvement limit reachedzrestart limit reachedzoptimal solution found infeasible unboundedzinfeasible or unboundedc,|jtSN)executableExtension scip_pathr!s r" defaultPathzSCIP_CMD.defaultPathps'' 22r$c8|j|jSTrue if the solver is available executablerr,s r" availablezSCIP_CMD.availablestyy))r$c|j|jstd|jz|j|jddd\}}}|j |g}|j |jd|j d|jvr!|jd|jdd |jvr!|jd |jd d |jvr!|jd |jd d |jvr0t|jd dkDrtjd|js"tj|jdg}|j|j|jd|g|js|jdd|jvr |jd|jdgt|j }|D]Z}|j#drt%|} |j|| g3d|vrt%|} |d| z }|j|\|jdd|dg|jddg|jdd|dg|jddgt'|d5} | j)dj+|ddd|j-} t/j0|| | | | j3t4jj7|std |jz|j9|\} } i}|j;D]+}| j=|jd!||j<-|j?||jA| |jC|||| S#1swYxYw)""Solve a well formulated lp problemPuLP: cannot execute lpsolsetN limits/time=r limits/gap=rlimits/absgap=r limits/nodes=rrzUSCIP can only run with a single thread - use FSCIP_CMD for a parallel version of SCIP' does not allow a problem to be relaxed-s-qr-l-=z-czread ""optimizezwrite solution "quitw stdoutstderrPuLP: Error while executing )"r2rrcreate_tmp_filesnamewriteLPrappend optionsDictintwarningswarnrextendriterr startswithnextopenwritejoinget_piper check_callcloseosexistsreadsol variablesgetassignVarsVals assignStatusdelete_tmp_files)r!r8tmpLptmpSol tmpOptions file_optionscommandroptionargument options_filepipestatusvalues finalValsvs r" actualSolvezSCIP_CMD.actualSolvewstyy)!"9DII"EF F$($9$9"''4PU$V!vz 5"$ >> %   ,t~~.> ? @ t'' '   +d.>.>x.H-I J K t'' '   .1A1A(1K0L M N )) )   -0@0@0L/M N O (( (S1A1A)1L-MPQ-Q MMg xx MMTYYK'NO Ptyy!j)*xx NN4 (( ( NND$"2"29"=> ?t||$ ,F  %=12f$#G}H(n,F##F+ , ugQ/01j)* 0:;<f~& *c " 8l   tyy6 7 8}}gd4@   JJLww~~f%!"@499"LM Mf-  8A & 1663 7Iaff  8 )$  eVZ8 + 8 8s !PP ct|5} |j}|jd}|ddk(sJt|dk(sJ t jj|djtj}i} |j}|jd}|ddk(sJt|dk(sJt|dj|D]'} |j}t|d||d<)tj}||fcdddS#t$rt dwxYw#t$r||fcYcdddSwxYw#t d |xYw#1swYyxYw) zRead a SCIP solution filez: rzsolution statusrzCan't get SCIP solver status: robjective valueNzCan't read SCIP solver output: )r[readlinesplitlen Exceptionrr SCIP_STATUSESrestripr LpStatusUndefinedfloatLpStatusOptimal)filenameflinecompsrrrss r"rczSCIP_CMD.readsols(^% "q Qzz| 4(Qx#44445zQ&++//a )"="=FF &zz| 4(Qx#44445zQ&eAhnn&'  VV JJLE',U1XF58$ V..F6>K% "% " Q%(Fth&OPP Q &v~%1% "% ", &V),KD8*TUUA% "% "s`E6;D, AE6 AE#E6)$E# E6,EE6 E E6E  E6#E33E66E? NTFTNNNNNNN)__name__ __module__ __qualname____doc__rPr r rLpStatusNotSolvedrLpStatusInfeasibleLpStatusUnboundedr}NO_SOLUTION_STATUSESr-r3rv staticmethodrcr$r"rr*s & D  ( V..#55'99$-$?$?$-$?$?'99 ) ; ;&66"+"="=.7.I.I!*!#>M$ $$#### 3*IV'"'"r$rc6eZdZdZdZ ddZejejdZ ejejejhZ dZ dZdZeded eefd Zeded eefd Zeded eeeeffd Zed Zy) FSCIP_CMDzDThe multi-threaded FiberSCIP version of the SCIP optimization solverNc Ftj|||||||||| | |  y)a+ :param bool msg: if False, no log is shown :param bool mip: if False, assume LP even if integer variables :param list options: list of additional options to pass to solver :param bool keepFiles: if True, files are saved in the current directory and not deleted after solving :param str path: path to the solver binary :param float timeLimit: maximum time for solver (in seconds) :param float gapRel: relative gap tolerance for the solver to stop (in fraction) :param float gapAbs: absolute gap tolerance for the solver to stop :param int maxNodes: max number of nodes during branching. Stops the solving when reached. :param int threads: sets the maximum number of threads :param str logPath: path to the log file rNr) r!rrrrrrrrrrrs r"r zFSCIP_CMD.__init__r#r$)z No SolutionzFinal Solutionc,|jtSr))r* fscip_pathr,s r"r-zFSCIP_CMD.defaultPath)s'' 33r$c8|j|jSr/r1r,s r"r3zFSCIP_CMD.available,r4r$c|j|jstd|jz|j|jdddd\}}}}|j |g}|j |jd|j d|jvr!|jd |jdd |jvr!|jd |jd d |jvr!|jd |jd |js"tj|jdg}|jdg}|j|j|j||j||jd|g|jd|g|js|jdd|jvr |jd|jdgd|jvr!|jd|jdgt|j} | D]r} | j!drt#| } |j| | g3d| v} d| vrt#| } | d| z } | r|j| b|j| t|j%|t'|d5} | j)dj+|dddt'|d5}|j)dj+|ddd|j-}t/j0|||||j3t4jj7|std|jz|j9|\}}i}|j;D]+}|j=|jd||j<-|j?||jA||jC|||||S#1swYr?zNoPreprocessingInLC = TRUEr@z-fsolrArrBrz-sthrC/rDrHrIrJrMrN)"r2rrrOrPrQrrRrSrrUrVrWrrXrrYrZ silent_remover[r\r]r^rr_r`rarbrcrdrerfrgrh)r!r8rirjrk tmpParamsrlfile_parametersrmrrnrois_file_optionsrpparameters_filerqrrrsrtrus r"rvzFSCIP_CMD.actualSolve0styy)!"9DII"EF F/3/D/D GGT5%0 ,vz9 5"$ >> %   ,t~~.> ? @ t'' '   +d.>.>x.H-I J K t'' '   .1A1A(1K0L M N )) )   -0@0@0L/M N Oxx MMTYYK'NO P%';<tyy!y!uj)*()xx NN4 (( ( NND$"2"29"=> ? (( ( NNFt'7'7 'B&CE Ft||$ 3F  %=12#&-f$#G}H(n,F" ''/#**62% 3* 6" *c " 8l   tyy6 7 8 )S ! >_  ! !$))O"< = >}}gd4@   JJLww~~f%!"@499"LM Mf-  8A & 1663 7Iaff  8 )$  eVZC - 8 8 > >s+!Q !QQQstringreturnc`tjjD] \}}||vs |cSyr))rFSCIP_STATUSESitems)r fscip_status pulp_statuss r" parse_statuszFSCIP_CMD.parse_statuss8)2)A)A)G)G)I # %L+v%"" #r$c|jd}t|dk7ry|\}}|dk7ry|j} t|}|S#t$rYywxYw)N:rrx)rzr{r~r ValueError)rfieldslabel objectives r"parse_objectivezFSCIP_CMD.parse_objectiveslc" v;! !y % %OO%  i(I  s A AAc|j}t|dkry|dd\}} t|}||fS#t$rYywxYw)Nr)rzr{rr)rrrPvalues r"parse_variablezFSCIP_CMD.parse_variablesX v;?Raj e %LEU{  s 8 AAct|5}|j}tj|}|t d||tj vr |ifcdddS|j}tj |}|t d|i}|D]1}tj|}|t d||\} } | || <3||fcdddS#1swYyxYw)zRead a FSCIP solution fileNzCan't get FSCIP solver status: z"Can't get FSCIP solver objective: z Can't read FSCIP solver output: )r[ryrrrrrr) rfile status_linerrobjective_linerrd variable_linevariablerPrs r"rczFSCIP_CMD.readsols(^ %t--/K++K8F~%(G &WXX777rz % %"]]_N!11.AI %88JK +-I!% ( $33MB#):=:KL' e"' $ (9$= % % %sA C!A0CC$r)rrrrrPr r rrrrrrr-r3rvrstrr rTrrrr rrcrr$r"rrsN D  ( V!22#33N $$#### 4*TlSXc] " s xc5j0A'B   % %r$rceZdZdZdZ ddla d fd ZdZdZ dZ dZ d Z d Z xZS#e$r d Z d Z YxZSwxYw)SCIP_PYa The SCIP Optimization Suite (via its python interface) The SCIP internals are available after calling solve as: - each variable in variable.solverVar - each constraint in constraint.solverConstraint - the model in problem.solverModel rNc 8t |||||||||| |  y)a :param bool mip: if False, assume LP even if integer variables :param bool msg: if False, no log is shown :param list options: list of additional options to pass to solver :param float timeLimit: maximum time for solver (in seconds) :param float gapRel: relative gap tolerance for the solver to stop (in fraction) :param float gapAbs: absolute gap tolerance for the solver to stop :param int maxNodes: max number of nodes during branching. Stops the solving when reached. :param str logPath: path to the log file :param int threads: sets the maximum number of threads :param bool warmStart: if True, the solver will use the current value of variables as a start ) rrrrrrrrr warmStartN)superr ) r!rrrrrrrrrr __class__s r"r zSCIP_PY.__init__s52 G #!#  r$cxd|_|jj}tjtj tj tjtjtjtjtjtjtjtjtjtjtjtjd}d}||}||vr |jj}|jD]}||j|_ |jjD]-}|jj|j ||_/|tjk(r!|j%|tj&n0tj}|j%|tj(|S|j%||S#|j%|tj*Y|SxYw)NT)optimalr'r& inforunbd timelimit userinterrupt nodelimittotalnodelimitstallnodelimitgaplimitmemlimitsollimit bestsollimit restartlimitr%)rrrrrrrr) resolveOK solverModel getStatusr rrrrr getBestSol _variables solverVarvarValue constraintsrsgetSlacksolverConstraintslackrgLpSolutionOptimalLpSolutionIntegerFeasibleLpSolutionNoSolutionFound) r!r8solutionStatusscip_to_pulp_status possible_solution_found_statusesrrsolutionr constraints r"findSolutionValueszSCIP_PY.findSolutionValuessBL^^557N$44&88'::&88&88!*!>+B+B&77, (!:!:: 0K0KL!*!:!: 0S0STM'MQOOFI,O,OPMs C-H"H9cy)r0Trr,s r"r3zSCIP_PY.availableSsr$c8|jjy)zSolves the problem with scipN)rrFr!r8s r" callSolverzSCIP_PY.callSolverXs NN # # %r$c  tj|j|_|jt j k(r|jjn|jj|js|jj|j&|jjd|jd|jvr)|jjd|jdd|jvr)|jjd|jdd|jvr)|jjd|jdd |jvr(|jj|jd d |jvr>t|jd d kDr#t!j"d |jd |j$s"t!j"|jdt'|j(}|D]F}d|vr|j+dd\}}n |t-|}}|jj||Ht j.dt j0dt j2di}|j5D]l}|jj7|j||j8|j:|j<|j>jA|d|_!nt jDtFjHt jJtFjLt jNtFjPi}|jRjUD]n\}} |jjW|| jtjYd| jUD| jZ || _.p|jjAddr|jj_} |j5D]@}|j`|jjc| |jB|j`B|jje| yy)zU Takes the pulp lp model and translates it into a scip model Nz limits/timerz limits/gaprz limits/absgaprz limits/nodesrrrz The solver z" can only run with a single threadr?rDr)maxsplitBCIrN)rPvtypelbubobjc3@K|]\}}||jzywr))r).0r coefficients r" z+SCIP_PY.buildSolverModel..s'& 5+((*<*<<&s)consrPrF)3rModelrPrsenser LpMaximize setMaximize setMinimizer hideOutputrsetParamrS setLogfilerTrUrVrrXrrzrZLpBinary LpContinuous LpIntegerrdaddVarcatlowBoundupBoundrrerLpConstraintLEoperatorleLpConstraintGEgeLpConstraintEQeqrraddConsquicksumconstantrcreatePartialSolr setSolValaddSol) r!r8rrnrPrcategory_to_vtypevarsense_to_operatorrss r"buildSolverModelzSCIP_PY.buildSolverModel\s"ZZ0BNxx9///**,**, 88))+~~)'' t~~F4+++'' d6F6Fx6PQ4+++''9I9I(9STT---''8H8H8TUD,,,))$*:*:9*EFD,,,T5E5Ei5P1QTU1U !$)),NO88  +RST4<<(G! 5&="(,,sQ,"?KD%"($w-%D''e4  5""C&&##S!  ||~  " 5 5+CGG4||{{ ((c2 !6!  (((++(((++(((++!  %'NN$8$8$:  j.0nn.D.D<*:+;+;< &9C9I9I9K&$,,, /E / + ##K7NN335<<>QC||/00CMM3<<PQ%%a( 8r$c|j||j||j|}|jD] }d|_ |j j D] }d|_ |S)z Solve a well formulated lp problem creates a scip model, variables and constraints and attaches them to the lp model which it then solves F)rrrrmodifiedrrs)r!r8rrrs r"rvzSCIP_PY.actualSolvess  ! !" % OOB !44R8NMM *$)! * nn335 , &+ # ,! !r$c4td|jd)z Solve a well formulated lp problem uses the old solver and modifies the rhs of the modified constraints The z$ solver does not implement resolvingrrPrs r" actualResolvezSCIP_PY.actualResolves""tyyk!EF r$cy)r0Frr,s r"r3zSCIP_PY.availablesr$c4td|jd)r6rz solver is not availablerrs r"rvzSCIP_PY.actualSolves!D 3K"LM Mr$) TTNNNNNNNF)rrrrrP pyscipoptrr rr3rrrvr ImportError __classcell__)rs@r"rrs~ Ds $ L; z   &Z )x " I N  N Ns5 AAr)rrasysrUcorerrrrr typingr r r r r+rrSCIPrFSCIPrrr$r"r#sf6 EE..   @"|@"F^% ^%B hr$