ó
UxTc           @   sT   d  d l  m Z d  d l m Z d  d l Z d   Z e d k rP e   d GHn  d S(   i˙˙˙˙(   t   division(   t   wignerNc          C   s?  t  j d  }  t j d  } t  j d d d d d d  } |  j d d d d d d  } t | |  d k  sv t  t | j   |  d k  s t  t j d  } t  j d d d	 d d d
  } |  j d d d	 d d d
  } t | |  d k  s˙ t  t | j   |  d k  s!t  t j d  d } t  j d d d d	 d
 d  } |  j d d d d	 d
 d  } t | |  d k  st  t | j   |  d k  sŽt  t  j d d d d d d  } d t j d  } |  j d d d d d d  } t | |  d k  st  t | j   |  d k  s;t  d  S(   Ni   g      đ?g      @i   i   iţ˙˙˙gíľ ÷Ć°>g      ,@i   iý˙˙˙g      .@g    `dé@g      Ŕi   i   i
   i    i   iú˙˙˙g      @@g      @g      :@g   BŇEAgÉ?g$I$I˛?g¸ŹfMç2?gUUUUUU%@gUUUUUU%ŔgůÍó\ă>(	   R   t   wigner3j_fastt   matht   sqrtt   wigner3jt   computet   abst   AssertionErrort	   get_value(   t   wigner_fastt   resultt   wt   value(    (    sH   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/wigner/tst_wigner3j.pyt   tst1   s*    """t   __main__t   OK(   t
   __future__R    t   scitbxR   R   R   t   __name__(    (    (    sH   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/wigner/tst_wigner3j.pyt   <module>   s   	