
    Џkh2                     *   d dl Z d dlZ d dlZd dlmZ d dlmZmZm	Z	m
Z
mZ d dlZd dlmZ d dlmZ ddlmZ ej&                  j(                  Z eddd	g
      gZe j,                  dk(  xr ej.                  dk  Z edd       G d d             Z edd       G d d             Zd Zd Zd Zd Zd Zd Z d Z!d Z"d Z#d Z$d Z% edd      d        Z&d  Z'd! Z( eddgd"#      d$        Z) eddgd"#      d%        Z* edd&      d'        Z+d( Z, edd)*      d+        Z-d, Z. edd-*       G d. d/             Z/d0 Z0d1 Z1d2 Z2d3 Z3d4 Z4d5 Z5d6 Z6d7 Z7d8 Z8d9 Z9d: Z:d; Z;d< Z<d= Z=d> Z>d? Z?d@ Z@dA ZAdB ZBdC ZCdD ZDdE ZEdF ZFdG ZGdH ZHdI ZIdJ ZJdK ZKdL ZLdM ZMdN ZNdO ZOdP ZPdQ ZQdR ZRdS ZSdT ZTdU ZUdV ZVdW ZWdX ZXdY ZYdZ ZZd[ Z[d\ Z\d] Z]d^ Z^d_ Z_d` Z`da Zadb Zbdc Zcdd Zdde Zedf Zfdg Zgdh Zhdi Zidj Zjdk Zkdl Zldm Zmdn Zndo Zodp Zpdq Zqdr Zrds Zsdt Ztdu Zudv Zvdw Zwdx Zxdy Zydz Zzd{ Z{ edd|      d}        Z|d~ Z} edd*       G d d             Z~ ed      ej&                  j                  dej                   ej                  g      d               Zy)    N)suppress_warnings)is_torchxp_assert_equalxp_assert_closeassert_array_almost_equalassert_almost_equal)raises   )typesTcupyz	jax.numpy)cpu_only
exceptionsnt2z test internal numpy-only helpersnp_onlyreasonc                   .    e Zd ZdZd Zd Zd Zd Zd Zy)Test_measurements_statszndimage._measurements._stats() is a utility used by other functions.

        Since internal ndimage/_measurements.py code is NumPy-only,
        so is this this test class.
    c           	         g d}g d}ddg}dD ]  }t        j                  |      j                  |      }t        j                  |      j                  |      }t        j                  j                  |||      \  }}t        rdt         j                  ini }t        |t        j                  ddgfi |       t        |t        j                  d	d
g              y )Nr   r
         r   r   r
   r
   r   r
   )   )r   r   labelsindexdtyper         ?       @
nparrayreshapendimage_measurements_statsIS_WINDOWS_AND_NP1int64r   asarray	selfxpxr   r   shpcountssums	dtype_args	            a/var/www/teggl/fontify/venv/lib/python3.12/site-packages/scipy/ndimage/tests/test_measurements.pytest_azTest_measurements_stats.test_a#   s    A! 	:C##C(AXXf%--c2F"0077& 8 /LFD 0B"((+rIFBJJ1v$C$CDD"**c3Z"89	:    c           	         g d}g d}ddg}dD ]  }t        j                  |      j                  |      }t        j                  |      j                  |      }t        j                  j                  |||      \  }}t        rdt         j                  ini }t        |t        j                  ddgfi |       t        |t        j                  d	d
g              y )Nr   r   r   	   r:   r   r:   r   r   r    r   r!   r"   r#   r-   s	            r5   test_bzTest_measurements_stats.test_b1   s     A! 	:C##C(AXXf%--c2F"0077& 8 /LFD 0B"((+rIFBJJ1v$C$CDD"**c3Z"89	:r7   c           	         g d}g d}ddg}dD ]  }t        j                  |      j                  |      }t        j                  |      j                  |      }t        j                  j                  |||d      \  }}}t        rdt         j                  ini }	t        |t        j                  d	d	gfi |	       t        |t        j                  d
dg             t        |t        j                  ddg              y )Nr   r   r   r
   r   Tr   r   centeredr    r   r!   r"         ?r#   
r.   r/   r0   r   r   r1   r2   r3   centersr4   s
             r5   test_a_centeredz'Test_measurements_stats.test_a_centeredA       A! 		=C##C(AXXf%--c2F$+$9$9$@$@& %A %>!FD' 0B"((+rIFBJJ1v$C$CDD"**c3Z"89GRZZc
%;<		=r7   c           	         g d}g d}ddg}dD ]  }t        j                  |      j                  |      }t        j                  |      j                  |      }t        j                  j                  |||d      \  }}}t        rdt         j                  ini }	t        |t        j                  d	d	gfi |	       t        |t        j                  d
dg             t        |t        j                  ddg              y )Nr   r9   r   r:   r   Tr=   r    r   r!   r"   r?   r#   r@   s
             r5   test_b_centeredz'Test_measurements_stats.test_b_centeredP   rC   r7   c           	         g d}g d}ddg}dD ]  }t        j                  |      j                  |      }t        j                  |      j                  |      }t        j                  j                  |||d      \  }}}t        rdt         j                  ini }	t        |t        j                  d	d	gfi |	       t        |t        j                  d
dg             t        |t        j                  ddg              y )Nr   )        rG         "@rH   rG   rH   r   Tr=   r    r   r!   r"   r?   r#   r@   s
             r5   test_nonint_labelsz*Test_measurements_stats.test_nonint_labels_   s    %c
! 		=C##C(AXXf%--c2F$+$9$9$@$@& %A %>!FD' 0B"((+rIFBJJ1v$C$CDD"**c3Z"89GRZZc
%;<		=r7   N)	__name__
__module____qualname____doc__r6   r;   rB   rE   rI    r7   r5   r   r      s     
:: ===r7   r   c                       e Zd ZdZd Zy)Test_measurements_selectzEndimage._measurements._select() is a utility used by other functions.c                    g d}g dddgfg dddgfg ddd	gfg}|D ]  \  }}t         j                  j                  |||
      }t        |      dk(  sJ t         j                  j                  |||d      }t        |      dk(  sJ t	        |d   ddg       t         j                  j                  |||d      }t        |      dk(  sJ t	        |d   ddg       t         j                  j                  |||dd      }t        |      dk(  sJ t	        |d   ddg       t	        |d   ddg       |d   j
                  j                  dk(  sJ t         j                  j                  |||dd      }t        |      dk(  sJ t	        |d   ddg       t	        |d   ddg       |d   j
                  j                  dk(  rJ  y )N)r   r
   r   r   r   r   r
   r9   r:   )rG   rG         @rR   rG   rR   r   T)r   r   find_maxr   )r   r   find_minr   )r   r   rT   find_min_positions   i)r   r   rS   find_max_positions)r'   r(   _selectlenr   r    kind)r.   r/   r0   casesr   r   results          r5   
