IR-e, @ddlZddlZddlZddlmZddlmZddl m Z ddl m Z m Z ddlmZddlmZe gdZd Zd Zd Zd Zd ZdZdZdZdZejdeddedddfeddeddeddfedgedddfedgeddddfgdZejddejzdejzdejzdejzgdZ dS)N) assert_equal)units)Time)aggregate_downsamplereduceat) TimeSeries)AstropyUserWarning)2016-03-22T12:30:31z2016-03-22T12:30:322016-03-22T12:30:332016-03-22T12:30:342016-03-22T12:30:35c tjtjdgd}ttjdgdtj}t ||ttjdtjddddtj}t |tjgdttjdgdtj}t |tjgdt ttjdtjddddtjttjdtjddddtjdS)N)rr)g?@@g@)g?rrrg @rrg@) npaddrarangesumrmeanarraynanmean) add_output sum_output mean_outputnanmean_outputs Hlib/python3.11/site-packages/astropy/timeseries/tests/test_downsample.py test_reduceatr%s^1/G/G/GHHJ")A,,(@(@(@"&IIJZ(((29Q<<1ccc):BGDDKbh';';';<<===bill,D,D,DbgNNN*L*L*L!M!MNNN1ry||CCaC0"'::1ry||CCaC0"*==ctjtd5tdddddS#1swxYwYdS)Nz"time_series should be a TimeSeriesmatch)pytestraises TypeErrorrr&r$test_timeseries_invalidr.*s y(L M M M##T"""##################s 9==ctjtd5tt dddddS#1swxYwYdS)Nz3'time_bin_size' should be a Quantity or a TimeDeltar(r) time_bin_size)r*r+r,rrr-r&r$test_time_bin_invalidr1/s O   << Z\\;;;;<<<<<<<<<<<<<<<<<r&cttgdgdg}t|ddg}t|jjddgdS) Nr3r4r5r r )r=rA2016-03-22T12:30:35.000)rr9rrr@rE)r: down_auto_ends r$test_time_bin_end_autorLVsq ???*;C5 I I IB) 13HIM"' "$=>r&cttgdgdg}t|ddg}t|jjddgdS) Nr3r4r5r r )r@r>rA)rr9rrr=rE)r:down_start_arrays r$test_time_bin_start_arrayrOdsr ???*;C5 I I IB , /1FG', "$=>r&c4ttgdgdg}t|d}t|jjddgt t}d|d z ||d zfD].}t|t| }t ||ksJ/dS) Nr3r4r5r)n_binsr>rArr)r=rQ)rr9rrr=rElen)r: down_nbinsn_timesrQs r$ test_nbinsrUts ???*;C5 I I IB&b333J!& "$=> *ooGgk7GaK8**)"ZPVWWW :')))))**r&c  ttgdgdg}ttgdtjzgdg}tjdkrdtjz}d}n,dtjz}|jddtjzz }t|||}tj |j gd |zt|j jtgd t|djjt#jgdt|d |z|}tj |j gd |zt|j jtgd t|djjt#jgdt|d|z|}tj |j ddg|zt|j jtddgt|djjt#jd dgt|d|z|}tj |j ddg|zt|j jtddgt|djjt#jd dgt|d|z|}tj |j ddg|zt|j jtddg|djjdksJt|djt#jddgt|gd|z|} tj | j gd|zt| j jtgdt| djjt#jgdt|tddgtddg} tj | j ddgtjzt| j jtddgt| djjt#jddgt-jt0d5t|tddgtd d!g} t| djt#jd dgddddS#1swxYwYdS)"Nr3r4r5lrg !?r)r0r=)rrrrr)r>2016-03-22T12:30:32.000rArBrJr)rrr)r>rArJ)rrrrr>rBrrJrctrg@)rrr)r>rArB)rrrrWr?zaOverlapping bins should be avoided since they can lead to double-counting of data during binning.r(rAr z2016-03-22T12:30:36.000)rr9rCcountsysmaxsizerDr6nsriscloser0rr=rErr7rrunitnamesecondr*warnsr ) r:ts_units time_bin_incrr=down_1down_2down_3down_4 down_unitsdown_uneven_binsdown_time_arraydown_overlap_binss r$test_downsamplerls  ???*;C5 I I IB 89#H  {UAC "QS( a!$h. ! -FIf"OOOm$CDDD"         !&(A(ABBB ! !m+NFIf"III $=>>>"         !&(;(;<<< ! !m+NFIf"QF]$:;;;" ')B CDD!&!Q(8(8999 ! !m+NFIf"QF]$:;;;" ')B CDD!&!Q(8(8999%M 1.JIj&A(>???!& ')B CDD c?  $ , , , ,C%rxc ';';<<<, )))m3NI,iii-.GHHH',         !#&+0"(9992E2EFFF+ 68QRSS46OPQQO Io+aVah->???&+ ')B CDD%*/1a&1A1ABBB  B    D D1 !:    D t99N++ + + + + tO$) * **** w{{}}~a(((L,D  ssuu      * * * IaLBsGAJ & & & & + * & &r&diff_from_base 2dcdtjz}ddlm}t ddd}||z}||z}|||}|||}||ksJdS) ae A test on time precision limit supported by downsample(). It is related to an implementation details: that time comparison (and sorting indirectly) is done with relative time for computational efficiency. The relative time converted has a slight loss of precision, which worsens as the gap between a time and the base time increases, e.g., when downsampling a timeseries that combines current observation with archival data years back. This test is to document the acceptable precision limit. see also: https://github.com/astropy/astropy/pull/13069#issuecomment-1093069184 ir)_to_relative_longdoublez1980-01-01T12:30:31.000rEtdb)formatscaleN)rCr\astropy.timeseries.downsampleryr)rtprecision_limitryt_baset2t3r_t2r_t3s r$test_time_precision_limitr!s"ADjOEEEEEE +F% H H HF . B o B " "2v . .D " "2v . .D $;;;;;;r&)!rZnumpyrr* numpy.testingrastropyrrC astropy.timerr}rrastropy.timeseries.sampledrastropy.utils.exceptionsr r9r%r.r1r;rHrLrOrUrlmark parametrizersyearrr-r&r$rsK  &&&&&&HHHHHHHH111111777777 T     ### <<<!!!(       ***"|D|D|D~( BQBABB. ABBBQBAaC9 QC*RaR.$/ QC*SSqS/40    (q16z2;QV S16\Jr&