ó
ª™xTc           @   s„   d  d l  m Z d  d l m Z d  d l m Z d  d l Z d  d l Z d  d l m	 Z	 d „  Z
 d „  Z e d k r€ e ƒ  n  d S(	   iÿÿÿÿ(   t   division(   t   scaling(   t   flexN(   t   approx_equalc          C   sé  t  j d g ƒ }  t  j d g ƒ } t  j d g ƒ } t  j t g ƒ } t j d |  d | d | d | d d	 ƒ } t | j d d ƒ d ƒ s“ t ‚ t | j	 d d ƒ d
 d d ƒsº t ‚ t | j
 d d ƒ d | j	 d d ƒ ƒ së t ‚ t | j
 d d ƒ d | j d d ƒ ƒ st ‚ d } t j ƒ  } x4 t | ƒ D]& } | j	 d d ƒ | j d d ƒ q;Wt j ƒ  } | | GHt j ƒ  } x$ t | ƒ D] } | j
 d d ƒ q“Wt j ƒ  } | | GHt j d |  d | d | d | d d	 ƒ }	 t | j d d ƒ d ƒ st ‚ t |	 j	 d d ƒ d d d ƒs5t ‚ t |	 j
 d d ƒ d |	 j	 d d ƒ ƒ sft ‚ t |	 j
 d d ƒ d |	 j d d ƒ ƒ s—t ‚ d } t  j d g | ƒ }  t  j d g | ƒ } t  j d g | ƒ } t  j t g | ƒ } t j d |  d | d | d | d d	 ƒ }
 t j ƒ  } x8 t d ƒ D]* } |
 j	 d d ƒ } |
 j d d ƒ } q8Wt j ƒ  } | | GHt j d |  d | d | d | d d	 ƒ }
 t j ƒ  } x& t d ƒ D] } |
 j
 d d ƒ } q¾Wt j ƒ  } | | GHd } d } xç t  j t d ƒ ƒ d D]Ì } | j	 | | ƒ } | j	 | | | ƒ } | j | | ƒ } | | | } |	 j	 | | ƒ } |	 j	 | | | ƒ } |	 j | | ƒ } | | | } t | | d d ƒsÆt ‚ t | | d d ƒst ‚ qWd  S(   Ni   g333333ó?g      à?t   e_obst   e_calct   centrict   d_star_cubedt   widthgš™™™™™¹?gT^-wÑ¿t   epsg-Cëâ6?i    i † g‘ÔBÉäTæ¿id   g�íµ ÷Æ°>i   g      4@gü©ñÒMbP?(   R   t   doublet   boolt   FalseR   t   sigmaa_estimatorR   t   sum_weightst   AssertionErrort   targett   target_and_gradientt   dtargett   timet   xranget   range(   t   eot   ect   dscR   t   tmp_at   Nt   startt   iit   endt   tmp_ct	   tmp_larget   trialt   at   ht   dt   fat   fa1t   gat   ga1t   fct   fc1t   gct   gc1(    (    sF   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/tst_sigmaa.pyt
   tst_sigmaa   sš    		!'11				!'11						 c           C   s   t  ƒ  d GHd  S(   Nt   OK(   R,   (    (    (    sF   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/tst_sigmaa.pyt   runp   s    t   __main__(   t
   __future__R    t   mmtbxR   t   cctbx.array_familyR   t   mmtbx.scalingR   t   libtbx.test_utilsR   R,   R.   t   __name__(    (    (    sF   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/tst_sigmaa.pyt   <module>   s   	h	