test_basicz#Test_measurements_select.test_basics   s   Aq6"Aq6"!C:.

 # 	/MFE**22& 3 /Fv;!###**22& 3 >Fv;!###F1I1v.**22& 3 >Fv;!###F1I1v.**22&#' 3 )F v;!###F1I1v.F1I1v.!9??''3...**22&#' 3 )F v;!###F1I1v.F1I1v.!9??''3...3	/r7   N)rJ   rK   rL   rM   r^   rN   r7   r5   rP   rP   o   s
    O /r7   rP   c                 r    | j                  d      }t        j                  |      \  }}|dk(  sJ |dk(  sJ y )NrN   r
   )onesr'   labelr/   dataoutns       r5   test_label01rf      s7    772;D]]4 FC!8O86M6r7   c                 r    | j                  d      }t        j                  |      \  }}|dk(  sJ |dk(  sJ y )NrN   r   )zerosr'   ra   rb   s       r5   test_label02ri      s7    88B<D]]4 FC!8O86M6r7   c                     | j                  dg      }t        j                  |      \  }}t        || j	                  dg             |dk(  sJ y )Nr
   r`   r'   ra   r   r,   rb   s       r5   test_label03rl      sB    77A3<D]]4 FCc2::qc?36M6r7   c                     | j                  dg      }t        j                  |      \  }}t        || j	                  dg             |dk(  sJ y )Nr
   r   )rh   r'   ra   r   r,   rb   s       r5   test_label04rn      sB    88QC=D]]4 FCc2::qc?36M6r7   c                     | j                  dg      }t        j                  |      \  }}t        || j	                  g d             |dk(  sJ y )N   )r
   r
   r
   r
   r
   r
   rk   rb   s       r5   test_label05rq      sA    77A3<D]]4 FCc2::o#>?6M6r7   c                     | j                  g d      }t        j                  |      \  }}t        || j                  g d             |dk(  sJ y )N)r
   r   r
   r
   r   r
   r
   r   r   r   r   rV   rV   r,   r'   ra   r   rb   s       r5   test_label06ru      sB    ::()D]]4 FCc2::.@#AB6M6r7   c                     | j                  g dg dg dg dg dg dg      }t        j                  |      \  }}t        || j                  g dg dg dg dg dg dg             |dk(  sJ y )Nr   r   r   r   r   r   r   rt   rb   s       r5   test_label07rx      st    ::))))))+ ,D ]]4 FCc2::%7%7%7%7%7%7%9$: ; 6M6r7   c                     | j                  g dg dg dg dg dg dg      }t        j                  |      \  }}t        || j                  g dg dg dg dg dg d	g             |d
k(  sJ y )Nr
   r   r   r   r   r   r   r   r
   r
   r   r   r   r   r
   r
   r
   r   r
   r
   r   r   r   r   r   r   r   r
   r
   r   r   r   r   r   r   r   r   r   r   r   r   r   rV   rV   r   r   r   r   r   r   r   r   r   r   r   rt   rb   s       r5   test_label08r      sw    ::))))))+ ,D ]]4 FCc2::/A/A/A/A/A/A/C $D E 6M6r7   c                 2   | j                  g dg dg dg dg dg dg      }t        j                  dd      }| j                  |      }t        j                  ||      \  }}t	        || j                  g dg dg dg d	g d	g d
g             |dk(  sJ y )Nrz   r{   r|   r}   r~   r   r   r   )r   r   r   r   r   r   )r   r   r   rV   rV   r   rV   r,   r'   generate_binary_structurera   r   r/   rc   structrd   re   s        r5   test_label09r      s    ::))))))+ ,D ..q!4FZZF]]4(FCc2::/A/A/A/A/A/A/C $D E 6M6r7   c           
         | j                  g dg dg dg dg      }t        j                  dd      }| j                  |      }t        j                  ||      \  }}t	        || j                  g dg dg dg dg             |dk(  sJ y )Nrw   )r   r
   r
   r   r
   r   )r   r
   r
   r
   r
   r   r   r
   r   r   s        r5   test_label10r      s    ::))))+ ,D ..q!4FZZF]]4(FCc2::/A/A/A/A/C $D E 6M6r7   c           
         t         D ]  }t        | |      }| j                  g dg dg dg dg dg dg|      }t        j                  |      \  }}g dg dg dg d	g d	g d
g}| j                  |      }t        ||       |dk(  rJ  y )Nrz   r{   r|   r}   r~   r    r   r   r   r   r   r   getattrr,   r'   ra   r   r/   typer    rc   rd   re   expecteds          r5   test_label11r     s     D!zz------/
 7<  = t$Q&&&&&&( ::h'!#x0Avv#r7   z inplace output is numpy-specificc           
         t         D ]  }t        | |      }| j                  g dg dg dg dg dg dg|      }t        j                  ||      }g dg dg d	g d
g d
g dg}| j                  |      }t        ||       |dk(  rJ  y )Nrz   r{   r|   r}   r~   r   outputr   r   r   r   r   r   )r/   r   r    rc   re   r   s         r5   test_label11_inplacer     s     D!zz------/
 7<  = MM$t,&&&&&&( ::h'!$1Avv#r7   c           	         t         D ]}  }t        | |      }| j                  g dg dg dg dg dg|      }t        j                  |      \  }}g dg dg dg dg dg}| j                  |      }t        ||       |dk(  r}J  y )N)r   r   r   r   r
   r
   )r   r   r   r   r   r
   )r   r   r
   r   r
   r
   )r   r   r
   r
   r
   r
   r~   r   r
   r   r   s          r5   test_label12r   0  s     D!zz-----	/ 7<	  =
 t$Q&&&&&	(
 ::h'!#x0Avvr7   c                    t         D ]w  }t        | |      }| j                  g dg dg dg dg|      }t        j                  |      \  }}g dg dg dg dg}| j                  |      }t        ||       |dk(  rwJ  y )N)r
   r   r
   r
   r
   r   r
   r
   r
   r   r
   )r
   r
   r
   r   r
   r
   r
   r   r
   r
   r
   )r
   r   r   r   r   r   r   r   r   r   r
   )r
   r
   r
   r
   r
   r
   r
   r
   r
   r
   r
   r   r
   r   r   s          r5   test_label13r   C  s     D!zz<<<<> !&	  '
 t$Q55557 ::h'!#x0Avvr7   zoutput=dtype is numpy-specific)r   r   r   c                    | j                  dg      }t        D ]q  }t        | |      }| j                  dg|      }t	        j
                  ||      }t        || j                  |j                  |j                               |dk(  rqJ  y )Nrp   r   r   r
   )	r`   r   r   rh   r'   ra   r   shaper    r/   rc   tr    r   re   s         r5   test_label_output_typedr   U  s     77A3<D A1#U+MM$v.!&"$''&,,fll'"K	MAvvr7   c                    | j                  dg      }t        D ]j  }t        | |      }t        j                  ||      \  }}t        || j                  |j                  |j                               |j                  |k(  rjJ  y )Nrp   r   r   )r`   r   r   r'   ra   r   r   r    r   s         r5   test_label_output_dtyper   b  sv     77A3<D !AMM$u5	!&"$''&,,fll'"K	M||q   !r7   z!in-place output is numpy-specificc                     | j                  dg      }t        D ]C  }t        | |      }| j                  dg|      }t	        t
        t        j                  ||       E y )Nrp   
   r   r   )r`   r   r   rh   assert_raises
ValueErrorr'   ra   )r/   rc   r   r    r   s        r5   test_label_output_wrong_sizer   n  sS    77A3<D FA2$e,j'--fEFr7   c           
         t        j                  t        j                  j	                  t        j                  j                  t              dd            }t        j                  t        j                  j	                  t        j                  j                  t              dd            }t        j                  t        j                  j	                  t        j                  j                  t              dd            }|j                  d      }|j                  d      }|j                  d      }| j                  |      }| j                  |      }| j                  |      }d}t        |j                  d         D ]k  }||d d d d f   }t        |j                  d         D ]B  }||d d d d f   }t        t        j                  ||      d   ||d d d d f   d	       |d
z  }D m y )Nrc   zlabel_inputs.txtzlabel_strels.txtzlabel_results.txt)   r   )r   rV   rV   r   Fcheck_dtyper
   )r$   loadtxtospathjoindirname__file__r&   r,   ranger   r   r'   ra   )	r/   rc   strelsresultsrrW   djss	            r5   test_label_structuring_elementsr   w  s~   ::bggll277??$-/ 0DZZ
!6+=? @Fjj
!6+>@ AG<<
#D^^J'Fooj)G::dDZZFjj!G	A4::a=! AqMv||A' 	Aq!QwAGMM!Q/2GAq!G4DRWXFA	r7   z2`cupyx.scipy.ndimage` does not have `find_objects`r   c                 >    d fd	}t        j                  t         j                        t        j                  d      k7  r]d}t         j                  j	                  t        j
                  |            j                  |      } j                  |      } ||       y y )Nc                     | |kD  }t        |j                        }t        j                  ||      }	j	                  |      }t        j
                  ||      \  }}t        j                  |      }y )N)rZ   r   r'   r   r,   ra   find_objects)
imgthreshsizemaskrankr   laco_r/   s
            r5   SEztest_ticket_742.<locals>.SE  s`    V|4::224>F#t%'B  $r7   rW   )rV     r   )gffffff?r   )r$   r    intprandomrandprodr&   r,   )r/   r   r   as   `   r5   test_ticket_742r     sj    % 
