ó
T™xTc           @   sW   d  d l  m Z d „  Z d „  Z d „  Z e d k rS e e ƒ e e ƒ d GHn  d S(   iÿÿÿÿ(   t   divisionc         C   s   t  t t t t |  ƒ ƒ ƒ S(   N(   t   tuplet   mapt   intt   round(   t   a(    (    sN   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/miller/tst_map_to_asu_isym.pyt   intify   s    c         C   sÍ  d d l  m } | j |  j ƒ  ƒ } g  } g  } x�| D]‡} t } x© t |  j ƒ  ƒ D]• } xŒ t |  j ƒ  ƒ D]x }	 |  d | |	 ƒ }
 t | |
 j	 ƒ  ƒ } | j
 | ƒ rp | j | ƒ | rÑ | j |	 ƒ n | j |	 ƒ t } Pqp qp WqW W| rü q8 n |  j ƒ  st ‚ x­ t |  j ƒ  ƒ D]™ } x� t |  j ƒ  ƒ D]| }	 |  d | |	 ƒ }
 g  | D] } | ^ qZ} t | |
 j	 ƒ  ƒ } | j
 | ƒ r;| j | ƒ | j |	 ƒ t } Pq;q;Wq"Wq8 W| | f S(   Niÿÿÿÿ(   t   sgtbxi    (   t   cctbxR   t   reciprocal_space_asut   typet   Falset   ranget   f_invt   n_smxR   t   rt	   is_insidet   appendt   Truet
   is_centrict   AssertionError(   t   sgt   miR   t   asut   isym_t   mi_t   hklt   foundt   i_invt   i_smxt   rt_mxt   hkl_t   ht   _hklt   mhkl_(    (    sN   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/miller/tst_map_to_asu_isym.pyt   reference_map   s>    c         C   s=  d d l  m } d d l m } d d l m } | j ƒ  } | j ƒ  } d d  l } d } xY t	 | ƒ D]K } g  t	 d ƒ D] } | j
 d d ƒ ^ qz }	 | j |	 ƒ | j d	 ƒ qg W| j d
 ƒ j ƒ  }
 | j |
 ƒ } t | | ƒ \ } } | | j ƒ  |  | | ƒ x. t	 | ƒ D]  } | | | | k st ‚ qWd  S(   Niÿÿÿÿ(   R   (   t   map_to_asu_isym(   t   flexiè  i   iöÿÿÿi
   i    iÃ   (   R   R   t   cctbx.millerR$   t   cctbx.array_familyR%   t   miller_indexR   t   randomR   t   randintR   t   space_group_symbolst   hallt   space_groupR#   R
   R   (   t   anomalous_flagR   R$   R%   R   t   iR)   t   nhklt   jR   t
   spacegroupR   R   R   (    (    sN   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/miller/tst_map_to_asu_isym.pyt   tst_map_to_asu_isym-   s"    +t   __main__t   OKN(   t
   __future__R    R   R#   R3   t   __name__R   R   (    (    (    sN   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/miller/tst_map_to_asu_isym.pyt   <module>   s   		'	

