ó
_™xTc           @   sg   d  d l  m Z d  d l Z d  d l m Z d  d l Z d d „ Z d „  Z e	 d k rc e ƒ  n  d S(   iÿÿÿÿ(   t   divisionN(   t   flexi    c         C   sõ  t  j t d ƒ ƒ } t j j | |  d ƒ } | j d ƒ } t j j | |  d ƒ } t j j | |  d ƒ } t  j | ƒ } t  j | | ƒ | | } d } d }	 d }
 d } d } d } x’ t t	 | ƒ ƒ D]~ } | j t	 | ƒ ƒ } | j t	 | ƒ ƒ } t  j | ƒ } | | 7} |	 | | 7}	 t  j | ƒ } |
 | 7}
 | | | 7} qÓ W| | } |	 | | | }	 t j
 |	 ƒ }	 |
 | }
 | | |
 |
 } t j
 | ƒ } t j d	 | ƒ |	 d
 k  sÎt ‚ t j d	 |
 ƒ | d
 k  sñt ‚ d  S(   Ni   i   id   i   i   i    g     @�@g      Y@g      #@i   (   R   t   doublet   ranget   scitbxt   matht   non_parametric_bootstrapt   drawt   smooth_bootstrapt   meant   intt   sqrtt   fabst   AssertionError(   t   seedt   obs_orit   npb_drawt   obst   npbt   sbst   mean_tt   var_tt   mean_of_meant   var_of_meant   mean_sbst   std_sbst   n_samplet   sizet	   iterationt   samplet   sample_2t   single_meant   tmp(    (    sL   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_resample.pyt   test_resample   s<    



#c          C   s*   x t  d ƒ D] }  t |  ƒ q Wd GHd  S(   Ni
   t   OK(   R   R!   (   t   ii(    (    sL   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_resample.pyt   run2   s    t   __main__(
   t
   __future__R    t   scitbx.mathR   t   scitbx.array_familyR   R   R!   R$   t   __name__(    (    (    sL   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_resample.pyt   <module>   s   ,	