xxBHHSM)IINN2775>*2259JJqM
1 *r7   c                     t        j                  d      }d|ddddf<   d|ddddf<   d|d<   d|d<   d|d<   | j                  |      }t        j                  || j                  d	            d   dk(  sJ y)
z/Github issue #3025 - improper merging of labels)<   i@  r
   N  i  )$   r   )#   i  )r   i  )rV   rV   )r$   rh   r,   r'   ra   r`   )r/   r   s     r5   test_gh_issue_3025r     s~    
AAa#gJAagJAgJAgJAgJ


1A==BGGFO,Q/1444r7   z-cupyx.scipy.ndimage does not have find_objectc                   H    e Zd Zd Zd Zd Zd Zd Zd Zd Z	d Z
d	 Zd
 Zy)TestFindObjectsc                    t         j                  j                  dd      }|j                  |      }t	        j
                  |dkD        \  }}|j                  |j                  |j                  fv sJ t	        j                  |       y )Nr   r?   )
r$   r   r   r,   r'   ra   r    int32r+   r   )r.   r/   
test_arrayra   no_featuress        r5   test_label_default_dtypez(TestFindObjects.test_label_default_dtype  si    YY^^B+
ZZ
+
$]]:+;<{{{rxx2222U#r7   c                 x    |j                  g |j                        }t        j                  |      }|dgk(  sJ y )Nr   rN   )r`   r+   r'   r   r.   r/   rc   rd   s       r5   test_find_objects01z#TestFindObjects.test_find_objects01  s6    wwrw*""4(rd{{r7   c                 v    |j                  g |j                        }t        j                  |      }|g k(  sJ y )Nr   rh   r+   r'   r   r   s       r5   test_find_objects02z#TestFindObjects.test_find_objects02  s4    xx"((x+""4(byyr7   c                     |j                  dg|j                        }t        j                  |      }|t	        ddd       fgk(  sJ y )Nr
   r   r   r`   r+   r'   r   slicer   s       r5   test_find_objects03z#TestFindObjects.test_find_objects03  G    wws"((w+""4(aD)+,,,,r7   c                 x    |j                  dg|j                        }t        j                  |      }|g k(  sJ y )Nr
   r   r   r   s       r5   test_find_objects04z#TestFindObjects.test_find_objects04  s6    xx288x,""4(byyr7   c                     |j                  dg|j                        }t        j                  |      }|t	        ddd       fgk(  sJ y )Nrp   r   r   r   r   s       r5   test_find_objects05z#TestFindObjects.test_find_objects05  r   r7   c                     |j                  g d      }t        j                  |      }|t        ddd       ft        ddd       ft        ddd       fgk(  sJ y )Nrs   r   r
   r   r   rp   r   r,   r'   r   r   r   s       r5   test_find_objects06z#TestFindObjects.test_find_objects06  sd    zz,-""4(aD)+aD)+aD)+- - 	- -r7   c           	          |j                  g dg dg dg dg dg dg      }t        j                  |      }|g k(  sJ y )Nrw   )r,   r'   r   r   s       r5   test_find_objects07z#TestFindObjects.test_find_objects07  sE    zz------/ 0 ""4(byyr7   c           
      J   |j                  g dg dg dg dg dg dg      }t        j                  |      }|t        ddd       t        ddd       ft        ddd       t        d	d
d       ft        dd
d       t        dd	d       ft        d
dd       t        dd
d       fgk(  sJ y )Nrz   r   r   r   r   r   r
   rV   r   rp   r   r   r   s       r5   test_find_objects08z#TestFindObjects.test_find_objects08  s    zz------/ 0 ""4(aD)5At+<=!!Q-uQ4/@A!!Q-uQ4/@A!!Q-uQ4/@AC C 	C Cr7   c           
         |j                  g dg dg dg dg dg dg      }t        j                  |      }|t        ddd       t        ddd       ft        ddd       t        d	d
d       fd t        d
dd       t        dd
d       fgk(  sJ y )Nrz   r   r   rw   r   r   r
   rV   r   rp   r   r   r   s       r5   test_find_objects09z#TestFindObjects.test_find_objects09  s    zz------/ 0 ""4(aD)5At+<=!!Q-uQ4/@A!!Q-uQ4/@AC C 	C Cr7   N)rJ   rK   rL   r   r   r   r   r   r   r   r   r   r   rN   r7   r5   r   r     s6    $---CCr7   r   c           	      2   | j                  g dg dg dg dg dg dg      }t        j                  |d      }g d}t        |j	                               |k(  sJ |D ci c]  }|| j                  ||k(         }}t        j                  |d      }|j	                         |j	                         k(  sJ |j	                         D ]G  }t        ||         t        ||         k(  sJ t        ||   ||         D ]  \  }}t        ||        I y	c c}w )
z Test dictionary keys and entriesrz   r   r   rw   r   r   )ignore_value)r
   r   r   N)	r,   r'   value_indiceslistkeysnonzerorZ   zipr   )	r/   rc   vi	true_keysktruevikeyvtrue_vs	            r5   test_value_indices01r    s   ::))))))+ ,D 
		t!	4BI	?i'''09:1aDAI&&:F:			t!	4B779%%%wwy '2c7|s6#;////RWfSk2 	'IAvAv&	''	 ;s   !Dc                     | j                  d| j                        }d}t        t        |      5  t	        j
                  |       ddd       y# 1 sw Y   yxY w)zTest input checking)rp   r   r   z(Parameter 'arr' must be an integer array)matchN)rh   float32r   r   r'   r   )r/   rc   msgs      r5   test_value_indices02r    sJ    88F"**8-D
4C	z	- $d#$ $ $s   AAc                 `   dD ]#  }t        j                  ddgz  ddgz  z   ddgz  z   t         j                        }t        j                  ||      }t        j                  |      }| j                  |      }t        j                  |      }t        |j                               t        |      k(  sJ |D cg c]  }t        |       c}D ]U  }| j                  ||k(        }t        ||         t        |      k(  sJ t        ||   |      D ]  \  }}	t        ||	        W & yc c}w )z2Test different input array shapes, from 1-D to 4-D))r   )   r   )rV   rV   r   )rV   rV   r   r      r
   r   rV   r   N)r$   r,   r   r&   uniquer'   r   r   r   intr   rZ   r   r   )
r/   r   r   trueKeysr   r0   r   trueNdxviktrue_viks
             r5   test_value_indices03r  "  s   : /JJA3r1#vb!f,RXX>JJq% 99Q<JJqM""1%BGGI$x.000"*+Q#a&+ 	/Ajja(Gr!u:W---!$RUG!4 /XX./	// ,s   :D+c                     t         D ]=  }t        | |      }| j                  g |      }t        j                  |      }|dk(  r=J  y )Nr   r   )r   r   r,   r'   sumr/   r   r    inputr   s        r5   
test_sum01r  3  sI     D!

2U
+U#{{	r7   c                     t         D ]?  }t        | |      }| j                  ddg|      }t        j                  |      }|dk(  r?J  y )Nr   r   r   )r   r   rh   r'   r  r  s        r5   
test_sum02r  ;  sM     D!!Qu-U#{{	r7   c                     t         D ]S  }t        | |      }| j                  g |      }t        j                  |      }t        || j                  d      d       U y )Nr   r!   Fcheck_0d)r   r   r`   r'   r  r   r,   r  s        r5   
test_sum03r  C  sS     ED!%(U#FBJJsOeD	Er7   c                     t         D ]U  }t        | |      }| j                  ddg|      }t        j                  |      }t        || j                  d      d       W y )Nr
   r   r         @Fr  r   r   r,   r'   r  r   r  s        r5   
test_sum04r   K  sW     ED!

Aq6
/U#FBJJsOeD	Er7   c                     t         D ]Y  }t        | |      }| j                  ddgddgg|      }t        j                  |      }t        || j                  d      d       [ y )	Nr
   r   rV   r   r   g      $@Fr  r  r  s        r5   
test_sum05r"  S  sa     FD!

QFQF+5
9U#FBJJt$4uE	Fr7   c                     t        j                  g t              }| j                  |      }t        D ]?  }t	        | |      }| j                  g |      }t        j                  ||      }|dk(  r?J  y )Nr   r   r   r$   r,   boolr   r   r'   r  r/   r   r   r    r  r   s         r5   
test_sum06r(  [  si    ZZ$'FZZF D!

2U
+U62{{	r7   c                     t        j                  ddgt              }| j                  |      }t        D ]A  }t        | |      }| j                  ddg|      }t        j                  ||      }|dk(  rAJ  y )Nr   r   r   r$  )	r$   r`   r&  r,   r   r   rh   r'   r  r'  s         r5   
test_sum07r*  e  sq    WWaV4(FZZF D!!Qu-U62{{	r7   c                     t        j                  ddgt              }| j                  |      }t        D ]A  }t	        | |      }| j                  ddg|      }t        j                  ||      }|dk(  rAJ  y )Nr
   r   r   r   r$  r%  r'  s         r5   
test_sum08r,  o  sq    ZZAd+FZZF D!

Aq6
/U62{{	r7   c                 (   t        j                  ddgt              }| j                  |      }t        D ][  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }t        || j                  d      d	
       ] y )Nr
   r   r   r   rV   r   r$        @Fr  )r$   r,   r&  r   r   r'   r  r   r'  s         r5   
test_sum09r/  y  s    ZZAd+FZZF ED!

QFQF+5
9U62FBJJsOeD	Er7   c                 ,   t        j                  ddgt              }t        j                  ddgddggt              }| j                  |      }| j                  |      }t        j                  ||      }t        || j                  d      d	
       y Nr
   r   r   r   rV   r   r$         @Fr  )r$   r,   r&  r'   r  r   r/   r   r  r   s       r5   
test_sum10r4    ss    ZZAd+FJJAA't4EZZFJJuE[[v.F

3%@r7   c                    | j                  ddg| j                        }t        D ]\  }t        | |      }| j                  ddgddgg|      }t	        j
                  ||d      }t        || j                  d      d	       ^ y )
Nr
   r   r   rV   r   r         @Fr  )r,   int8r   r   r'   r  r   r'  s         r5   
test_sum11r8    s    ZZAbggZ.F ED!

QFQF+5
9U6#$&FBJJsOeDEr7   c           	      6   | j                  ddgddgg| j                        }t        D ]m  }t        | |      }| j                  ddgddgg|      }t	        j
                  ||| j                  g d            }t        || j                  g d             o y 	Nr
   r   r   r   rV   r      r   r   )r.  rG   g      @)r,   r7  r   r   r'   r  r   r'  s         r5   
test_sum12r=    s    ZZ!Q!Q(Z8F GD!

QFQF+5
9U6I9NO!&"**_*EF	Gr7   c           	         | j                  ddgddgg| j                        }t        D ]  }t        | |      }| j                  ddgddgg|      }t	        j
                  ||| j                  g d            }t	        j                  ||| j                  g d            }| j                  ||k(        sJ t        || j                  g d              y r:  )	r,   r7  r   r   r'   r  
sum_labelsallr   )r/   r   r   r    r  
output_sumoutput_labelss          r5   test_sum_labelsrC    s    ZZ!Q!Q(Z8F ND!

QFQF+5
9[[vRZZ	=RS
**&

9(=? vvjM1222!-O1LMNr7   c                 (   t        j                  ddgt              }| j                  |      }t        D ][  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }t        || j                  d      d	
       ] y r1  )r$   r,   r&  r   r   r'   meanr   r'  s         r5   test_mean01rF    s    ZZAd+FZZF ED!

QFQF+5
9eF3FBJJsOeD	Er7   c                 ,   t        j                  ddgt              }t        j                  ddgddggt              }| j                  |      }| j                  |      }t        j                  ||      }t        || j                  d      d	
       y Nr
   r   r   r   rV   r   r$  r!   Fr  )r$   r,   r&  r'   rE  r   r3  s       r5   test_mean02rI    ss    ZZAd+FJJAA't4EZZFJJuE\\%/F

3%@r7   c                     | j                  ddg      }t        D ]\  }t        | |      }| j                  ddgddgg|      }t        j                  ||d      }t        || j                  d      d	       ^ y )
Nr
   r   rV   r   r   r   r  Fr  )r,   r   r   r'   rE  r   r'  s         r5   test_mean03rK    sy    ZZAF ED!

QFQF+5
9eF$%'FBJJsOeDEr7   c           
         | j                  ddgddgg| j                        }t        j                  d      5  t        D ]z  }t        | |      }| j                  ddgddgg|      }t        j                  ||| j                  g d      	      }|d
   dk(  sJ |d   dk(  sJ | j                  |d         rzJ  	 d d d        y # 1 sw Y   y xY w)Nr
   r   r   r   ignorer@  rV   r;  r   r   r.        @)	r,   r7  r$   errstater   r   r'   rE  isnanr'  s         r5   test_mean04rR    s    ZZ!Q!Q(Z8F		" 
' 		'DB%EJJAA/uJ=E\\%(*

9(=?F !9###!9###88F1I&&&		'
' 
' 
's   B C;CCc                 (   t        j                  ddgt              }| j                  |      }t        D ][  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }t        || j                  d      d	
       ] y rH  )r$   r,   r&  r   r   r'   minimumr   r'  s         r5   test_minimum01rU        ZZAd+FZZF ED!

QFQF+5
9v6FBJJsOeD	Er7   c                 ,   t        j                  ddgt              }t        j                  ddgddggt              }| j                  |      }| j                  |      }t        j                  ||      }t        || j                  d      d	       y 
Nr
   r   r   r   r   r$  r!   Fr  )r$   r,   r&  r'   rT  r   r3  s       r5   test_minimum02rY    ss    ZZAd+FJJAA't4EZZFJJuE__U62F

3%@r7   c                     | j                  ddg      }t        D ]\  }t        | |      }| j                  ddgddgg|      }t        j                  ||d      }t        || j                  d      d	       ^ y )
Nr
   r   rV   r   r   r   r2  Fr  )r,   r   r   r'   rT  r   r'  s         r5   test_minimum03r[    sy    ZZAF ED!

QFQF+5
9v'(*FBJJsOeDEr7   c           	         | j                  ddgddgg      }t        D ]m  }t        | |      }| j                  ddgddgg|      }t        j                  ||| j                  g d            }t        || j                  g d             o y )	Nr
   r   rV   r   r   r   rV   r<  r   )r2  r.  rG   )r,   r   r   r'   rT  r   r'  s         r5   test_minimum04r^        ZZ!Q!Q()F GD!

QFQF+5
9v')zz)'<>!&"**_*EFGr7   c                 (   t        j                  ddgt              }| j                  |      }t        D ][  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }t        || j                  d      d	
       ] y )Nr
   r   r   r   rV   r   r$  r  Fr  )r$   r,   r&  r   r   r'   maximumr   r'  s         r5   test_maximum01rb    rV  r7   c                 ,   t        j                  ddgt              }t        j                  ddgddggt              }| j                  |      }| j                  |      }t        j                  ||      }t        || j                  d      d	       y rX  )r$   r,   r&  r'   ra  r   r3  s       r5   test_maximum02rd    ss    ZZAd+FJJAA't4EZZFJJuE__U62F

3%@r7   c                     | j                  ddg      }t        D ]\  }t        | |      }| j                  ddgddgg|      }t        j                  ||d      }t        || j                  d      d	       ^ y )
Nr
   r   rV   r   r   r   r.  Fr  )r,   r   r   r'   ra  r   r'  s         r5   test_maximum03rf    sy    ZZAF ED!

QFQF+5
9v'(*FBJJsOeDEr7   c           	         | j                  ddgddgg      }t        D ]m  }t        | |      }| j                  ddgddgg|      }t        j                  ||| j                  g d            }t        || j                  g d             o y )	Nr
   r   rV   r   r   r]  r   )r  r.  rG   )r,   r   r   r'   ra  r   r'  s         r5   test_maximum04rh     r_  r7   c                 ^    | j                  g d      }t        j                  |      dk(  sJ y )N)r   r   )r,   r'   ra  )r/   r0   s     r5   test_maximum05rl  *  s'    


< A??1###r7   c                    | j                  g dg dg dg dg      }| j                  g dg dg dg dg      }t        j                  ||| j                  g d      	      }t        || j                  g d
             y )Nr
   r   r   r
   rp   rV   r   r   r   r   r   r   r:   rV   r   r   r
   r
   r   r   r   r   r   r   rV   rV   r   r   )r
   r   rV   r   )rO  r.  r6  )r,   r'   medianr   r/   r   r   r   s       r5   test_median01rw  0  sp    


L   " 	#A ZZ%%%' (F ^^AfBJJy4IJFfbjj&ABr7   c                     | j                  g dg dg dg dg      }t        j                  |      }t        || j                  d      d       y )Nrn  ro  rp  rq  r!   Fr  r,   r'   ru  r   r/   r   r   s      r5   test_median02r{  =  sD    


L   " 	#A ^^AF

3%@r7   c                     | j                  g dg dg dg dg      }| j                  g dg dg dg dg      }t        j                  ||      }t        || j                  d	      d
       y )Nrn  ro  rp  rq  rr  rs  rt  r$  r  Fr  ry  rv  s       r5   test_median03r}  F  sf    


L   " 	#A ZZ%%%' (F ^^Af-F

3%@r7   c                    t        j                  ddgt              }| j                  |      }t        j                  || j                  d      | j                  dg            }t        || j                  dg             y )Nr
   r   r   r   r!   )r$   r,   r&  r'   ru  r`   r   rz  s      r5   test_median_gh12836_boolr  S  s\    


