o 5ήc!@sfdZddlmZddlmZddlmZmZmZddl Z ddl m Z dgZ d d Zddefd dZdS) z&Byrd-Omojokun Trust-Region SQP method.)eye) projections)modified_dogleg projected_cgbox_intersectionsN)normequality_constrained_sqpcCst|\}t|S)N)npshapespeye)xnrb/tmp/pip-target-vg8gfxp4/lib/python/scipy/optimize/_trustregion_constr/equality_constrained_sqp.pydefault_scaling s rcN CsLd}d}d}d}d}d}d}d}d}d}d}t|\}| d ur(t|tj } |d ur3t|tj}t|}| }| }|}|} |}!|}"||}#t|"| \}$}%}&|%|  }'|||'}(t| |"j|'tj})t |!d krst|!tjnd }*d d d d }+d },|| ||,|)|*|||+s"t |"|&|!|||| ||}-|(|-| }.t |!}/t |d tj |-d }0| |-}1||-}2t|(|.|$|&|/|0|1|2\}3}+|-|3}4d|(|4|4| j|4}5|"|4|!}6t|!t|6}7td|7}7|}8|5d kr|5d||7}9t||9}|5 ||7}:||t|!};||#|4}<||<\}=}>|=|t|>}?|;|?}@|@|:}A|A|krt|-|t|3kr|&|> }Bt|4|B| |\}C}D}E||#|4|D|B}F||F\}G}H|G|t|H}I|;|I}J|J|:}K|Er|K|kr|F}<|G}=|H}>|K}A|A|krt|t|4|}n<|A|krt|t|4|}n-|A|krd|d|A}L|Lt|4}M|M||kr||9}n|M||kr|M}n||9}|A|kr|<}|=|>}}!||\} }"||}#t|"| \}$}%}&|%|  }'|||'}(d },t| |"j|'tj})t |!d krt|!tjnd }*n|8}d},|| ||,|)|*|||+r|| fS)auSolve nonlinear equality-constrained problem using trust-region SQP. Solve optimization problem: minimize fun(x) subject to: constr(x) = 0 using Byrd-Omojokun Trust-Region SQP method described in [1]_. Several implementation details are based on [2]_ and [3]_, p. 549. References ---------- .. [1] Lalee, Marucha, Jorge Nocedal, and Todd Plantenga. "On the implementation of an algorithm for large-scale equality constrained optimization." SIAM Journal on Optimization 8.3 (1998): 682-706. .. [2] Byrd, Richard H., Mary E. Hribar, and Jorge Nocedal. "An interior point algorithm for large-scale nonlinear programming." SIAM Journal on Optimization 9.4 (1999): 877-900. .. [3] Nocedal, Jorge, and Stephen J. Wright. "Numerical optimization" Second Edition (2006). g333333?g?g:0yE>g@g@g?g?g?NrF)niter stop_cond hits_boundarygؗҜs