o
    i}                     @   s  d dl Z d dl mZ d dlZd dlmZmZ d dlmZm	Z	m
Z
mZ d dlmZ d dlmZ d dlmZ d dlm  mZ G dd dZG d	d
 d
ZG dd dZG dd dZG dd dZG dd dZG dd dZG dd dZG dd dZ	dAddZdd Zdd  Z d!d" Z!d#d$ Z"d%d& Z#d'd( Z$d)d* Z%d+d, Z&d-d. Z'd/d0 Z(d1d2 Z)d3d4 Z*d5d6 Z+d7d8 Z,d9d: Z-d;d< Z.d=d> Z/d?d@ Z0dS )B    N)raises)arraysqrt)assert_array_almost_equalassert_equalassert_almost_equalassert_allclose)	integrategammac                   @   s,   e Zd Zdd Zdd Zdd Zdd Zd	S )
	TestChebyc                 C   s   t d}t d}tjdd t d}t d}t d}t d}W d    n1 s/w   Y  t|jdgd	 t|jddgd	 t|jg d
d	 t|jg dd	 t|jg dd	 t|jg dd	 d S )Nr      ignoreall               )r   r   )r   r   r   )r   r   r   r   )r   r   r   r   r   )orthchebycnperrstater   c)selfZC0ZC1ZC2ZC3ZC4ZC5 r    ]/home/kim/smarthome/.venv/lib/python3.10/site-packages/scipy/special/tests/test_orthogonal.pytest_chebyc   s   




zTestCheby.test_chebycc                 C      t d}t d}t d}t d}t d}t d}t|jdgd t|jddgd t|jg dd t|jg d	d t|jg d
d t|jg dd d S )Nr   r   r   r   r   r   r   )r   r   )r   r   r   r   )r   r   r   r   r   )r   r   r   r   r   r   )r   chebysr   r   )r   ZS0ZS1ZS2ZS3ZS4ZS5r    r    r!   test_chebys       





zTestCheby.test_chebysc                 C   r#   )Nr   r   r   r   r   r   r   )r   r   r$   )r   r   r   r   )   r   ir   r   )   r   r   r   r   )r   chebytr   r   )r   ZT0ZT1ZT2ZT3ZT4ZT5r    r    r!   test_chebyt.   r'   zTestCheby.test_chebytc                 C   s   t d}t d}t d}t d}t d}t d}t|jdgd t|jddgd t|jg dd t|jg d	d t|jg d
d t|jg dd d S )Nr   r   r   r   r   r   r   )r   r   r$   )r(   r   r   r   )r)   r   r   r   )    r   ir      r   )r   chebyur   r   )r   ZU0ZU1ZU2ZU3ZU4ZU5r    r    r!   test_chebyu<   r'   zTestCheby.test_chebyuN)__name__
__module____qualname__r"   r&   r,   r1   r    r    r    r!   r      s
    r   c                   @   s6   e Zd Zdd Zejdddgdd Zdd	 Zd
S )TestGegenbauerc                 C   s  dt j  d }t |dkrd}td|}td|}td|}td|}td|}td|}t|jtdgd	 t|jtd| dgd	 t|jtd| |d  d| gd	 t|jtdt	|d dd
| |d  dgd d t|jtdt	|d ddt	|d dd| |d  gd d t|jtdt	|d ddt	|d ddt	|d dgd d d S )Nr         ?r   皙ɿr   r   r   r   r   ig      @   r-   g      @r*      g      .@)
r   randomanyr   
gegenbauerr   r   r   scZpoch)r   aZCa0ZCa1ZCa2ZCa3ZCa4ZCa5r    r    r!   test_gegenbauerM   sB   &&&&zTestGegenbauer.test_gegenbauerr>   r   r   c                 C   s$   t d|}t|jtdg d S )Nr         ?)r   r<   r   r   r   Zasarray)r   r>   ZCn0r    r    r!   test_n_zero_gh8888b   s   z!TestGegenbauer.test_n_zero_gh8888c                 C   s   d}t jt|d tdtj W d    n1 sw   Y  t jt|d tdd W d    n1 s8w   Y  t jt|d tdtj  W d    d S 1 sXw   Y  d S )Nz*`alpha` must be a finite number greater...)matchr   r         r   )pytestr   
ValueErrorr   r<   r   naninf)r   messager    r    r!   test_valid_alphah   s   "zTestGegenbauer.test_valid_alphaN)	r2   r3   r4   r?   rD   markZparametrizerA   rI   r    r    r    r!   r5   K   s
    