Aq6&A


1A^^Abggdm2::qc?KFfbjj#&78r7   c                     | j                  ddg| j                        }t        j                  || j	                  d      | j                  dg            }t        || j                  dg             y )NA   F   r   r  r
   r   g     P@)r,   r7  r'   ru  r`   r   rz  s      r5   test_median_no_int_overflowr  [  sU    


B8277
+A^^Abggdm2::qc?KFfbjj$&89r7   c                 j   t        j                  d      5  t        D ]r  }t        | |      }| j	                  g |      }t               5 }|j                  t        d       t        j                  |      }d d d        | j                        rrJ  	 d d d        y # 1 sw Y   (xY w# 1 sw Y   y xY wNrM  rN  r   zMean of empty slice)r$   rP  r   r   r,   r   filterRuntimeWarningr'   variancerQ  r/   r   r    r  supr   s         r5   test_variance01r  b  s    		" $ 	$DB%EJJrJ/E"$ 1

>+@A ))%01 88F###	$$ $1 1	$ $)   3B)
,B6B)B)B&"B))B2c                     t         D ]T  }t        | |      }| j                  dg|      }t        j                  |      }t        || j                  d      d       V y Nr
   r   rG   Fr  r   r   r,   r'   r  r   r  s        r5   test_variance02r  m  sW     ED!

