ó
N™xTc           @   sL   d  d l  m Z d „  Z e d k rH d  d l Z e d e j d ƒ n  d S(   iÿÿÿÿ(   t   divisionc         C   sð  t  |  ƒ d k s t ‚ d d l m } d d l m } | j df dg dh di dj dk g ƒ } | j d d d d d d g ƒ } | d | d | ƒ } g  | D] } t | ƒ ^ q– d g d d  g d d! g d  g g  g  g k sÞ t ‚ | j dl dm dn do dp dq dr ds dt du dv dw dx dy dz d{ d| d} d~ d d€ d� d‚ dƒ g ƒ } | j d^ d_ d^ d_ d^ d_ d^ d^ d` d_ d^ d^ d` d^ da da da da da da da da da da g ƒ } | d | d | ƒ } | j	 ƒ  | j	 ƒ  k s¿t ‚ t
 | d ƒ d db dc dd g k sçt ‚ de GHd  S(„   Ni    iÿÿÿÿ(   t   build_simple_two_way_bond_sets(   t   flexgHáz®§9@g5^ºI¢E@gq=
×£PE@g}?5^º	8@gw¾Ÿ/íE@gB`åÐ"«E@g?5^ºI7@g33333CF@g…ëQ¸¥E@g+‡95@gÑ"Ûù~F@gÍÌÌÌÌ¼D@gú~j¼tó8@gö(\�Â•G@g{®GáÊB@g?5^ºIL3@gj¼t“xG@g�•C‹lçI@t   St   Ct   Nt   CUt   ZNt   CAt
   sites_cartt   elementsi   i   i   gÙÎ÷SãÀg¸…ëQ¸Ž?g}?5^ºIú?g¸…ëQ¸ø¿gÛù~j¼tç¿gœÄ °rh@g‹lçû©ñâ?g)\�Âõ(@g-²�ï§Æã?g%�•C‹ð?gTã¥›Ä ü?gé&1¬º?g#Ûù~j¼@g¸…ëQ¸®¿gh‘í|?5æ?gffffffø?gÃõ(\�Âõ¿g¶óýÔxé@gú~j¼t“˜¿gé&1¬ÀgÓMbX9Ä?gh‘í|?5þ¿gš™™™™™™¿gú~j¼t“˜?g¤p=
×£
ÀgZd;ßO�ï¿gã¥›Ä °î¿g=
×£p=Àgü©ñÒMb€?gºI+‡ö¿g®Gázâ?gøSã¥›Ä
Àgü©ñÒMb�?g²�ï§ÆK@g-²�ï§Æ
Àgyé&1¬ì?gTã¥›Ä ò?g-²�ï§Æë¿g“V-ò?g‡ÙÎ÷Sñ¿g{®Gáz”?gú~j¼t“@gF¶óýÔx	@g333333ï¿g)\�Âõ(ô?gö(\�Âõ	@g…ëQ¸é?gÅ °rh‘õ?g®Gáz®@gªñÒMb¨¿g
×£p=
·¿g¨ÆK7‰Aà?gZd;ßO�ë?g?5^ºIÀgžï§ÆK7ñ¿g!°rh‘í¼?g{®GázÀg¶óýÔxéÖ?gyé&1¬ì¿g–C‹lçûÀs    Cs    Ns    Os    Hi   i   i   t   OK(   gHáz®§9@g5^ºI¢E@gq=
×£PE@(   g}?5^º	8@gw¾Ÿ/íE@gB`åÐ"«E@(   g?5^ºI7@g33333CF@g…ëQ¸¥E@(   g+‡95@gÑ"Ûù~F@gÍÌÌÌÌ¼D@(   gú~j¼tó8@gö(\�Â•G@g{®GáÊB@(   g?5^ºIL3@gj¼t“xG@g�•C‹lçI@(   gÙÎ÷SãÀg¸…ëQ¸Ž?g}?5^ºIú?(   g¸…ëQ¸ø¿g¸…ëQ¸Ž?g}?5^ºIú?(   gÛù~j¼tç¿g¸…ëQ¸Ž?gœÄ °rh@(   g‹lçû©ñâ?g¸…ëQ¸Ž?g)\�Âõ(@(   g-²�ï§Æã?g¸…ëQ¸Ž?g%�•C‹ð?(   gTã¥›Ä ü?g¸…ëQ¸Ž?gé&1¬º?(   g#Ûù~j¼@g¸…ëQ¸®¿gh‘í|?5æ?(   gffffffø?g¸…ëQ¸Ž?gÃõ(\�Âõ¿(   g¶óýÔxé@gú~j¼t“˜¿gé&1¬À(   gÓMbX9Ä?g¸…ëQ¸Ž?gh‘í|?5þ¿(   gš™™™™™™¿gú~j¼t“˜?g¤p=
×£
À(   gZd;ßO�ï¿g¸…ëQ¸Ž?gã¥›Ä °î¿(   g=
×£p=Àgü©ñÒMb€?gºI+‡ö¿(   gÛù~j¼tç¿g¸…ëQ¸Ž?g®Gázâ?(   gøSã¥›Ä
Àgü©ñÒMb�?g²�ï§ÆK@(   g-²�ï§Æ
Àgyé&1¬ì?gTã¥›Ä ò?(   g-²�ï§Æ
Àg-²�ï§Æë¿g“V-ò?(   g‡ÙÎ÷Sñ¿g{®Gáz”?gú~j¼t“@(   gF¶óýÔx	@g333333ï¿g)\�Âõ(ô?(   gö(\�Âõ	@g…ëQ¸é?gÅ °rh‘õ?(   g®Gáz®@gªñÒMb¨¿g
×£p=
·¿(   g¨ÆK7‰Aà?gZd;ßO�ë?g?5^ºIÀ(   gžï§ÆK7ñ¿g!°rh‘í¼?g{®GázÀ(   g¶óýÔxéÖ?gyé&1¬ì¿g–C‹lçûÀ(   t   lent   AssertionErrort)   cctbx.crystal.distance_based_connectivityR   t   scitbx.array_familyR   t   vec3_doublet
   std_stringt   sortedt   sizet   list(   t   argsR   R   R	   R
   t	   bond_listt   bt   bonds(    (    s[   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/crystal/tst_distance_based_connectivity.pyt   run   s^    !O$-(t   __main__NR   i   (   t
   __future__R    R   t   __name__t   syst   argv(    (    (    s[   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/crystal/tst_distance_based_connectivity.pyt   <module>   s   	6