a
    hÝEbgD  ã                   @   st  d Z ddlmZ ddlZddlm  mZ ddl	m
Z
 ddlmZmZmZmZ e dg¡d Ze ddg¡d Ze g d¢¡d	 Ze g d
¢¡d Ze g d¢¡d Ze g d¢¡d Ze g d¢¡d ZeeeeeeegZ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!G d$d%„ d%ƒZ"G d&d'„ d'ƒZ#dS )(zTests for laguerre module.

é    )ÚreduceN©Úpolyval)Úassert_almost_equalÚassert_raisesÚassert_equalÚassert_é   éÿÿÿÿ)é   éüÿÿÿr	   r   )é   iîÿÿÿé	   r
   r   )é   i ÿÿÿéH   iðÿÿÿr	   r   )éx   i¨ýÿÿiX  i8ÿÿÿé   r
   r   )éÐ  i ïÿÿi  i öÿÿiÂ  iÜÿÿÿr	   r   c                 C   s   t j| dd�S )Ng�íµ ÷Æ°>)Ztol)ÚlagÚlagtrim©Úx© r   úd/home/ja/django-apps/lartica_env/lib/python3.9/site-packages/numpy/polynomial/tests/test_laguerre.pyÚtrim   s    r   c                   @   s,   e Zd Zdd„ Zdd„ Zdd„ Zdd„ Zd	S )
ÚTestConstantsc                 C   s   t tjddgƒ d S )Nr   r	   )r   r   Z	lagdomain©Úselfr   r   r   Útest_lagdomain   s    zTestConstants.test_lagdomainc                 C   s   t tjdgƒ d S )Nr   )r   r   Zlagzeror   r   r   r   Útest_lagzero!   s    zTestConstants.test_lagzeroc                 C   s   t tjdgƒ d S ©Nr	   )r   r   Zlagoner   r   r   r   Útest_lagone$   s    zTestConstants.test_lagonec                 C   s   t tjddgƒ d S )Nr	   r
   )r   r   Zlagxr   r   r   r   Ú	test_lagx'   s    zTestConstants.test_lagxN)Ú__name__Ú
__module__Ú__qualname__r   r   r!   r"   r   r   r   r   r      s   r   c                   @   sJ   e Zd Ze ddd¡Zdd„ Zdd„ Zdd	„ Zd
d„ Z	dd„ Z
dd„ ZdS )ÚTestArithmeticéýÿÿÿé   éd   c                 C   sž   t dƒD ]�}t dƒD ]‚}d|› d|› �}t t||ƒd ¡}||  d7  < ||  d7  < t dg| dg dg| dg ¡}tt|ƒt|ƒ|d� qqd S ©Né   úAt i=ú, j=r	   r   ©Úerr_msg)ÚrangeÚnpÚzerosÚmaxr   Úlagaddr   r   ©r   ÚiÚjÚmsgÚtgtÚresr   r   r   Útest_lagadd.   s    $zTestArithmetic.test_lagaddc                 C   sž   t dƒD ]�}t dƒD ]‚}d|› d|› �}t t||ƒd ¡}||  d7  < ||  d8  < t dg| dg dg| dg ¡}tt|ƒt|ƒ|d� qqd S r*   )r0   r1   r2   r3   r   Zlagsubr   r   r5   r   r   r   Útest_lagsub8   s    $zTestArithmetic.test_lagsubc                 C   s†   t t dg¡dgƒ t t dg¡ddgƒ tddƒD ]L}dg| dg }dg|d  | d| d |d  g }tt |¡|ƒ q4d S )Nr   r	   r
   r+   r   )r   r   Zlagmulxr0   r   )r   r6   Zserr9   r   r   r   Útest_lagmulxB   s    (zTestArithmetic.test_lagmulxc           
      C   s²   t dƒD ]¤}dg| dg }t | j|¡}t dƒD ]x}d|› d|› �}dg| dg }t | j|¡}t ||¡}t | j|¡}	tt|ƒ|| d k|ƒ t|	|| |d� q2qd S )Nr+   r   r	   r,   r-   r.   )r0   r   Úlagvalr   Úlagmulr   Úlenr   )
r   r6   Zpol1Zval1r7   r8   Zpol2Zval2Zpol3Zval3r   r   r   Útest_lagmulJ   s    zTestArithmetic.test_lagmulc           
      C   s–   t dƒD ]ˆ}t dƒD ]z}d|› d|› �}dg| dg }dg| dg }t ||¡}t ||¡\}}t t ||¡|¡}	tt|	ƒt|ƒ|d� qqd S )Nr+   r,   r-   r   r	   r.   )r0   r   r4   Zlagdivr?   r   r   )