A3e
,!!%(FBJJsOeD	Er7   c                     t         D ]U  }t        | |      }| j                  ddg|      }t        j                  |      }t        || j                  d      d       W y Nr
   rV   r   r!   Fr  r  r  s        r5   test_variance03r  u  sY     ED!

Aq6
/!!%(FBJJsOeD	Er7   c                     t        j                  ddgt              }| j                  |      }t        j                  |      }t        || j                  d      d       y )Nr
   r   r   g      ?Fr  )r$   r,   r&  r'   r  r   r/   r  r   s      r5   test_variance04r  }  sI    JJ1vT*EJJuEe$F

4 05Ar7   c                     | j                  g d      }t        D ]W  }t        | |      }| j                  g d|      }t        j                  ||d      }t        || j                  d      d       Y y N)r   r   rV   )r
   rV   r<  r   r   r!   Fr  )r,   r   r   r'   r  r   r'  s         r5   test_variance05r    sg    ZZ	"F ED!

9E
2!!%3FBJJsOeDEr7   c           
      d   | j                  g d      }t        j                  d      5  t        D ]h  }t	        | |      }| j                  g d|      }t        j                  ||| j                  g d            }t        || j                  g d             j 	 d d d        y # 1 sw Y   y xY wN)r   r   rV   rV   r   rM  rN  )r
   rV   r<  r   r<  r   )r   rV   r   )r!   r!   rG   )r,   r$   rP  r   r   r'   r  r   r'  s         r5   test_variance06r    s    ZZ(F		" K 	KDB%EJJ/uJ=E%%eVRZZ	5JKF%fbjj.IJ		KK K K   A2B&&B/c                 j   t        j                  d      5  t        D ]r  }t        | |      }| j	                  g |      }t               5 }|j                  t        d       t        j                  |      }d d d        | j                        rrJ  	 d d d        y # 1 sw Y   (xY w# 1 sw Y   y xY wr  )r$   rP  r   r   r,   r   r  r  r'   standard_deviationrQ  r  s         r5   test_standard_deviation01r    s    		" $ 	$DB%EJJrJ/E"$ ;

>+@A 33E:; 88F###	$$ $; ;	$ $r  c                     t         D ]T  }t        | |      }| j                  dg|      }t        j                  |      }t        || j                  d      d       V y r  r   r   r,   r'   r  r   r  s        r5   test_standard_deviation02r    sW     ED!

A3e
,++E2FBJJsOeD	Er7   c                     t         D ]U  }t        | |      }| j                  ddg|      }t        j                  |      }t        || j                  d      d       W y r  r  r  s        r5   test_standard_deviation03r    sY     ED!

