ó
2™xTc           @   se   d  d l  m Z d  d l m Z d  d l m Z d „  Z d „  Z d „  Z e	 d k ra e ƒ  n  d S(	   iÿÿÿÿ(   t   division(   t   henke(   t   approx_equalc         C   sI   |  j  | ƒ } t | j ƒ  | ƒ s* t ‚ t | j ƒ  | ƒ sE t ‚ d  S(   N(   t   at_evR   t   fpt   AssertionErrort   fdp(   t   tablet   energyR   R   t   fp_fdp(    (    sI   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_henke.pyt   verify   s    c          C   s�  t  j d ƒ }  |  j ƒ  d k s' t ‚ |  j ƒ  d k s? t ‚ |  j d ƒ } | j ƒ  s` t ‚ | j ƒ  sr t ‚ | j ƒ  s„ t ‚ d d l	 m
 } t | j ƒ  |  j | d ƒ j ƒ  ƒ sÂ t ‚ t | j ƒ  |  j | d ƒ j ƒ  ƒ sð t ‚ t | j ƒ  |  j d | d ƒ j ƒ  ƒ s"t ‚ t | j ƒ  |  j d | d ƒ j ƒ  ƒ sTt ‚ | j ƒ  } | j | j ƒ  k s{t ‚ | j | j ƒ  k s–t ‚ t |  d d	 d
 ƒ t |  d d d ƒ t |  d d d ƒ d } xŸ t  j ƒ  D]‘ }  | d 7} | d k r|  j ƒ  d k s@t ‚ n' | d k r@|  j ƒ  d k s@t ‚ n  t  j |  j ƒ  ƒ } | j ƒ  |  j ƒ  k sât ‚ qâW| d k s‰t ‚ d  S(   Nt   SIt   Sii   i   iÿÿÿÿ(   t   factor_kev_angstromiè  g      $@g    €‡ÃÀgeüûŒ@gÍÌÌÌÌL=@g˜†á#b*@g,-#õžÊ×?g     LÝ@g¯%äƒž,@gçdeêd—?i    i   t   Hi\   t   Ug´<îÎê#Àg ^KÈ=›?(   R   R   t   labelR   t   atomic_numbert   at_angstromt   is_valid_fpt   is_valid_fdpt   is_validt   cctbxR   R   R   t   at_kevR   R   t
   as_complext   realt   imagR
   t   table_iterator(   t   tt   fR   t   ct   nt   u(    (    sI   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_henke.pyt   exercise
   s8    ..22
"c           C   s   t  ƒ  d GHd  S(   Nt   OK(   R!   (    (    (    sI   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_henke.pyt   run(   s    t   __main__N(
   t
   __future__R    t   cctbx.eltbxR   t   libtbx.test_utilsR   R
   R!   R#   t   __name__(    (    (    sI   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_henke.pyt   <module>   s   			