ghvdZddlZddlmZmZmZmZmZddlm Z m Z m Z m Z m Z mZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9ddl:m;Z;mZ>m?Z?m@Z@mAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZHmIZImJZJmKZKmLZLmMZMmNZNmOZOmPZPmQZQmRZRe efZSeefZTe!efZUe"efZVe#e$fZWe#e$e%fZXdZYdZZd Z[d efd Z\d Z]efd Z^dZ_efdZ`efdZadefdZbdefdZcefdZddefdZedZfefdZgefdZhefdZiefdZjefdZkefdZlefdZmefdZnefdZoefdZpefd Zqefd!Zrefd"Zsefd#Ztefd$Zud%Zvefd&Zwefd'Zxefd(Zyefd)Zzd*Z{efd+Z|efd,Z}efd-Z~efd.Zefd/Zefd0Zefd1Zefd2Zefd3Zefd4Zefd5Zefd6Zefd7Zefd8Zefd9Zed:Zed;Zd<Zefd=Zefd>Zefd?Zefd@ZefdAZefdBZdJdDZdJdEZdJdFZdJdGZedHkr6 ddlmcmcmZej~Z~ejzZzdS#eef$redIYdSwxYwdS)Kz- Low-level functions for complex arithmetic. N)MPZMPZ_ZEROMPZ_ONEMPZ_TWOBACKEND)1 round_floor round_ceiling round_downround_up round_nearest round_fastbitcountbctable normalize normalize1reciprocal_rndrshiftlshift giant_steps negative_rndto_strto_fixed from_man_exp from_floatto_floatfrom_intto_intfzerofoneftwofhalffinffninffnanfnonempf_absmpf_posmpf_negmpf_addmpf_submpf_mulmpf_div mpf_mul_int mpf_shiftmpf_sqrt mpf_hypot mpf_rdiv_int mpf_floormpf_ceilmpf_nintmpf_fracmpf_signmpf_hash ComplexResult)mpf_pimpf_expmpf_log mpf_cos_sin mpf_cosh_sinhmpf_tan mpf_pow_int mpf_log_hypotmpf_cos_sin_pimpf_phimpf_cosmpf_sin mpf_cos_pi mpf_sin_pimpf_atan mpf_atan2mpf_coshmpf_sinhmpf_tanhmpf_asinmpf_acos mpf_acosh mpf_nthroot mpf_fibonaccic<|\}}|tvrdS|tvrdSdS)z2Check if either real or imaginary part is infiniteTF)_infszreims O/var/www/html/ai-engine/env/lib/python3.11/site-packages/mpmath/libmp/libmpc.py mpc_is_infrY)s* FB U{{44 U{{44 5c<|\}}|tvrdS|tvrdSdS)z9Check if either real or imaginary part is infinite or nanTF) _infs_nanrTs rX mpc_is_infnanr]0s* FB Ytt Ytt 5rZc |\}}t||}|dr$|dztt||fi|zdzS|dzt||fi|zdzS)Nrz - jz + )rr))rUdpskwargsrVrWrss rX mpc_to_strrc7sv FB CB !u<EzF72;;>>v>>>DDEzF2s55f555;;rZFch|\}}tt|||t|||SN)complexr)rUstrictrndrVrWs rXmpc_to_complexri?s3 FB 8B,,hr63.G.G H HHrZcVtjdkr\|\}}t|tjjt|zz}|dtjjzz}t |S tt|dS#t$rt|cYSwxYw)N)rlT)rg) sys version_infor8 hash_infoimagwidthinthashri OverflowError)rUrVrWhs rXmpc_hashrvCs 6!!B RLL3=- < < CM'' (1vv  q66677 7   77NNN s.B B('B(c2|\}}|t|||fSrer)rUprecrhrVrWs rX mpc_conjugater{Ps" FB wr4%% %%rZc|tkSre)mpc_zero)rUs rXmpc_is_nonzeror~Ts =rZc^|\}}|\}}t||||t||||fSrer*rUwrzrhabcds rXmpc_addrW; DAq DAq 1as # #WQ4%=%= ==rZc4|\}}t|||||fSrer)rUxrzrhrrs rX mpc_add_mpfr\$ DAq 1as # #Q &&rZc^|\}}|\}}t||||t||||fSrer+rs rXmpc_subr`rrZc4|\}}t|||||fSrer)rUprzrhrrs rX mpc_sub_mpfrerrZcP|\}}t|||t|||fSre)r(rUrzrhrrs rXmpc_posri. DAq 1dC '!T3"7"7 77rZcP|\}}t|||t|||fSrerxrs rXmpc_negrmrrZcL|\}}t||t||fSre)r/)rUnrrs rX mpc_shiftrqs& DAq Q??IaOO ++rZc0|\}}t||||S)zEAbsolute value of a complex number, |a+bi|. Returns an mpf value.)r1rs rXmpc_absrus! DAq Q4 % %%rZc0|\}}t||||S)z3Argument of a complex number. Returns an mpf value.)rIrs rXmpc_argr{s DAq Q4 % %%rZcP|\}}t|||t|||fSre)r3rs rX mpc_floorrs. DAq Qc " "Ias$;$; ;;rZcP|\}}t|||t|||fSre)r4rs rXmpc_ceilr. DAq AtS ! !8AtS#9#9 99rZcP|\}}t|||t|||fSre)r5rs rXmpc_nintrrrZcP|\}}t|||t|||fSre)r6rs rXmpc_fracrrrZc|\}}|\}}t||}t||} t||} t||} t|| ||} t| | ||} | | fS)z Complex multiplication. Returns the real and imaginary part of (a+bi)*(c+di), rounded to the specified precision. The rounding mode applies to the real and imaginary parts separately. )r,r+r*)rUrrzrhrrrrrqrsrVrWs rXmpc_mulrsy DAq DAq1 A1 A1 A1 A AtS ! !B AtS ! !B r6MrZc|\}}t||}t||}t||||}t||||}t|d} || fSNr)r,r+r/) rUrzrhrrrrrrVrWs rX mpc_squarersa DAq! A! A!T3A AtS ! !B 1aB r6MrZc\|\}}t||||}t||||}||fSre)r,rUrrzrhrrrVrWs rX mpc_mul_mpfrs: DAq AtS ! !B AtS ! !B r6MrZcv|\}}tt||||}t||||}||fS)zB Multiply the mpc value z by I*x where x is an mpf value. )r)r,)rUrrzrhrrrVrWs rXmpc_mul_imag_mpfrsD DAq AtS)) * *B AtS ! !B r6MrZc\|\}}t||||}t||||}||fSre)r.)rUrrzrhrrrVrWs rX mpc_mul_intrs: DAq Q4 % %B Q4 % %B r6MrZcv|\}}|\}}|dz}tt||t|||} tt||t|||} tt||t|||} t| | ||t| | ||fSN )r*r,r+r-) rUrrzrhrrrrwpmagtus rXmpc_divrs DAq DAq B '!Q--A 3 3C! gallB//A! gallB//A 1Sc " "GAc$s$;$; ;;rZc\|\}}t||||}t||||}||fS)zCalculate z/p where p is real)r-rs rX mpc_div_mpfrs: DAq AtS ! !B AtS ! !B r6MrZc|\}}tt||t|||dz}t||||}tt||||}||fS)zCalculate 1/z efficientlyrr*r,r-r))rUrzrhrrmrVrWs rXmpc_reciprocalrsg DAq! WQq\\$r'22A AtS ! !B AtS)) * *B r6MrZc|\}}tt||t|||dz}tt|||||}ttt|||||}||fS)z)Calculate p/z where p is real efficientlyrr) rrUrzrhrrrrVrWs rX mpc_mpf_divrsw DAq! WQq\\4733A 1q$ , ,B 1&&4 5 5B r6MrZcd}d}|r:|dzr||z||zz ||z||zz}}|dz}||z||zz d|z|z}}|dz}|:||fS)zgComplex integer power: computes (a+b*I)**n exactly for nonnegative n (a and b must be Python ints).rrrl)rrrwrewims rXcomplex_int_powrs C C  q5 1us1u}c!ec!emC FAsQqSy!A#a%1 a  8OrZc |dtkrt||d||Sttt ||dz||dz||S)Nrrr)r mpc_pow_mpfmpc_exprmpc_log)rUrrzrhs rXmpc_powrsWtu}}1adD#... 771d2g..47;;T3 G GGrZc |\}}}}|dkrt|d|z||zz||S|dkr+t||dz}t|d|z|z||Sttt ||dz||dz||S)Nrr) mpc_pow_intmpc_sqrtrrr) rUrrzrhpsignpmanpexppbcsqrtzs rXrrsE4s qyy1rEkT4Z8$DDD rzzDG$$52+"4dC@@@ ;wq$r'22AtBw??s K KKrZc |\}}|tkrt||||tfS|tkrmt||||}|dz}|dkr |tfS|dkr t|fS|dkrt|tfS|dkrtt|fS|dkrtS|dkrt |||S|dkrt |||S|dkrt |||S|dkr$t t|| |dz||S|\}}} } |\} } } }|r| }| r| } | | z }t|}||t| |zz}|dkrr|dkr||z}| } n| | z} | } t|| |\}}t|t|| z||}t|t|| z||}||fSttt||dz||dz||S) Nrrrlrkri'r)rr@r)mpc_onerrrrabsmaxrrrrrrr)rUrrzrhrrvasignamanaexpabcbsignbmanbexpbbcdeabs_de exact_sizerVrWs rXrrsL DAqEzz1as++U22Ezz 1dC ( ( Q 66e8O !VV!8O !VV1::u$ $ !VV'!**$ $Avvg~Avvgas+++AvvjD#...Bww~as3331uu^KA2tAv$>$>cJJJE4sE4s dUd dUd B WWFFSc]]*+JE 66 RKDDD rcNDD tQ//B "c!D&kk4 5 5 "c!D&kk4 5 52v ;wq$r'22AtBw??s K KKrZc|\}}|tkrX|tkr||fS|dr'tt|||}t|fSt|||}|tfS|dz}|dsutt ||f|||}t |d} t| ||}t |d} t| |} t || ||}ntt ||f|||}t |d} t| ||}t |d} t| |} t || ||}|drt|}t|}||fS)zComplex square root (principal branch). We have sqrt(a+bi) = sqrt((r+a)/2) + b/sqrt(2*(r+a))*i where r = abs(a+bi), when a+bi is not a negative real number.rrr)rr0r)r*rr/r-r+) rUrzrhrrrWrVrrrrrs rXrr's DAqEzz ::q6M Q4 '!**dC00B2; !T3''B;  bB Q4 WaVR((!R 0 0 a   as # # aOO a__ Q4 % % GQFB''B / / a   as # # aOO a__ Q4 % % Q4 BB r6MrZczd}tt||||zz }tt||||zz } |d|zzd|z z}|j}|j} t t|}t t| } n#t $r~t ||}t ||}t |} td| |} t||f| tf|\}} t|}t| } YnwxYwd} |} |}t||| zD]}t|| |dz \}}t||dz | z|z |z }t||dz | z|z |z }||z||zz||zz }t|||z }t|||z }||z||zz|z }| |z||zz|z }||z|z}||z|z}||dz t||| z zz|z}||dz t| || z zz|z} |} || fS)N2y?g?rr)rrrrealrprrtrr2rrrrrr)rrrrzstarta1b1rrVrWfnnthextraprevpextra1rre2im2r4apbprecimcrebimbs rXmpc_nthroot_fixedrKs E VAtag~ & & ' 'B VAtag~ & & ' 'B  "r'\SU # V V R\\ R\\  b%  b%  a[[1b%(("bC<77B BZZ BZZ E E F U + +"2r1Q3//SS1Q3+/F233S1Q3+/F233#gCQZ 0 Atax  Atax Cx"s("q(sSy28#)axBaxBQqS&QuW---- 1QqS&QuW---- 1 r6Ms ABBD'&D'cV|\}}|ddkr&|tkrt||||}|tfS|dkr|dkrtS|dkrt||f||S|dkrt t||f||St ||f| |dzt |}t t|||S|dkrtd|dzz}|\} } } } |\} }}}t||f|}|d |dzd kr||d |dz|krgt||}t||}t||||\}}d}t|| |z ||}t|| |z ||}||fSt|}|dzdz}td||}t||f|tf||\}}t|d|d|d|d ||}t|d|d|d|d ||}||fS) zu Complex n-th root. Use Newton method as in the real case when it is faster, otherwise use z**(1/n) rrlrrrg333333?rirk)rrPrrr mpc_nthrootrrrrrrrrr2rr)rUrrzrhrrrVinverseprec2rrrrrrrrpfafbfrWrrrs rXrrrs[ DAqtqyyQ%ZZ AtS ) )E{1uu 66N 66Aq64-- - 777QFD#66 6q!fqb$q&.2EFFws333BwwC4"9%&&!"tT3!"tT3 aUD ! ! b6BrF?S bfr"vo&<&<!U##B!U##B&r2q%88FBEb5&,s;;Bb5&,s;;Br6M !B GbLE q"e $ $C aVc5\5# 6 6FB 2a5"Q%A1tS 9 9B 2a5"Q%A1tS 9 9B r6MrZc&t|d||S)z Complex cubic root. rk)rrUrzrhs rXmpc_cbrtrs q!T3 ' ''rZc0|\}}|tkrt|||S|tkrt|||tfSt||dz|}t||dz|\}}t||||}t||||} || fS)av Complex exponential function. We use the direct formula exp(a+bi) = exp(a) * (cos(b) + sin(b)*i) for the computation. This formula is very nice because it is pefectly stable; since we just do real multiplications, the only numerical errors that can creep in are single-ulp rounding errors. The formula is efficient since mpmath's real exp is quite fast and since we can compute cos and sin simultaneously. It is no problem if a and b are large; if the implementations of exp/cos/sin are accurate and efficient for all real numbers, then so is this function for all complex numbers. r)rr=r;r,) rUrzrhrrrrrrVrWs rXrrs DAqEzz1dC(((Ezzq$$$e++ !T!VS ! !C q$q&# & &DAq as # #B as # #B r6MrZcht|d|d||}t|||}||fS)Nrr)rArrys rXrrs7 qtQqT4 - -B D#  B r6MrZcX|\}}|tkrt|||tfS|tkrt|||tfS|dz}t||\}}t ||\}} t ||||} t || ||} | t | fS)aSComplex cosine. The formula used is cos(a+bi) = cos(a)*cosh(b) - sin(a)*sinh(b)*i. The same comments apply as for the complex exp: only real multiplications are pewrormed, so no cancellation errors are possible. The formula is also efficient since we can compute both pairs (cos, sin) and (cosh, sinh) in single stwps.)rrDrJr=r>r,r) rUrzrhrrrrrchshrVrWs rXmpc_cosrs DAqEzzq$$$e++Ezz4%%u,, B q"  DAq 1b ! !FB Bc " "B Bc " "B wr{{?rZc>|\}}|tkrt|||tfS|tkrtt|||fS|dz}t||\}}t ||\}} t ||||} t || ||} | | fS)zComplex sine. We have sin(a+bi) = sin(a)*cosh(b) + cos(a)*sinh(b)*i. See the docstring for mpc_cos for additional comments.r)rrErKr=r>r,rs rXmpc_sinrs DAqEzzq$$$e++Ezzhq$,,,, B q"  DAq 1b ! !FB Bc " "B Bc " "B r6MrZc|\}}|\}}}}|\} } } } |tkrt|||tfS|tkrtt|||fS|dz} t|d}t|d}t || \}}t || \}}t ||| }t||||}t||||}||fS)zcComplex tangent. Computed as tan(a+bi) = sin(2a)/M + sinh(2b)/M*i where M = cos(2a) + cosh(2b).r)rr?rLr/r=r>r*r-)rUrzrhrrrrrrrrrrrrrrrrrVrWs rXmpc_tanrs DAqE4sE4sEzz'!T3//66Ezz%!T3!7!777 B!QA!QA q"  DAq 1b ! !FB !R  C Cs # #B S$ $ $B r6MrZc|\}}|tkrt|||tfSt|t|dz|dz}|tkrt |||tfS|dz}t ||\}}t ||\}} t||||} t|| ||} | t| fSNrr)rrFr,r:rJrBr>r)rs rX mpc_cos_pirs DAqEzz!T3''..6$q&>>46**AEzz4%%u,, B !R DAq 1b ! !FB Bc " "B Bc " "B wr{{?rZc|\}}|tkrt|||tfSt|t|dz|dz}|tkrtt |||fS|dz}t ||\}}t ||\}} t||||} t|| ||} | | fSr)rrGr,r:rKrBr>rs rX mpc_sin_pir s DAqEzz!T3''..6$q&>>46**AEzzhq$,,,, B !R DAq 1b ! !FB Bc " "B Bc " "B r6MrZc|\}}|tkr&t|||\}}|tft|ffS|tkr&t|||\}}|tf|tffS|dz} t|| \}}t|| \}}t||||} t||||} t||||} t||||} | t | f| | ffS)Nr)rr>r=r,r))rUrzrhrrrrrrrcrecimsresims rX mpc_cos_sinr&s DAqEzzq$,,BE{UBK''Ezz1dC((15zAu:%% B q"  DAq 1b ! !FB !Rs # #C !Rs # #C !Rs # #C !Rs # #C  c **rZc(|\}}|tkr&t|||\}}|tf|tffSt|t|dz|dz}|tkr&t |||\}}|tft|ffS|dz} t|| \}}t || \}}t||||} t||||} t||||} t||||} | t | f| | ffSr)rrBr,r:r>r))rUrzrhrrrrrrrr"r#r$r%s rXmpc_cos_sin_pir(%s! DAqEzzas++15zAu:%%6$q&>>46**AEzzq$,,BE{UBK'' B !R DAq 1b ! !FB !Rs # #C !Rs # #C !Rs # #C !Rs # #C  c **rZcL|\}}t|t|f||S)z:Complex hyperbolic cosine. Computed as cosh(z) = cos(z*i).)rr)rs rXmpc_coshr*7s' DAq Awqzz?D# . ..rZc@|\}}t||f||\}}||fS)z;Complex hyperbolic sine. Computed as sinh(z) = -i*sin(z*i).)rrs rXmpc_sinhr,<- DAq Aq64 % %DAq a4KrZc@|\}}t||f||\}}||fS)z>Complex hyperbolic tangent. Computed as tanh(z) = -i*tan(z*i).)rrs rXmpc_tanhr/Br-rZc|\}}|dz}tt||t|f}tt|||f}t ||}t ||} t || ||\}}tt |dt |df} | dtkrt|r| dtf} | S)Nrrrr) r*r r)r+rrr/r%rYr) rUrzrhrrrryl1l2rs rXmpc_atanr4Is DAq Bagajj(AaaA BB BB 2r4 % %DAq !B  )Ab//1A tt|| 1 | qT5M HrZg:pΈ?g?c 0 |\}}|dz}|tkrttt||}|ds6|dkrt |||tfSt |||tfS|drdt ||}tt|||} |dkr|t| fStt|d| fSt|||} |dkr t| fSt ||}t|dt| fSdx} } |drt|}d} |drt|}d} tt||}tt||} t| ||} t|||}tt| ||d}t|||}t|||}tt||ds*|dkrt ||}nt ||}nt|||}|dst|t| | ||} t|||}tt|t| |||d}|dkr.ttt!|||||}ntt|t!||||}nt|t| | ||} t|t||||}tt| ||d}t|t!|||}|dkr tt||||}ntt||||}tt"||dst|t| | ||}t|drBt|||}t|||}tt|||d}n0t|||}tt|||d}t|t|t||}t%ttt|t!|||||}nRt!tt|||t||}t%t||||}| r4|dkrtt |||}nt|}| s|dkrt|}| r|dkrt|}t'|d|d|d|d||}t'|d|d|d|d||}||fS)a& complex acos for n = 0, asin for n = 1 The algorithm is described in T.E. Hull, T.F. Fairgrieve and P.T.P. Tang 'Implementing the Complex Arcsine and Arcosine Functions using Exception Handling', ACM Trans. on Math. Software Vol. 23 (1997), p299 The complex acos and asin can be defined as acos(z) = acos(beta) - I*sign(a)* log(alpha + sqrt(alpha**2 -1)) asin(z) = asin(beta) + I*sign(a)* log(alpha + sqrt(alpha**2 -1)) where z = a + I*b alpha = (1/2)*(r + s); beta = (1/2)*(r - s) = a/alpha r = sqrt((a+1)**2 + y**2); s = sqrt((a-1)**2 + y**2) These expressions are rewritten in different ways in different regions, delimited by two crossovers alpha_crossover and beta_crossover, and by abs(a) <= 1, in order to improve the numerical accuracy. rrrrrlrk)rr+r r'rNrMr:rOr)r/r*r1r-r,beta_crossoverrHr0alpha_crossoverr<r)rUrzrhrrrrampirrrrrralphabetab2rVAxrc1c2Am1rWs rX acos_asinrA_s" DAq BEzz T71::r * *!u 9Avv4--u444--u44t 9D#&&gajj$4466wqzz>)"9R#4#455q88as++66 !8Oc**B$R,,gajj88EEt AJJt AJJ q"  B q"  B"aA"aA gaB'' , ,E 1eR D 1b  B >4 , ,Q /&6 66$##BB$##BBUAr " "!u 6 GAr2..33A2r""A72wq!R'8'8"==rBBBAvvghr2&6&62>>CCga"b)9)92>>CC GAr2..33AGAr2..33A71a,,b11BHR,,b11BAvvgb!R00"55gaR00"55 ?E2 . .q 11 RB++R 0 0 2;;q> 5B##BR$$BGBB//44CCB##BGBB//44C S'%r22B 7 7 WT73R0@0@"#E#ErJJB O OggeUB77rBBB G G WUB++R 0 0  66R,,BBB Q!VV R[[ a R[[ 2a5"Q%A1tS 9 9B 2a5"Q%A1tS 9 9B r6MrZc&t|||dSNrrAr s rXmpc_acosrE Qc1 % %%rZc&t|||dSrrDr s rXmpc_asinrHrFrZct|\}}t|t|f||\}}t||fSre)rHr)rs rX mpc_asinhrJs: DAq a_dC 0 0DAq 1::q=rZct|||\}}|ds |tkrt||fS|t|fSrC)rErr)rs rX mpc_acoshrLsM AtS ! !DAqtqEzzqzz1}'!**}rZcD|dz}t|t|}tt||}t||}t||}t t|||d}|dt krt |rt|df}|S)Nrrrr)rrrrrr%rYr)rUrzrhrrrrs rX mpc_atanhrNs B7BABA2A2A'!Q##R((A tt|| 1 | AaDM HrZc|\}}|tkrt|||tfStt|d|dzt|d|dz}||zdz}t |}t t |dt|}t|tf||} t||} t| | |} t| | |} t| |||} | S)Nrlrkrr) rrQrrrCr*r/r&rrrrr) rUrzrhrVrWsizerrrrrs rX mpc_fibonaccirQ s FB U{{b$,,e44 s2a5A;RU2a5[!1!1 2 2D r B A !Q++AE Ar""A1bA1bA1bAAq$$$A HrZfctrer9rrzrhs rXmpf_expjrV rZc`|\}}|tkrt|||S|tkr%tt|||tfStt||dz}t||dz\}}t ||||}t ||||}||fSr)rr=r;r)r,)rUrzrhrVrWeyrrs rXmpc_expjrZs FB U{{2tS))) U{{wr{{D#..55 d2g & &B r47 # #DAq Qc " "B Qc " "B r6MrZctrerTrUs rX mpf_expjpir\(rWrZc|\}}|tkrt|||S|\}}}}|dz} |r| td||zz } tt t | || }|tkrt |||tfSt ||dz} t||dz\} } t | | ||}t | | ||}||fS)Nrr)rrBrr)r,r:r;) rUrzrhrVrWsignmanexpbcrrYrrs rX mpc_expjpirb+s FB U{{b$,,,D#sB bB  c!SVnn R,, - -B U{{r4%%u,, T"W  B "d2g & &DAq Qc " "B Qc " "B r6MrZsagez&Warning: Sage imports in libmpc failed)rR)__doc__rmbackendrrrrrlibmpfr r r r r rrrrrrrrrrrrrrrrrrr r!r"r#r$r%r&r'r(r)r*r+r,r-r.r/r0r1r2r3r4r5r6r7r8r9 libelefunr:r;r<r=r>r?r@rArBrCrDrErFrGrHrIrJrKrLrMrNrOrPrQrr}mpc_twompc_halfrSr\rYr]rcrirvr{r~rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr r&r(r*r,r/r4r6r7rArErHrJrLrNrQrVrZr\rbsage.libs.mpmath.ext_libmplibsmpmath ext_libmp_lbmp ImportErrorAttributeErrorprintrrZrXrrs  ==============                                                                                                       + %< + 5> u  5$  <<<# IIII    *&&&&'>>>> !+''''j>>>> *''''$8888j8888,,,$&&&& $&&&& &<<<<%::::%::::%:::: '$'!+ &0!+ ' < < < <!+!+!+   'HHHH !+LLLL!+&L&L&L&LP%""""H%%%N!+''''R%(((( $6$ $($ $&'    '    (++++"!+++++$%//// % %%    &F##*S//FFFP%&&&&%&&&&& &&      *         $ f8222222222222->  (888 67777778 sH H65H6