r5   c                   @   s   e Zd Zdd Zdd ZdS )TestHermitec                 C   s   t d}t d}t d}t d}t d}t d}t|jdgd t|jddgd t|jg dd t|jg d	d t|jg d
d t|jg dd d S )Nr   r   r   r   r   r   r   )r   r   r   )r(   r   r-   r   )r)   r   ir      rL   )r.   r   i`r   x   r   )r   hermiter   r   )r   H0H1H2H3H4H5r    r    r!   test_hermitet   r'   zTestHermite.test_hermitec                 C   s2  t dtd dg}td}td}td}td}td}td}td|}td|td }	td|d }
td|dtd  }td|d	 }td|d	td  }t|j|jd
 t|j|	jd
 t|j|
jd
 t|j|jd
 t|j|jd
 t|j|jd
 d S )Nr@   r   r   r   r   r   r          @      @r   )r   poly1dr   r   hermitenormrN   r   r   )r   psubrO   rP   rQ   rR   rS   rT   Zhe0Zhe1Zhe2Zhe3Zhe4Zhe5r    r    r!   test_hermitenorm   s&   





zTestHermite.test_hermitenormN)r2   r3   r4   rU   r[   r    r    r    r!   rK   s   s    rK   c                   @      e Zd Zdd ZdS )TestShLegendrec                 C     t ddg}td}td}td}td}td}td}td|}td|}	td|}
td|}td|}td|}t|j|jd t|j|	jd t|j|
jd t|j|jd t|j|jd	 t|j|jd	 d S 
Nr   r$   r   r   r   r   r   r   rL   )r   rX   r   sh_legendrelegendrer   r   )r   rZ   ZPs0ZPs1ZPs2ZPs3ZPs4ZPs5Zpse0Zpse1Zpse2Zpse3Zpse4Zpse5r    r    r!   test_sh_legendre   &   





zTestShLegendre.test_sh_legendreN)r2   r3   r4   rb   r    r    r    r!   r]          r]   c                   @   r\   )TestShChebytc                 C   r^   r_   )r   rX   r   	sh_chebytr+   r   r   )r   rZ   ZTs0ZTs1ZTs2ZTs3ZTs4ZTs5Ztse0Ztse1Ztse2Ztse3Ztse4Ztse5r    r    r!   test_sh_chebyt   rc   zTestShChebyt.test_sh_chebytN)r2   r3   r4   rg   r    r    r    r!   re      rd   re   c                   @   r\   )TestShChebyuc                 C   s  t ddg}td}td}td}td}td}td}td|}td|}	td|}
td|}td|}td|}t|j|jd t|j|	jd t|j|
jd t|j|jd t|j|jd	 t|j|jd
 d S )Nr   r$   r   r   r   r   r   r   rL   r8   )r   rX   r   	sh_chebyur0   r   r   )r   rZ   ZUs0ZUs1ZUs2ZUs3ZUs4ZUs5Zuse0Zuse1Zuse2Zuse3Zuse4Zuse5r    r    r!   test_sh_chebyu   rc   zTestShChebyu.test_sh_chebyuN)r2   r3   r4   rj   r    r    r    r!   rh      rd   rh   c                   @   r\   )TestShJacobic                 C   s  dd }t ddg}dt j  }|d dt j   }td||}td||}td||}td||}td||}	td	||}
td|| |d ||d| }td|| |d ||d| }td|| |d ||d| }td|| |d ||d| }td|| |d ||d| }td	|| |d ||d	| }t|j|jd
 t|j|jd
 t|j|jd
 t|j|jd
 t|	j|jd
 t|
j|jd
 d S )Nc                 S   s(   t | d t | |  t d|  |  S )Nr   r   r
   )npr    r    r!   conv   s   (z)TestShJacobi.test_sh_jacobi.<locals>.convr   r$   r   r   r   r   r   r   )r   rX   r:   r   Z	sh_jacobiZjacobir   r   )r   rn   rZ   qrm   ZG0ZG1ZG2ZG3ZG4ZG5Zge0Zge1Zge2Zge3Zge4Zge5r    r    r!   test_sh_jacobi   s,   $$$$$$zTestShJacobi.test_sh_jacobiN)r2   r3   r4   rp   r    r    r    r!   rk      rd   rk   c                   @   r\   )TestCallc                 C   s   g }t dD ]}|dd dt|d  D  qtjdd! |D ]}t|}t|dt|j	d|d	 q#W d    d S 1 sDw   Y  d S )
Nr   c                 S   s   g | ]}|  qS r    )strip).0xr    r    r!   
<listcomp>      z&TestCall.test_call.<locals>.<listcomp>aW  
                orth.jacobi(%(n)d,0.3,0.9)
                orth.sh_jacobi(%(n)d,0.3,0.9)
                orth.genlaguerre(%(n)d,0.3)
                orth.laguerre(%(n)d)
                orth.hermite(%(n)d)
                orth.hermitenorm(%(n)d)
                orth.gegenbauer(%(n)d,0.3)
                orth.chebyt(%(n)d)
                orth.chebyu(%(n)d)
                orth.chebyc(%(n)d)
                orth.chebys(%(n)d)
                orth.sh_chebyt(%(n)d)
                orth.sh_chebyu(%(n)d)
                orth.legendre(%(n)d)
                orth.sh_legendre(%(n)d)
                )rl   r   r   g)\(?)err_msg)
rangeextenddictsplitr   r   evalr   rX   Zcoef)r   Zpolyrl   Zpstrrm   r    r    r!   	test_call  s    

"zTestCall.test_callN)r2   r3   r4   r}   r    r    r    r!   rq      rd   rq   c                   @   r\   )TestGenlaguerrec                 C   sx   t tjdddddd t tjdddddd t tjddddtdd	g t tjddddtdd
g d S )Nr   F)Zmonicr   rV   T       r$   r   r   )r   r   Zgenlaguerrer   rX   )r   r    r    r!   test_regression  s    $zTestGenlaguerre.test_regressionN)r2   r3   r4   r   r    r    r    r!   r~     rd   r~   V瞯<vIh%,=c                    s  | |d\}}	}
t j|t dd}||d d t jf |}t ||	 |j}dt |  }|d d t jf | | }t|t 	||| t|	
 |
|| dd  t fdd||}t  ||	}d	|d k rpd	n|d d
 }t|d ||d d S )NTlong)dtyper   c                 S   s   | d d| d   |  d S )Nr   r   r    rt   r    r    r!   f6  s   zverify_gauss_quad.<locals>.fc                    s    | |  S Nr    r   r   weight_funcr    r!   <lambda>8      z#verify_gauss_quad.<locals>.<lambda>gư>
   r   rtol)r   Zaranger   ZnewaxisdotTr   Zdiagonalr   eyesumr	   quadZvdot)Z	root_funcZ	eval_funcr   r>   bNr   atolrt   wmurl   vZvvZvdZresIZresGr    r   r!   verify_gauss_quad&  s   r   c               	   C   s  dd } dd }dd }t }|| 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
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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|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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
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 || 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#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$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dd || d&d'|d&d'|d&d'd	d
ddd td(d)d)\}}td(\}}t||d*d* t||d*d* td(d+d+\}}td(d,\}}	t||d*d* t||	d*d* tddd-d.\}
}tddd-d/\}}}t|
|d*d* t||d*d* t|dd-d0d\}}t|||d1 tt	tjd2dd tt	tjd3dd tt	tjd-d4d tt	tjd-dd4 tt	tjd-d4d4 d S )5Nc                        fddS )Nc                       t |  |S r   )r=   roots_jacobirl   r   r>   r   r    r!   r   ?  r   z/test_roots_jacobi.<locals>.rf.<locals>.<lambda>r    r   r    r   r!   rf>     ztest_roots_jacobi.<locals>.rfc                    r   )Nc                    r   r   )r=   Zeval_jacobirl   rt   r   r    r!   r   A  r   z/test_roots_jacobi.<locals>.ef.<locals>.<lambda>r    r   r    r   r!   ef@  r   ztest_roots_jacobi.<locals>.efc                    r   )Nc                    s   d|    d|    S )Nr   r    r   r   r    r!   r   C  s    z/test_roots_jacobi.<locals>.wf.<locals>.<lambda>r    r   r    r   r!   wfB  r   ztest_roots_jacobi.<locals>.wfrC               r@   r      -q=r   d   dy=r6   7QE=-=r   vIh%L=?r   vIh%<=g7QU==
ףp=2@L;@g^~z=̌G@r7   g     @   g     @绽|=g      @g     @@g     @@      @皙9Q@r/           +=rW   g      @r   FTr$   r   r   ffffff
@r   )
r   r=   r   roots_legendrer   roots_gegenbauerr	   r   assert_raisesrE   )r   r   r   vgqZxjZwjZxlZwlZxcwcrt   r   yr   mmuImuI_errr    r    r!   test_roots_jacobi=  s   $$((((($(($ ((((((         $  r   c               	   C   s  dd } dd }dd }t }|| 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
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	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|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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dd || dd|dd|ddd	d
dd d || dd|dd|ddd	d
ddd tdd!dd"\}}tdd!dd#\}}}t||d$d$ t||d$d$ t|d!dd%d\}	}
t||	|
d& tttjd%dd tttjd'dd tttjd!dd tttjd!dd( tttjd!d)d( d S )*Nc                    r   )Nc                    r   r   )r=   roots_sh_jacobir   r   r    r!   r     r   z2test_roots_sh_jacobi.<locals>.rf.<locals>.<lambda>r    r   r    r   r!   r     r   z test_roots_sh_jacobi.<locals>.rfc                    r   )Nc                    r   r   )r=   Zeval_sh_jacobir   r   r    r!   r     r   z2test_roots_sh_jacobi.<locals>.ef.<locals>.<lambda>r    r   r    r   r!   r     r   z test_roots_sh_jacobi.<locals>.efc                    r   )Nc                    s   d|     | d   S )Nr@   r    r   r   r    r!   r     s    z2test_roots_sh_jacobi.<locals>.wf.<locals>.<lambda>r    r   r    r   r!   r     r   z test_roots_sh_jacobi.<locals>.wfrC   g      ?r   r@   r   r   r   r   r   r   r6   r   r   r   r   r   r   r   r   r   皙?r   r   r   g(#=r   r   FTr   r   r   r   r$   r   )r   r=   r   r   r	   r   r   rE   )r   r   r   r   rt   r   r   r   r   r   r   r    r    r!   test_roots_sh_jacobi  sX   $$(($(($(($$(((((r   c            
   	   C   s<  t j} t j}tdj}t| ||tj tjd t| ||tj tjddd t| ||tj tjddd t dd\}}t dd\}}}t	||d	d	 t	||d	d	 t
|tj tj\}}	t	|||	d
 t dd\}}t dd\}}}t	||d	d	 t	||d	d	 t	t||d	d	 ttt jd ttt jd d S )Nr   r   r   r   r   r   FTr   r      r   r   )r=   Zroots_hermiteZeval_hermiter   rN   r   r   r   rG   r   r	   r   r   r   rE   
rootfevalfweightfrt   r   r   r   r   r   r   r    r    r!   test_roots_hermite  s&   r   c                     s   dd  d f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 S )Nc                 S   s   t | |jf}t jd t d|d   |dd d f< | dkrftd| |dd d f  |dd d f< td| D ],}td| | ||d d d f  t|d | ||d d d f   ||d d f< q9|S )N      пrC   r   r   r   rV   r@   )r   zerossizepiexpr   rx   )rl   nodesHkr    r    r!   hermite_recursion  s   ((Vz1test_roots_hermite_asy.<locals>.hermite_recursionr   r   c                    sX   t | \}} | d |}t|dd d f t| || tt|ttj|| d S )Nr   r$   )r   Z_roots_hermite_asyr   r   r   r   r   r   )r   r   r   rt   r   r   r   r    r!   test  s    z$test_roots_hermite_asy.<locals>.test   r   r      i,  i-  i  i  i  i  i  i  )r   r   r    )r   r    r   r!   test_roots_hermite_asy  s   
r   c            
   	   C   s   t j} t j}tdj}t| ||tj tjd t| ||tj tjddd t| ||tj tjddd t dd\}}t dd\}}}t	||d	d	 t	||d	d	 t
|tj tj\}}	t	|||	d
 ttt jd ttt jd d S )Nr   r   r   r   r   r   FTr   r   r   r   )r=   Zroots_hermitenormZeval_hermitenormr   rY   r   r   r   rG   r   r	   r   r   rE   r   r    r    r!   test_roots_hermitenorm  s   r   c               	   C   s  dd } dd }dd }t }|| 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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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	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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dd	ddd || d|d|ddd	ddd || dtj|ddd	d
 || dtj|ddd	d || dtj|ddd	ddd td
dd\}}td
dd\}}}t||dd t||dd t|dd d\}	}
t||	|
d! tttjdd tttjd"d tttjd#d$ d S )%Nc                        fddS )Nc                       t |  |S r   )r=   r   r   r>   r    r!   r   #      z6test_roots_gegenbauer.<locals>.rootf.<locals>.<lambda>r    r   r    r   r!   r   "     z$test_roots_gegenbauer.<locals>.rootfc                    r   )Nc                    r   r   )r=   Zeval_gegenbauerr   r   r    r!   r   %  r   z6test_roots_gegenbauer.<locals>.evalf.<locals>.<lambda>r    r   r    r   r!   r   $  r   z$test_roots_gegenbauer.<locals>.evalfc                    r   )Nc                    s   d| d   d  S )Nr   r   r6   r    r   r   r    r!   r   '  rv   z8test_roots_gegenbauer.<locals>.weightf.<locals>.<lambda>r    r   r    r   r!   r   &  r   z&test_roots_gegenbauer.<locals>.weightfr   r   r@   r   r   r   r   r   r   皙?r   r   r   2      g     Pe@g-A=   g     @r   r   r   FTr   r$   r   r   r   r   )	r   r=   eval_chebytr   r   r	   r   r   rE   r   r   r   r   rt   r   r   r   r   r   r   r    r    r!   test_roots_gegenbauer!  sV   """"""""""""""""""""""" r   c               	   C      t dj} ttjtj| ddd ttjtj| ddd ttjtj| ddddd tdd\}}tdd	\}}}t||d
d
 t||d
d
 t	| dd\}}t|||d t
ttjd t
ttjd d S )Nr   r   r@   r   r   r   r   FTr   r$   r   r   r   r   )r   r+   r   r   r=   roots_chebytr   r   r	   r   r   rE   r   rt   r   r   r   r   r   r   r    r    r!   test_roots_chebyta     r   c                  C   sP   t d\} }| d d | dd  }}t||d d d   t| d d d S )N   r   r8   r$   r   )r=   r   r   )rt   r   posnegr    r    r!   test_chebyt_symmetrys  s   r   c                  C      t dj} ttjtj| ddd ttjtj| ddd ttjtj| ddd tdd\}}tdd\}}}t||dd t||dd t	| d	d
\}}t|||d t
ttjd t
ttjd d S )Nr   r   r@   r   r   FTr   r$   r   r   r   r   )r   r0   r   r   r=   Zroots_chebyuZeval_chebyur   r	   r   r   rE   r   r    r    r!   test_roots_chebyuy     r   c               	   C   r   )Nr   r   rV   r   r   r   r   FTr   r   r   r   r   r   )r   r   r   r   r=   Zroots_chebycZeval_chebycr   r	   r   r   rE   r   r    r    r!   test_roots_chebyc  r   r   c                  C   r   )Nr   r   rV   r   r   FTr   r   r   r   r   r   )r   r%   r   r   r=   Zroots_chebysZeval_chebysr   r	   r   r   rE   r   r    r    r!   test_roots_chebys  r   r   c               	   C      t dj} ttjtj| ddd ttjtj| ddd ttjtj| ddddd tdd\}}tdd	\}}}t||d
d
 t||d
d
 t	| dd\}}t|||d t
ttjd t
ttjd d S Nr   r   r@   r   r   r   r   FTr   r   r   r   r   )r   rf   r   r   r=   Zroots_sh_chebytZeval_sh_chebytr   r	   r   r   rE   r   r    r    r!   test_roots_sh_chebyt     r   c               	   C   r   r   )r   ri   r   r   r=   Zroots_sh_chebyuZeval_sh_chebyur   r	   r   r   rE   r   r    r    r!   test_roots_sh_chebyu  r   r   c               	   C   s   t dj} ttjtj| ddd ttjtj| ddddd ttjtj| ddddd tdd	\}}tdd
\}}}t||dd t||dd t	| dd\}}t|||d t
ttjd t
ttjd d S )Nr   r   r@   r   r   r   r   r   FTr   r$   r   r   r   r   )r   ra   r   r   r=   r   Zeval_legendrer   r	   r   r   rE   r   r    r    r!   test_roots_legendre      r   c               	   C   s   t dj} ttjtj| ddd ttjtj| ddddd ttjtj| ddddd tdd	\}}tdd
\}}}t||dd t||dd t	| dd\}}t|||d t
ttjd t
ttjd d S )Nr   r   r@   r   r   r   r   r   FTr   r   r   r   r   )r   r`   r   r   r=   Zroots_sh_legendreZeval_sh_legendrer   r	   r   r   rE   r   r    r    r!   test_roots_sh_legendre  r   r   c               	   C   s   t dj} ttjtj| dtjd ttjtj| dtjddd ttjtj| dtjddd tdd\}}tdd	\}}}t	||d
d
 t	||d