Aq6
/++E2FBJJsOeD	Er7   c                     t        j                  ddgt              }| j                  |      }t        j                  |      }t        || j                  d      d       y )Nr
   r   r   r?   Fr  )r$   r,   r&  r'   r  r   r  s      r5   test_standard_deviation04r    sH    JJ1vT*EJJuE''.F

3%@r7   c                     | j                  g d      }t        D ]W  }t        | |      }| j                  g d|      }t        j                  ||d      }t        || j                  d      d       Y y r  )r,   r   r   r'   r  r   r'  s         r5   test_standard_deviation05r    sg    ZZ	"F ED!

9E
2++E61=FBJJsOeD	Er7   c           
      d   | j                  g d      }t        j                  d      5  t        D ]h  }t	        | |      }| j                  g d|      }t        j                  ||| j                  g d            }t        || j                  g d             j 	 d d d        y # 1 sw Y   y xY wr  )r,   r$   rP  r   r   r'   r  r   r'  s         r5   test_standard_deviation06r    s    ZZ(F		" K 	KDB%EJJ/uJ=E//vrzz)4F &fbjj.IJ	KK K Kr  c           	         | j                  dg      }t        j                  d      5  t        D ]  }t	        |       r|dk(  rt        j                  d       t        | |      }| j                  dg|      }t        j                  ||| j                  dg            }t        || j                  dg              	 d d d        y # 1 sw Y   y xY w)	Nr
   rM  rN  uint8z8value cannot be converted to type uint8 without overflowgl `yr   r   )r,   r$   rP  r   r   pytestxfailr   r'   r  r   r'  s         r5   test_standard_deviation07r    s    ZZ_F		" ? 	?D| 0 1B%EJJ}EJ:E//vrzz1#OF%fbjj!o>	?? ? ?s   BCCc                     t        j                  ddgt              }| j                  |      }t        D ]E  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }|dk(  rEJ  y )	Nr
   r   r   r   rV   r   r$  r   r   )r$   r,   r&  r   r   r'   minimum_positionr'  s         r5   test_minimum_position01r    s~    ZZAd+FZZF  D!

QFQF+5
9))%?	 r7   c                     t         D ]F  }t        | |      }| j                  g dg dg dg|      }t        j                  |      }|dk(  rFJ  y Nrp   r   r   rp   rV   r   r   r   r
   rp   r
   r
   r   r
   r   )r   r   r,   r'   r  r  s        r5   test_minimum_position02r    `      D!

L((*16  8 ))%0 r7   c                     t        j                  g dg dg dgt              }| j                  |      }t        j                  |      }|dk(  sJ y r  r$   r,   r&  r'   r  r  s      r5   test_minimum_position03r    N    JJ$$&-13E JJuE%%e,FVr7   c                     t        j                  g dg dg dgt              }| j                  |      }t        j                  |      }|dk(  sJ y )Nr  )rV   r   r
   r   r  r   r  r  r  s      r5   test_minimum_position04r    r  r7   c                     | j                  g d      }t        D ]G  }t        | |      }| j                  g dg dg dg|      }t        j                  ||      }|dk(  rGJ  y )Nr
   r   r   r   r  r  )r
   rp   r   rV   r   )r   r   r,   r   r   r'   r  r'  s         r5   test_minimum_position05r    p    ZZ%F  D!

L((*16  8 ))%8 r7   c                     | j                  g d      }t        D ]H  }t        | |      }| j                  g dg dg dg|      }t        j                  ||d      }|dk(  rHJ  y )Nr
   r   rV   r   r  r  r  r   r   r   r
   r  r'  s         r5   test_minimum_position06r    r    ZZ%F  D!

L((*16  8 ))%; r7   c           	         | j                  g d      }t        D ]f  }t        | |      }| j                  g dg dg dg|      }t        j                  ||| j                  ddg            }|d   d	k(  sJ |d
   dk(  rfJ  y )Nr  r  r  r  r   r   rV   r   r  r
   r  r  r'  s         r5   test_minimum_position07r        ZZ%F #D!

L((*16  8 ))%*,**aV*<>ayF"""ayF"""#r7   c                     t        j                  ddgt              }| j                  |      }t        D ]E  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }|dk(  rEJ  y )	Nr
   r   r   r   rV   r   r$  r
   r   )r$   r,   r&  r   r   r'   maximum_positionr'  s         r5   test_maximum_position01r  %  s    ZZAd+FZZF  D!

QFQF+5
9))%179 r7   c                     t         D ]F  }t        | |      }| j                  g dg dg dg|      }t        j                  |      }|dk(  rFJ  y )Nr  rV   r   r<  r   r  r   r  )r   r   r,   r'   r  r  s        r5   test_maximum_position02r  0  r  r7   c                     t        j                  g dg dg dgt              }| j                  |      }t        j                  |      }|dk(  sJ y )Nr  r  r  r   r  )r$   r,   r&  r'   r  r  s      r5   test_maximum_position03r  :  r  r7   c                     | j                  g d      }t        D ]G  }t        | |      }| j                  g dg dg dg|      }t        j                  ||      }|dk(  rGJ  y )Nr  r  r  r  r   r
   r
   r,   r   r   r'   r  r'  s         r5   test_maximum_position04r  C  r  r7   c                     | j                  g d      }t        D ]H  }t        | |      }| j                  g dg dg dg|      }t        j                  ||d      }|dk(  rHJ  y )Nr  r  r  r  r   r
   r  r  r'  s         r5   test_maximum_position05r  N  r  r7   c           	         | j                  g d      }t        D ]f  }t        | |      }| j                  g dg dg dg|      }t        j                  ||| j                  ddg            }|d   d	k(  sJ |d   d
k(  rfJ  y )Nr  r  r  r  r   r
   r   r   r  r  r  r'  s         r5   test_maximum_position06r  Y  r  r7   c           	         | j                  g d      }t        D ]f  }t        | |      }| j                  g dg dg dg|      }t        j                  ||| j                  ddg            }|d   d	k(  sJ |d
   dk(  rfJ  y )N)r!   rO  rG         @r  r  r  r   r!   r  r   r  r
   )r   rV   r  r'  s         r5   test_maximum_position07r  f  s    ZZ,-F #D!

L((*16  8 ))%*,**c3Z*@BayF"""ayF"""#r7   c                    t        j                  ddgt              }| j                  |      }t        D ]  }t	        | |      }| j                  ddgddgg|      }t        j                  ||      }t        j                  ||      }t        j                  ||      }t        j                  ||      }t        j                  ||      }	|||||	fk(  rJ  y )Nr
   r   r   r   rV   r   r$  )r$   r,   r&  r   r   r'   extremarT  ra  r  r  
r/   r   r   r    r  output1output2output3output4output5s
             r5   test_extrema01r  t  s    ZZAd+FZZF 
?D!

