Edit File: test_hermite_e.cpython-37.pyc
B ��Fd�H � @ s d Z ddlmZmZmZ ddlZddlm m Z ddl m Z ddl mZmZmZmZmZmZ e�dg�Ze�ddg�Ze�dddg�Ze�ddddg�Ze�d dd ddg�Ze�dd dd ddg�Ze�d dddd ddg�Ze�ddddddddg�Z e�dddddddddg �Z e�ddddddddddg �Z eeeeeeee e e g Z dd� Z G dd � 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 d-d.� d.e�Z*e+d/k�r e� dS )0z Tests for hermite_e 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��i\����� i��i� i���iz i���c C s t j| dd�S )Ng���ư>)�tol)�herme� hermetrim)�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_hermite_e.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 � hermedomain)�selfr r r �test_hermedomain! s z TestConstants.test_hermedomainc C s t tjdg� d S )Nr )r r � hermezero)r r r r �test_hermezero$ s z TestConstants.test_hermezeroc C s t tjdg� d S )Nr )r r �hermeone)r r r r � test_hermeone' s zTestConstants.test_hermeonec C s t tjddg� d S )Nr r )r r �hermex)r r r r � test_hermex* s zTestConstants.test_hermexN)�__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 )N� z At i=%d, j=%dr r )�err_msg)�range�np�zeros�maxr �hermeaddr r )r �i�j�msg�tgt�resr r r � test_hermeadd1 s $z TestArithmetic.test_hermeaddc 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+ )r, r- r. r/ r �hermesubr r )r r1 r2 r3 r4 r5 r r r � test_hermesub; s $z TestArithmetic.test_hermesubc C sz t t�dg�dg� t t�dg�ddg� xJtdd�D ]<}dg| dg }dg|d |ddg }t t�|�|� q6W d S )Nr r r* )r r � hermemulxr, )r r1 �serr4 r r r �test_hermemulxE s z TestArithmetic.test_hermemulxc 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+ )r, r �hermevalr �hermemulr �lenr ) r r1 �pol1�val1r2 r3 �pol2�val2�pol3�val3r r r � test_hermemulM s z TestArithmetic.test_hermemulc 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+ )r, r r0 �hermedivr= r r ) r r1 r2 r3 �ci�cjr4 �quo�remr5 r r r � test_hermediv[ s z TestArithmetic.test_hermedivN) r% r& r' r- �linspacer r6 r8 r; rE rK 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 ddg�Z d d � Z d d � Zdd� Zdd� Zdd� ZdS )�TestEvaluationg @g @g @zi,j->ijz i,j,k->ijk)r r* � r 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>w s z0TestEvaluation.test_hermeval.<locals>.<listcomp>� zAt i=%d)r+ r rN ) r r r< �sizer- rL �Helistr, r r. �shape)r �yr1 r3 r4 r5 �dimsr )r r � test_hermevalq s z TestEvaluation.test_hermevalc C s� | j \}}}| j\}}}tttj||d d� | j� || }t�||| j�}t||� t� d�} t�| | | j�}t |j dk� d S )NrN )rN r ) r rV r � ValueErrorr � hermeval2d�c2dr r- �onesr rU ) r �x1�x2�x3�y1�y2�y3r4 r5 �zr r r �test_hermeval2d� s z TestEvaluation.test_hermeval2dc C s� | j \}}}| j\}}}tttj|||d d� | j� || | }t�|||| j�}t||� t� d�} t�| | | | j�}t |j dk� d S )NrN )rN r ) r rV r rY r � hermeval3d�c3dr r- r\ r rU ) r r] r^ r_ r` ra rb r4 r5 rc r r r �test_hermeval3d� s z TestEvaluation.test_hermeval3dc C sl | j \}}}| j\}}}t�d||�}t�||| j�}t||� t�d�} t�| | | j�}t |j dk� d S )Nzi,j->ij)rN r )rN r rN r ) r rV r- �einsumr � hermegrid2dr[ r r\ r rU ) r r] r^ r_ r` ra rb r4 r5 rc r r r �test_hermegrid2d� s z TestEvaluation.test_hermegrid2dc C sr | j \}}}| j\}}}t�d|||�}t�|||| j�}t||� t�d�} t�| | | | j�}t |j dk� d S )Nz i,j,k->ijk)rN r )rN r rN r rN r ) r rV r- rh r � hermegrid3drf r r\ r rU ) r r] r^ r_ r` ra rb r4 r5 rc r r r �test_hermegrid3d� s z TestEvaluation.test_hermegrid3dN)r% r& r' r- �array�c1drh r[ rf �randomr r rV rX rd rg rj rl r r r r rM g s rM 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 rN r* )�m�k)rq rr �lbnd)rq rr �scl)rq ) r rY r �hermeintr, r � poly2herme� herme2polyr r< �list) r r1 rr r5 rt �polr4 Zhermepolru r2 r r r � test_hermeint� sz zTestIntegral.test_hermeintc 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 ru )rO rP r r r rQ s z3TestIntegral.test_hermeint_axis.<locals>.<listcomp>r )�axisc S s g | ]}t �|��qS r )r ru )rO rP r r r rQ s r c S s g | ]}t j|d d��qS )r )rr )r ru )rO rP r r r rQ # s r )rr r| )r- ro �vstack�Tr ru r )r r[ r4 r5 r r r �test_hermeint_axis s z TestIntegral.test_hermeint_axisN)r% r& r' rz r r r r r rp � s Prp 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 )rq rN )rq rt ) r rY r �hermederr, r r ru r )r r1 r4 r5 r2 r r r � test_hermeder* s" z TestDerivative.test_hermederc 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� )rO rP r r r rQ H s z5TestDerivative.test_hermeder_axis.<locals>.<listcomp>r )r| c S s g | ]}t �|��qS r )r r� )rO rP r r r rQ L s r )r- ro r} r~ r r� r )r r[ r4 r5 r r r �test_hermeder_axisD s z!TestDerivative.test_hermeder_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* rN 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 .rN r* � )r rN r{ ) r- �aranger � hermevanderr rU r, r r< rm )r r �vr1 �coefr r r �test_hermevanderU s zTestVander.test_hermevanderc C sx | j \}}}tj�d�}t�||ddg�}t�|||�}t�||j�}t||� t�|g|gddg�}t |j dk� d S )N)rN r r rN )r r* r� ) r r- ro r � hermevander2drZ �dot�flatr r rU )r r] r^ r_ rP �vanr4 r5 r r r �test_hermevander2df s z TestVander.test_hermevander2dc C s� | j \}}}tj�d�}t�|||dddg�}t�||||�}t�||j�}t||� t�|g|g|gdddg�}t |j dk� d S )N)rN r r{ r rN r )r r* � ) r r- ro r � hermevander3dre r� r� r r rU )r r] r^ r_ rP r� r4 r5 r r r �test_hermevander3ds s z TestVander.test_hermevander3dN) r% r& r' r- ro 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 rN r )r r r r �f� s z$TestFitting.test_hermefit.<locals>.fc S s | d | d d S )Nr{ rN r r )r r r r �f2� s z%TestFitting.test_hermefit.<locals>.f2r r r rN )�wr� r r{ r* y �?y � �)r rY r �hermefit� TypeErrorr- rL r r> r r< rm r~ � zeros_like�copy)r r� r� r rV �coef3�coef4�coef2dr� �yw�wcoef3�wcoef2d�coef1�coef2r r r � test_hermefit� sp " & zTestFitting.test_hermefitN)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 rY r �hermecompanion)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 )r, r r r� rU )r r1 r� r r r �test_dimensions� s z TestCompanion.test_dimensionsc C s t t�ddg�d dk� d S )Nr rN )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�� t�dtj �}t|� � |� d S )Nr) �c r rN ) r � hermegaussr� r- r� r~ �sqrt�diagonalr �eye�pi�sum)r r r� r� �vv�vdr4 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 rN r ) r �hermefromrootsr r r, r- �cosrL r� r< r r> rw )r r5 r1 �rootsry r4 r r r �test_hermefromroots� s * z TestMisc.test_hermefromrootsc 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 rN r* )r r � hermerootsr, r- rL r� r )r r1 r4 r5 r r r �test_hermeroots� s zTestMisc.test_hermerootsc C sf ddddg}t ttj|d� tt�|�|d d� � tt�|d�|d d� � tt�|d�dg� d S )NrN r r r r )r rY r r r )r r� r r r �test_hermetrim s zTestMisc.test_hermetrimc C s t t�dd�ddg� d S )Nr r{ )r r � hermeline)r r r r �test_hermeline s zTestMisc.test_hermelinec C s6 x0t d�D ]$}tt�dg| dg �t| � q W d S )NrR r r )r, r r rw rT )r r1 r r r �test_herme2poly s zTestMisc.test_herme2polyc C s6 x0t d�D ]$}tt�t| �dg| dg � q W d S )NrR r r )r, r r rv rT )r r1 r r r �test_poly2herme s zTestMisc.test_poly2hermec C s8 t �ddd�}t �d|d �}t�|�}t||� d S )N���r* � g �rN )r- rL �expr � hermeweightr )r r r4 r5 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 �numpyr- Znumpy.polynomial.hermite_e� polynomial� hermite_er Znumpy.polynomial.polynomialr Z numpy.testingr r r r r r rm ZHe0ZHe1ZHe2ZHe3ZHe4ZHe5ZHe6ZHe7ZHe8ZHe9rT r r r( rM rp r� r� r� r� r� r� r% r r r r �<module> s: 9^c)0M3