r   r6   r7   r8   ÚciÚcjr9   ZquoÚremr:   r   r   r   Útest_lagdivX   s    zTestArithmetic.test_lagdivc                 C   s|   t dƒD ]n}t dƒD ]`}d|› d|› �}t |d ¡}ttj|g| t dg¡ƒ}t ||¡}tt	|ƒt	|ƒ|d� qqd S )Nr+   r,   r-   r	   r.   )
r0   r1   Úaranger   r   r?   ÚarrayZlagpowr   r   )r   r6   r7   r8   Úcr9   r:   r   r   r   Útest_lagpowc   s    zTestArithmetic.test_lagpowN)r#   r$   r%   r1   Úlinspacer   r;   r<   r=   rA   rE   rI   r   r   r   r   r&   +   s   

r&   c                   @   s‚   e Zd Ze g d¢¡Ze dee¡Ze deee¡Zej	 	d¡d d Z
ee
g d¢ƒZdd	„ Zd
d„ Zdd„ Zdd„ Zdd„ ZdS )ÚTestEvaluation)g      "@g      ,Àg      @úi,j->ijú
i,j,k->ijk©r(   r+   r   r	   )ç      ð?g       @g      @c                    sæ   t t g dg¡jdƒ t dd¡‰ ‡ fdd„tD ƒ}tdƒD ]<}d|› �}|| }t ˆ dg| dg ¡}t|||d� q<td	ƒD ]^}d
g| }t 	|¡‰ t t ˆ dg¡j
|ƒ t t ˆ ddg¡j
|ƒ t t ˆ g d¢¡j
|ƒ q‚d S )Nr	   r   r
   c                    s   g | ]}t ˆ |ƒ‘qS r   r   ©Ú.0rH   r   r   r   Ú
<listcomp>}   ó    z.TestEvaluation.test_lagval.<locals>.<listcomp>é   r,   r.   r(   r   )r	   r   r   )r   r   r>   Úsizer1   rJ   ÚLlistr0   r   r2   Úshape)r   Úyr6   r8   r9   r:   Údimsr   r   r   Útest_lagvalw   s    


zTestEvaluation.test_lagvalc           
      C   s‚   | j \}}}| j\}}}tttj||d d… | jƒ || }t ||| j¡}t||ƒ t 	d¡}	t |	|	| j¡}t
|jdkƒ d S ©Nr   ©r   r(   )r   rX   r   Ú
ValueErrorr   Úlagval2dÚc2dr   r1   Úonesr   rW   ©
r   Úx1Úx2Úx3Úy1Úy2Zy3r9   r:   Úzr   r   r   Útest_lagval2dŒ   s    

zTestEvaluation.test_lagval2dc           
      C   sŒ   | j \}}}| j\}}}tttj|||d d… | jƒ || | }t |||| j¡}t||ƒ t 	d¡}	t |	|	|	| j¡}t
|jdkƒ d S r[   )r   rX   r   r]   r   Úlagval3dÚc3dr   r1   r`   r   rW   ra   r   r   r   Útest_lagval3d�   s    

zTestEvaluation.test_lagval3dc           
      C   sl   | j \}}}| j\}}}t d||¡}t ||| j¡}t||ƒ t d¡}	t |	|	| j¡}t	|j
dkƒ d S )NrL   r\   )r   r(   r   r(   )r   rX   r1   Úeinsumr   Z	laggrid2dr_   r   r`   r   rW   ra   r   r   r   Útest_laggrid2d®   s    

zTestEvaluation.test_laggrid2dc           
      C   sr   | j \}}}| j\}}}t d|||¡}t |||| j¡}t||ƒ t d¡}	t |	|	|	| j¡}t	|j
dkƒ d S )NrM   r\   )r   r(   r   r(   r   r(   )r   rX   r1   rl   r   Z	laggrid3drj   r   r`   r   rW   ra   r   r   r   Útest_laggrid3d¼   s    

zTestEvaluation.test_laggrid3dN)r#   r$   r%   r1   rG   Zc1drl   r_   rj   Úrandomr   r   rX   rZ   rh   rk   rm   rn   r   r   r   r   rK   m   s   rK   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚTestIntegralc           
   	   C   s2  t ttjdgdƒ t ttjdgdƒ t ttjdgdddgƒ t ttjdgdgd� t ttjdgdgd� t ttjdgdd� tdd	ƒD ]8}dg|d  dg }tjdg||d
