ó
]xTc           @   sS   d  d l  m Z d  d l m Z d  d l m Z d   Z e d k rO e   n  d S(   i˙˙˙˙(   t   division(   t   fit_quadratic_function(   t   approx_equalc          C   s  g  }  g  } d } d } d } d } d } g  } x t  d  D] } x t  d  D]y }	 |  j | d  | j |	 d  | j | | d | d	 | |	 d | d	 d	 | | d | |	 d |  qP Wq= Wt j |  | |  }
 t |
 j | d
 d st  t |
 j | d
 d s"t  t |
 j | d
 d s@t  t |
 j	 | d
 d s^t  t |
 j
 | d
 d s|t  d GHd  S(   Ng      đ?g      @i   i   g      ŕżi
   i   g      @g       @t   epsgńhăľřä>t   OK(   t   ranget   appendt   fqft   fit_quadratic_function_2d_dataR   t   x1mt   AssertionErrort   x2mt   at   bt   c(   t   x1_obst   x2_obsR	   R   R   R   R   t   y_obst   iit   jjt   fit(    (    sZ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_fit_quadratic_function.pyt   test_fit   s,    ',t   __main__N(	   t
   __future__R    t   scitbx.mathR   R   t   libtbx.test_utilsR   R   t   __name__(    (    (    sZ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_fit_quadratic_function.pyt   <module>   s
   	