ó
²™xTc           @   sÀ   d  d l  m Z d  d l Z d  d l m Z d  d l Z d  d l Z d „  Z e	 d k r¼ x` d d d d d g D]F Z
 e d e
 d d e
 d d e
 d d e
 d d e
 d ƒ Z e GHqo Wn  d S(   iÿÿÿÿ(   t   divisionN(   t   flexc         C   s·   t  j ƒ  } t  j ƒ  } xW t d |  ƒ D]F } | |  d } | t j | | ƒ } | j | ƒ | j | ƒ q( Wt j ƒ  }	 t j j	 | | | | ƒ }
 t
 |
 ƒ d t j ƒ  |	 f S(   Ni   g      ð?s   Time: %6.4f(   R   t   doublet   ranget   matht   expt   appendt   timet   mmtbxt   bulk_solventt   fit_k_exp_b_to_k_totalt   list(   t   nreft   kt   bt   k_startt   b_startt   datat   sst   ss_t   yt   t0t   r(    (    sX   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/bulk_solvent/tst_fit_kexpb_to_ktotal.pyt   run   s    t   __main__i'  i   i
   id   i   i    iZ   iöÿÿÿR   R   R   i   R   i   R   i   (   i'  i   i
   i   i
   (   i'  i
   id   i
   id   (   i'  i   i
   i   i
   (   i'  i    i    i    i    (   i'  iZ   iöÿÿÿiöÿÿÿi
   (   t
   __future__R    t   mmtbx.bulk_solventR   t   scitbx.array_familyR   R   R   R   t   __name__t   iR   (    (    (    sX   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/bulk_solvent/tst_fit_kexpb_to_ktotal.pyt   <module>   s   	;