ftj<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/ddl0m1Z1m2Z2dZ3dZ4dZ5dS)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)initxfcinitfiltcg}g}g}g}d}tjtj|tj|f}tjtj|dz}t}|tj|nd}tj|}||z }tj|}trt |dksJ||kr|dksJ|dksJ||dzksJ| |kr|dksJt tj|sJd|cxkr | cxkrdksnJd| cxkr dcxkr| ksnJd|ksJd|ksJd|cxkr|ksnJ|tj|||fksJtt||cxkr|ksnJt||||||||| | ||||||\} }!}"}#}$}%}&}'tjtj tj |d|}(tjtj|tj|(f})tjtj|(}*tjtj|tj|(f}+t|!|||"|%|#| |(|)|*|+ },|'tkr|'}-t!|*d|,|(d|,||}.|+dd|.f}|*|.}|)dd|.f}|(|.}/t#||-||&|||/|trm|&|ksJtj||kr!t%tj|rJtj|stj|rJ||||/|&|||||-f S| }0| }1tj |tj dt+|!|%}2d}3d }4d}5tj dtj d | zd }6d |z}7t,}-t/|7D]}8t1|||!|2|"|1|%|0|#|$ }2t3|2|!|"|%|#|$\}!}"}%}#}$}'|'t4kr|'}-nt t7t9|#ddd|fd dt9|1zk}9t;|%d||%|z |$}:||jn||ddd|f<t;|!|dd|ftj|!|d|f|dfjz |$j|dd|df<tA||!dd|f |1|:};t|1tC|;}<||=|2|>zz}?|?dt|2dz|0zk }@|3s|@r|1dz}1|1|6|0zkr|0}1n;|#dd|f|;z}t7t9||#dd|fz ||<t7t9|%|d|#dd|f%|d|#ddd|fzz d |d|<tj&|}A||At9d|zkr|%|A}|!dd|Af}|"|A}/nxtO|||||\}}tj#tj$d|}/|&dz }&tQ||||||/||| tS|/|||||,|(|*|+||) \},}(}*}+})tU|d||&|1|||/| |%||2|"|zz||2|/zzz }BtW|B|?|}4tY|1|<|| | | |4}1|1|6|0zkr|0}1|Bdk}Ct[|C|;|1|0|#|$}5t]|5||2|/|;||!|"|%|#|$ \}#}$}%}!}"}'|'t4kr|'}-nt_||&|/||| |}'|'tkr|'}-n|3p |@p |4dkp|5du}D|Do|9 }E|Do|9otG|1|<|0k}F|Eryt t7t9|#ddd|fd dt9|1zks-tj0t7t9|#ddd|fd d }G|1dz }Htc|G|||!|2|"|H|%|$ };|#dd|f|;z}t7t9||#dd|fz ||<t7t9|%|d|#dd|f%|d|#ddd|fzz d |d|<tj&|}A||At9d|zkr|%|A}|!dd|Af}|"|A}/nxtO|||||\}}tj#tj$d|}/|&dz }&tQ||||||/||| tS|/|||||,|(|*|+||) \},}(}*}+})tU|d||&|1|||/| t]|G||2|/|;||!|"|%|#|$ \}#}$}%}!}"}'|'t4kr|'}-n>t_||&|/||| |}'|'tkr|'}-n|Fr|0|kr td}-ntGd|0ztg|0|}1tg|0|}0tj |tj |2t+|!|%}2ti|||&|%||0|#dd|f|"||!dd|f|2 t3|2|!|"|%|#|$\}!}"}%}#}$}'|'t4kr|'}-n?|r;||#dd|f|%||&|8|"||!dd|f}I|Ir tj}-n|#dd|f|;z}|-tdkr|3rtC||#dd|fz d|zkr|&|krtO|||||\}}tj#tj$d|}/|&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 ^/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/_lib/pyprima/cobyla/cobylb.pycobylbrs0 E E EGF "'!**bgfoo.//A Xbgajj1n % %FG#.RWT]]]AFgfooO&GwqzzH =6{{a&((Vq[[[81}}}}A%%%%FQJJJ.2;q>>"""""D$$$$D$$$$1$$$$$$6&&&&A&&&&&&&&&&DyyyyG||||G%%%%v%%%%%%|rx~~&(1CCCCC;''7<<<udGUV?@%PU?Dfg7W7W3ueUG UAq% P P P8ntd8n'<<TE\AQRF66400E<%dFFEJJEs "" I ")QsCFFH6?xQUW\^_abdjlprvx{~B6C6C 2CtVT7+++%VRa!LLG,&&GvG!Gx47G  5%5 MLMS5F5F#5M l N s8Ic!!!YhY,6G,H,Hq#Q#Q#QUVYbchYiYiUi#ijjN : (9SIXI5F+G+Ga"P"P"PWXYYYI1WF 4vtT64QUVVAAAAxK 1$A' !c!!!X+6F2F(G(GHHF8  (199Xq3I3IQQQ[!))(A66QQQ  \9JJ4L*M*MST!V!V!VF9H9  &!!AayIdVm4444G1Q%VQtDD 6ryF3344a!UAueUFGTTT6>udGUV?@%PU?Dfg7W7W3ueUG VR1eV L L L5>y&RVX]_`bcekmqswy|C6D6D 2CtVT7+++%VRa!LLG,&&  f}}&c 6#v#6#677Ef%%C:grz$8M8M'N'NOOD 662tH~sC8 2r4PX>[abcbcbcembm[nooI ) AAAxK1A   SH%% & && 8 8 KKVQt<< 6ryF++,, a!UAueUFGLLL.6udGQPQSXZ_afhmouw~..+ueUG V^VRaE6JJJ 5%=%-T71C1CT J JD aaagA d A QQQW F $KE 64Q5&999 aE5%$ FFc  |}|}|}|}| } tjtj|dtj|df} d} |tj|nd} tj|d} tj|d}trc| dksJ|dksJ|dksJtj|d| krtj|d|dzksJtj|tj|zrJtj||dzkr=t|dtj|ztj|zkrJtj||dzkr6ttj|tj|zrJtj|d|krtj|d|dzksJtj | sJttj t|ddd|fddksJtj| d|krtj| d|ksJtj |  sJt|ddd|f| | sJ||kr|dksJt}t|dzD]y}t!|||||| \}}}}} }|t"krnOt%|d|||z | }||jn|| ddd| f<t%|| dd|ftj|| d|f|dfjz | j| dd| df<t+| |dd|f ||}t-|| }||tj tjd|dd|ft%|| zz }|dkr|dksnBt|t1d|z|z t2}t5||||krn{trZ||kr|dksJ|||zzdks@|t"ks5|dkr/tj||dkrtj |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@2sf[[]]F 99;;D 99;;D ((**C 99;;D "'#q//2761#5#5677A D#.RWT]]]AFgfa((OwsAH(!####1}}}}axxxxwvq!!_449K9KxZ[|9[9[9[[HV$$r{6':'::??AAAAAwt}}1 ,,D1rx~~- D0A0AAABB-,Cwt}}1 ,,S$"+VZJ[J[9[5\5\,,\wsA(**rwsA(Q,/N/N/NN{3##%%%%%26#c!!!YhY,/00q999A=>>>>>wtQ8++a0@0@H0L0L0LL{4  $$&&&&&SIXI&d33333||a' D"hl####.8vtTSV9=/?/?+dCt $ $ $ E D(Od8n4d ; ;"&"2!!!WfW* ((:!;'&():";h]KKM"NPTVVVW !!!VWW*  1vaaak**E1 5 5A,,h"&1fQQQ[6IGTUWXMM6Y)Z)Z"["[[ vzz E4R&[617;;<< D$ % % 1 1 E 2\t||q(v %))T5F-F-F1FF!;!;a!?!? TZH[H[!?!?\ Krc"trGtj|dksJtj|dtj|ksJtj| d}tj| d}t|}t |}t |d|zkrM||krGtjtj|d|z |d|zktj}||z |z }nd}tr|dksJ|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>s3wt}}!!!!wvq!!RWT]]2222 66' " " "D 66' " " "D t99D t99D 4#* $++rz$**T1d 9JTVTZ[[[D[E ! Avvvv 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......@@@@@@@@@@@@88888888888888&&&&&&%%%%%%MMMMMMMMMMMMMMMM1111111111!!!!!!""""""........3333333333))))))))&&&&&&&&))))))))QGQGQGjqqqh$ $ $ $ $ r