U e{j @sddlZddlmZddlZddlZddlmZmZm Z ddl m Z ddl mZmZddlmZdZeeeeZddd gZd Zd d gZd d gZd ZeeeeddZdZeeeeddZd*ddZd+ddZddZ GdddZ!GdddZ"ej#$deddZ%ej#$deeeeeeddZ&ej#$deeeeeed d!Z'ej#$d"d#d id$d id%d igd&d'Z(d(d)Z)dS),N) nullcontext)assert_allcloseassert_array_almost_equal_nulpassert_array_equal)units)convolve convolve_fft)AstropyUserWarning)z>f4zf8zrrrrr'r(r? answer_key answer_dictZ result_dictkrrrtest_uniform_3sF  zTestConvolve1D.test_uniform_3c Cs>tjdddgdd}tjdddgdd}t|dt||dt||||||d }tjdd d gddtjdd d gddtjdd d gddtjd d dgddtjd d d gddtjd d d gddtjd d dgddtjdd d gddd}|rd} nd} |dkr| d7} n| d7} t||| W5QRXW5QRXdS)z Test that the different modes are producing the correct results using a uniform, non-unity kernel with three elements r+r"rBr1r2?r6r7rrA?rFrG)sum sum_zeros sum_nozerosaveragerJrK average_zerosaverage_nozerosrYrVr _wrap_zerosNrC) r>rrrrr'r(r?rQrPrrrtest_halfity_3s>    zTestConvolve1D.test_halfity_3c Cstjdtjdgdd}tjdddgdd}t|dt||||||d}W5QRX|rjt|dsjtt|}t|dddgd S) Test that a unit kernel with three elements returns the same array (except when boundary is None). This version includes a NaN value in the original array. r+rBr1r2r"r6rrMN r#r8nanrrisnanr< nan_to_numr)r>rrrrr'r(r?rrrtest_unity_3_withnans    z#TestConvolve1D.test_unity_3_withnanr+rBr1r2r"r)invaloutvalc Csl|}tjdgdd}t|dt||||||d} W5QRX|rTt| dsTtt| } t| |dS)r_r+r1r2r6rrMN)r#r8rrrbr<rcr)) r>rrrrrfrgr'r(r?rrrtest_unity_1_withnan:s    z#TestConvolve1D.test_unity_1_withnanc Cstjdtjdgdd}tjdddgdd}t|dt||||||d}W5QRX|rjt|dsjttjdddgddtjdddgddtjdddgddtjdddgddtjdd dgddtjdddgddtjd d d gddtjd d d gddtjd d d gddtjdd dgddtjdd dgddtjd d d gddd }t|D](} d| krT|| dd || d<qT|rd} nd} |dkr| d7} n| d7} |dkr| d7} t |} || | } |dkr| drd|rdnd| d<t || | dS)z Test that the different modes are producing the correct results using a uniform kernel with three elements. This version includes a NaN value in the original array. r+rBr1r2r6rrMr,rArGrLrTrUrF) rVrXrWZsum_nozeros_interpnanrYrJrKaverage_wrap_interpnanr[Zaverage_nozeros_interpnanrZZaverage_zeros_interpnanrV _interpnanrYr r\r]r r N) r#r8rarrrbr<rNrOisfiniter)) r>rrrrr'r(r?rQkeyrPposnsZanswerrrrtest_uniform_3_withnan[sZ         z%TestConvolve1D.test_uniform_3_withnancCsdtjdtjdgdd}tdddg}tjj|dddgd}t||dd tjd }t|dd d gdS) Nr+rBr1r2rMrmaskr r rr fill_valuerL)r#r8rama masked_arrayrr))r>r8kernelrvresultrrrtest_nan_interpolatesz#TestConvolve1D.test_nan_interpolatecCsLtjdtjdgdd}tdddg}t||dddd}t|d d dgdS) Nr+rBr1r2rMr rrrrFrGr#r8rarr)r>r8rwrxrrr test_nan_fillszTestConvolve1D.test_nan_fillcCsLtjdtjdgdd}tdddg}t||dddd}t|dddgdS) Nr+rBr1r2rMr rrg?rzr{rrrtest_nan_fill_twosz TestConvolve1D.test_nan_fill_twocCstjdddgdd}tdddg}tjj|dddgd}t||d d d }t|d d dgt||d d d }t||tjdddgdd}tdddg}tjj|dddgd}t||d d d }t|dd dgt||d d d }t||dS)zF Check whether convolve_fft works with masked arrays. r+rArBr1r2rMrrpr r")rrsrTrLrUN)r#r8rurvrr)r)r>r8rwrvrxZconvolve_resultZ masked_kernelrrrtest_masked_arrays0 z TestConvolve1D.test_masked_arraycCs8dddg}dddg}t||tjd}t|dddgdS)zP Check if convolve_fft works when passing a normalize function. rMrLrt)rN)rr#r%r)r{rrrtest_normalize_functions  z&TestConvolve1D.test_normalize_functionc Csvtdddg}tdg}d}|d|dg}|D]<} || } t|| |||d} |rbt| |q4t| || q4dS) zp Check that if normalize_kernel is False then the normalization tolerance is respected. rMrLrtr+g-C6?r!)rrZnormalization_zero_tolN)r#r8rr)) r>rrrrr8Z base_kernelZnormalization_rtolZ norm_errorerrrwrxrrrtest_normalization_is_respecteds   z.TestConvolve1D.test_normalization_is_respectedN)$__name__ __module__ __qualname__rmark parametrize option_namesoptionsr@rDrErSr^option_names_preserve_nanoptions_preserve_nanrer#r8rainfinputsoutputsrN itertoolsproductBOUNDARY_OPTIONSNANTREATMENT_OPTIONSZoptions_unity1withnanrhroryr|r}r~rrrrrrr*`sT        1  /      J r*c@seZdZejeeddZejeeddZ ejeeddZ eje e ddZ eje e d d Zd d Zd dZejdeddZdS)TestConvolve2Dc Cstjdddgdddgddd ggd d }tjdggd d }t|d <t||d $t||||||d}t||W5QRXW5QRXdS)zD Test that a 1x1 unit kernel returns the same array r+rArBr,r-r.r/r0"@r1r2r6r7rNrCr=rrrtest_unity_1x1_none s   z"TestConvolve2D.test_unity_1x1_nonec Cstjdddgdddgddd ggd d }tjd d d gd dd gd d d ggd d }t|d <t||d$t||||||d}t||W5QRXW5QRXdS)zl Test that a 3x3 unit kernel returns the same array (except when boundary is None). r+rArBr,r-r.r/r0rr1r2r"r6r7rNrCr=rrrtest_unity_3x39s&  zTestConvolve2D.test_unity_3x3c Cstjdddgdddgdddggdd}tjdddgdddgdddggdd}t|d8t||dt|||||r~tjnd ||d }tjd d d gd d d gd d d ggdd}tjdd dgdd dgdddggddtjd d d gd d d gd d d ggddd} | d|| d<| dd | d<| dd | d<| d| d<|rFd} nd} |dkr^| d7} n|dkrp| d7} | | } t|| W5QRXW5QRXdS)zu Test that the different modes are producing the correct results using a 3x3 uniform kernel. r"rBr+rAr1r2r6r7r)rrrsrrr,r.rr-rVrJrVrYrJrKaverage_withzeros sum_withzerosr r\r _withzerosN)r#r8rr rrar)) r>rrrrr'r(r?wrQrParrrtest_uniform_3x3UsZ       zTestConvolve2D.test_uniform_3x3c Cstjdddgdtjdgdddggd d }tjd d d gd dd gd d d ggd d }t|d t||||||d }W5QRX|rt|dstt|}t|}t||dS)z Test that a 3x3 unit kernel returns the same array (except when boundary is None). This version includes a NaN value in the original array. r+rArBr,r.r/r0rr1r2r"r6rrMrMNr`rdrrrtest_unity_3x3_withnans,     z%TestConvolve2D.test_unity_3x3_withnanc CsRtjdddgdtjdgdddggdd}tjdddgdddgdddggdd}t|d.t|||||rv|dkrvtjnd ||d }W5QRX|rt|d sttjdd dgd d d gdd dggdd}tjdddgdddgdddggdd} tjdddgddd gdddggddtjdddgdddgdddggddd} | d| | d<| d|| d<| dd | d<| dd| d<| dd| d<| dd | d<| d| d<| ddd | d<| d| d<| ddd | d<|rd} nd} |dkr| d 7} n|d!kr| d"7} |dkr$| d#7} | | tt|} t || || d$S)%z Test that the different modes are producing the correct results using a 3x3 uniform kernel. This version includes a NaN value in the original array. r"rBr+rAr1r2r6r r)rrrsrrrr-r0r,r.rrrVrYZaverage_interpnanrJrirKrZaverage_withzeros_interpnanr Z sum_interpnanZsum_withzeros_interpnanZsum_wrap_interpnanr r\r rrjN) r#r8rarrrbr<whererlr)) r>rrrrr'r(r?Zw_nZw_zrQrPrnrrrtest_uniform_3x3_withnansv         z'TestConvolve2D.test_uniform_3x3_withnanc Cs<tttf$tjdddgtd}t||W5QRXdS)z_ Test that convolve_fft raises an exception if a too-large array is passed in. ir2N)rrr MemoryErrorr#emptycomplexr)r>Zarrrrr test_big_failszTestConvolve2D.test_big_failc Csd}tj|dd}d|ddddd f<tjd dd}d|d <tjtd d Dt||dddd}t|j|t||dddd}|jdkstW5QRXt||dddd}t|jt |t |jt||dddd}|jdkstdS)z Test that convolve_fft pads to _next_fast_lengths and does not expand all dimensions to length of longest side (#11242/#10047). )rMiir1r2r+riXibi,i0)rMr)rrtrtz=psf_pad was set to False, which overrides the boundary='fill'rTF)Z return_fftpsf_padfft_pad))rM)rMi ))rLrr)rLrrN) r#Zzerosrrr rrshaper<r8)r>rZimgrwZimg_fftrrr test_paddingsF zTestConvolve2D.test_paddingrc Cstjdddgdddgdddggdd}tjdddgdddgdddggdd}t|d t|||d d d }W5QRX|d krt|tjdddgdddgdddggddn>|dkrt|tjdddgdddgdddggddntddS)Nr"r,r+rArBfloatr2gr6r F)rrrNr gggr g r.r-gzInvalid boundary specification)r#r8rrr)r)r>rr'r(r?rrrtest_non_normalized_kernel@s<& z)TestConvolve2D.test_non_normalized_kernelN)rrrrrrrrrrrrrrrrrrrrrrrrs      9  #  ] " rrc Cstjdddgdd}tjdddgdd}t|d t|||d d }W5QRX|d krtt|tjd ddgdddn$|dkrt|tjdddgddddS)zT Make sure that asymmetric convolution functions go the right direction rBr"r+r r2rMrLrtr6F)rrrr.g$@rArr!r rr-N)r#r8rrr)rr'r(r?rrrtest_asymmetric_kernelas r)rrrrr3c Csdddddddg}dddg}tj||d }tj||d }d |j_d |j_t|d t||||||d W5QRXttj||d |kstttj||d |kstd S)zF Test that convolve_fft works correctly when inputs are lists r+r,r-r.r/r0r4r5r2Fr6rN)r#r8flags writeablerrallr<) rrrrr3r8rwr'r(rrrtest_input_unmodifiedus"   rc Csdddtjdddg}dddg}tj||d}tj||d}d |j_d |j_|} |} t|d t||||||d W5QRXt|} t|} | } | }t | | || kst t | |||kst t t|| st t t|| s t d S) z> Test that convolve_fft doesn't modify the input data r+r,r-r/r0r4r5r2Fr6rN) r#rar8rrcopyrrrbrr<)rrrrr3r8rwr'r(Zx_copyZy_copyZ array_is_nanZ kernel_is_nanZ array_not_nanZkernel_not_nanrrrtest_input_unmodified_with_nans2    r error_kwargrrrc Cstjdddgdddgddd ggd d }tjdggd d }t||d d |ksRttt}t||fdd i|W5QRXt|j dt | ddkstdS)Nr+rArBr,r-r.r/r0rr r2r r6rzWith boundary='wrap', rz cannot be enabled.) r#r8rrr<rrrstrvaluerNrO)rr'r(rrrr%test_convolve_fft_boundary_wrap_errors&  rc Csdtjdddgdddgddd ggd d }tjdggd d }tjtd d t||ddW5QRXdS)Nr+rArBr,r-r.r/r0rr r2z@The 'extend' option is not implemented for fft-based convolutionrextendr6)r#r8rrNotImplementedErrorrr&rrr'test_convolve_fft_boundary_extend_errors&r)N)NF)*r contextlibrnumpyr#rZ numpy.testingrrrZastropyrr9Zastropy.convolution.convolverrZastropy.utils.exceptionsr Z VALID_DTYPESrNrZVALID_DTYPE_MATRIXrrZNORMALIZE_OPTIONSZPRESERVE_NAN_OPTIONSrrrrrr r)r*rrrrrrrrrrrrs       BD    +