Q(h]7jddlZddlZddlmZej ejZdZ dZ dZ y)N)get_arrays_tolc L|rNt|tsJt|tjr|jdk(sJt j j|sJt|tjr|j|jk(sJt|tjr|j|jk(sJt|tsJt|tsJt||}tj||ksJtj|| k\sJtj|r|dkDsJtj|d}tj|d}t||||||\}} t| | fd||||\} } t!| t!| k\r|n| } |r[tj|| ksJtj| |ksJtj"j%| d|zksJ| S)a' Maximize approximately the absolute value of a quadratic function subject to bound constraints in a trust region. This function solves approximately .. math:: \max_{s \in \mathbb{R}^n} \quad \bigg\lvert c + g^{\mathsf{T}} s + \frac{1}{2} s^{\mathsf{T}} H s \bigg\rvert \quad \text{s.t.} \quad \left\{ \begin{array}{l} l \le s \le u,\\ \lVert s \rVert \le \Delta, \end{array} \right. by maximizing the objective function along the constrained Cauchy direction. Parameters ---------- const : float Constant :math:`c` as shown above. grad : `numpy.ndarray`, shape (n,) Gradient :math:`g` as shown above. curv : callable Curvature of :math:`H` along any vector. ``curv(s) -> float`` returns :math:`s^{\mathsf{T}} H s`. xl : `numpy.ndarray`, shape (n,) Lower bounds :math:`l` as shown above. xu : `numpy.ndarray`, shape (n,) Upper bounds :math:`u` as shown above. delta : float Trust-region radius :math:`\Delta` as shown above. debug : bool Whether to make debugging tests during the execution. Returns ------- `numpy.ndarray`, shape (n,) Approximate solution :math:`s`. Notes ----- This function is described as the first alternative in Section 6.5 of [1]_. It is assumed that the origin is feasible with respect to the bound constraints and that `delta` is finite and positive. References ---------- .. [1] T. M. Ragonneau. *Model-Based Derivative-Free Optimization Methods and Software*. PhD thesis, Department of Applied Mathematics, The Hong Kong Polytechnic University, Hong Kong, China, 2022. URL: https://theses.lib.polyu.edu.hk/handle/200/12294. c| S)N)xcurvs c/var/www/skyplay_api_hub/venv/lib/python3.12/site-packages/scipy/_lib/cobyqa/subsolvers/geometry.pyz!cauchy_geometry..[s 47(皙?) isinstancefloatnpndarrayndiminspect signaturebindshapeboolrallisfiniteminimummaximum _cauchy_geomabslinalgnorm) constgradr xlxudeltadebugtolstep1q_val1step2q_val2steps ` r cauchy_geometryr. st %'''$ + Q>>  &++D111"bjj)bhh$**.DDD"bjj)bhh$**.DDD%'''%&&&R$vvbCi   vvbSDj!!!{{5!eck11 B B B B !dBE5IME6      ME6K3v;.5ED vvbDj!!!vvdbj!!!yy~~d#cEk111 Krc|rt|tsJt|tjr|jdk(sJt j |j|sJt|tjr+|jdk(r|jd|jk(sJt|tjr|j|jk(sJt|tjr|j|jk(sJt|tsJt|tsJt||}tj||ksJtj|| k\sJtj|r|dkDsJtj|d}tj|d}tj |} |} tj"j%|d} |tj& kD|j(t* |zkDz} |tj& kD|j(t*|zkz} |tj&k|j(t*|zkDz}|tj&k|j(t* |zkz}tj,tj.|| j| |j(| z }tj0|dtj& }tj.tj2||jd}tj,tj.|| j| |j(| z }tj0|dtj&}tj.tj2||jd}tj,tj.||j||j(|z }tj0|dtj& }tj.tj2||jd}tj,tj.||j||j(|z }tj0|dtj&}tj.tj2||jd}t5|jdD]S}| |t*|zkDrt1|| |z d}n't1t7||||d}t7t1||||d}||dd|fz}||dd|f}|dk\r |t* |zks|dkr|t* |zkDrt1| |z d}ntj&}|dk\r |t*|zkDs|dkr|t*|zkrt7| |z d}ntj& }t7||}t1| |}|||zzd|d zz|zz}|||zzd|d zz|zz}||kr/|||zzd|d zz|zz} t9| t9|kDr|}| }||kDr/|||zzd|d zz|zz}!t9|!t9|kDr|}|!}t9|t9|k\r float`` returns :math:`s^{\mathsf{T}} H s`. xpt : `numpy.ndarray`, shape (n, npt) Points defining the straight lines. The straight lines considered are the ones passing through the origin and the points in `xpt`. xl : `numpy.ndarray`, shape (n,) Lower bounds :math:`l` as shown above. xu : `numpy.ndarray`, shape (n,) Upper bounds :math:`u` as shown above. delta : float Trust-region radius :math:`\Delta` as shown above. debug : bool Whether to make debugging tests during the execution. Returns ------- `numpy.ndarray`, shape (n,) Approximate solution :math:`s`. Notes ----- This function is described as the second alternative in Section 6.5 of [1]_. It is assumed that the origin is feasible with respect to the bound constraints and that `delta` is finite and positive. References ---------- .. [1] T. M. Ragonneau. *Model-Based Derivative-Free Optimization Methods and Software*. PhD thesis, Department of Applied Mathematics, The Hong Kong Polytechnic University, Hong Kong, China, 2022. URL: https://theses.lib.polyu.edu.hk/handle/200/12294. rrrr)axis)r0initialN?@r)rrrrrrrrrsizerrrrrr zeros_liker r!infTTINY atleast_2d broadcast_tomax atleast_1drangeminrclip)"r"r#r xptr$r%r&r'r(r-q_vals_normi_xl_posi_xl_negi_xu_posi_xu_neg alpha_xl_pos alpha_xl_neg alpha_xu_neg alpha_xu_poskalpha_tr alpha_bd_pos alpha_bd_neg grad_step curv_stepalpha_quad_posalpha_quad_neg alpha_pos alpha_neg q_val_pos q_val_negq_val_quad_posq_val_quad_negs" r spider_geometryrYjsx %'''$ + Q>>  &++D111 sBJJ 'A  !  )  *"bjj)bhh$**.DDD"bjj)bhh$**.DDD%'''%&&&R$vvbCi   vvbSDj!!!{{5!eck11 B B B B == D E YY^^Ca^ (FbffW $!34HbffW !23HRVV r 12HRVV  23H== HNN+H5hGL66,Q@L??2==#> ! ML== HNN+H5hGL66,Q?L??2==#> ! ML== HNN+H5hGL66,Q@L??2==#> ! ML== HNN+H5hGL66,Q?L??2==#> ! ML 399Q< J !9te| #56!9,c2H 3|A Q@#F 3|A Q@#F 3q!t9$ QTO  TEI--CTEI-- )i!7=NVVN  D9,,CD9,, )i!7=N ffWN,/  <0 I ) )C)S.,@9,L L  I ) )C)S.,@9,L L  I % 9,-++i78  >"S^3* * I % 9,-++i78  >"S^3* * y>S^ +IU0K779s1a4y0"b9DE ^c)n ,Y#e*1L779s1a4y0"b9DEUJX vvbDj!!!vvdbj!!!yy~~d#cEk111 Krc|dk|dkDz}|dkD|dkz}tj|} ||| |<||| |<tjj| |kDr||z} tjj|| } tj|dz| | | | zz } | t t | zkDrt| | z d} nn`| || z| | <| | |kz}| | |kDz}tj|stj|sn||| |<||| |<| ||zz} || z}|dk\rXtjj| }|t |zkDrt||z d}nd}|| }|t |zkrt| |z d}ntj}|tj kD| t |zkz}|tjk| t |zkDz}tj||| |z tj}tj||| |z tj}t||}t|||}tj|| z||}|||zzd|dzz|zz}ntj|}|}|r[tj||ksJtj||ksJtjj|d|zksJ||fS)zM Same as `bound_constrained_cauchy_step` without the absolute value. rr3)r1r2r) rr5r r!sqrtr8rr;anyr6r>r?r)r"r#r r$r%r&r'fixed_xlfixed_xu cauchy_stepworkingg_norm delta_reducedmurOrBrLrP alpha_quadi_xli_xualpha_xlalpha_xualpha_bdalphar-rAs r rr8s! STCZ(HSTCZ(H--%KxLKxLK yy~~k"U*X%YY^^DM2FGGs ['2['5JJJMs=111/5#%W #5K +"23H+"23H66(#BFF8,<$&xLK !$&xLK !(X"5 66G).{"IC , D5L 56>3/HH% uy( (iZ)3S9JJbffW tby!89RVV dRi 7866"T([%66G66"T([%66Gx*Hj(3wwu{*B3 ))C%*,ros9"rxx\~K\Kr