ó
m™xTc           @   sx   d  d l  m Z d  d l m Z d  d l Z d  d l m Z d  d l m Z d „  Z	 d „  Z
 e d k rt e
 ƒ  n  d S(	   iÿÿÿÿ(   t   division(   t   linalgN(   t   flex(   t   approx_equalc         C   s�   t  j d d ƒ j } xt t |  ƒ D]f } d } | t j | ƒ ƒ j ƒ  } t  j j | ƒ } t	 | j
 ƒ  t | d t ƒd d ƒs" t ‚ q" Wd  S(	   Ni   i   i   i   t   reverset   epsgê-�™—q=(   i   i   i   i   (   R   t   random_normal_matrix_generatort!   symmetric_matrix_with_eigenvaluest   xrangeR   t   doublet   matrix_packed_u_as_symmetrict   eigensystemt   real_symmetricR   t   valuest   sortedt   Truet   AssertionError(   t   n_trialst   symm_matt   it   lambdast   at   eigen(    (    sQ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/matrix/tests/tst_householder.pyt    exercise_random_symmetric_matrix   s    c           C   s   t  d ƒ d GHd  S(   Ni2   t   OK(   R   (    (    (    sQ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/matrix/tests/tst_householder.pyt   run   s    
t   __main__(   t
   __future__R    t   scitbxR   t   scitbx.linalg.eigensystemt   scitbx.array_familyR   t   libtbx.test_utilsR   R   R   t   __name__(    (    (    sQ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/matrix/tests/tst_householder.pyt   <module>   s   	
	