(ph?`SSKrSSKrSSKJs Jr SSKJr SSK J r J r J r SSSSSSSS.Sjr SrSSSSSSSS.S jrg) N) _RichResult)xp_signxp_copyxp_take_along_axis)argsxatolxrtolfatolfrtolmaxitercallbackc$^[XXEXgX5n U uppEpgp[R"XU4U5n U up ppmU unnU unnTRUTR [R 5TR S9nSunnTRU5nUcSUR-OUnUcSUR-OUnUc UROUnUTRTRU5TRU55-nUcA[R"UR5[R"UR5- OUn[UUUUSSSXEXgUUUS9n/SQnSnU4S jnU4S jnU4S jnU4S jn[R "UXXX?UUUUUUTS 9$)aFind the root of an elementwise function using Chandrupatla's algorithm. For each element of the output of `func`, `chandrupatla` seeks the scalar root that makes the element 0. This function allows for `a`, `b`, and the output of `func` to be of any broadcastable shapes. Parameters ---------- func : callable The function whose root is desired. The signature must be:: func(x: ndarray, *args) -> ndarray where each element of ``x`` is a finite real and ``args`` is a tuple, which may contain an arbitrary number of components of any type(s). ``func`` must be an elementwise function: each element ``func(x)[i]`` must equal ``func(x[i])`` for all indices ``i``. `_chandrupatla` seeks an array ``x`` such that ``func(x)`` is an array of zeros. a, b : array_like The lower and upper bounds of the root of the function. Must be broadcastable with one another. args : tuple, optional Additional positional arguments to be passed to `func`. xatol, xrtol, fatol, frtol : float, optional Absolute and relative tolerances on the root and function value. See Notes for details. maxiter : int, optional The maximum number of iterations of the algorithm to perform. The default is the maximum possible number of bisections within the (normal) floating point numbers of the relevant dtype. callback : callable, optional An optional user-supplied function to be called before the first iteration and after each iteration. Called as ``callback(res)``, where ``res`` is a ``_RichResult`` similar to that returned by `_chandrupatla` (but containing the current iterate's values of all variables). If `callback` raises a ``StopIteration``, the algorithm will terminate immediately and `_chandrupatla` will return a result. Returns ------- res : _RichResult An instance of `scipy._lib._util._RichResult` with the following attributes. The descriptions are written as though the values will be scalars; however, if `func` returns an array, the outputs will be arrays of the same shape. x : float The root of the function, if the algorithm terminated successfully. nfev : int The number of times the function was called to find the root. nit : int The number of iterations of Chandrupatla's algorithm performed. status : int An integer representing the exit status of the algorithm. ``0`` : The algorithm converged to the specified tolerances. ``-1`` : The algorithm encountered an invalid bracket. ``-2`` : The maximum number of iterations was reached. ``-3`` : A non-finite value was encountered. ``-4`` : Iteration was terminated by `callback`. ``1`` : The algorithm is proceeding normally (in `callback` only). success : bool ``True`` when the algorithm terminated successfully (status ``0``). fun : float The value of `func` evaluated at `x`. xl, xr : float The lower and upper ends of the bracket. fl, fr : float The function value at the lower and upper ends of the bracket. Notes ----- Implemented based on Chandrupatla's original paper [1]_. If ``xl`` and ``xr`` are the left and right ends of the bracket, ``xmin = xl if abs(func(xl)) <= abs(func(xr)) else xr``, and ``fmin0 = min(func(a), func(b))``, then the algorithm is considered to have converged when ``abs(xr - xl) < xatol + abs(xmin) * xrtol`` or ``fun(xmin) <= fatol + abs(fmin0) * frtol``. This is equivalent to the termination condition described in [1]_ with ``xrtol = 4e-10``, ``xatol = 1e-5``, and ``fatol = frtol = 0``. The default values are ``xatol = 4*tiny``, ``xrtol = 4*eps``, ``frtol = 0``, and ``fatol = tiny``, where ``eps`` and ``tiny`` are the precision and smallest normal number of the result ``dtype`` of function inputs and outputs. References ---------- .. [1] Chandrupatla, Tirupathi R. "A new hybrid quadratic/bisection algorithm for finding the zero of a nonlinear function without using derivatives". Advances in Engineering Software, 28(3), 145-149. https://doi.org/10.1016/s0965-9978(96)00051-8 See Also -------- brentq, brenth, ridder, bisect, newton Examples -------- >>> from scipy import optimize >>> def f(x, c): ... return x**3 - 2*x - c >>> c = 5 >>> res = optimize._chandrupatla._chandrupatla(f, 0, 3, args=(c,)) >>> res.x 2.0945514818937463 >>> c = [3, 4, 5] >>> res = optimize._chandrupatla._chandrupatla(f, 0, 3, args=(c,)) >>> res.x array([1.8932892 , 2. , 2.09455148]) dtype)rN?)x1f1x2f2x3f3tr r r r nitnfevstatus) rr)xxmin)funfminrrrrxlrflr)xrr)frrclURURURUR- --nU$N)rrr)workr s O/var/www/html/venv/lib/python3.13/site-packages/scipy/optimize/_chandrupatla.py pre_func_eval$_chandrupatla..pre_func_evals+ GGdff$'' 12 2c>TRURSS9TRURSS9sUlUlTR U5TR UR 5:HnU)nURUUR UsURU'URU'URUUR UsURU'URU'XsUlUlg)NT)copy)asarrayrrrrsignrr)r fr.jnjxps r/post_func_eval%_chandrupatla..post_func_evalsJJtwwTJ:JJtwwTJ:  GGAJ"''$''* *R!%TWWQZ DGGAJ#'772;   TWWR[r2c>TRUR5TRUR5:nTRXRUR 5UlTRXRUR5UlTRURTRS9nTRUR5URUR-:*n[RURU'SX!'[UR5[UR5:HU)-nTR!TR"UR R$S9nX3[R&sUR U'URU'URU'SX!'TR)UR5TR)UR 5-)nTR+UR5TR+UR5-nXE-U)-nX3[R,sUR U'URU'URU'SX!'TRUR UR- 5UlTRUR 5UR0-UR2-UlUR.UR4:n[RURU'SX!'U$)NrT)absrrwhererrr!r# zeros_likeboolr r eim _ECONVERGEDrrr5nanr _ESIGNERRisfiniteisnan _EVALUEERRdxr r tol)r.istopNaN x_nonfinitef_nanr:s r/check_termination(_chandrupatla..check_terminations( FF477ObffTWWo -HHQ1 HHQ1 }}TWWBGG}4 FF499 djj!8 8 ATWW !1 1dU :jjtyyj758s}}2 ! diilDKKN DGG,r{{477/CCD !BHHTWW$55  TE )58s~~2 ! diilDKKN&&477*+66$))$tzz1DJJ> GGdhh  A r2cN>URUR- URUR- - n[R"SSS9 UR UR - URUR - - nSSS5 URUR- URUR- - nST RSU- 5- W:UT RU5:-nUR UUR UURUX44upVpxT RUT RS55n XUU- - U-Xv- - X-Xu- - U-Xg- - - X'SUR-UR- n T RXSU - 5Ulg!,(df  GN=f)Nignore)divideinvalidr)rrrnperrstaterrrsqrt full_liker5rJrIclipr) r.xi1phi1alphar8f1jf2jf3jalphajrtlr:s r/post_termination_check-_chandrupatla..post_termination_checksdww TWWtww%67 [[( ;GGdgg%$''DGG*;USUSUSUS4up#pEUSUS:nTRXbU5US'TRXcU5US'TRXdU5US'TRXeU5US'U$Nr'r*r)r+r?resshaper'r*r)r+rKr:s r/customize_result'_chandrupatla..customize_resultsTCIs4y#d)C ID !HHQB'D HHQB'D HHQB'D HHQB'D  r2r:)_chandrupatla_ivrB _initializerZr5 _EINPROGRESSint32finfosmallest_normalepsminimumr>mathlog2maxr_loop) funcabrr r r r r rrjtempxsfsrkrrrrrrrrrsr.res_work_pairsr0r;rPrdrlr:s @r/ _chandrupatlar sh 4u  >344 ;;t  w :: 1uq % 1uq % 1uq % 1uq: ;D MM$**bii 0 0BFF4!84D4Dvvbhhtn%%t);CDD'l  !Wq[HI IHX$6$6788 uU7 DDr2dc^![XXVXxX5n U uppVpxpXU4n [R"X U5n U up pnnm!U upnUunnnT!RSUS9SnT!R UT!R[R 5T!R S9nSunnUcT!RU5ROUnUcT!RU5ROUnUcT!RU5ROUnUc/[R"T!RU5R5OUnT!RXU45T!RUUU45pT!RU SS9n[U USS9upn[UUSS9unnn[U5n[!S0SU_SU_S U_S U_S U_S U_S U_SU_SU_SU_SU_SU_SU_SU_SU_SU_6n/SQnU!4SjnSnU!4SjnSnU!4Sjn [R""UXXUUUUUUU UT!S9$)aFind the minimizer of an elementwise function. For each element of the output of `func`, `_chandrupatla_minimize` seeks the scalar minimizer that minimizes the element. This function allows for `x1`, `x2`, `x3`, and the elements of `args` to be arrays of any broadcastable shapes. Parameters ---------- func : callable The function whose minimizer is desired. The signature must be:: func(x: ndarray, *args) -> ndarray where each element of ``x`` is a finite real and ``args`` is a tuple, which may contain an arbitrary number of arrays that are broadcastable with `x`. ``func`` must be an elementwise function: each element ``func(x)[i]`` must equal ``func(x[i])`` for all indices ``i``. `_chandrupatla` seeks an array ``x`` such that ``func(x)`` is an array of minima. x1, x2, x3 : array_like The abscissae of a standard scalar minimization bracket. A bracket is valid if ``x1 < x2 < x3`` and ``func(x1) > func(x2) <= func(x3)``. Must be broadcastable with one another and `args`. args : tuple, optional Additional positional arguments to be passed to `func`. Must be arrays broadcastable with `x1`, `x2`, and `x3`. If the callable to be differentiated requires arguments that are not broadcastable with `x`, wrap that callable with `func` such that `func` accepts only `x` and broadcastable arrays. xatol, xrtol, fatol, frtol : float, optional Absolute and relative tolerances on the minimizer and function value. See Notes for details. maxiter : int, optional The maximum number of iterations of the algorithm to perform. callback : callable, optional An optional user-supplied function to be called before the first iteration and after each iteration. Called as ``callback(res)``, where ``res`` is a ``_RichResult`` similar to that returned by `_chandrupatla_minimize` (but containing the current iterate's values of all variables). If `callback` raises a ``StopIteration``, the algorithm will terminate immediately and `_chandrupatla_minimize` will return a result. Returns ------- res : _RichResult An instance of `scipy._lib._util._RichResult` with the following attributes. (The descriptions are written as though the values will be scalars; however, if `func` returns an array, the outputs will be arrays of the same shape.) success : bool ``True`` when the algorithm terminated successfully (status ``0``). status : int An integer representing the exit status of the algorithm. ``0`` : The algorithm converged to the specified tolerances. ``-1`` : The algorithm encountered an invalid bracket. ``-2`` : The maximum number of iterations was reached. ``-3`` : A non-finite value was encountered. ``-4`` : Iteration was terminated by `callback`. ``1`` : The algorithm is proceeding normally (in `callback` only). x : float The minimizer of the function, if the algorithm terminated successfully. fun : float The value of `func` evaluated at `x`. nfev : int The number of points at which `func` was evaluated. nit : int The number of iterations of the algorithm that were performed. xl, xm, xr : float The final three-point bracket. fl, fm, fr : float The function value at the bracket points. Notes ----- Implemented based on Chandrupatla's original paper [1]_. If ``x1 < x2 < x3`` are the points of the bracket and ``f1 > f2 <= f3`` are the values of ``func`` at those points, then the algorithm is considered to have converged when ``x3 - x1 <= abs(x2)*xrtol + xatol`` or ``(f1 - 2*f2 + f3)/2 <= abs(f2)*frtol + fatol``. Note that first of these differs from the termination conditions described in [1]_. The default values of `xrtol` is the square root of the precision of the appropriate dtype, and ``xatol = fatol = frtol`` is the smallest normal number of the appropriate dtype. References ---------- .. [1] Chandrupatla, Tirupathi R. (1998). "An efficient quadratic fit-sectioning algorithm for minimization without derivatives". Computer Methods in Applied Mechanics and Engineering, 152 (1-2), 211-217. https://doi.org/10.1016/S0045-7825(97)00190-4 See Also -------- golden, brent, bounded Examples -------- >>> from scipy.optimize._chandrupatla import _chandrupatla_minimize >>> def f(x, args=1): ... return (x - args)**2 >>> res = _chandrupatla_minimize(f, -5, 0, 5) >>> res.x 1.0 >>> c = [1, 1.5, 2] >>> res = _chandrupatla_minimize(f, -5, 0, 5, args=(c,)) >>> res.x array([1. , 1.5, 2. ]) gw?rr)rr)axisrrrrrrphir r r r rrrq0r) r)r r)r"rr$r%r&)xmr)r*rr()fmr)r+rc>URUR- nURUR- nXRUR- -nX R UR- -nX3U-- nSXPRUR- -UR-UR--nT R X`R- 5ST R U5-:nXgnT R XgURU- 5URU:*n URUU [X'U 5URUU --X'URSUR- U--n XU'X`lU $)Nrr) rrrrrrr>rxtolrr) r.x21x32ABCq1rKxir8r r:s r/r0-_chandrupatla_minimize..pre_func_evalsKgggg 77TWW$ % 77TWW$ % QK Aww()DGG3dgg= > FF2< 3#4 4 U FF254771:% &$))A, 6 1 q 2TYYq\!_ DD GGq488|s* *!r2cP[XR- 5[URUR- 5:HnXURUURUURU4upEpgXURUUR UUR U4uppX:n XLXsX|'X'U )n XlXXLX4uX\'X'Xl'X'U)n X URU URU URU 4upnnXURU UR U UR U 4unnnnUU:n XUU sX'UU 'U )n UU UU XUU 4uUU 'UU 'UU 'UU 'XVUsURU'URU'URU'XU sURU'UR U'UR U'UUUsURU 'URU 'URU 'UUUsURU 'UR U 'UR U 'gr-)rrrrrrr)r r7r.rKrx1ix2ix3ifif1if2if3ir8nixnix1nix2nix3nifnif1nif2nif3nis r/r;._chandrupatla_minimize..post_func_evals- AK GDGGdgg,=$> >D$''!*dggaj$''!*ED$''!*dggaj$''!*D H B),)E&R !twwr{DGGBK M4 !twwr{DGGBK LT4 $J63q6a B-1!Wd1gsvs1v-M*Qa$q'47-0s* DGGAJ -0s* DGGAJ 04dD- TWWR[$''"+04dD- TWWR[$''"+r2cD>TRUR[S9nURUR:URUR :-nTR TR sURU'URU'S[RsX'URU'TRURUR-UR-UR-UR-UR -5nX1-)nTR TR sURU'URU'S[RsX'URU'TRURUR- 5TRURUR- 5:nURUnURUURU'X@RU'URUnUR UURU'X@R U'TRUR5UR-UR -UlTRURUR- 5SUR"-:*nTRUR5UR$-UR&-nX RSUR-- UR -SU-:*-nX!)-nS[R(sX'URU'U$)NrTr)r@rrArrrrDrrBrErrFrrHr>r r rr r rC)r.rLrKfiniter~ftolr:s r/rP1_chandrupatla_minimize..check_terminationsb}}TWWD}1ggDGGdgg$5 6!# DGGAJ"& QTWWTWW_TWW4TWW.post_termination_checks r2c>USUSUSUS4up#pEUSUS:nTRXcU5US'TRXbU5US'TRXeU5US'TRXdU5US'U$rgrhris r/rl0_chandrupatla_minimize..customize_resultsTCIs4y#d)C IT "HHQB'D HHQB'D HHQB'D HHQB'D  r2rn)rorBrpr5rZrqrrrsrtrwrYrustackargsortrrrrz)"r{rrrrr r r r r rrjrr~rrkrrrrrrrrrKrr.rr0r;rPrdrlr:s" @r/_chandrupatla_minimizersoj 4u  rsR 66(FF')DAtdgTE<68t!%Ts$(Tr2