ó
1™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   sasaki(   t   approx_equalc         C   sI   |  j  | ƒ } t | j ƒ  | ƒ s* t ‚ t | j ƒ  | ƒ sE t ‚ d  S(   N(   t   at_angstromR   t   fpt   AssertionErrort   fdp(   t   tablet   wave_lengthR   R   t   fp_fdp(    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_sasaki.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 ƒ t  j d ƒ }  |  j ƒ  d k söt ‚ |  j ƒ  d k st ‚ t |  d d d ƒ 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 smt ‚ qmW| d# k st ‚ d  S(%   Nt   SIt   Sii   i   iÿÿÿÿ(   t   factor_kev_angstromiè  gš™™™™™¹?gMóŽSt$—¿gü©ñÒMbP?g…ëQ¸@g·bÙ=yØ?gO@aÃÓð?gµ7øÂdê@g?5^ºIÌÀg	ù g³j@t   agt   Agi/   g©ÐDØð	@g9´Èv¾_!ÀgŸ«­Ø_6,@gQkšwœ"@g¥½Á&ƒ0ÀgÏfÕçj‹+@gj¼t“˜@gyX¨5Íû=Àgé·¯%@g¸@‚âÇ˜@g’\þCú%@ÀgcîZB>h	@i    i   t   LiiR   t   U(   R   R   t   labelR   t   atomic_numberR   t   is_valid_fpt   is_valid_fdpt   is_validt   cctbxR   R   R   t   at_kevR   t   at_evt
   as_complext   realt   imagR
   t   table_iterator(   t   tt   fR   t   ct   nt   u(    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_sasaki.pyt   exercise
   sF    ..22
"c           C   s   t  ƒ  d GHd  S(   Nt   OK(   R#   (    (    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_sasaki.pyt   run/   s    t   __main__N(
   t
   __future__R    t   cctbx.eltbxR   t   libtbx.test_utilsR   R
   R#   R%   t   __name__(    (    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/tests/tst_sasaki.pyt   <module>   s   		%	