QFQF+5
9//%7//%7//%7**528:**528:7GWg>>>>
?r7   c                    | j                  ddg      }t        D ]  }t        | |      }| j                  ddgddgg|      }t        j                  ||d      }t        j
                  ||d      }t        j                  ||d      }t        j                  ||d      }t        j                  ||d      }	|||||	fk(  rJ  y )Nr
   r   rV   r   r   r   )	r,   r   r   r'   r  rT  ra  r  r  r  s
             r5   test_extrema02r    s    ZZAF ?D!

QFQF+5
9//%()+//%()+//%()+**528C**528C7GWg>>>>?r7   c           	         | j                  ddgddgg      }t        D ]K  }t        |       r|dv rt        j                  d       t        | |      }| j                  ddgddgg|      }t        j                  ||| j                  g d      	      }t        j                  ||| j                  g d      	      }t        j                  ||| j                  g d      	      }t        j                  ||| j                  g d      	      }t        j                  ||| j                  g d      	      }	t        |d
   |       t        |d   |       |d   |k(  sJ |d   |	k(  rLJ  y )Nr
   r   rV   uint16uint32uint64/https://github.com/pytorch/pytorch/issues/58734r   r   r]  r   r   r,   r   r   r  r  r   r'   r  rT  ra  r  r  r   r  s
             r5   test_extrema03r    s[   ZZ!Q!Q()F %B<D$BB\\KLD!

QFQF+5
9//%)/(*

9(=? //%)/(*

9(=? //%(*

9(=?**52813I1FH **52813I1FH 	"'!*g6!'!*g6qzW$$$qzW$$$/%r7   c           	         | j                  g d      }t        D ]I  }t        |       r|dv rt        j                  d       t        | |      }| j                  g dg dg dg|      }t        j                  ||| j                  dd	g            }t        j                  ||| j                  dd	g            }t        j                  ||| j                  dd	g            }t        j                  ||| j                  dd	g            }t        j                  ||| j                  dd	g            }	t        |d
   |       t        |d   |       |d	   |k(  sJ |d   |	k(  rJJ  y )Nr  r  r  r  r  r  r   r
   r   r   rV   r  r  s
             r5   test_extrema04r    sQ   ZZ%F %B<D$BB\\KLD!

L((*16  8 //%QF1CD//%QF1CD//%QF1CD**5&+-::q!f+=?**5&+-::q!f+=?!'!*g6!'!*g6qzW$$$qzW$$$%%r7   c                     d}t         D ]C  }t        | |      }| j                  ddgddgg|      }t        j                  |      }||k(  rCJ  y )N)rG   rG   r
   r   r   r   r   r,   r'   center_of_massr/   r   r   r    r  r   s         r5   test_center_of_mass01r    _    H "D!

QFQF+5
9''.!!!	"r7   c                     d}t         D ]C  }t        | |      }| j                  ddgddgg|      }t        j                  |      }||k(  rCJ  y )Nr  r   r
   r   r  r  s         r5   test_center_of_mass02r    _    H "D!

QFQF+5
9''.!!!	"r7   c                     d}t         D ]C  }t        | |      }| j                  ddgddgg|      }t        j                  |      }||k(  rCJ  y )Nr  r   r
   r   r  r  s         r5   test_center_of_mass03r    r  r7   c                     d}t         D ]C  }t        | |      }| j                  ddgddgg|      }t        j                  |      }||k(  rCJ  y )Nr  r   r
   r   r  r  s         r5   test_center_of_mass04r    r  r7   c                     d}t         D ]C  }t        | |      }| j                  ddgddgg|      }t        j                  |      }||k(  rCJ  y )Nr?   r?   r
   r   r  r  s         r5   test_center_of_mass05r	    r  r7   c                     d}t        j                  ddgddggt              }| j                  |      }t        j                  |      }||k(  sJ y )Nr  r
   r   rV   r   )r$   r,   r&  r'   r  r/   r   r  r   s       r5   test_center_of_mass06r    sP    HJJAA't4EJJuE##E*FXr7   c                     | j                  ddg      }d}t        j                   ddgddggt              }| j                  |      }t        j                  ||      }||k(  sJ y )Nr
   r   r?   rG   r   rV   r   r,   r$   r&  r'   r  r/   r   r   r  r   s        r5   test_center_of_mass07r    sd    ZZAFHJJAA't4EJJuE##E62FXr7   c                     | j                  ddg      }d}t        j                   ddgddggt              }| j                  |      }t        j                  ||d      }||k(  sJ y )Nr
   r   r?   r!   rp   rV   r   r  r  s        r5   test_center_of_mass08r    sf    ZZAFHJJAA't4EJJuE##E615FXr7   c                 l   | j                  d      }| j                  ddg| j                        }t        j                   ddgddggt              }| j                  |      }t	        j
                  ||| j                  ddg            }t        | j                  |      | j                  |             y )Nr  r  r  r   r
   r   )r,   float64r$   r&  r'   r  r   r  s        r5   test_center_of_mass09r    s    ZZFzz:z2"**zEHJJAA't4EJJuE##E62::q!f3EFFBJJv&

