Ë
    qwj8	  ã                   óö   — d Z ddlZddlZddlmZmZ ddlmZ d„ Z	dd„Z
 eej                  «      ej                  j                  dg d¢«      ej                  j                  d	g d
¢«      d„ «       «       «       Zy)zTests for spline filtering.é    N)Úassert_almost_equalÚmake_xp_test_case)Úndimagec                 ó2   — dgdgddgddgg d¢g d¢dœ}||    S )z0Knot values to the right of a B-spline's center.é   é   é   )éæ   éL   r   )éB   é   r   ©r   r   é   é   r	   é   © )ÚorderÚknot_valuess     úe/var/www/html/newmanjeet/manjet/venv/lib/python3.12/site-packages/scipy/ndimage/tests/test_splines.pyÚget_spline_knot_valuesr   
   s4   € à�cØ�cØ˜!�fØ˜!�fÚ"Ú!ñ#€Kð �uÑÐó    c                 ó\  — t        |«      }t        j                  ||f«      }t        |«      D ]<  \  }}t        j                  ||«      }|dk(  r||||f<   Œ)|||||z
  f<   ||||z
  |f<   Œ> |d   dt        |dd «      z  z   }	|dk(  rd\  }
}n$|dk(  rd\  }
}n|d	k(  rd
\  }
}nt        d|› �«      ‚t        t        |«      dz
  «      D ]M  }t        ||dz   d «      D ]7  \  }}|||
||z  z   fxx   |z  cc<   || dz
  |
 dz
  ||z  z
  fxx   |z  cc<   Œ9 ŒO | j                  ||	z  «      S )z1Matrix to invert to find the spline coefficients.r   r   r   NÚmirror)r   r   Úreflect)r   r   ú	grid-wrap)éÿÿÿÿr   zunsupported mode )
r   ÚnpÚzerosÚ	enumerateÚarangeÚsumÚ
ValueErrorÚrangeÚlenÚasarray)ÚxpÚnr   Úmoder   ÚmatrixÚdiagÚ
knot_valueÚindicesÚknot_values_sumÚstartÚstepÚrowÚidxs                 r   Úmake_spline_knot_matrixr2      s|  € ä(¨Ó/€Kô �X‰X�q˜!�fÓ€FÜ% kÖ2ÑˆˆjÜ—)‘)˜D !Ó$ˆØ�1Š9Ø'1ˆF�7˜GÐ#Ò$à.8ˆF�7˜G d™NÐ*Ñ+Ø.8ˆF�7˜T‘> 7Ð*Ò+ð 3ð " !‘n q¬3¨{¸1¸2¨Ó+?Ñ'?Ñ?€OàˆxÒØ‰ˆ‰tØ	�Ò	Ø‰ˆ‰tØ	�Ò	Ø‰ˆ‰täÐ,¨T¨FÐ3Ó4Ð4ä”S˜Ó%¨Ñ)Ö*ˆÜ(¨°S¸1±W°XÐ)>Ö?‰OˆC�Ø�3˜  S¡Ñ(Ð(Ó)¨ZÑ7Ó)Ø�C�4˜!‘8˜e˜V a™Z¨$¨s©(Ñ2Ð2Ó3°zÑAÔ3ñ  @ð +ð
 �:‰:�f˜Ñ.Ó/Ð/r   r   r   r(   )r   r   r   c                 ó  — d}|j                  ||j                  ¬«      }t        j                  |d| |¬«      }t        j                  |d| |¬«      }t	        ||| |¬«      }t        |||z  «       t        |||j                  z  «       y )Néd   )Údtyper   )Úaxisr   r(   r   )r(   )ÚeyeÚfloat64r   Úspline_filter1dr2   r   ÚT)r   r(   r&   r'   r7   Úspline_filter_axis_0Úspline_filter_axis_1r)   s           r   Ú%test_spline_filter_vs_matrix_solutionr=   8   s‡   € ð 	€AØ
�&‰&�˜"Ÿ*™*ˆ&Ó
%€CÜ"×2Ñ2°3¸QÀeØ8<ô>Ðä"×2Ñ2°3¸QÀeØ8<ô>Ðä$ R¨¨E¸Ô=€FÜ˜Ð1°FÑ:Ô;Ü˜Ð1°F·H±HÑ<Õ=r   )r   )Ú__doc__ÚpytestÚnumpyr   Úscipy._lib._array_apir   r   Úscipyr   r   r2   r9   ÚmarkÚparametrizer=   r   r   r   Ú<module>rE      ss   ðÙ !Û ã ß Hå ò	ó0ñD �7×*Ñ*Ó+Ø‡�×Ñ˜Ò"4Ó5Ø‡�×Ñ˜Ò!CÓDñ	>ó Eó 6ó ,ñ	>r   