
    tKg                     8    d dl mZmZ d dlmZ d Zd Zd Zd Zy)    )assert_array_equalassert_array_almost_equal)padec                      t        dgd      \  } }t        | j                  dg       t        |j                  dg       t        dgdd      \  } }t        | j                  dg       t        |j                  dg       y )N      ?r   )r   r   c)numpdenomps     e/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/scipy/interpolate/tests/test_pade.pytest_pade_trivialr      sg    q>LD&tvvu%vxx#'q!$LD&tvvu%vxx#'    c                  x   g d} t        | d      \  }}t        |j                  g d       t        |j                  dg       t        | d      \  }}t        |j                  g d       t        |j                  ddg       t        | d      \  }}t        |j                  d	dg       t        |j                  g d
       t        | d      \  }}t        |j                  dg       t        |j                  g d       t        | dd      \  }}t        |j                  g d       t        |j                  dg       t        | dd      \  }}t        |j                  g d       t        |j                  ddg       t        | dd      \  }}t        |j                  d	dg       t        |j                  g d
       t        | dd      \  }}t        |j                  dg       t        |j                  g d       t        | dd      \  }}t        |j                  g d       t        |j                  dg       t        | dd      \  }}t        |j                  ddg       t        |j                  ddg       t        | dd      \  }}t        |j                  dg       t        |j                  g d       y )N)r   r         ?UUUUUU?r   )r   r   r   r   r      )r   gUUUUUU?r   gUUUUUUտ   gUUUUUU?)r   gUUUUUUr      )gUUUUUUſr         r   )r   r   r   r   g      )r   r   r   )r   r   r   )anr	   r
   s      r   test_pade_4term_expr      s!    
 BA;LD&dff&<=fhh.A;LD&dff&9:fhh6A;LD&dffucl3fhh(<=A;LD&dffse,fhh(@A Aq>LD&dff&<=fhh.Aq>LD&dff&9:fhh6Aq>LD&dffucl3fhh(<=Aq>LD&dffse,fhh(@A Aq>LD&dffo6fhh.Aq>LD&dffucl3fhh6Aq>LD&dffse,fhh(:;r   c                  B   g d} g d}t        dt        |             D ]  }t        dt        |       |z
        D ]b  }t        | ||      \  }}t        |||      \  }}t        |j                  |j                         t        |j                  |j                         d  y )N)r   r   r      )r   g       @g      @g      @r   )rangelenr   r   r   )an_intan_fltijnump_int
denomp_intnump_flt
denomp_flts           r   test_pade_intsr#   B   s    F!F 1c&k"q#f+/*A $(1#5 Hj#'1#5 Hj xzz8::6z||Z\\: + #r   c            	         d} d| |  | j                         z  | j                         | dz  z  | | j                         dz  z  z   | dz   | j                         z  d| | j                         z  dz  z  z
  | | j                         dz  z  z
  g}t        |dd      \  }}t        |j                  | | j                         z   dg       t        |j                  | j                         dg       t        |dd      \  }}t        |j                  | dz  d| z  | j                         z   dg       t        |j                  | | j                         z   dg       t        |dd      \  }}t        |j                  | dz  | | j                         z  z   | j                         dz  z   d| | j                         z   z  dg       t        |j                  | j                         dz  | d| j                         z  z   dg       y )Ny?333333?r   r   r   r   )	conjugater   r   r   )xr   r	   r
   s       r   test_pade_complexr'   T   s    	A
q1"Q[[]"AKKM1a4$81akkmQ>N;O$OqD'!++-
!Qq{{}_q$8"8
81akkmQ>N;O
O
QB Aq>LD&dffq1;;='8#&>?fhh(<=Aq>LD&dffq!tQqS1;;=-@#&FGfhhQ[[]):C(@AAq>LD&	
A!++-	!++-"2	2Aq1;;=7H4I3O fhh)91q;NPS(TUr   N)	numpy.testingr   r   scipy.interpolater   r   r   r#   r'    r   r   <module>r+      s!    I "(1<h;$Vr   