0^j<dZddlZddlmZddlmZmZmZm Z ddl m Z m Z m Z mZmZddlmZddlmZdd lmZmZmZmZmZmZdd lmZmZmZdd lm Z dd l!m"Z"dd l#m$Z$m%Z%ddl&m'Z'm(Z(m)Z)ddl*m+Z+m,Z,ddl-m.Z.m/Z/ddl0m1Z1m2Z2dZ3dZ4dZ5y)z This module performs the major calculations of COBYLA. Translated from Zaikun Zhang's modern-Fortran reference implementation in PRIMA. Dedicated to late Professor M. J. D. Powell FRS (1936--2015). Python translation by Nickolai Belakovski. N)checkbreak_con)REALMAXEPS DEBUGGING MIN_MAXFILT) INFO_DEFAULT MAXTR_REACHEDDAMAGING_ROUNDINGSMALL_TR_RADIUSCALLBACK_TERMINATE)evaluate)savehist)isinvmatprodinprodnormprimasum primapow2)fmsgretmsgrhomsg)redrat)redrho)savefiltselectx) updatepolefindpole updatexfc) setdrop_trgeostep)trstlptrrad)initxfcinitfiltcrg}g}g}g}d}tjtj|tj|f}tjtj|dz}t}|tj|nd}tj|}||z }tj|}trt |dksJ||k\r|dk\sJ|dk\sJ||dzk\sJ| |k\r|dkDsJt tj|sJd|cxkr | cxkrdksJJd| cxkr dcxkr| ksJJd|ksJd|ksJd|cxkr|ksJJ|tj|||fk(sJtt||cxkr|ksJJt||||||||| | ||||||\} }!}"}#}$}%}&}'tjtjtj|d|}(tjtj|tj|(f})tjtj|(}*tjtj|tj|(f}+t|!|||"|%|#| |(|)|*|+ },|'tk7r|'}-t!|*d|,|(d|,||}.|+dd|.f}|*|.}|)dd|.f}|(|.}/t#||-||&|||/|trk|&|ksJtj||k(rt%tj&|rJtj&|stj(|rJ||||/|&|||||-f S| }0| }1tj|tjdt+|!|%}2d}3d }4d}5tjdtjd | zd }6d |z}7t,}-t/|7D]}8t1|||!|2|"|1|%|0|#|$ }2t3|2|!|"|%|#|$\}!}"}%}#}$}'|'t4k(r|'}-ndt t7t9|#ddd|fd dt9|1zk}9t;|%d||%|z |$}:| |j<n||ddd|f<t;|!|dd|ftj>|!|d|f|dfj<z |$j<|dd|df<tA||!dd|f |1|:};t|1tC|;}<||=|2|>zz}?|?dt|2dz|0zkD }@|3s@r|1dz}1|1|6|0zkr|0}1n|#dd|f|;z}t7t9||#dd|fz ||<t7t9|jK|d|#dd|fjK|d|#ddd|fzz d |d|tjL|}A||At9d|zkr|%A}|!dd|Af}|"|A}/nntO|||||\}}tjFtjHd|}/|&dz }&tQ||||||/||| tS|/|||||,|(|*|+||) \},}(}*}+})tU|d||&|1|||/| |%||2|"|zz||2|/zzz }BtW|B|?|}4tY|1|<|| | | |4}1|1|6|0zkr|0}1BdkD}Ct[|C|;|1|0|#|$}5t]|5||2|/|;||!|"|%|#|$ \}#}$}%}!}"}'|'t4k(r|'}-n>t_||&|/||| |}'|'tk7r|'}-n|3xs@xs |4dkxs|5du}D|Dxr|9 }EDxr|9xrtG|1|<|0k}FErt t7t9|#ddd|fd dt9|1zkstj`t7t9|#ddd|fd d }G|1dz }Htc|G|||!|2|"|H|%|$ };|#dd|f|;z}t7t9||#dd|fz ||<t7t9|jK|d|#dd|fjK|d|#ddd|fzz d |d|tjL|}A||At9d|zkr|%A}|!dd|Af}|"|A}/nntO|||||\}}tjFtjHd|}/|&dz }&tQ||||||/||| tS|/|||||,|(|*|+||) \},}(}*}+})tU|d||&|1|||/| t]G||2|/|;||!|"|%|#|$ \}#}$}%}!}"}'|'t4k(r|'}-nt_||&|/||| |}'|'tk7r|'}-nFr|0|krtd}-ntGd|0ztg|0|}1tg|0|}0tj|tj|2t+|!|%}2ti|||&|%||0|#dd|f|"||!dd|f|2 t3|2|!|"|%|#|$\}!}"}%}#}$}'|'t4k(r|'}-n2|ss||#dd|f|%||&|8|"||!dd|f}I|Istj}-n|#dd|f;z}|-tdk(r|3rtC||#dd|fz d|zkDr|&|krtO|||||\}}tjFtjHd|}/|&dz }&tQ||||||/||| tS|/|||||,|(|*|+||) \},}(}*}+})tU|d||&|0|||/| t!|*d|,|(d|,tG|2||}.|+dd|.f}|*|.}|)dd|.f}|(|.}/t#||-||&|||/|||||/|&|||||-f S)zE This subroutine performs the actual computations of COBYLA. COBYLArNrrg@@Fg?g? axisg?gư>g-C6?z Trust regionGeometry?gMbP?)6npzerossizerrabsallisfiniteshapeminrr%minimummaximumr&r rranyisnanisposinffcratior rangegetcpenrr rrrTtiler#rrmaxappendreshapeargminrrrrrr$r!r rargmaxr"r rrr )Jcalcfciprintmaxfiltmaxfunamatbvecctolcweighteta1eta2ftargetgamma1gamma2rhobegrhoendconstrfxmaxhistcallbackxhistfhistchistconhistsolverAdistsqcpenminm_lconnum_constraintsm_nlconnum_vars evaluatedconmatcvalsimsimifvalnfsubinfocfiltconfiltffiltxfiltnfiltinfokoptcstrvrhodeltacpenshortdratiojdrop_trgamma3maxtrtr adequate_geogddnormprerefprerecpreremtrfailjactrem ximproved bad_trstep improve_geo reduce_rho jdrop_geodelbar terminatesJ p/var/www/ramen.bs-engineer-server.com/venv/lib/python3.12/site-packages/scipy/_external/pyprima/cobyla/cobylb.pycobylbrs  E E EGF "''!*bggfo./A XXbggaj1n %FG#.RWWT]AFggfoO&GwwqzH6{a&(Vq[881}}A%%%FQJ..2;;q>"""D$D$1$$$$$6&A&&&&&&DyyG||G%v%%%%%|rxx~&(1CCCC;'7udGUV?@%PU?Dfg7W3ueUG UAq% P8ntd8n'<<TE\AQRF6640E<%dFFEJEs " I ")QsCFH6?xQUW\^_abdjlprvx{~B6C 2CtVT7++%VRa!LG,&GvG!Gx47G  5%5 MLMS5F#5M l s8Ic!YhY,6G,Hq#QUVYbchYiUi#ij: (9SIXI5F+Ga"PWXYI1WF 4vtT64QUVAAxK 1$A' !c!X+6F2F(GHF8  (199Xq3IQ[!))(A6Q  \9JJ4L*MST!VF9H  &!AayIdVm44G1Q%VQtD 6ryyF34a!UAueUFGT6>udGUV?@%PU?Dfg7W3ueUG VR1eV L5>y&RVX]_`bcekmqswy|C6D 2CtVT7++%VRa!LG,& f}&c 6#v#67Ef%C::grzz$8M'NOD 662tH~sC8 2r4PX>[abcembm[noI)y \B AxK1A   SH%% && 8 KVQt< 6ryyF+, a!UAueUFGL.6udGQPQSXZ_afhmouw~.+ueUG V^VRaE6J 5%=%-T71CT JD agA d A QW F $KE 64Q5&9 aE5%$ FFc  |j}|j}|j}|j}| j} tjtj|dtj|df} d} |tj|nd} tj|d} tj|d}tr?| dk\sJ|dk\sJ|dkDsJtj|d| k(rtj|d|dzk(sJtj |tj |zjrJtj||dzk(r:t|dtj |ztj|zkrJtj||dzk(r4ttj |tj|zrJtj|d|k(rtj|d|dzk(sJtj|jsJttjt|ddd|fddkDsJtj| d|k(rtj| d|k(sJtj| jsJt|ddd|f| | sJ||k\r|dkDsJt}t|dzD]J}t!|||||| \}}}}} }|t"k(rn't%|d|||z | }| |j&n|| ddd| f<t%|| dd|ftj(|| d|f|dfj&z | j&| dd| df<t+| |dd|f ||}t-|| }||tjtj.d|dd|ft%|| zz }|dkDr|dksn5t|t1d|z|z t2}t5||||k(sKntrU||k\r|dkDsJ|zzdkDs>|t"k(s5|dk\r0tj6||dkDrtj|rJ|S)z This function gets the penalty parameter CPEN so that PREREM = PREREF + CPEN * PREREC > 0. See the discussions around equation (9) of the COBYLA paper. rrNr,)copyr1r2r3rr<isneginfr;r=r6r5rCr4rr r?rr rrArBr#rrDr8rrr:)rLrMrirzrjryrmrxrkrlraitolrdrergruiterrrrrs rr@r@2s[[]F 99;D 99;D ((*C 99;D "''#q/27761#567A D#.RWWT]AFggfa(OwwsAH!###1}}axxwwvq!_49KxZ[|9[[[HHV$r{{6'::??AAAwwt}1 ,D1rxx~- D0AAAB CCwwt}1 ,S$"++VZJ[9[5\\\wwsA(*rwwsA(Q,/NNN{{3##%%%266#c!YhY,/0q9A=>>>wwtQ8+a0@H0LLL{{4 $$&&&SIXI&d333|a'' D"hl##.8vtTSV9=/?+dCt $ $  D(Od8n4d ;"&"2!WfW* ((:!;''&():";h]KMM"NPTVVWVW !VW*  1vak**E1 5A,h"&&1fQ[6IGTUWXM6Y)Z"[[ vz 4R&[617;< D$ % 1 G#Tt|q((v %)T5F-F1FF!;a!? TZH[ \\ KrctrHtj|dk\sJtj|dtj|k(sJtj| d}tj| d}t|}t |}t |d|zkrQ||krLtjtj |d|z |d|zktj}||z |z }nd}tr|dk\sJ|S)z This function calculates the ratio between the "typical change" of F and that of CONSTR. See equations (12)-(13) in Section 3 of the COBYLA paper for the definition of the ratio. rr,r0r)whereinitial)rr1r3r8rCr;r:inf)rirmcmincmaxfminfmaxdenomrs rr>r>swwt}!!!wwvq!RWWT]222 666' "D 666' "D t9D t9D 4#* $+rzz$*T1d 9JTVTZTZ[D[E ! Av v Hr)6__doc__numpyr1common.checkbreakr common.constsrrrr common.infosr r r r r common.evaluatercommon.historyr common.linalgrrrrrrcommon.messagerrr common.ratior common.redrhorcommon.selectxrrupdaterrr geometryr!r" trustregionr#r$ initializer%r&rr@r>rrrs`.@@88&%MM11!".33)&)QGjqh$ r