Edit File: test_legendre.cpython-37.pyc
B ��FdVF � @ s$ d Z ddlmZmZmZ ddlZddlm m Z ddl m Z ddl mZmZmZmZmZmZ e�dg�Ze�ddg�Ze�dddg�d Ze�dd dd g�d Ze�ddd dd g�d Ze�ddddddg�d Ze�dddddddg�d Ze�ddddddddg�d Z e�d dddd dd dd g �d Z e�dddd dd!dd"dd#g �d Z eeeeeeee e e g Z 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 d0d1� d1e�Z&G d2d3� d3e�Z'G d4d5� d5e�Z(G d6d7� d7e�Z)G d8d9� d9e�Z*e+d:k�r e� dS );z Tests for legendre module. � )�division�absolute_import�print_functionN)�polyval)�TestCase�assert_almost_equal� assert_raises� assert_equal�assert_�run_module_suite� ���� � ���� i���# � � i�����? ����i i����� � i���i; iK���i� i���i i��i# � i��ibF it���i{/ c C s t j| dd�S )Ng���ư>)�tol)�leg�legtrim)�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_legendre.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 � legdomain)�selfr r r �test_legdomain! s z TestConstants.test_legdomainc C s t tjdg� d S )Nr )r r �legzero)r$ r r r � test_legzero$ s zTestConstants.test_legzeroc C s t tjdg� d S )Nr )r r �legone)r$ r r r � test_legone' s zTestConstants.test_legonec C s t tjddg� d S )Nr r )r r �legx)r$ r r r � test_legx* s zTestConstants.test_legxN)�__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 )�TestArithmeticr r �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 )Nr z At i=%d, j=%dr r )�err_msg)�range�np�zeros�maxr �legaddr r! )r$ �i�j�msg�tgt�resr r r � test_legadd1 s $zTestArithmetic.test_legaddc 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 )r1 )r2 r3 r4 r5 r �legsubr r! )r$ r7 r8 r9 r: r; r r r � test_legsub; s $zTestArithmetic.test_legsubc C s� t t�dg�dg� t t�dg�ddg� xbtdd�D ]T}d| d }dg| dg }dg|d || d|d | g }t t�|�|� q6W d S )Nr r r r )r r �legmulxr2 )r$ r7 �tmp�serr: r r r � test_legmulxE s $zTestArithmetic.test_legmulxc 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)r1 )r2 r �legvalr �legmulr �lenr ) r$ r7 �pol1�val1r8 r9 �pol2�val2�pol3�val3r r r � test_legmulN s zTestArithmetic.test_legmulc 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 )r1 )r2 r r6 �legdivrD r r! ) r$ r7 r8 r9 �ci�cjr: �quo�remr; r r r � test_legdiv\ s zTestArithmetic.test_legdivN) r, r- r. r3 �linspacer r< r> rB rL rR r r r r r/ . s r/ c @ s� e Zd Ze�dddg�Ze�dee�Ze�deee�Zej � d�d d Z e e dddg�Z d d � Z d d � Zd d� Zdd� Zdd� ZdS )�TestEvaluationg @zi,j->ijz i,j,k->ijk)r r r r 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>x s z.TestEvaluation.test_legval.<locals>.<listcomp>� zAt i=%d)r1 r r ) r r rC �sizer3 rS �Llistr2 r r4 �shape)r$ �yr7 r9 r: r; �dimsr )r r � test_legvalr s zTestEvaluation.test_legvalc C s� | j \}}}| j\}}}tttj||d d� | j� || }t�||| j�}t||� t� d�} t�| | | j�}t |j dk� d S )Nr )r r ) r r\ r � ValueErrorr �legval2d�c2dr r3 �onesr r[ ) r$ �x1�x2�x3�y1�y2�y3r: r; �zr r r � test_legval2d� s z TestEvaluation.test_legval2dc C s� | j \}}}| j\}}}tttj|||d d� | j� || | }t�|||| j�}t||� t� d�} t�| | | | j�}t |j dk� d S )Nr )r r ) r r\ r r_ r �legval3d�c3dr r3 rb r r[ ) r$ rc rd re rf rg rh r: r; ri r r r � test_legval3d� s z TestEvaluation.test_legval3dc 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 r\ r3 �einsumr � leggrid2dra r rb r r[ ) r$ rc rd re rf rg rh r: r; ri r r r �test_leggrid2d� s z TestEvaluation.test_leggrid2dc 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 r\ r3 rn r � leggrid3drl r rb r r[ ) r$ rc rd re rf rg rh r: r; ri r r r �test_leggrid3d� s z TestEvaluation.test_leggrid3dN)r, r- r. r3 �array�c1drn ra rl �randomr r r\ r^ rj rm rp rr r r r r rT h s rT 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)rw rx �lbnd)rw rx �scl)rw ) r r_ r �legintr2 r �poly2leg�leg2polyr! rC �list) r$ r7 rx r; rz �polr: Zlegpolr{ r8 r r r � test_legint� sz zTestIntegral.test_legintc 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 r{ )rU rV r r r rW s z1TestIntegral.test_legint_axis.<locals>.<listcomp>r )�axisc S s g | ]}t �|��qS r )r r{ )rU rV r r r rW s r c S s g | ]}t j|d d��qS )r )rx )r r{ )rU rV r r r rW $ s r )rx r� )r3 ru �vstack�Tr r{ r )r$ ra r: r; r r r �test_legint_axis s z TestIntegral.test_legint_axisN)r, r- r. r� r� r r r r rv � s Prv 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 )rw r )rw rz ) r r_ r �legderr2 r r! r{ r )r$ r7 r: r; r8 r r r � test_legder+ s zTestDerivative.test_legderc 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� )rU rV r r r rW H s z3TestDerivative.test_legder_axis.<locals>.<listcomp>r )r� c S s g | ]}t �|��qS r )r r� )rU rV r r r rW L s r )r3 ru r� r� r r� r )r$ ra r: r; r r r �test_legder_axisD s z TestDerivative.test_legder_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� ) r3 �aranger � legvanderr r[ r2 r rC rs )r$ r �vr7 �coefr r r �test_legvanderU s zTestVander.test_legvanderc 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 r3 ru r � legvander2dr` �dot�flatr r r[ )r$ rc rd re rV �vanr: r; r r r �test_legvander2df s zTestVander.test_legvander2dc 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 r3 ru r � legvander3drk r� r� r r r[ )r$ rc rd re rV r� r: r; r r r �test_legvander3ds s zTestVander.test_legvander3dN) r, r- r. r3 ru r r� r� r� r r r r r� Q s r� c @ s e Zd Zdd� ZdS )� TestFittingc C s& dd� }dd� }t ttjdgdgd� t ttjdggdgd� t ttjg dgd� t ttjdgdgggd� t ttjddgdgd� t ttjdgddgd� t ttjdgdgddggd � t ttjdgdgdddgd � t ttjdgdgdg� t ttjdgdgddd g� 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�||dd d ddg�}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 � ddddg}t t�||d�ddg� t t�||ddg�ddg� t�dd�}||�}t�||d �} t t� || �|� t�||ddd g�} t t� || �|� t | | � d S )Nc S s | | d | d S )Nr r r )r r r r �f� s z"TestFitting.test_legfit.<locals>.fc S s | d | d d S )Nr� r r r )r r r r �f2� s z#TestFitting.test_legfit.<locals>.f2r r r r )�wr� r r� r y �?y � �)r r_ r �legfit� TypeErrorr3 rS r rE r rC rs r� � zeros_like�copy)r$ r� r� r r\ �coef3�coef4�coef2dr� �yw�wcoef3�wcoef2d�coef1�coef2r r r � test_legfit� sp " & zTestFitting.test_legfitN)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 � legcompanion)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 )r2 r r r� r[ )r$ r7 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 )Nr0 �c r g @) r �leggaussr� r3 r� r� �sqrt�diagonalr �eye�sum)r$ r r� r� �vv�vdr: 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 � legfromrootsr r! r2 r3 �cosrS �pirC r rE r} )r$ r; r7 �rootsr r: r r r �test_legfromroots� s * zTestMisc.test_legfromrootsc 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 g �r r )r r �legrootsr2 r3 rS r� r! )r$ r7 r: r; r r r � test_legroots� s zTestMisc.test_legrootsc 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_legtrim s zTestMisc.test_legtrimc C s t t�dd�ddg� d S )Nr r� )r r �legline)r$ r r r � test_legline s zTestMisc.test_leglinec C s6 x0t d�D ]$}tt�dg| dg �t| � q W d S )NrX r r )r2 r r r} rZ )r$ r7 r r r � test_leg2poly s zTestMisc.test_leg2polyc C s6 x0t d�D ]$}tt�t| �dg| dg � q W d S )NrX r r )r2 r r r| rZ )r$ r7 r r r � test_poly2leg s zTestMisc.test_poly2legc C s* t �ddd�}d}t�|�}t||� d S )Nr r � g �?)r3 rS r � legweightr )r$ r r: r; 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 �numpyr3 Znumpy.polynomial.legendre� polynomial�legendrer Znumpy.polynomial.polynomialr Z numpy.testingr r r r r r rs �L0�L1�L2�L3�L4�L5�L6ZL7�L8ZL9rZ r! r" r/ rT rv r� r� r� r� r� r� r, r r r r �<module> s: ":^c(0M3