ó
\xTc           @   s[   d  d l  m Z d  d l m Z d  d l m Z d d  Z e d k rW e   d GHn  d S(	   i’’’’(   t   division(   t   math(   t   flexi   c         C   s/  t  j |   d d } t j |  } xW t d  D]I } d d d t  j |   t  j t |    d d } | j |  q3 W| j   } | j   } | j	   } x* | D]" } t
 | d  d k  s« t  q« Wx* | D]" }	 t
 |	 d  d	 k  sŲ t  qŲ Wx* | D]" }
 t
 |
 d
  d k  st  qWd  S(   Ng        g      š?i  g¹?g       @gü©ńŅMbP?gÉ?i   g-Cėā6?gF%uk?g|®Gįz¤?gPč“Nk?(   R   t   doubleR   t   multivariate_momentst   ranget   random_doublet   updatet   vcv_upper_triangle_packedt   meant   variancet   abst   AssertionError(   t   Nt   weightst   mvot   iit   tmpt   vcvR	   t   vart   mt   vt   c(    (    sR   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tst_multivariate_moments.pyt   tst   s    6  t   __main__t   OKN(   t
   __future__R    t   scitbxR   t   scitbx.array_familyR   R   t   __name__(    (    (    sR   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tst_multivariate_moments.pyt   <module>   s   