bi dZddlZddlmZddlmZddlmZddl m Z ddl m Z dZ ddZdefdZ dd Z dd Z dd ed ed edefdZy)a+ Copyright, the CVXPY authors Licensed under the Apache License, Version 2.0 (the "License"); you may not use this file except in compliance with the License. You may obtain a copy of the License at http://www.apache.org/licenses/LICENSE-2.0 Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License. N)Minimize)failure_solutionc tjjtd|j|j Dcgc] }| c}zScc}w)zEvaluates lazy constraints.r)problemsproblemProblemr constraints_lazy_constraints)rcs ]/home/cdr/jupyterlab/.venv/lib/python3.12/site-packages/cvxpy/reductions/solvers/bisection.py_lower_problemr sH    # # QK   0I0I"J113"J J LL"JsA returnctj5tjdd|j|dddy#1swYyxYw)Nignorez.*DPP.*)messagesolver)warningscatch_warningsfilterwarningssolve)rrs r _solver sA  "%*= V $%%%s *AAcf|duxs,|jtjtjfvS)N)statuss INFEASIBLEINFEASIBLE_INACCURATE)rs r _infeasibler's6 d? Jgnn121H1H1JJJc||jrdnd}||jrdnd}|j}|j}t|D]}|sr||_t |} t | |t | r|}|dz}6| jtjvrd}n"tjd| jz|ss||_t |} t | |t | rd}nF| jtjvr|}|dz}tjd| jz|s|s||fcStjd) z Finds an interval for bisection.rTSolver failed with status %srzOUnable to find suitable interval for bisection; your problem may be unbounded..) is_nonneg is_nonposrangevaluer rrrrSOLUTION_PRESENTerror SolverError) rtrlowhigh max_itersinfeasible_low feasible_high_lowereds r _find_bisection_intervalr41sb {;;=ab |KKMqq[[]NKKMM 9 AG$W-G 7F #7# 1#5#55 $ ''2W^^CEEAG$W-G 76 *7#!%1#5#55q''2W^^CEE m9 =@    . //rc d} d} t| D]*} ||ksJ| ||z |kr| ||fcS||zdz } |r8| | zdk(r0td| |fztd| |fztd| | fz| |_t|}t ||t |r!|r| | zdk(rtd | z|| }|j tjvr9|r"| | zdk(rtd | |jfz|j} || }|r td tjd |j ztjd )z&Bisect `problem` on the parameter `t`.N@rz!(iteration %d) lower bound: %0.6fz!(iteration %d) upper bound: %0.6fz"(iteration %d) query point: %0.6f rz%(iteration %d) query was infeasible. z'(iteration %d) query was feasible. %s) z Aborting; the solver failed ... r$zMax iters hit during bisection.) r'printr(r rrrrr)solutionr*r+)rrr,r-r. tighten_lowertighten_higherepsverboser/ verbose_freqsolniquery_ptr3s r _bisectrB`sL D 9 Ad{{   3 d? "$J#% q<'1, 5C@ A 5D A B 6!XF G )wv& w 1|+q0>BC)C ^^q11 11|+q0@'**+,-##D!(+D9:##.?A A=A@   = >>rr<r=r/max_iters_interval_searchc  t|ds td|j\}} } } |rtdt |zt |} t | |t | r&|r tdttjS|| |r tdt|| ||||\}}|rtd|ztd|zt||| ||| | ||| \} }}||zdz | _ |rtd ||fz| S) aBisection on a one-parameter family of DCP problems. Bisects on a one-parameter family of DCP problems emitted by `Dqcp2Dcp`. Parameters ------ problem : Problem problem emitted by Dqcp2Dcp solver : Solver solver to use for bisection low : float lower bound for bisection (optional) high : float upper bound for bisection (optional) eps : float terminate bisection when width of interval is < eps verbose : bool whether to print verbose output related to the bisection max_iters : int the maximum number of iterations to run bisection Returns ------- A Solution object. _bisection_dataz3`bisect` only accepts problems emitted by Dqcp2Dcp.zr ******************************************************************************** Preparing to bisect problem %s zProblem is infeasible.z"Finding interval for bisection ...zinitial lower bound: %0.6fzinitial upper bound: %0.6f r7zBisection completed, with lower bound %0.6f and upper bound %0.7f ******************************************************************************** ) hasattr ValueErrorrEr8r rrrrrr4rBopt_val)rrr-r.r<r=r/rC feas_problemr,r:r; lowered_feasr?s r bisectrKs46 7- .NOO5<5L5L2L!]N 46DW6MN O",/L < <  * + -- {dl  6 7,Wad-FH T *S01 ,t34C}n Wi!OD#t$J#%DL 9d   Kr)rN)NNNd)ư>FrL)NNNrMFrLrL)__doc__r cvxpy.errorr*cvxpy.problemsrcvxpy.settingssettingsrcvxpy.problems.objectivercvxpy.reductions.solutionrr rboolrr4rBfloatintrKrr rYs!-6L%JDJFJ'*,/`03&?RZ_BE=5=RV==