ó
]xTc           @   sX   d  d l  m Z d  d l m Z d  d l m Z d   Z e d k rT e   d GHn  d S(   i’’’’(   t   division(   t   dmatrix(   t   mathc          C   sU  d d d d d g }  d } t  d d  } xJ t d	 d
  D]9 } t |  | d | j d d	 |   | k  s: t  q: Wt  d d  } xJ t d	 d
  D]9 } t |  | d | j d d	 |   | k  s t  q Wt  d d  } xJ t d	 d
  D]9 } t |  | d | j d d	 |   | k  sņ t  qņ Wxm t d  D]_ } d t j |  t j |  d } t t  d |  j d d d  |  | k  s<t  q<Wt t  d t j  j d d d   | k  sŅt  t t  d t j d  j d d d   | k  s	t  t t  d t j d  j d d d  t j	 d  d  | k  sQt  d  S(   Ng1µ„ņśā?gw;Sč¼ä?g t//ĄŪ?gXä×±ĮČ?g pvk«?g-Cėā6?i   g      š?iž’’’i   i   i
   g       @i   i   i    g      ą?g      @g      @(
   R   t   ranget   abst   djmnt   AssertionErrort   smatht   cost   sint   pit   sqrt(   t   expectt   epst   d2t   nt   d4t   d7t   iit   expt(    (    sK   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_dmatrix.pyt   tst_dmatrix   s"    777%837t   __main__t   OKN(	   t
   __future__R    t   scitbx.mathR   t   stdlibR   R   R   t   __name__(    (    (    sK   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/math/tests/tst_dmatrix.pyt   <module>   s   	