o i"*@sdZddlZddlZddlZddlZddlmZm Z m Z m Z m Z ddl mZddlmZgZddZGd d d Z    dddZdS)zTrust-region optimization.N)_check_unknown_options_status_messageOptimizeResult_prepare_scalar_function_call_callback_maybe_halt)HessianUpdateStrategy) FD_METHODScs.dgdur dfSfdd}|fS)Nrcs*dd7<t|g|RS)Nrr)npcopy)xZ wrapper_argsargsfunctionZncallsj/var/www/html/eduruby.in/lip-sync/lip-sync-env/lib/python3.10/site-packages/scipy/optimize/_trustregion.pyfunction_wrappersz(_wrap_function..function_wrapperr)rrrrr r_wrap_functions rc@sjeZdZdZdddZddZeddZed d Zed d Z d dZ eddZ ddZ ddZ dS)BaseQuadraticSubproblemaQ Base/abstract class defining the quadratic model for trust-region minimization. Child classes must implement the ``solve`` method. Values of the objective function, Jacobian and Hessian (if provided) at the current iterate ``x`` are evaluated on demand and then stored as attributes ``fun``, ``jac``, ``hess``. NcCsF||_d|_d|_d|_d|_d|_d|_||_||_||_ ||_ dSN) _x_f_g_h_g_magZ _cauchy_pointZ _newton_point_fun_jac_hess_hessp)selfr funjachesshessprrr__init__(s z BaseQuadraticSubproblem.__init__cCs*|jt|j|dt|||S)Ng?)r r dotr!r#rprrr__call__5s*z BaseQuadraticSubproblem.__call__cC|jdur ||j|_|jS)z1Value of objective function at current iteration.N)rrrrrrrr 8 zBaseQuadraticSubproblem.funcCr))z=Value of Jacobian of objective function at current iteration.N)rrrr*rrrr!?r+zBaseQuadraticSubproblem.jaccCr))zr@rrrrrs       rr?@@333333?-C6?FTc# slt||dur td|dur|durtd|dur tdd| kr-dks2tdtd|dkr:td|dkrBtd ||krJtd t|}t|||||d j}j}t |rhj }nt |rmn|t vsvt |t rd}fd d }ntdt||\}}| durt|d} d}|}|}| r|g}||||||}d}|j| kr@z ||\}}Wn tjjyd}Yn{w||}||}||||||}|j|j}|j|}|dkrd}nX||}|dkr|d9}n|dkr|rtd||}|| kr |}|}| r|t||d7}t||jd} t|| r*n|j| kr3d}n || kr;d}n|j| kstdtdddf}!| r|dkrYt|!|n t|!|tdtd|jtd|tdjtdjtdj |dt||dk||j|j!jjj |d||!|d }"|dur|j |"d<| r||"d<|"S) a Minimization of scalar function of one or more variables using a trust-region algorithm. Options for the trust-region algorithm are: initial_trust_radius : float Initial trust radius. max_trust_radius : float Never propose steps that are longer than this value. eta : float Trust region related acceptance stringency for proposed steps. gtol : float Gradient norm must be less than `gtol` before successful termination. maxiter : int Maximum number of iterations to perform. disp : bool If True, print convergence message. inexact : bool Accuracy to solve subproblems. If True requires less nonlinear iterations, but more vector products. Only effective for method trust-krylov. This function is called by the `minimize` function. It is not supposed to be called directly. Nz7Jacobian is currently required for trust-region methodsz_Either the Hessian or the Hessian-vector product is currently required for trust-region methodszBA subproblem solving strategy is required for trust-region methodsrg?zinvalid acceptance stringencyz%the max trust radius must be positivez)the initial trust radius must be positivez?the initial trust radius must be less than the max trust radius)r!r"rcs||Sr)r"r%)r r'rZsfrrr#sz%_minimize_trust_region..hesspr/g?r)r r successmaxiterz:A bad approximation caused failure to predict improvement.z3A linalg error occurred, such as a non-psd Hessian.z# Current function value: %fz Iterations: %dz! Function evaluations: %dz! Gradient evaluations: %dz Hessian evaluations: %d) r rMstatusr r!nfevZnjevnhevnitmessager"allvecs)"r ValueError Exceptionr Zasarrayflattenrr Zgradcallabler"r isinstancerrlenr.r@r-Z LinAlgErrorminappendr rrrprintwarningswarnRuntimeWarningrPZngevrQr!)#r Zx0rr!r"r#Z subproblemZinitial_trust_radiusZmax_trust_radiusetaZgtolrNZdispZ return_allcallbackZinexactZunknown_optionsZnhesspZwarnflagr8r rTmkr'Z hits_boundaryZpredicted_valueZ x_proposedZ m_proposedZactual_reductionZpredicted_reductionrhoZintermediate_resultZstatus_messagesresultrrJr_minimize_trust_regionvs          ;     rg)rNNNNrFrGrHrINFFNT)rDr2r^numpyr Z scipy.linalgr, _optimizerrrrrZ'scipy.optimize._hessian_update_strategyrZ(scipy.optimize._differentiable_functionsr __all__rrrgrrrrs"  X