
    ^Mh-                        d dl Z d dlmZ d dlZd dlZd dlZd dlm	Z	m
Z
 d Zd Zd Zd Zd Zej                            dg d	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7ej         ej        d8z  z   d9fd:ej        ej        d8z  z
  d9f          d;             Ze j        j        d<         d=         d>         Zed?k    o
ej        d@k    ZedAk    o
ej        dBk    Zej                            epe dCD          ej                            dg dE          dF                         ZdS )G    N)assert_equalassert_allclosec                  v    t          t          j        dd          t          j        dz  dz  dz
  d           d S )N         -q=rtol)r   sczetanppi     ]/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/scipy/special/tests/test_zeta.py	test_zetar   
   s5    BGAaLL"%(1*q.u======r   c                  ^    g d} g d}t          t          j        |           |d           d S )N)g g?gH.?	   2   K   )g^ңge
gp#	6gj-s`?gXk   <g    @;r	   r
   )r   r   zetac)xdesireds     r   
test_zetacr      sE     	'&&A  G BHQKKu555555r   c                  x   t          j        t          j                  dk    sJ t          j        t          j        t          j                             sJ t          j        d          dk    sJ t          j        d          t          j        k    sJ t          t          j        g d          d           d S )Nr   g      g      ?)ii)r   r   r   infisnanr   r   r   r   test_zetac_special_casesr!      s    8BFq    8BHbfW%%&&&&&8A;;$8C==BF""""///**B/////r   c                     t          j        t          j        t           j                            sJ t          j        t           j                  dk    sJ t          j        d          dk    sJ t          t          j        g d          d           t          t          j        d          t           j        dz  dz  d           t          t          j        d	          t           j        d	z  d
z  d           d S )Nr   r   g      )r   iiir   r   r	   r
      Z   )	r   r    r   r   nanr   r   r   r   r   r   r   test_riemann_zeta_special_casesr'   &   s    8BGBFOO$$$$$726??a71:: ...//333BGAJJq
7777BGAJJq%888888r   c                  X    d} d}t          t          j        |           |dd           d S )Ng    @pgkr   gvIh%,=)atolr   r   r   r   )sr   s     r    test_riemann_zeta_avoid_overflowr,   2   s1    A)GBGAJJae<<<<<<r   zz, desired, rtol)y܍	@ޓq@yӻQN?%Kh?vIh%<=)yz&E@ayf|~?xL{r-   )ynN9 S;@yV(?eyr-   )yK"Bnyy[r6nr-   )y?e5@MtN@y̌q]?AaO?r-   )yRp{fӰ@|/PyIU ?ž?r-   )y g
S@yϜ<@]@r-   )y@0Ry
)اAr-   )yv_)?0PZ?ygPq i?{^A޿r-   )y2?cB@y]'??r-   )y}d?wOP@y+j<?.ѿr-   )yvQ_?m5@y]hN?>ä?r-   )yi6?{V@@ypjXۦ`4?r-   )yv_)?0PZygPq i?{^A?r-   )y2?cBy]'?ƿr-   )y}d?wOPy+j<?.?r-   )yvQ_?m5y]hN?>är-   )yi6?{V@ypjXۦ`4ǿr-   )y       V瞯<ygOkWe9>/fHz>)y      V瞯<yqU-$yhb<r.   )y      V瞯<y|XZ_45[r.   )y       V瞯<ycQ*9G9,c<r.   )y       V瞯ҼyOkWe9>/f<r.   )y      V瞯ҼyQqU-$yhbr.   )y      V瞯ҼyrXZ_45[<r.   )y       V瞯ҼyZQ*9G9,cr.   )y      V瞯<yʿ)Cp r.   )y        V瞯<y       ۍмr.   )y?V瞯<yHG|Hd5ڼr.   )y      ?V瞯<y;	5r-   )y      V瞯Ҽyʿ)Cp <r.   )y       V瞯Ҽy      `ۍ<r.   )y?V瞯ҼyHG|0t5<r.   )y      ?V瞯Ҽy;	3=r-   )y;ziH@ESJ@y     ?,莕e<r-   )y8-CnGI@ESJ@y     ?ՌRP<r-   )y;ziH@)|QD@y?	Idr-   )y8-CnGI@)|QD@y?(r-   )y;ziH@yuU@y?}6 Gwr-   )y8-CnGI@yuU@y?r-   )y      ?    .Ayz?ی^p@r.   )y_+_hA@    .Ayf    ?
lR=r-   )y8-CnGI@    .Ay     ?Z <r-   )y      ?    .yz?ی^pr.   )y_+_hA@    .yf    ?
lRr-   )y8-CnGI@    .y     ?Z ļr-   yX\mrh`<+@y              ?r-   yX\mr%0c                 N    t          t          j        |           ||           d S Nr
   r*   zr   r   s      r   test_riemann_zeta_complexr3   8   s(    b BGAJJd333333r   	Compilerszc++namegcclinuxclangdarwinz/Underflow may not be avoided on other platforms)reason))yF5,kh`<+@yqdݻIZ-a=)ypmh`<+@y:q#Nxm Vkr-   )yF5,k%0y61;Jjr;   )ypm%0yC@Ҷ*Ѣ˺dr;   c                 N    t          t          j        |           ||           d S r0   r*   r1   s      r   )test_riemann_zeta_complex_avoid_underflowr=     s'    6 BGAJJd333333r   )scipyscipy.specialspecialr   sysnumpyr   pytestnumpy.testingr   r   r   r   r!   r'   r,   markparametrizer   r3   
__config__CONFIGcpp_compilerplatform	gcc_linuxclang_macOSskipifr=   r   r   r   <module>rN      sS          



      7 7 7 7 7 7 7 7> > >6 6 60 0 0	9 	9 	9= = = L|	}LB	CLH	ILN	OLV	WL\	]Lb	cLh	iLp	qLv	wL|	}LB	CLH	ILN	OLT	ULZ	[L`	aLf	gLn 	NoLp 	NqLr 	PsLt 	NuLv 	NwLx 	OyLz 	N{L| 	N}L@ 	NALB 	7CLD 	MELF 	NGLH 	NILJ 	8KLL 	MMLN 	NOLR	SLX	YL^	_Ld	eLj	kLp	qLx 	EyLz	{L@	ALF 	EGLH	ILN	OLT 2bfWRVBY5FOULV 2RVBF2I5ENWLO O`4 4aO O`4
 &{3E:6BE!=clg&=	g%B#,(*B 		!k"<       ,4 4- 	 44 4 4r   