�}t|ddgƒ q†td	ƒD ]n}|d }dg| dg }|gdg|  d| g }t |¡}tj|d|gd
�}t |¡}tt	|ƒt	|ƒƒ qÈtd	ƒD ]N}|d }dg| dg }t |¡}tj|d|gdd�}tt 
d|¡|ƒ �q@td	ƒD ]r}|d }dg| dg }|gdg|  d| g }t |¡}tj|d|gdd�}t |¡}tt	|ƒt	|ƒƒ �q˜td	ƒD ]r}tdd	ƒD ]`}	dg| dg }|d d … }t|	ƒD ]}tj|dd�}�qJtj||	d�}tt	|ƒt	|ƒƒ �q"�qtd	ƒD ]€}tdd	ƒD ]n}	dg| dg }|d d … }t|	ƒD ]}tj|d|gd
�}�qÆtj||	tt|	ƒƒd
�}tt	|ƒt	|ƒƒ �qž�q�td	ƒD ]„}tdd	ƒD ]r}	dg| dg }|d d … }t|	ƒD ]}tj|d|gdd�}�qPtj||	tt|	ƒƒdd�}tt	|ƒt	|ƒƒ �q(�qtd	ƒD ]„}tdd	ƒD ]r}	dg| dg }|d d … }t|	ƒD ]}tj|d|gdd�}�qÞtj||	tt|	ƒƒdd�}tt	|ƒt	|ƒƒ �q¶�q¨d S )Nr   ç      à?r
   r	   )Úlbnd)Úscl©Úaxisr   r+   )ÚmÚk)rv   rw   rr   )rv   rw   rs   ©rv   )r   Ú	TypeErrorr   Úlagintr]   r0   r   Úpoly2lagÚlag2polyr   r>   Úlist)
r   r6   rw   r:   rs   Úpolr9   Zlagpolrz   r7   r   r   r   Útest_lagintÍ   s€    




zTestIntegral.test_lagintc                 C   sš   t j d¡}t  dd„ |jD ƒ¡j}tj|dd�}t||ƒ t  dd„ |D ƒ¡}tj|dd�}t||ƒ t  dd„ |D ƒ¡}tj|d	dd
�}t||ƒ d S )N©r(   é   c                 S   s   g | ]}t  |¡‘qS r   ©r   rz   rP   r   r   r   rR   $  rS   z1TestIntegral.test_lagint_axis.<locals>.<listcomp>r   rt   c                 S   s   g | ]}t  |¡‘qS r   r‚   rP   r   r   r   rR   (  rS   r	   c                 S   s   g | ]}t j|d d�‘qS )r(   )rw   r‚   rP   r   r   r   rR   ,  rS   r(   )rw   ru   )r1   ro   ÚvstackÚTr   rz   r   ©r   r_   r9   r:   r   r   r   Útest_lagint_axis   s    

zTestIntegral.test_lagint_axisN)r#   r$   r%   r   r†   r   r   r   r   rp   Ë   s   Srp   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚTestDerivativec                 C   s  t ttjdgdƒ t ttjdgdƒ tdƒD ]4}dg| dg }tj|dd�}tt|ƒt|ƒƒ q,tdƒD ]N}tddƒD ]>}dg| dg }tjtj||d�|d�}t	t|ƒt|ƒƒ qxqjtdƒD ]R}tddƒD ]B}dg| dg }tjtj||dd�|dd�}t	t|ƒt|ƒƒ qÐqÂd S )	Nr   rq   r
   r+   r	   rx   r   )rv   rs   )
r   ry   r   Úlagderr]   r0   r   r   rz   r   )r   r6   r9   r:   r7   r   r   r   Útest_lagder3  s     zTestDerivative.test_lagderc                 C   sl   t j d¡}t  dd„ |jD ƒ¡j}tj|dd�}t||ƒ t  dd„ |D ƒ¡}tj|dd�}t||ƒ d S )Nr€   c                 S   s   g | ]}t  |¡‘qS r   ©r   rˆ   rP   r   r   r   rR   P  rS   z3TestDerivative.test_lagder_axis.<locals>.<listcomp>r   rt   c                 S   s   g | ]}t  |¡‘qS r   rŠ   rP   r   r   r   rR   T  rS   r	   )r1   ro   rƒ   r„   r   rˆ   r   r…   r   r   r   Útest_lagder_axisL  s    
zTestDerivative.test_lagder_axisN)r#   r$   r%   r‰   r‹   r   r   r   r   r‡   1  s   r‡   c                   @   s8   e Zd Zej d¡d d Zdd„ Zdd„ Zdd	„ Zd
S )Ú
TestVanderrN   r   r	   c                 C   sÎ   t  d¡}t |d¡}t|jdkƒ tdƒD ].}dg| dg }t|d|f t ||¡ƒ q,t  	ddgddgdd	gg¡}t |d¡}t|jd
kƒ tdƒD ].}dg| dg }t|d|f t ||¡ƒ qšd S )Nr(   r€   r�   r   r	   .r   r+   r   )r(   r   r�   )
r1   rF   r   Ú	lagvanderr   rW   r0   r   r>   rG   )r   r   Úvr6   Úcoefr   r   r   Útest_lagvander]  s    
zTestVander.test_lagvanderc                 C   sx   | j \}}}tj d¡}t ||ddg¡}t |||¡}t ||j¡}t||ƒ t |g|gddg¡}t	|j
dkƒ d S )Nr\   r	   r   )r	   r+   r   )r   r1   ro   r   Zlagvander2dr^   ÚdotÚflatr   r   rW   ©r   rb   rc   rd   rH   Zvanr9   r:   r   r   r   Útest_lagvander2dn  s    
zTestVander.test_lagvander2dc                 C   s€   | j \}}}tj d¡}t |||g d¢¡}t ||||¡}t ||j¡}t||ƒ t |g|g|gg d¢¡}t	|j
dkƒ d S )N)r   r(   r�   )r	   r   r(   )r	   r+   r   )r   r1   ro   r   Zlagvander3dri   r‘   r’   r   r   rW   r“   r   r   r   Útest_lagvander3d{  s    
zTestVander.test_lagvander3dN)	r#   r$   r%   r1   ro   r   r�   r”   r•   r   r   r   r   rŒ   Y  s   rŒ   c                   @   s   e Zd Zdd„ ZdS )ÚTestFittingc                 C   sf  dd„ }t ttjdgdgdƒ t ttjdggdgdƒ t ttjg dgdƒ t ttjdgdgggdƒ t ttjddgdgdƒ t ttjdgddgdƒ t ttjdgdgddggd� t ttjdgdgdddgd� t ttjdgdgdgƒ t ttjdgdgg d¢ƒ t ttjdgdgg ƒ t dd¡}||ƒ}t ||d	¡}tt|ƒd
ƒ t	t 
||¡|ƒ t ||g d¢¡}tt|ƒd
ƒ t	t 
||¡|ƒ t ||d
¡}tt|ƒdƒ t	t 
||¡|ƒ t ||g d¢¡}tt|ƒdƒ t	t 
||¡|ƒ t |t ||g¡jd	¡}t	|t ||g¡jƒ t |t ||g¡jg d¢¡}t	|t ||g¡jƒ t |¡}| ¡ }d|dd d…< d|dd d…< tj||d	|d�}	t	|	|ƒ tj||g d¢|d�}	t	|	|ƒ tj|t ||g¡jd	|d�}
t	|
t ||g¡jƒ tj|t ||g¡jg d¢|d�}
t	|
t ||g¡jƒ g d¢}t	t ||d¡ddgƒ t	t ||ddg¡ddgƒ d S )Nc                 S   s   | | d  | d  S )Nr	   r   r   r   r   r   r   ÚfŒ  s    z"TestFitting.test_lagfit.<locals>.fr	   r
   r   r   )Úw)r   r
   r   r(   r�   )r   r	   r   r(   r+   )r   r	   r   r(   r�   )r	   y              ð?r
   y       €      ð¿)r   r]   r   Zlagfitry   r1   rJ   r   r@   r   r>   rG   r„   Z
zeros_likeÚcopy)r   r—   r   rX   Zcoef3Zcoef4Zcoef2dr˜   ZywZwcoef3Zwcoef2dr   r   r   Útest_lagfit‹  sZ    


"zTestFitting.test_lagfitN)r#   r$   r%   rš   r   r   r   r   r–   ‰  s   r–   c                   @   s$   e Zd Zdd„ Zdd„ Zdd„ ZdS )ÚTestCompanionc                 C   s"   t ttjg ƒ t ttjdgƒ d S r    )r   r]   r   Úlagcompanionr   r   r   r   Útest_raisesÉ  s    zTestCompanion.test_raisesc                 C   s<   t ddƒD ],}dg| dg }tt |¡j||fkƒ q
d S )Nr	   r+   r   )r0   r   r   rœ   rW   )r   r6   r�   r   r   r   Útest_dimensionsÍ  s    zTestCompanion.test_dimensionsc                 C   s   t t ddg¡d dkƒ d S )Nr	   r   )r   r   g      ø?)r   r   rœ   r   r   r   r   Útest_linear_rootÒ  s    zTestCompanion.test_linear_rootN)r#   r$   r%   r�   rž   rŸ   r   r   r   r   r›   Ç  s   r›   c                   @   s   e Zd Zdd„ ZdS )Ú	TestGaussc                 C   s|   t  d¡\}}t  |d¡}t |j| |¡}dt | ¡ ¡ }|d d …d f | | }t|t 	d¡ƒ d}t| 
¡ |ƒ d S )Nr)   éc   r	   rO   )r   Zlaggaussr�   r1   r‘   r„   ÚsqrtZdiagonalr   ÚeyeÚsum)r   r   r˜   rŽ   ÚvvZvdr9   r   r   r   Útest_100Ø  s    zTestGauss.test_100N)r#   r$   r%   r¦   r   r   r   r   r    Ö  s   r    c                   @   sD   e Zd Zdd„ Zdd„ Zdd„ Zdd„ Zd	d
„ Zdd„ Zdd„ Z	dS )ÚTestMiscc              	   C   s¤   t  g ¡}tt|ƒdgƒ tddƒD ]z}t t tj dd| d ¡dd d… ¡}t  |¡}t  	||¡}d}t
t|ƒ|d kƒ tt  |¡d dƒ t||ƒ q$d S )Nr	   r+   r   r   r
   )r   Úlagfromrootsr   r   r0   r1   ÚcosrJ   Úpir>   r   r@   r|   )r   r:   r6   Úrootsr~   r9   r   r   r   Útest_lagfromrootsë  s    
*
zTestMisc.test_lagfromrootsc                 C   sl   t t dg¡g ƒ t t ddg¡dgƒ tddƒD ]4}t dd|¡}t t |¡¡}t t|ƒt|ƒƒ q2d S )Nr	   r   r   r+   r(   )r   r   Zlagrootsr0   r1   rJ   r¨   r   )r   r6   r9   r:   r   r   r   Útest_lagroots÷  s    zTestMisc.test_lagrootsc                 C   sb   g d¢}t ttj|dƒ tt |¡|d d… ƒ tt |d¡|d d… ƒ tt |d¡dgƒ d S )N)r   r
   r	   r   r
   r	   r'   r   r   )r   r]   r   r   r   )r   r�   r   r   r   Útest_lagtrimÿ  s
    zTestMisc.test_lagtrimc                 C   s   t t dd¡ddgƒ d S )Nr(   r�   rT   r   )r   r   Zlagliner   r   r   r   Útest_lagline
  s    zTestMisc.test_laglinec                 C   s2   t dƒD ]$}tt dg| dg ¡t| ƒ qd S ©NrT   r   r	   )r0   r   r   r|   rV   ©r   r6   r   r   r   Útest_lag2poly  s    zTestMisc.test_lag2polyc                 C   s2   t dƒD ]$}tt t| ¡dg| dg ƒ qd S r°   )r0   r   r   r{   rV   r±   r   r   r   Útest_poly2lag  s    zTestMisc.test_poly2lagc                 C   s2   t  ddd¡}t  | ¡}t |¡}t||ƒ d S )Nr   é
   é   )r1   rJ   Úexpr   Z	lagweightr   )r   r   r9   r:   r   r   r   Útest_weight  s    
zTestMisc.test_weightN)
r#   r$   r%   r¬   r­   r®   r¯   r²   r³   r·   r   r   r   r   r§   é  s   r§   )$Ú__doc__Ú	functoolsr   Únumpyr1   Znumpy.polynomial.laguerreZ
polynomialZlaguerrer   Znumpy.polynomial.polynomialr   Znumpy.testingr   r   r   r   rG   ZL0ZL1ZL2ZL3ZL4ZL5ZL6rV   r   r   r&   rK   rp   r‡   rŒ   r–   r›   r    r§   r   r   r   r   Ú<module>   s0   B^f(0>