a ߙfb @sddlZddlmZddlZddlZddlmZmZm Z ddl m Z m Z ddl mZddlmZdZeeeeZgdZd Zd d gZd d gZd Zeeeed d ZdZeeeed d Zd)ddZd*ddZddZ GdddZ!GdddZ"ej#$deddZ%ej#$deeeeeeddZ&ej#$deeeeeedd Z'ej#$d!d"d id#d id$d igd%d&Z(d'd(Z)dS)+N) nullcontext)assert_allcloseassert_array_equalassert_array_almost_equal_nulp) convolve_fftconvolve)AstropyUserWarning)units)z>f4zf8z?rG)r>r>r>rJrWrHrY)sum sum_zeros sum_nozerosaveragerMrN average_zerosaverage_nozerosr^r[r _wrap_zerosNr@) r:rrrrr%r&r;rTrSrrrtest_halfity_3s2  zTestConvolve1D.test_halfity_3c Cstjdtjdgdd}tjgddd}t|d$t||||||d}Wdn1s\0Y|r|t|d s|Jt|}t|gd dS) 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*r?r0r1rBr5rNrPrDr#r7nanrrisnan nan_to_numr'r:rrrrr%r&r;rrrtest_unity_3_withnans  $ z#TestConvolve1D.test_unity_3_withnanr*r?r0r1rDr)invaloutvalc Cs|}tjdgdd}t|d$t||||||d} Wdn1sH0Y|rht| dshJt| } t| |dS)rdr*r0r1r5rNrP)r#r7rrrgrhr') r:rrrrrkrlr%r&r;rrrtest_unity_1_withnans  $ z#TestConvolve1D.test_unity_1_withnanc Cstjdtjdgdd}tjgddd}t|d$t||||||d}Wdn1s\0Y|r|t|d s|Jtjgd ddtjgd ddtjgd ddtjgd ddtjgd ddtjgd ddtjgd ddtjgdddtjgdddtjgd ddtjgdddtjgdddd }t|D](} d| vrN|| dd|| d<qN|rd} nd} |dkr| d7} n| d7} |dkr| d7} t|} || | } |dkr| d rd|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*r?r0r1rEr5rNrPrFr=rIrJ)rOrOrOrZrGrX) r[r]r\Zsum_nozeros_interpnanr^rMrNaverage_wrap_interpnanr`Zaverage_nozeros_interpnanr_Zaverage_zeros_interpnanr[r> _interpnanr^r rarbr r ) r#r7rfrrrgrQrRisfiniter') r:rrrrr%r&r;rTkeyrSposnsZanswerrrrtest_uniform_3_withnan!sN  $      z%TestConvolve1D.test_uniform_3_withnancCs^tjdtjdgdd}tgd}tjj|gdd}t||dd tjd }t|gd dS) Nr*r?r0r1rPrPrPrrPrmaskr r rr fill_valuerPrO)r#r7rfma masked_arrayrr')r:r7kernelr~resultrrrtest_nan_interpolateesz#TestConvolve1D.test_nan_interpolatecCsHtjdtjdgdd}tgd}t||dddd}t|gd dS) Nr*r?r0r1rur rryrGr#r7rfrr'r:r7rrrrr test_nan_fillrszTestConvolve1D.test_nan_fillcCsHtjdtjdgdd}tgd}t||dddd}t|gd dS) Nr*r?r0r1rur rPry)r*竪?rrrrrrtest_nan_fill_twosz TestConvolve1D.test_nan_fill_twocCstjgddd}tgd}tjj|gdd}t||ddd }t|gd t||ddd }t||tjgddd}tgd}tjj|gdd}t||ddd }t|gd t||ddd }t||d S) zF Check whether convolve_fft works with masked arrays. r=r0r1rurvrwr r")rrz)rWrOrY)rPrOrPN)r#r7r}r~rr'r)r:r7rr~rZconvolve_resultZ masked_kernelrrrtest_masked_arrays, z TestConvolve1D.test_masked_arraycCs2gd}gd}t||tjd}t|gddS)zP Check if convolve_fft works when passing a normalize function. r{)r|r|r|)r)r|N)rr#r$r'rrrrtest_normalize_functionsz&TestConvolve1D.test_normalize_functionc Csttgd}tdg}d}|d|dg}|D]<} || } t|| |||d} |r`t| |q2t| || q2dS)zp Check that if normalize_kernel is False then the normalization tolerance is respected. r{r*g-C6?r!)rrZnormalization_zero_tolN)r#r7rr') r:rrrrr7Z base_kernelZnormalization_rtolZ norm_errorerrrrrrrtest_normalization_is_respecteds   z.TestConvolve1D.test_normalization_is_respectedN)$__name__ __module__ __qualname__rmark parametrize option_namesoptionsr<rArCrVrcoption_names_preserve_nanoptions_preserve_nanrjr#r7rfinfZinputsZoutputsrQ itertoolsproductBOUNDARY_OPTIONSNANTREATMENT_OPTIONSZoptions_unity1withnanrmrtrrrrrrrrrrr(SsJ        )  )      C  r(c@seZdZejeeddZejeeddZ ejeeddZ eje e ddZ eje e d d Zd d Zd dZejdeddZdS)TestConvolve2Dc Cstjgdgdgdgdd}tjdggdd}t|dZt||d.t||||||d }t||Wd n1s~0YWd n1s0Yd S) zD Test that a 1x1 unit kernel returns the same array r=r+r,r-r.r/"@r0r1r*r5r6rNr@r9rrrtest_unity_1x1_nones  z"TestConvolve2D.test_unity_1x1_nonec Cstjgdgdgdgdd}tjgdgdgdgdd}t|dZt||d .t||||||d }t||Wd n1s0YWd n1s0Yd S) zl Test that a 3x3 unit kernel returns the same array (except when boundary is None). r=rrr0r1r"r"r"rBr5r6rNr@r9rrrtest_unity_3x3s(   zTestConvolve2D.test_unity_3x3c Cstjgdgdgdgdd}tjgdgdgdgdd}t|dFt||dt|||||rrtjnd ||d }tjgd gd gd gdd}tjgd gdgdgddtjgdgdgdgddd} | d|| d<| dd| d<| dd| d<| d| d<|r(d} nd} |dkr@| d7} n|dkrR| d7} | | } t|| Wdn1sz0YWdn1s0YdS)zu Test that the different modes are producing the correct results using a 3x3 uniform kernel. r"r"r?)r*r"r"r"r>r"r0r1rEr5r6r)rrrzrrr+r-r+r-rr-rFr?r-r,r?r?r>r-r-r-r[rMr[r^rMrrNaverage_withzeros sum_withzerosr rar _withzerosN)r#r7rr rrfr') r:rrrrr%r&r;wrTrSarrrtest_uniform_3x3sh          zTestConvolve2D.test_uniform_3x3c Cstjgddtjdggdgdd}tjgdgdgdgdd}t|d $t||||||d }Wd n1sx0Y|rt|d sJt|}t|}t||d S) 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=r+r-rr0r1rrBr5rNrPrPrerirrrtest_unity_3x3_withnan<s.    $  z%TestConvolve2D.test_unity_3x3_withnanc CsDtjgddtjdggdgdd}tjgdgdgdgdd}t|d8t|||||rl|d krltjnd ||d }Wd n1s0Y|rt|d sJtjgdgdgdgdd}tjgdgdgdgdd} tjgdgdgdgddtjgdgdgdgddd} | d| | d<| d|| d<| dd| d<| dd| d<| dd| d<| dd| d <| d| d!<| dd"d| d#<| d| d$<| dd"d| d%<|rd} nd} |d&kr| d'7} n|d(kr| d)7} |d kr| 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. rr*r"rr0r1rEr5r r)rrrzrrNr)r?r,r?)r,r/r,rrrFrrrrr[r^Zaverage_interpnanrMr/rnrrNrZaverage_withzeros_interpnanr Z sum_interpnanZsum_withzeros_interpnanZsum_wrap_interpnanr rar rro) r#r7rfrrrgwhererqr') r:rrrrr%r&r;Zw_nZw_zrTrSrsrrrtest_uniform_3x3_withnan[s    $          z'TestConvolve2D.test_uniform_3x3_withnancCsNtttf,tjgdtd}t||Wdn1s@0YdS)zN Test that convolve_fft raises an exception if a too-large array is passed in.)rrr1N)rrr MemoryErrorr#emptycomplexr)r:Zarrrrr test_big_failszTestConvolve2D.test_big_failcCsd}tj|dd}d|ddddd f<tjd dd}d|d <tjtd d Nt||dddd}t|j|t||dddd}|jdvsJWdn1s0Yt||dddd}t|jt|t|jt||dddd}|jdvsJdS)z Test that convolve_fft pads to _next_fast_lengths and does not expand all dimensions to length of longest side (#11242/#10047). )rPiir0r1r*riXibi,i0)rPr)rr|r|z=psf_pad was set to False, which overrides the boundary='fill'rTF)Z return_fftpsf_padfft_pad))rP)rPi N))rOrr)rOrr) r#Zzerosrrrrrshaper7)r:rZimgrZimg_fftrrr test_paddings  ,zTestConvolve2D.test_paddingrcCstjgdgdgdgdd}tjgdgdgdgdd}t|d"t|||d d d }Wdn1sr0Y|d vrt|tjgd gdgdgddn8|dkrt|tjgdgdgdgddntddS)N)r"r"r+)r*r>r")r"r?r"floatr1)r*r*)rr"rr5r F)rrrNr )r*gr>)r*r"g)grrr )r"g r-)r,r")r>r?rzInvalid boundary specification)r#r7rrr'r)r:rr%r&r;rrrtest_non_normalized_kernels>    $  z)TestConvolve2D.test_non_normalized_kernelN)rrrrrrrrrrrrrrrrrrrrrrrrs      3    V  rrcCstjgddd}tjgddd}t|d t|||dd}Wdn1sT0Y|d vrt|tjgd d dd n"|d krt|tjgdd dd dS)zT Make sure that asymmetric convolution functions go the right direction )r?r"r*r r1r{r5F)rrNr)r-$@r>rr!r )rrr,)r#r7rrr)rr%r&r;rrrtest_asymmetric_kernels .r)rrrrr2c Csgd}gd}tj||d}tj||d}d|j_d|j_t|d$t||||||d} Wdn1sp0Yttj||d|ksJttj||d|ksJdS)zF Test that convolve_fft works correctly when inputs are lists r)r3r1Fr5rN)r#r7flags writeablerrall) rrrrr2r7rr%r&r;rrrtest_input_unmodifieds   $rc Cs"dddtjdddg}gd}tj||d}tj||d}d|j_d|j_|} |} t|d $t||||||d } Wd n1s0Yt|} t|} | }| }t | |||ksJt | |||ksJt t|| sJt t|| sJd S) z> Test that convolve_fft doesn't modify the input data r*r+r,r.r/r3r1Fr5rN) r#rfr7rrcopyrrrgr)rrrrr2r7rr%r&Zx_copyZy_copyr;Z array_is_nanZ kernel_is_nanZ array_not_nanZkernel_not_nanrrrtest_input_unmodified_with_nan$s(   $  r error_kwargrrrcCstjgdgdgdgdd}tjdggdd}t||dd|ksLJtt&}t||fd di|Wdn1s0Yt|jd t | d d ksJdS) Nr=rrr r1r*r r5rzWith boundary='wrap', rz cannot be enabled.) r#r7rrrrrstrvaluerQrR)rr%r&rrrr%test_convolve_fft_boundary_wrap_errorOs  4rcCstjgdgdgdgdd}tjdggdd}tt}t||ddWdn1s`0Yt|jd ks|JdS) Nr=rrr r1r*extendr5z@The 'extend' option is not implemented for fft-based convolution)r#r7rrNotImplementedErrorrrr)r%r&rrrr'test_convolve_fft_boundary_extend_error]s  ,r)N)NF)*r contextlibrrZnumpyr#Z numpy.testingrrrZastropy.convolution.convolverrZastropy.utils.exceptionsrZastropyr r8Z VALID_DTYPESrQrZVALID_DTYPE_MATRIXrrZNORMALIZE_OPTIONSZPRESERVE_NAN_OPTIONSrrrrrr r'r(rrrrrrrrrrrrst          $