L im ddlZddlZddlZddlmZmZmZmZmZm Z m Z ddl m Z m Z mZmZddlZddlmZmZmZmZddlmZmZddlmZddlZddlZddlmZddlm cm!Z"ddZ#d Z$Gd d Z%Gd d Z&GddZ'GddZ(y)N)sqrtcossinarctanexplogpi)assert_assert_allcloseassert_array_lessassert_almost_equal)quaddblquadtplquadnquad)erferfc)LowLevelCallable) sine_ctypescJ|\}}t|||d| t||yy)Nr)atolrtol)r r ) value_and_err tabled_valueerror_tolerancevalueerrs i/mnt/ssd/data/python-lab/Trading/venv/lib/python3.12/site-packages/scipy/integrate/tests/test_quadpack.py assert_quadrs.JE3Ecallbacks rtest_ctypes_variantsz#TestCtypesQuad.test_ctypes_variantsCs%h&,oovH&h&,llFNN6??4S&,oo7&h&,oo7&h&,llFNN6??4SU&h&,llFOOE5%6en !7!* ?GAt'-H'' j$!R@Xq" 5a 8#>  ?!- =GAt{"T1b 1! 4c: j$a<  =r N)__name__ __module__ __qualname__r?rErHr_r rr.r. s%4G 2"=r r.c$eZdZdZdZdZdZy)TestMultivariateCtypesQuadctj}tjtjf}dD]}t||g|}t |||y)N)_multivariate_typical_multivariate_indefinite_multivariate_sin)r%r:rPr,setattr)r;r)r*r(r>s rr?z'TestMultivariateCtypesQuad.setup_methodisK//LL&//2* &D(wBBD D$ % &r cPtt|jdtddy)Nrg?ۗ?)rrrgr rDs rrEz'TestMultivariateCtypesQuad.test_typicalqs D33QHE* ,r cbtt|jdtjdy)Nrox?)rrrhnpinfrDs rtest_indefinitez*TestMultivariateCtypesQuad.test_indefinitevs"D66266B. 0r c>fd}tt|dddy)Nc>|tjdddzSrG)rri)yr;s r threadsafetyzBTestMultivariateCtypesQuad.test_threadsafety..threadsafety}s"tD22Aq9!<< .myfuncs"qs1SV8|$R' 'r rrlrn)rrr r;rs rrEzTestQuad.test_typicals (DB13IJr cTd}tt|dtjdy)Nc4t|  t|zSr")rrr}s rrz(TestQuad.test_indefinite..myfuncsG8CF? "r rrprrrqrrrs rrszTestQuad.test_indefinites #DBFF+-GHr c d}tt|ddddgdtdz tdztd z y) Ncfd|cxkrdkr t|Sd|cxkrdkrnyt| Sy)Nr@@g)rrrs rrz&TestQuad.test_singular..myfuncs;1{s{1v SA2wr r rr)pointsrBgg)rrrrrs r test_singularzTestQuad.test_singulars?  DBSz:CL3t9,s4y8 :r c d}d}tt|dddd|dt|z|t|zz |t dzzd |d zzz y) Nc$t||dz zSNrBrr}as rrz2TestQuad.test_sine_weighted_finite..myfuncsq!A#w< r glSS%@rrBrargsweightwvariirm)rrrrr)r;romes rtest_sine_weighted_finitez"TestQuad.test_sine_weighted_finitesZ DABu3GC[SX-c#c(l:US!V^L Nr c |d}d}d}tt|dtj|d|||dz|dzzz y)Nc t| |zSr"rrs rrz4TestQuad.test_sine_weighted_infinite..myfuncsr!t9 r @@rrrrmrr;rrrs rtest_sine_weighted_infinitez$TestQuad.test_sine_weighted_infinitesD  DBFF5sKAQ' )r c ~d}d}d}tt|tj d|d|||dz|dzzz y)Nct||zSr"rrs rrz6TestQuad.test_cosine_weighted_infinite..myfuncsqs8Or rgffffff@rrrrmrrs rtest_cosine_weighted_infinitez&TestQuad.test_cosine_weighted_infinitesF  D"&&!!ELq!tc1f}% 'r c ~d}d}tt|dd|ddttdd| zzdzdz z y) Nc dd|zd| zzz SNrBrmrcrs rrz2TestQuad.test_algebraic_log_weight..myfuncsac!qb'k? "r ?rBalg)rrrm)rrr r)r;rrs rtest_algebraic_log_weightz"TestQuad.test_algebraic_log_weightsJ # DQQu*,tQqA2wYNQ.// 1r c d}d}dtdzdtd| zdzd| zdzz zz td |d zzz td |zz d| zdzz }tt|d d dd d |dy)Nc.d| z|dz dzd| zzz S)NrNrBrmrrcrs rrz-TestQuad.test_cauchypv_weight..myfuncs&!9qsQhsaRy01 1r g?g\`@?rg]`@?rrBrNrmrrAcauchyrg#fT>r)rrrr)r;rr tabledValues rtest_cauchypv_weightzTestQuad.test_cauchypv_weights 2 #CH,#Cr2#)A+(F$GGHsQqSz*+sAv'qb A ' DACsK 9r cd}t|dtjd\}}t|tjdd\}}t|| t ||y)Nc:|tj| |zzSr"rqr)r}pqs rfz&TestQuad.test_b_less_than_a..fsrvvqbd|# #r rrmrrrrqrrr maxr;rval_1err_1val_2err_2s rtest_b_less_than_azTestQuad.test_b_less_than_asN $Aq"&&v6 uArvvqv6 uvCu,=>r cd}t|tj tjd\}}t|tjtj d\}}t|| t ||y)Ncrtj|dz dz |z tjd|zz S)NrmrN)rqrr)r}ss rrz(TestQuad.test_b_less_than_a_2..fs0661a4%!)a-(2772a4=8 8r )rmrrrrs rtest_b_less_than_a_2zTestQuad.test_b_less_than_a_2sZ 9AwT: uArvvwT: uvCu,=>r cd}t|dddd\}}t|dddd\}}t|| t||y)NcyN?rcrs rrz(TestQuad.test_b_less_than_a_3..fr rrBrrr)rrr)rr rrs rtest_b_less_than_a_3zTestQuad.test_b_less_than_a_3sJ Aq!E? uAq!E? uvCu,=>r cd}t|ddddd}t|ddddd}t|d|d}t|d|d |y) Ncyrrcrs rrz2TestQuad.test_b_less_than_a_full_output..frr rrBrrT)rr full_outputr)rrr )r;rres_1res_2rs rtest_b_less_than_a_full_outputz'TestQuad.test_b_less_than_a_full_outputs[ Q1UTJQ1UTJ%(E!H%a58)#6r c bd}d\}}tt|||ddd|dz|dzz zy)Nc ||zSr"rcrvr}s rsimpfuncz/TestQuad.test_double_integral..simpfuncs Q3Jr rrNc|Sr"rcrs rz/TestQuad.test_double_integral..sar c d|zSNrmrcrs rrz/TestQuad.test_double_integral..s 1Q3r g?rrrr;rrbs rtest_double_integralzTestQuad.test_double_integrals= 1GHaKGQVAsF]+ -r c Pd}d}d}d}tt|dd|||dy) Nc||z|z|zSr"rc)x0x1t0t1s rr>z,TestQuad.test_double_integral2..funcs7R<"$ $r c|Sr"rcrs rgz)TestQuad.test_double_integral2..gsHr c d|zSrrcrs rhz)TestQuad.test_double_integral2..hs q5Lr )rBrmrBrmrg$@r)r;r>rrrs rtest_double_integral2zTestQuad.test_double_integral2s/ %  GD!Q148Fr c <d}tt|dddddy)Nc||zdzdzSrrcrrs rr>z,TestQuad.test_double_integral3..funcs7Q;? "r rBrmg@r)r;r>s rtest_double_integral3zTestQuad.test_double_integral3s #GD!Q1-b1r z,x_lower, x_upper, y_lower, y_upper, expectedrrrBrmc @d}tt||||||dy)Nc@tj|dz |dzz Srrr}rvs rrz1TestQuad.test_double_integral_improper..fDs 6616'AF*+ +r g+i)+`>rr)r;x_lowerx_uppery_lowery_upperexpectedrs rtest_double_integral_improperz&TestQuad.test_double_integral_improper s(t ,  Aw' :   r c ld}d\}}tt|||dddddd|d z|d zz zy) Nc||z|z|zSr"rc)rrvr}ts rrz/TestQuad.test_triple_integral..simpfuncOsaCE19 r rc|Sr"rcrs rrz/TestQuad.test_triple_integral..Tsar c d|zSrrcrs rrz/TestQuad.test_triple_integral..Ts 1Q3r c ||z Sr"rcrs rrz/TestQuad.test_triple_integral..Us Qr c ||zSr"rcrs rrz/TestQuad.test_triple_integral..Us QUr )rNgUUUUUU@rrrrs rtest_triple_integralzTestQuad.test_triple_integralMsK 1GHa'.0B!#3C0  2r z>x_lower, x_upper, y_lower, y_upper, z_lower, z_upper, expectedrrc Dd}tt||||||||dy)NcLtj|dz |dzz |dzz Srr)r}rvrs rrz1TestQuad.test_triple_integral_improper..fs)6616'AF*Q!V34 4r g+i)+p>rr) r;rrrrz_lowerz_upperrrs rtest_triple_integral_improperz&TestQuad.test_triple_integral_improperYs,v 5  Aw'7G L   r cd}tjt|dtjdz dddsJdtjdddd }t|fd di|}td fd d i|}td fd d i|}tj j |d|dd|dzztj j |d|dd|dzzt|ddt|ddcxk(rdk(sJJ|ddd|dk(sJ|ddd|dk(sJ|dddd|ddk(sJt|ddt|ddcxk(rdk(sJJ|dddd|ddk(sJy)Nc2tjd|zSN?rrs rtfuncz$TestQuad.test_complex..tfuncs66"Q$< r rrmT) complex_funcy??rrB)rrrrrrcXtjtjd|zSr )rqrealrrs rrz'TestQuad.test_complex.. rwwrvvbd|4r FcXtjtjd|zSr )rqimagrrs rrz'TestQuad.test_complex..rr r rrrlstr)rqallcloserr rrtesting assert_equallen)r;rkwargsres_cres_rres_is r test_complexzTestQuad.test_complexs {{2557>qA rvvd!1.U8884"'4"' a%(Ra[*@A a%(Ra[*@A58F#$E!"I;!;;;;;Qx"eAh...Qx"eAh...Qx"5)U1Xe_<<<58F#$E!"I;!;;;;;Qx"5)U1Xe_<<ffWa"&&!RUUQY#a&1*%= >ffWa"&&!RUUQY3q6A:!2C%D EffWb266'1beeaiCFQJ$q'3I&J KffWa"&&"beeaiCFQJ$q'3I&J K26625519 -26625519tAw#6 7 26625519tAw#6 726625519Q1 #= >BFFBEEAIQ!$< = BFFBEEAIQ!$< =RVVRUUQY3q6A:!2C%D EBFFBEEAI#a&1*Q1G$H IBFFBEEAI#a&1*Q1G$H IffWbffrvvgrvvruu 5c2 5l  m5l   2 [[ [[HK ffWa"&&!bffWa"%%E2Ba1G HK ffWb266'1rvvgqee! #d1g - / K ffWa"&&"rvvgqee! #d1g - /K ffWa"&&!bffWbee! #d1g - /K ffWb266'2wee! #tAw!| 4 6K ffWb266'1rvvgree! #tAw!| 4 6!K $ffWa"&&"rvvgree! #tAw!| 4 6%K ,ffWb266'2wee! #tAw!| 4 6-K 4ffWb266'1rvvgqee! #Q!(9T!W'D E G5K <ffWb266'2wee! #A tAw!|'D E G=K DffWb266'1rvvgree! #A tAw!|'D E GEK LffWa"&&"rvvgree! #A tAw!|'D E GMK TffWa"&&!bffWbee! #Q!(9T!W'D E GUK \ffWa"&&"rvvgqee! #Q!(9T!W'D E G]K dffWa"&&!bffWaee! #s1vz 2 4eK hffWa"&&!bffWaee! #s1vz 2 4iK lffWa"&&!bffWaee! #s1vz 2 4mK tffWa"&&!bffWaee! #A q'8 9 ;uK xffWa"&&!bffWaee! #A q'8 9 ;yK |ffWa"&&!bffWaee! #A q'8 9 ;}K DffWa"&&!bffWaee! #A q'8 9 ;EK L2661bffruu/?1.D EMK R2661bffee! #d1g - /SK V2661bffee! #d1g - /WK Z2661bffee! #d1g - /[K b2661bffee! #tAw!| 4 6cK f2661bffee! #tAw!| 4 6gK j2661bffee! #tAw!| 4 6kK r2661bffee! #tAw!| 4 6sK zBFFArvvee! #s1vz 2 4{K ~BFFArvvee! #s1vz 2 4K B2662rvvee! #s1vz 2 4CK JRVVQee! #A q'8 9 ;KK NBFFBee! #A q'8 9 ;OK RBFFBee! #A q'8 9 ;SK ZRVVRee! #A q'8 9 ;[K bBFFBee! #Q!(9T!W'D E GcK j2662rvvee! #A tAw!|'D E GkK rBFFArvvee! #A tAw!|'D E GsK zBFFArvvee! #A tAw!|'D E G{K BRVVQee! #Q!(9T!W'D E GCK JBFFBee! #Q!(9T!W'D E GKK RffWbffrvvgrvvw UUu  SK N^ _N` (=r rzceZdZejj ddZejj ddZdZdZ dZ dZ dZ d Z d Zd Zy ) TestNQuadrAcd}d}t|ddgddgddgddgg|iiigd }t|ddd t|dd dkDxr |dd d ky)Nc|dz||zz|dzz tj|z|d|zz dz d|zz dkDrdz}|Sdz}|S)Nrmr皙???rrBrqr)rrx2x3vals rfunc1z*TestNQuad.test_fixed_limits..func1"shq52b5=2q5(266":5s2v++d2g59ACCJABCCJr c0dd|dzdzd|dzzgiS)Nrr&rmr'r(rrcrs r opts_basicz/TestNQuad.test_fixed_limits..opts_basic's*s47{S04Q<?@A Ar rrBrgp= ף?g?g333333ÿT)optsrg 5m?nevalgjA)rrr )r;r-r/ress rtest_fixed_limitszTestNQuad.test_fixed_limits s   BEQFRGc2Yq B$b"b1tECH01B 1$?R)9C)?@r c  d d} fd} fd} fd} fd}d}d}d }d } t|||||gd |||| g } t| d y)Ng?c|||z|dzztj|zdz|||zz|z dkDrdz}|Sdz}|S)NrmrBrr))rrr*r+rrr,s rfunc2z-TestNQuad.test_variable_limits..func23sVb5Q;+a/beb1,A5CJ345CJr c|dz|ztj||z|zzdzzdz |dz|ztj||z|zzdzzdzgSNrmrB)rqr)rr*r+rrscales rlim0z,TestNQuad.test_variable_limits..lim08soRURZ"&&*R-*::Q>?!CRURZ"&&*R-*::Q>?!CE Er cH||z||zzzdz ||z||zzzdzgSrrc)r*r+rrr9s rlim1z,TestNQuad.test_variable_limits..lim1<s?RURU]+a/RURU]+a/1 1r cT||dz|dzzzzdz ||dz|dzzzzdzgSNrmrrBrc)r+rrr9s rlim2z,TestNQuad.test_variable_limits..lim2@sIR"a%A+-.2R"a%A+-.24 4r c0||zzdz ||zzdzgSrrc)rrr9s rlim3z,TestNQuad.test_variable_limits..lim3Ds*R"W%)5BG+.opts0GsrBrEzl+ +r ciSr"rc)r*r+rrs ropts1z-TestNQuad.test_variable_limits..opts1JIr ciSr"rc)r+rrs ropts2z-TestNQuad.test_variable_limits..opts2MrFr ciSr"rc)rrs ropts3z-TestNQuad.test_variable_limits..opts3PrFr r)rr0g9@)rr) r;r6r:r<r?rArCrErHrJr2r9s @rtest_variable_limitszTestNQuad.test_variable_limits/si   E 1 4 B ,   ED$d3&u57C+,r cJd}tt|ddgddggiigdy)Ncyrrcrs rrz9TestNQuad.test_square_separate_ranges_and_opts..fXrr rrBr0rrr)r;rs r$test_square_separate_ranges_and_optsz.TestNQuad.test_square_separate_ranges_and_optsWs-  E!r1gAw/r2h?Er cNd}ddg}i}tt|||g||gdy)Ncyrrcrs rrz8TestNQuad.test_square_aliased_ranges_and_opts..f^rr rrBrNrrO)r;rropts r#test_square_aliased_ranges_and_optsz-TestNQuad.test_square_aliased_ranges_and_opts]s2 GE!aV3*5s;r cbd}d}d}d}d}||g}||g}tt|||dy)Ncyrrcrs rrz.ffrr cyN)rrBrcrs r fn_range0zDTestNQuad.test_square_separate_fn_ranges_and_opts..fn_range0ir cyrYrcrs r fn_range1zDTestNQuad.test_square_separate_fn_ranges_and_opts..fn_range1lr[r ciSr"rcrs rfn_opt0zBTestNQuad.test_square_separate_fn_ranges_and_opts..fn_opt0orFr ciSr"rcrs rfn_opt1zBTestNQuad.test_square_separate_fn_ranges_and_opts..fn_opt1rrFr rNrrO)r;rrZr]r_rarangesr0s r'test_square_separate_fn_ranges_and_optsz1TestNQuad.test_square_separate_fn_ranges_and_optsesD     Y'!E!V$/5r cVd}d}d}||g}||g}tt|||dy)Ncyrrcrs rrz;TestNQuad.test_square_aliased_fn_ranges_and_opts..fzrr cyrYrcrs rfn_rangezBTestNQuad.test_square_aliased_fn_ranges_and_opts..fn_range}r[r ciSr"rcrs rfn_optz@TestNQuad.test_square_aliased_fn_ranges_and_opts..fn_optrFr rNrrO)r;rrgrirbr0s r&test_square_aliased_fn_ranges_and_optsz0TestNQuad.test_square_aliased_fn_ranges_and_optsys:   H%E!V$/5r cd}t|dd\}}t|ddgg\}}t||t||y)Nc|dzdzSr8rcrs rr>z*TestNQuad.test_matching_quad..funcsa4!8Or rr)rb)rrr )r;r>r2reserrres2reserr2s rtest_matching_quadzTestNQuad.test_matching_quadsE 4A& VdQF84 gC&FG,r cd}t|dddd\}}t|ddgdg\}}t||t||y) Nc*|dz|dzz||zz dzSr>rcrs rfunc2dz/TestNQuad.test_matching_dblquad..func2ds"q52q5=27*Q. .r rmcy)Nrcrs rrz1TestNQuad.test_matching_dblquad..r cy)Nrrcrs rrz1TestNQuad.test_matching_dblquad..rwr rvr)rtrm)rrr )r;rsr2rmrnros rtest_matching_dblquadzTestNQuad.test_matching_dblquadsL /fb!\;G VfAw&89 gC&FG,r c d}t|ddddddd }t|tj tjgd dgd gd }t ||y) Ncb|dz||dzzz||zz dz|tj|zzSr>r))rrr*c0c1s rfunc3dz/TestNQuad.test_matching_tplquad..func3ds8q52A:%R/!3b266":oE Er rrmcy)Nrtrcrs rrz1TestNQuad.test_matching_tplquad..rwr cyrrcrs rrz1TestNQuad.test_matching_tplquad..rwr c$tj Sr"rqr rs rrz1TestNQuad.test_matching_tplquad..s BEE6r c"tjSr"rrs rrz1TestNQuad.test_matching_tplquad..s r rrrt)rrm)rrrqr r )r;r~r2rns rtest_matching_tplquadzTestNQuad.test_matching_tplquadsZ Ffb!\;)+=!#VvruuoAw@vNC&r cR tdddgddggddiy#t$rJwxYw)Nc ||zSr"rcrs rrz-TestNQuad.test_dict_as_opts..s q1ur rrBepsrelg-C6?rN)r TypeErrorrDs rtest_dict_as_optszTestNQuad.test_dict_as_optss:  $1v1v&6h=O P  5 s &N)r`rarbr8r fail_slowr3rKrPrUrcrjrpryrrrcr rr#r#sr [[1 A A [[1%-%-NF <6( 6--'r r#)g+i)+P>))r3rCnumpyrqrrrrrrr numpy.testingr r r r r8scipy.integraterrrr scipy.specialrrscipy._lib._ccallbackrr% ctypes.utilscipy._lib._ccallback_cr"scipy.integrate._test_multivariate integrate_test_multivariater$rr,r.rerzr#rcr rrs~ 666AA 99#2 /660B E=E=PBB4Z=Z=z HHr