Edit File: test_laguerre.cpython-37.pyc
B ��FdZC � @ s� d Z ddlmZmZmZ ddlZddlm m Z ddl m Z ddl mZmZmZmZmZmZ e�dg�d Ze�ddg�d Ze�dd dg�d Ze�d d d dg�d Ze�d ddddg�d Ze�ddddddg�d Ze�dddddddg�d ZeeeeeeegZ d d � Z G d d � d e�Z G d d!� d!e�Z G d"d#� d#e�Z G d$d%� d%e�Z!G d&d'� d'e�Z"G d(d)� d)e�Z#G d*d+� d+e�Z$G d,d-� d-e�Z%G d.d/� d/e�Z&G d0d1� d1e�Z'e(d2k�r�e� dS )3z Tests for laguerre module. � )�division�absolute_import�print_functionN)�polyval)�TestCase�assert_almost_equal� assert_raises� assert_equal�assert_�run_module_suite� ���� ���� i��� � i�����H i��x i����iX i8���� i� i �i i����i� i���c C s t j| dd�S )Ng���ư>)�tol)�lag�lagtrim)�x� r ��/builddir/build/BUILDROOT/alt-python37-numpy-1.13.3-8.el8.x86_64/opt/alt/python37/lib64/python3.7/site-packages/numpy/polynomial/tests/test_laguerre.py�trim s r c @ s, e Zd Zdd� Zdd� Zdd� Zdd� Zd S ) � TestConstantsc C s t tjddg� d S )Nr r )r r � lagdomain)�selfr r r �test_lagdomain s z TestConstants.test_lagdomainc C s t tjdg� d S )Nr )r r �lagzero)r r r r � test_lagzero! s zTestConstants.test_lagzeroc C s t tjdg� d S )Nr )r r �lagone)r r r r � test_lagone$ s zTestConstants.test_lagonec C s t tjddg� d S )Nr r )r r �lagx)r r r r � test_lagx' s zTestConstants.test_lagxN)�__name__� __module__� __qualname__r r" r$ r& r r r r r s r c @ sB e Zd Ze�ddd�Zdd� Zdd� Zdd � Zd d � Z d d � Z dS )�TestArithmetic���� �d c C s� x�t d�D ]�}x�t d�D ]~}d||f }t�t||�d �}|| d7 < || d7 < t�dg| dg dg| dg �}tt|�t|�|d� qW q W d S )N� z At i=%d, j=%dr r )�err_msg)�range�np�zeros�maxr �lagaddr r )r �i�j�msg�tgt�resr r r � test_lagadd. s $zTestArithmetic.test_lagaddc C s� x�t d�D ]�}x�t d�D ]~}d||f }t�t||�d �}|| d7 < || d8 < t�dg| dg dg| dg �}tt|�t|�|d� qW q W d S )Nr. z At i=%d, j=%dr r )r/ )r0 r1 r2 r3 r �lagsubr r )r r5 r6 r7 r8 r9 r r r � test_lagsub8 s $zTestArithmetic.test_lagsubc C s� t t�dg�dg� t t�dg�ddg� xZtdd�D ]L}dg| dg }dg|d | d| d |d g }tt�|�|� q6W d S )Nr r r r. r )r r �lagmulxr0 r )r r5 �serr8 r r r � test_lagmulxB s (zTestArithmetic.test_lagmulxc C s� x�t d�D ]�}dg| dg }t�| j|�}x�t d�D ]t}d||f }dg| dg }t�| j|�}t�||�}t�| j|�} tt|�|| d k|� t| || |d� q6W q W d S )Nr. r r z At i=%d, j=%d)r/ )r0 r �lagvalr �lagmulr �lenr ) r r5 �pol1�val1r6 r7 �pol2�val2�pol3�val3r r r � test_lagmulJ s zTestArithmetic.test_lagmulc C s� x�t d�D ]�}x�t d�D ]v}d||f }dg| dg }dg| dg }t�||�}t�||�\}}t�t�||�|�} tt| �t|�|d� qW q W d S )Nr. z At i=%d, j=%dr r )r/ )r0 r r4 �lagdivrA r r ) r r5 r6 r7 �ci�cjr8 �quo�remr9 r r r � test_lagdivX s zTestArithmetic.test_lagdivN) r' r( r) r1 �linspacer r: r< r? rI rO r r r r r* + s r* c @ s� e Zd Ze�dddg�Ze�dee�Ze�deee�Zej � d�d d Z e e d d d g�Z d d � Z dd� Zdd� Zdd� Zdd� ZdS )�TestEvaluationg "@g ,�g @zi,j->ijz i,j,k->ijk)r, r. r r g �?g @g @c s� t t�g dg�jd� t�dd�� � fdd�tD �}xFtd�D ]:}d| }|| }t�� dg| dg �}t|||d� q>W xltd �D ]`}d g| }t� |�� t t�� dg�j |� t t�� ddg�j |� t t�� dddg�j |� q�W d S ) Nr r r c s g | ]}t � |��qS r )r )�.0�c)r r r � <listcomp>t s z.TestEvaluation.test_lagval.<locals>.<listcomp>� zAt i=%d)r/ r, r ) r r r@ �sizer1 rP �Llistr0 r r2 �shape)r �yr5 r7 r8 r9 �dimsr )r r � test_lagvaln s zTestEvaluation.test_lagvalc C s� | j \}}}| j\}}}tttj||d d� | j� || }t�||| j�}t||� t� d�} t�| | | j�}t |j dk� d S )Nr )r r, ) r rY r � ValueErrorr �lagval2d�c2dr r1 �onesr rX ) r �x1�x2�x3�y1�y2�y3r8 r9 �zr r r � test_lagval2d� s z TestEvaluation.test_lagval2dc C s� | j \}}}| j\}}}tttj|||d d� | j� || | }t�|||| j�}t||� t� d�} t�| | | | j�}t |j dk� d S )Nr )r r, ) r rY r r\ r �lagval3d�c3dr r1 r_ r rX ) r r` ra rb rc rd re r8 r9 rf 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 dk� d S )Nzi,j->ij)r r, )r r, r r, ) r rY r1 �einsumr � laggrid2dr^ r r_ r rX ) r r` ra rb rc rd re r8 r9 rf 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 dk� d S )Nz i,j,k->ijk)r r, )r r, r r, r r, ) r rY r1 rk r � laggrid3dri r r_ r rX ) r r` ra rb rc rd re r8 r9 rf r r r �test_laggrid3d� s z TestEvaluation.test_laggrid3dN)r' r( r) r1 �array�c1drk r^ ri �randomr r rY r[ rg rj rm ro r r r r rQ d s rQ c @ s e Zd Zdd� Zdd� ZdS )� TestIntegralc C s2 t ttjdgd� t ttjdgd� t ttjdgdddg� xFtdd�D ]8}dg|d dg }tjdg||d�}t|ddg� qHW xztd�D ]n}|d }dg| dg }|gdg| d| g }t�|�}tj|d|gd�}t�|�}tt|�t|�� q�W xZtd�D ]N}|d }dg| dg }t�|�}tj|d|gdd�}tt� d|�|� �q W x~td�D ]r}|d }dg| dg }|gdg| d| g }t�|�}tj|d|gdd �}t�|�}tt|�t|�� �qfW x�td�D ]z}xrtdd�D ]d} dg| dg }|d d � }x t| �D ]}tj|dd �}�q W tj|| d �}tt|�t|�� �q�W �q�W x�td�D ]�}x�tdd�D ]r} dg| dg }|d d � }x$t| �D ]}tj|d|gd�}�q�W tj|| t t| ��d�}tt|�t|�� �q~W �qnW x�td�D ]�}x�tdd�D ]v} dg| dg }|d d � }x&t| �D ]}tj|d|gdd�}�q>W tj|| t t| ��dd�}tt|�t|�� �qW �qW x�td�D ]�}x�tdd�D ]v} dg| dg }|d d � }x&t| �D ]}tj|d|gdd �}�q�W tj|| t t| ��dd �}tt|�t|�� �q�W �q�W d S ) Nr g �?r r r r. )�m�k)rt ru �lbnd)rt ru �scl)rt ) r r\ r �lagintr0 r �poly2lag�lag2polyr r@ �list) r r5 ru r9 rw �polr8 Zlagpolrx r6 r r r � test_lagint� sz zTestIntegral.test_lagintc C s� t j�d�}t �dd� |jD ��j}tj|dd�}t||� t �dd� |D ��}tj|dd�}t||� t �dd� |D ��}tj|d dd �}t||� d S ) N)r, � c S s g | ]}t �|��qS r )r rx )rR rS r r r rT s z1TestIntegral.test_lagint_axis.<locals>.<listcomp>r )�axisc S s g | ]}t �|��qS r )r rx )rR rS r r r rT s r c S s g | ]}t j|d d��qS )r, )ru )r rx )rR rS r r r rT s r, )ru r )r1 rr �vstack�Tr rx r )r r^ r8 r9 r r r �test_lagint_axis s z TestIntegral.test_lagint_axisN)r' r( r) r} r� r r r r rs � s Prs c @ s e Zd Zdd� Zdd� ZdS )�TestDerivativec C s. t ttjdgd� t ttjdgd� x@td�D ]4}dg| dg }tj|dd�}tt|�t|�� q.W x^td�D ]R}xLtdd�D ]>}dg| dg }tjtj||d�|d�}tt|�t|�� q�W qpW xbtd�D ]V}xPtdd�D ]B}dg| dg }tjtj||dd�|dd�}tt|�t|�� q�W q�W d S ) Nr g �?r r. r )rt r )rt rw ) r r\ r �lagderr0 r r rx r )r r5 r8 r9 r6 r r r � test_lagder' s zTestDerivative.test_lagderc C sl t j�d�}t �dd� |jD ��j}tj|dd�}t||� t �dd� |D ��}tj|dd�}t||� d S )N)r, r~ c S s g | ]}t �|��qS r )r r� )rR rS r r r rT D s z3TestDerivative.test_lagder_axis.<locals>.<listcomp>r )r c S s g | ]}t �|��qS r )r r� )rR rS r r r rT H s r )r1 rr r� r� r r� r )r r^ r8 r9 r r r �test_lagder_axis@ s z TestDerivative.test_lagder_axisN)r' r( r) r� r� r r r r r� % s r� c @ s8 e Zd Zej�d�d d Zdd� Zdd� Zdd � Zd S ) � TestVander)r, r. r r c C s� t �d�}t�|d�}t|jdk� x:td�D ].}dg| dg }t|d|f t�||�� q.W t � ddgddgdd gg�}t�|d�}t|jd k� x:td�D ].}dg| dg }t|d|f t�||�� q�W d S ) Nr, )r, r~ r~ r r .r r. r )r, r r~ ) r1 �aranger � lagvanderr rX r0 r r@ rp )r r �vr5 �coefr r r �test_lagvanderQ s zTestVander.test_lagvanderc C sx | j \}}}tj�d�}t�||ddg�}t�|||�}t�||j�}t||� t�|g|gddg�}t |j dk� d S )N)r r, r r )r r. r ) r r1 rr r � lagvander2dr] �dot�flatr r rX )r r` ra rb rS �vanr8 r9 r r r �test_lagvander2db s zTestVander.test_lagvander2dc C s� | j \}}}tj�d�}t�|||dddg�}t�||||�}t�||j�}t||� t�|g|g|gdddg�}t |j dk� d S )N)r r, r~ r r r, )r r. r ) r r1 rr r � lagvander3drh r� r� r r rX )r r` ra rb rS r� r8 r9 r r r �test_lagvander3do s zTestVander.test_lagvander3dN) r' r( r) r1 rr r r� r� r� r r r r r� M s r� c @ s e Zd Zdd� ZdS )� TestFittingc C s� dd� }t ttjdgdgd� t ttjdggdgd� t ttjg dgd� t ttjdgdgggd� t ttjddgdgd� t ttjdgddgd� t ttjdgdgddggd� t ttjdgdgdddgd� t ttjdgdgdg� t ttjdgdgdddg� t ttjdgdgg � t�dd�}||�}t�||d �}tt|�d � t t� ||�|� t�||dddd g�}tt|�d � t t� ||�|� t�||d �}tt|�d � t t� ||�|� t�||dddd d g�}tt|�d � t t� ||�|� t�|t� ||g�j d �}t |t� ||g�j � t�|t� ||g�j dddd g�}t |t� ||g�j � t� |�}|�� }d|dd d�< d|dd d�< tj||d |d�} t | |� tj||dddd g|d�} t | |� tj|t� ||g�j d |d�} t | t� ||g�j � tj|t� ||g�j dddd g|d�} t | t� ||g�j � dd dd g}t t�||d�ddg� t t�||ddg�ddg� 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 )�wr r, r~ r. y �?y � �)r r\ r �lagfit� TypeErrorr1 rP r rB r r@ rp r� � zeros_like�copy) r r� r rY �coef3�coef4�coef2dr� �yw�wcoef3�wcoef2dr r r � test_lagfit sZ " & zTestFitting.test_lagfitN)r' r( r) r� r r r r r� } s r� c @ s$ e Zd Zdd� Zdd� Zdd� ZdS )� TestCompanionc C s" t ttjg � t ttjdg� d S )Nr )r r\ r � lagcompanion)r r r r � test_raises� s zTestCompanion.test_raisesc C s@ x:t dd�D ],}dg| dg }tt�|�j||fk� q W d S )Nr r. r )r0 r r r� rX )r r5 r� r r r �test_dimensions� s z TestCompanion.test_dimensionsc C s t t�ddg�d dk� 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 Zd Zdd� ZdS )� TestGaussc C s| t �d�\}}t �|d�}t�|j| |�}dt�|�� � }|d d �d f | | }t|t� d�� d}t|� � |� d S )Nr- �c r g �?) r �laggaussr� r1 r� r� �sqrt�diagonalr �eye�sum)r r r� r� �vv�vdr8 r r r �test_100� s zTestGauss.test_100N)r' r( r) r� r r r r r� � s r� c @ sD e Zd Zdd� Zdd� Zdd� Zdd� Zd d � Zd d � Zd d� Z dS )�TestMiscc C s� t �g �}tt|�dg� x�tdd�D ]z}t�t�tj dd| d �dd d� �}t �|�}t � ||�}d}t t |�|d k� tt � |�d d� t||� q&W d S )Nr r. r r r ) r � lagfromrootsr r r0 r1 �cosrP �pir@ r rB rz )r r9 r5 �rootsr| r8 r r r �test_lagfromroots� s * zTestMisc.test_lagfromrootsc C sp t t�dg�g � t t�ddg�dg� xBtdd�D ]4}t�dd|�}t�t�|��}t t|�t|�� q4W d S )Nr r r r. r, )r r �lagrootsr0 r1 rP r� r )r r5 r8 r9 r r r � test_lagroots� s zTestMisc.test_lagrootsc C sf ddddg}t ttj|d� tt�|�|d d� � tt�|d�|d d� � tt�|d�dg� d S )Nr r r r r+ )r r\ r r r )r r� r r r � test_lagtrim� s zTestMisc.test_lagtrimc C s t t�dd�ddg� d S )Nr, r~ rU r )r r �lagline)r r r r � test_lagline� s zTestMisc.test_laglinec C s6 x0t d�D ]$}tt�dg| dg �t| � q W d S )NrU r r )r0 r r rz rW )r r5 r r r � test_lag2poly s zTestMisc.test_lag2polyc C s6 x0t d�D ]$}tt�t| �dg| dg � q W d S )NrU r r )r0 r r ry rW )r r5 r r r � test_poly2lag s zTestMisc.test_poly2lagc C s2 t �ddd�}t �| �}t�|�}t||� d S )Nr � � )r1 rP �expr � lagweightr )r r r8 r9 r r r � test_weight s zTestMisc.test_weightN) r' r( r) r� r� r� r� r� r� r� r r r r r� � s r� �__main__))�__doc__� __future__r r r �numpyr1 Znumpy.polynomial.laguerre� polynomial�laguerrer Znumpy.polynomial.polynomialr Z numpy.testingr r r r r r rp ZL0ZL1ZL2ZL3�L4ZL5ZL6rW r r r* rQ rs r� r� r� r� r� r� r' r r r r �<module> s4 9^c(0>3