ó
�™xTc           @   s\   d  d l  m Z d  d l m Z d  d l m Z d „  Z d „  Z e d k rX e ƒ  n  d S(   iÿÿÿÿ(   t   division(   t   crystal(   t   millerc          C   sR  d }  d } d d l  m } d d l m } t j d, d	 ƒ } | j d- d. d/ d0 d1 d2 d3 d4 d5 f	 ƒ } | j d6 ƒ } | j d7 ƒ } t j	 | | d# t
 ƒ} t j | d$ | d% | ƒj | j ƒ  ƒ }	 |	 j d& d' d( t
 ƒ |	 j d) d* ƒ }	 |	 j d& d+ d( t
 ƒ t d' ƒ j ƒ  }
 t d+ ƒ j ƒ  } |
 |  k s<t ‚ | | k sNt ‚ d  S(8   Nsˆ      4 I 4
  1  0  0  0  1  0  0  0  1
  0  0  0
  0 -1  0  1  0  0  0  0  1
  0  0  0
  0  1  0 -1  0  0  0  0  1
  0  0  0
 -1  0  0  0 -1  0  0  0  1
  0  0  0
   0   0   4   0   0   4     0 1 0  1 80331.8  8648.0
   0   0  12   0   0  12     0 1 0  1104526.7 11623.3
   0   0  -4   0   0   4     0 2 0  1 44883.6  4527.6
   0   0  -8   0   0   8     0 2 0  1 41134.1  4431.9
   0   0   8   0   0   8     0 1 0  1 50401.8  5464.6
   0   0   8   0   0   8     0 1 0  1 53386.0  5564.1
   0   0  -8   0   0   8     0 2 0  1119801.4 12231.2
   0   0  12   0   0  12     0 1 0  1105312.9 12602.2
   0   0  16   0   0  16     0 1 0  1 14877.6  2161.5
sˆ      4 I 4
  1  0  0  0  1  0  0  0  1
  0  0  0
  0 -1  0  1  0  0  0  0  1
  0  0  0
  0  1  0 -1  0  0  0  0  1
  0  0  0
 -1  0  0  0 -1  0  0  0  1
  0  0  0
   0   0   4   0   0   4     0 1 0  1 80331.8  8648.0
   0   0  -4   0   0   4     0 2 0  1 44883.6  4527.6
   0   0   8   0   0   8     0 1 0  1 50401.8  5464.6
   0   0   8   0   0   8     0 1 0  1 53386.0  5564.1
   0   0  -8   0   0   8     0 2 0  1 41134.1  4431.9
   0   0  -8   0   0   8     0 2 0  1119801.4 12231.2
   0   0  12   0   0  12     0 1 0  1104526.7 11623.3
   0   0  12   0   0  12     0 1 0  1105312.9 12602.2
   0   0  16   0   0  16     0 1 0  1 14877.6  2161.5
iÿÿÿÿ(   t   observation_types(   t   flexgÛù~j¼|\@gé&1¬<@@g     €V@t   I4i    i   i   iüÿÿÿiøÿÿÿi   i   gÍÌÌÌ¼œó@g3333ë„ù@g3333sêå@g3333Ãä@gš™™™9œè@g    @ê@gffff–?ý@gffff¶ù@gÍÌÌÌÌÍ@g     äÀ@gffff¦³Æ@gš™™™™¯±@gffffæO±@gš™™™™Xµ@gš™™™¼µ@gš™™™™ãÇ@gš™™™�È@g     ã @t   anomalous_flagt   datat   sigmast	   file_names   unsorted.scat%   scale_intensities_for_scalepack_merget   by_valuet   asu_indicess
   sorted.sca(   gÛù~j¼|\@gÛù~j¼|\@gé&1¬<@@g     €V@g     €V@g     €V@(   i    i    i   (   i    i    i   (   i    i    iüÿÿÿ(   i    i    iøÿÿÿ(   i    i    i   (   i    i    i   (   i    i    iøÿÿÿ(   i    i    i   (   i    i    i   (	   gÍÌÌÌ¼œó@g3333ë„ù@g3333sêå@g3333Ãä@gš™™™9œè@g    @ê@gffff–?ý@gffff¶ù@gÍÌÌÌÌÍ@(	   g     äÀ@gffff¦³Æ@gš™™™™¯±@gffffæO±@gš™™™™Xµ@gš™™™¼µ@gš™™™™ãÇ@gš™™™�È@g     ã @(   t
   cctbx.xrayR   t   cctbx.array_familyR   R   t   symmetryt   miller_indext   doubleR   t   sett   Truet   arrayt   set_observation_typet	   intensityt   export_as_scalepack_unmergedt   sortt   opent   readt   AssertionError(   t   expected_unsorted_datat   expected_sorted_dataR   R   t   xst   miR   R   t   mst   i_obst   unsortedt   sorted(    (    sR   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/regression/tst_miller_sort_asu.pyt   exercise_sort   s:    c           C   s   t  ƒ  d GHd  S(   Nt   OK(   R$   (    (    (    sR   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/regression/tst_miller_sort_asu.pyt   exerciseL   s    t   __main__N(   t
   __future__R    t   cctbxR   R   R$   R&   t   __name__(    (    (    sR   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/regression/tst_miller_sort_asu.pyt   <module>   s   	F	