d
 t
| dtj\}}t	|||d tttjd tttjd d S )Nr   r   r   r   r   r   r   FTr   r   r   r   )r   Zlaguerrer   r   r=   Zroots_laguerreZeval_laguerrer   rG   r   r	   r   r   rE   r   r    r    r!   test_roots_laguerre  s    r   c               
   C   s  dd } dd }dd }t }|| d|d|ddtjd	 || d|d|ddtjd
dd || d|d|ddtjddd || d|d|ddtjd	 || d|d|ddtjd
dd || d|d|ddtjddd || d|d|ddtjd	 || d|d|ddtjd
dd || d|d|ddtjddd || d|d|ddtjd	 || d|d|ddtjd
dd || d|d|ddtjddd || d|d|ddtjd	 || d|d|ddtjd
dd || d|d|ddtjdddd td	dd\}}td	dd\}}}t||dd t||dd t|ddtj\}	}
t||	|
d tt	tjdd tt	tjdd tt	tjdd  d S )!Nc                    r   )Nc                    r   r   )r=   roots_genlaguerrer   r   r    r!   r     r   z7test_roots_genlaguerre.<locals>.rootf.<locals>.<lambda>r    r   r    r   r!   r     r   z%test_roots_genlaguerre.<locals>.rootfc                    r   )Nc                    r   r   )r=   Zeval_genlaguerrer   r   r    r!   r     r   z7test_roots_genlaguerre.<locals>.evalf.<locals>.<lambda>r    r   r    r   r!   r     r   z%test_roots_genlaguerre.<locals>.evalfc                    r   )Nc                    s   |   t |   S r   )r   r   r   r   r    r!   r     rv   z9test_roots_genlaguerre.<locals>.weightf.<locals>.<lambda>r    r   r    r   r!   r     r   z'test_roots_genlaguerre.<locals>.weightfrC   r   r   r   r   r   r   r   r   g+F=r   gN<=r   r   r   r   )r   r   r   FTrV   r   r   r   r   g)
r   r   rG   r=   r   r   r	   r   r   rE   r   r    r    r!   test_roots_genlaguerre
  s8    $$ $$ $$ $$ $&r  c                   C   s   t dd d S )NA   r   )r=   r+   r    r    r    r!   test_gh_67214  s   r  )r   r   )1rD   r   r   numpyr   r   r   Znumpy.testingr   r   r   r   Zscipyr	   Zscipy.specialZspecialr=   r   Zscipy.special._orthogonalZ_orthogonalr   r   r5   rK   r]   re   rh   rk   rq   r~   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r    r    r    r!   <module>   sL    <('	
`7@*