
    (ph                     `    S SK rS SKJrJr  S SKJr  S SKJrJ	r	J
r
  S rS rS rS rS	 rS
 rg)    N)assert_allcloseassert_)FuncData)gammagammalnloggammac            	      B   [         R                  " / SQ5      n U R                  5       n[         R                  " X5      u  pU SU-  -   R	                  5       n[         R
                  " U[        U5      45      R                  nS n[        XCSSSSS9R                  5         g )Ng     Xg      #g            ?g      #@g     X@              ?c                 @    [         R                  " [        U 5      5      $ )N)npexpr   zs    T/var/www/html/venv/lib/python3.13/site-packages/scipy/special/tests/test_loggamma.pyftest_identities1.<locals>.f   s    vvhqk""    r      +=rtolatol)
r   arraycopymeshgridflattenvstackr   Tr   checkxyr   datasetr   s        r   test_identities1r&      s~    
45A	A;;qDA	
RTAiiE!H&((G# QAE6<<>r   c            	      p   [         R                  " / SQ5      n U R                  5       n[         R                  " X5      u  pU SU-  -   R	                  5       n[         R
                  " U[         R                  " U5      [        U5      -   45      R                  nS n[        XCSSSSS9R                  5         g )Nr
   r   c                     [        U S-   5      $ )Nr   )r   r   s    r   r   test_identities2.<locals>.f   s    Ar   r   r   r   r   )r   r   r   r   r   r   logr   r    r   r!   r"   s        r   test_identities2r+      s    
45A	A;;qDA	
RTAiiBFF1I34577G QAE6<<>r   c            	         [         R                  [         R                  " SS5      * [         R                  " SS5      4   S-   n [         R                  " U [	        U 5      45      R
                  nS n[        X!SSSSS9R                  5         g )	N
   r   c                 V    [         R                  " U SS9n [        U 5      R                  $ )N
complex128)dtype)r   r   r   realr   s    r   r   )test_complex_dispatch_realpart.<locals>.f+   s"    HHQl+{r   r   r   r   r   )r   r_logspacer   r   r    r   r!   )r#   r%   r   s      r   test_complex_dispatch_realpartr6   $   sq     	r{{2s##R[[b%99:S@AiiGAJ(**G  QAE6<<>r   c            	      h   [         R                  " SS5      S-   n [         R                  " U [        U 5      45      R                  n[        [        USSSSS9R                  5         [        [        S5      [         R                  :H  5        [        [         R                  " [        S5      5      5        g )	Nr.   r-   r   r   r   r   r   )r   r5   r   r   r    r   r   r!   r   infisnan)r#   r%   s     r   test_real_dispatchr;   2   sx    
Cs"AiiGAJ(**GXw15u=CCEHQK266!"BHHXb\"#r   c                      [        [        SS5      5      n [        [        SS5      5      n[        XR                  5       SSS9  g )Ng333333               V瞯<r   r   )r   complexr   	conjugate)r   zbars     r   test_gh_6536rC   ;   s7    t$%AGD$'(DA~~'e!<r   c            	      @   [         R                  " SSS5      * n [         R                  " U  Vs/ s H  n[        US5      PM     sn5      n[         R                  " U  Vs/ s H  n[        US5      PM     sn5      n[	        X#R                  5       SSS9  g s  snf s  snf )	Ni,  id   r=   r>   r?   r   r   )r   r5   asarrayr@   r   rA   )r#   x0r   rB   s       r   test_branch_cutrH   A   sz    	S#s	##A


q1qGB$q12A::151Rwr4(156DA~~'e!< 25s   BB)numpyr   numpy.testingr   r   scipy.special._testutilsr   scipy.specialr   r   r   r&   r+   r6   r;   rC   rH    r   r   <module>rN      s2     2 - 2 2???$==r   