glddlZddlZddlZddlZddlmZmZmZmZddlm Z m Z ddl m Z ddl mZmZmZmZGdd eZeZGd d eZeZGd d eZy)N) LpSolver_CMDLpSolver subprocessPulpSolverError) scip_path fscip_path) constants)DictListOptionalTupleceZdZdZdZ 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 O/var/www/dash_apps/app1/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)executableExtensionrr!s r" defaultPathzSCIP_CMD.defaultPathms'' 22r$c8|j|jSTrue if the solver is available executablerr+s r" availablezSCIP_CMD.availableptyy))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+|dddt-j.|t0j2t0j4t6jj9|std |jz|j;|\} } i} |j=D]+}| j?|jd!| |j<-|jA| |jC| |jE|||| 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 )#r1rrcreate_tmp_filesnamewriteLPrappend optionsDictintwarningswarnrextendriterr startswithnextopenwritejoinr check_callsysrJrKosexistsreadsol variablesgetassignVarsVals assignStatusdelete_tmp_files)r!r7tmpLptmpSol tmpOptions file_optionscommandroptionargument options_filestatusvalues finalValsvs r" actualSolvezSCIP_CMD.actualSolvetstyy)!"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 8gcjjLww~~f%!"@499"LM Mf-  8A & 1663 7Iaff  8 )$  eVZ8  8 8s !O::Pc t|5} |j}|jd}|ddk(sJt|dk(sJ t jj|djtj}i}|t jvr ||fcdddS |j}|jd}|ddk(sJt|dk(sJt|dj|D]'} |j}t|d||d<)||fcdddS#t$rt dwxYw#t$rt d |wxYw#t d |xYw#1swYyxYw) zRead a SCIP solution filez: rzsolution statusr zCan't get SCIP solver status: rNobjective valuez!Can't get SCIP solver objective: zCan't read SCIP solver output: )rZreadlinesplitlen Exceptionrr SCIP_STATUSESrcstripr LpStatusUndefinedNO_SOLUTION_STATUSESfloat)filenameflinecompsrorps r"razSCIP_CMD.readsols(^$ "q Qzz| 4(Qx#44445zQ&++//a )"="=FF666v~!$ "$ "& Tzz| 4(Qx#44445zQ&eAhnn&'  VV JJLE',U1XF58$ V6>I$ "$ " Q%(Fth&OPP Q$ T%(I$&RSS TV),KD8*TUUE$ "$ "sTF;D; AF+AEF$E1,F;EFE..F1FFF  NTFTNNNNNNN)__name__ __module__ __qualname____doc__rOr r r|LpStatusNotSolvedLpStatusOptimalLpStatusInfeasibleLpStatusUnboundedrzr}r,r2rs staticmethodrar$r"rr's & D  ( V..#55'99$-$?$?$-$?$?'99 ) ; ;&66"+"="=.7.I.I!*!#>M$ $$#### 3*CJ&"&"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*r r+s r"r,zFSCIP_CMD.defaultPaths'' 33r$c8|j|jSr.r0r+s r"r2zFSCIP_CMD.available"r3r$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+|dddt-j.||jrt0j2nt,j4|jrt0j6nt,j4t8jj;|std|jz|j=|\}}i}|j?D]+}|jA|jd||j<-|jC||jE||jG|||||S#1swYnxYw#1swYFxYw)r5r6r7r8r9prmNr:rr;rr<rr=r>zNoPreprocessingInLC = TRUEr?z-fsolr@rrArz-sthrB/rCrGrHrIrLrM)$r1rrrNrOrPrrQrRrrTrUrVrrWrrXrY silent_removerZr[r\rr]r^rJDEVNULLrKr_r`rarbrcrdrerf)r!r7rgrhri tmpParamsrjfile_parametersrkrrlrmis_file_optionsrnparameters_filerorprqrrs r"rszFSCIP_CMD.actualSolve&styy)!"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"< = > !%3::z/A/A!%3::z/A/A ww~~f%!"@499"LM Mf-  8A & 1663 7Iaff  8 )$  eVZC + 8 8 > >s+!Q8 !R8RRstringreturnc`tjjD] \}}||vs |cSyr))rFSCIP_STATUSESitems)r fscip_status pulp_statuss r" parse_statuszFSCIP_CMD.parse_status{s8)2)A)A)G)G)I # %L+v%"" #r$c|jd}t|dk7ry|\}}|dk7ry|j} t|}|S#t$rYywxYw)N:r ru)rwrxr{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 )rwrxr~r)rrrOvalues 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: )rZrvrrrr}rr) rfile status_lineroobjective_linerrb variable_linevariablerOrs r"razFSCIP_CMD.readsols(^ %t--/K++K8F~%(G &WXX777rz % %"]]_N!11.AI %88JK +-I!% ( $33MB#):=:KL' e"' $ (9$= % % %sA C!A0CC$r)rrrrrOr r rrrrrr}r,r2rsrstrrrSrr~rrrrarr$r"rrsN D  ( V!22#33N $$#### 4*SjSXc] " 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 constraintsrpgetSlacksolverConstraintslackreLpSolutionOptimalLpSolutionIntegerFeasibleLpSolutionNoSolutionFound) r!r7solutionStatusscip_to_pulp_status possible_solution_found_statusesrosolutionr constraints r"findSolutionValueszSCIP_PY.findSolutionValues sBL^^557N$44&88'::&88&88!*!>+B+B&77, (!:!:: 0K0KL!*!:!: 0S0STM'MQOOFI,O,OPMs C-H"H9cy)r/Trr+s r"r2zSCIP_PY.availableHsr$c8|jjy)zSolves the problem with scipN)rrEr!r7s r" callSolverzSCIP_PY.callSolverMs 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>rCr )maxsplitBCIrM)rOvtypelbubobjc3@K|]\}}||jzywr))r).0r coefficients r" z+SCIP_PY.buildSolverModel..s'& 5+((*<*<<&s)consrOrF)3scipModelrOrsenser LpMaximize setMaximize setMinimizer hideOutputrsetParamrR setLogfilerSrTrUrrWrrwrYLpBinary LpContinuous LpIntegerrbaddVarcatlowBoundupBoundrrcrLpConstraintLEoperatorleLpConstraintGEgeLpConstraintEQeqrraddConsquicksumconstantrcreatePartialSolr setSolValaddSol) r!r7rrlrOrcategory_to_vtypevarsense_to_operatorrss r"buildSolverModelzSCIP_PY.buildSolverModelQs"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)rrrrmodifiedrrp)r!r7rrrs r"rszSCIP_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 resolvingrrOrs r" actualResolvezSCIP_PY.actualResolves""tyyk!EF r$cy)r/Frr+s r"r2zSCIP_PY.availablesr$c4td|jd)r5rz solver is not availablerrs r"rszSCIP_PY.actualSolves!D 3K"LM Mr$) TTNNNNNNNF)rrrrrO pyscipoptrr rr2rrrsr ImportError __classcell__)rs@r"rrs~ Ds $ L; z   &Z )x " I N  N Ns5 AAr)rr_r^rTcorerrrrrr r typingr r rrrSCIPrFSCIPrrr$r"r s_6 EE'..y"|y"x]% ]%@ hr$