8(<=r7   c                     | j                  d      }| j                  d      }t        j                  |ddd      }t	        ||       y )Nr   r   )r`   aranger'   	histogramr   r  s       r5   test_histogram01r    s<    wwr{HIIbMEuaR0Ffh/r7   c                     | j                  g d      }| j                  g d      }| j                  g d      }t        j                  |ddd|d      }t        ||       y )N)r
   r
   r
   r
   r   r   r   r   )r   r   r   r
   r
   )r
   r
   rV   r   rV   rV   rV   rV   r   r   rp   r
   r,   r'   r  r   r  s        r5   test_histogram02r     sR    ZZ01Fzz/*HJJ/0EuaAvq9Ffh/r7   zobject arraysc                    | j                  g d      }| j                  g d      }| j                  g d      }| j                  g d      }t        j                  |ddd|d      }t        |d   |       t        |d	   |       y )
N)r
   r   r
   r
   r   r   r   r   )r   r
   r   r
   r
   )r   r   r   rV   r   )r
   r
   rV   r   rV   rp   rV   rV   r   r   rp   r  r
   r  )r/   r   	expected1	expected2r  r   s         r5   test_histogram03r"  (  sr    ZZ01F

?+I

?+IJJ/0EuaAvv>FfQi3fQi3r7   c           	      :   | j                  g dg dg dg      }| j                  g dg dg dg      }t        j                  ||| j                  ddg            }t        || j                  d	d
g| j                               t        j
                  ||| j                  ddg            }t        || j                  ddg| j                               t        j                  ||| j                  ddg            }t        || j                  | j                  ddg| j                                     t        j                  ||| j                  ddg            }t        || j                  d	d
g| j                               t        j                  ||| j                  ddg            }t        || j                  ddg      d       t        j                  ||| j                  ddg            }t        || j                  ddg      d       y )N)rp   r   r   r   r   )r<  r:   r   r   r   )r   r   r   rV   rp   )r
   r
   r   r   r   )r   r   r   r   r   r
   r   r   rR   r.  r   rO  r!   rp   rV   Fr   r:   )r,   r'   rE  r   r  r  r  r   sqrtru  rT  ra  )	r/   r   lblrE  varstdmedminmaxs	            r5   test_stat_funcs_2dr+  5  s   


O_oFGA
**oH
IC<<#RZZA-?@DD"**c3Zrzz*BC


1S

Aq60B
CCCS#JbjjAB

$
$Qs"**aV:L
MCc2772::sCj

:+S#TU
..3bjj!Q.@
ACCS#JbjjAB
//!Crzz1a&/A
BCCQF+?
//!Crzz1a&/A
BCCQF+?r7   no watershed_ift on CuPyc                       e Zd Zd Zd Zd Zd Zd Zd Z e	dd	      d
        Z
 e	dd      d        Z e	dd      d        Zy)TestWatershedIftc                    |j                  g dg dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg      }t        j                  |||      }g d	g d
g d
g d
g d
g d
g d	g d	g}t        ||j                  |             y )Nr   r   r   r   r   r   r   r   r
   r
   r
   r
   r
   r   r   r
   r   r   r   r
   r   r   r   r   r   r   r   r   r   r   r   r   r
   r   r   r   r
   r
   r
   	structurer   r   r   r   r   r   r   r   r
   r
   r
   r
   r
   r   r,   r  r7  r'   watershed_iftr   r.   r/   rc   markersr7  rd   r   s          r5   test_watershed_ift01z%TestWatershedIft.test_watershed_ift01O  s    zz000000002 :<  C **433333335 =?GG  E **i'') *	 ##D'YG0+++++002 	"#rzz(';<r7   c                 n   |j                  g dg dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg dg dg|j                        }t        j                  ||      }g dg dg d	g d	g d	g dg dg dg}t        ||j                  |             y )
Nr0  r1  r2  r   r3  r4  r8  )r   r   r
   r
   r
   r   r   r9  r:  r.   r/   rc   r=  rd   r   s         r5   test_watershed_ift02z%TestWatershedIft.test_watershed_ift02n  s    zz000000002 :<  C **433333335 =?GG  E ##D'20-+++-002 	"#rzz(';<r7   c           
      \   |j                  g dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg dg|j                        }t        j                  ||      }g dg dg d	g d	g d	g dg dg}t        ||j                  |             y )
Nr0  r1  r   r
   r   r
   r   r
   r   r   r   r   r   r   rV   r   r   r   r   r   r   r   r   r   r8  )r   r   r   r   rV   r   r   r   r   r   rV   rV   rV   r   r:  r@  s         r5   test_watershed_ift03z%TestWatershedIft.test_watershed_ift03  s    zz00000002 :<  C **33333346 >@WW  F ##D'20.+++.02 	"#rzz(';<r7   c           
         |j                  g dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg      }t        j                  |||      }g d	g d
g d
g d
g d
g d
g d	g}t        ||j                  |             y )Nr0  r1  rC  r   rD  rE  r5  r6  r8  rF  r:  r<  s          r5   test_watershed_ift04z%TestWatershedIft.test_watershed_ift04  s    zz00000002 :<  C **33333346 $&77  , **i'') *	 ##D'YG0+++++02 	"#rzz(';<r7   c           
         |j                  g dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg dg|j                        }|j                  g dg dg dg      }t        j                  |||      }g d	g d
g d
g d
g d
g d
g d	g}t        ||j                  |             y )Nr0  r1  rC  r   )r   r   rV   r   r   r   r   rE  r5  r6  r8  )r   rV   rV   r   r   r   r   r:  r<  s          r5   test_watershed_ift05z%TestWatershedIft.test_watershed_ift05  s    zz00000002 :<  C **33333346 $&77  , JJ	 ) ) + ,	 ##D'YG0+++++02 	"#rzz(';<r7   c           	         |j                  g dg dg dg dg dg dg|j                        }|j                  g dg dg dg dg dg dg|j                        }|j                  g dg dg dg      }t        j                  |||      }g d	g d	g d	g d	g d
g d
g}t        ||j                  |             y )Nr2  r1  r0  r   r3  r4  r5  r6  r9  r8  r:  r<  s          r5   test_watershed_ift06z%TestWatershedIft.test_watershed_ift06  s    zz0000002
 :<  C **4333335
 =?GG  E **i'') *	 ##D'YG++++002 	"#rzz(';<r7   Tzinplace ops are numpy-specificr   c           	      |   d}t        j                  |t         j                        }|j                         }t        j                  g dg dg dg dg dg dgt         j                        |d<   |j	                  |      }|j	                  g dg dg dg dg dg dg|j
                        }|j                  ||j                        }|j                  }|j	                  g d	g d	g d	g      }t        j                  ||||
       g dg dg dg dg dg dg}t        ||j	                  |             y )N)r   r   r   r2  r1  r0  .r3  r4  r5  )r7  r   r9  r8  )r$   rh   r  	transposer,   r7  int16Tr'   r;  r   )r.   r/   r   rc   r=  rd   r7  r   s           r5   test_watershed_ift07z%TestWatershedIft.test_watershed_ift07  s"   xxRXX.~~JJ 5 5 5 5 5 5 7
 ?AhhHS	 zz$**4333335
 =?GG  E hhuBHHh-ee**i'') *	 	dGy%(	*++++002 	"#rzz(';<r7   r   r,  r   c                    |j                  ddgddgg|j                        }|j                  ddgddgg|j                        }t        j                  ||      }ddgddgg}t        ||j                  |             y )N   r   r   r
   )r,   r  r7  r'   r;  r   r@  s         r5   test_watershed_ift08z%TestWatershedIft.test_watershed_ift08  s     zzC8q6#*,))  5**q!f !f&-/WW  6##D'2FF!#rzz(';<r7   c                 L   |j                  |j                  |j                        j                  dgddgg|j                        }|j                  ddgddgg|j                        }t        j                  ||      }ddgddgg}t        ||j                  |      d       y )Nr   r   r
   Fr   )r,   iinfor  r*  r7  r'   r;  r   r@  s         r5   test_watershed_ift09z%TestWatershedIft.test_watershed_ift09$  s     zzBHHRYY/33Q7q6#*,))  5**q!f !f&-/WW  6##D'2FFRZZ1uEr7   N)rJ   rK   rL   r>  rA  rG  rI  rK  rM  skip_xp_backendsrR  rU  rX  rN   r7   r5   r.  r.  L  sz    =>=8=2=<=:=2 d+KL= M=@ f%?@	= A	= f%?@	F A	Fr7   r.  )r   dtc                     t         j                  j                  d      }d}|j                  d|d      j	                  |       }t        j                  |      }t        |j                               |k(  sJ y )N{   r<  r   )r   r
  )lowhighr   r   )	r$   r   default_rngintegersastyper'   r   rZ   r   )rZ  r/   rngmax_valimageval_idxs         r5   test_gh_19423rf  1  sh     ))


$CGLLQW8L<CC"CME##E*Gw||~')))r7   )r   os.pathnumpyr$   numpy.testingr   scipy._lib._array_apir   r   r   r   r   r  r	   r   scipy.ndimager'    r   markrY  
pytestmarkname__version__r*   r   rP   rf   ri   rl   rn   rq   ru   rx   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r  r  r  r  r   r"  r(  r*  r,  r/  r4  r8  r=  rC  rF  rI  rK  rR  rU  rY  r[  r^  rb  rd  rf  rh  rl  rw  r{  r}  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r	  r  r  r  r  r  r  r"  r+  r.  parametrizeintcuintcrf  rN   r7   r5   <module>rt     s   	   +   *  ;;// &+9NOP
WW_=#)=  $'IJO= O= KO=d $'IJ#/ #/ K#/L$"&* $'IJ K*&$ $F89;; $F89;!;! $'JKF LF, &!UV W$	5 &!PQTC TC RTCn',$/"EEFEAEG
NEAE'EAEGEAEG$
CA
A9:$EEBEK$EEAE	K
?    
#    
##? ?$%8%.""""">00 $7	4 8	4@. &!;<aF aF =aFH $23* 4  *r7   