ó
µ™xTc           @   sE   d  d l  m Z d „  Z d „  Z d „  Z e d k rA e ƒ  n  d S(   iÿÿÿÿ(   t   divisionc         C   sj  d d l  m } t | d ƒ t | ƒ k rD t | ƒ t | ƒ k sJ t ‚ |  j d ƒ } | j d d d d g ƒ | j | d d | j d	 d
 ƒ} | j | d d d g d d ƒ|  j	 | d | ƒ| j
 t t | ƒ ƒ ƒ | j t t | ƒ ƒ ƒ | j g  | D] } d | ^ qƒ | j g  | D] } d | ^ q(ƒ | j d ƒ | j d ƒ | j d ƒ | S(   Niÿÿÿÿ(   t   cmi    io   gš™™™™™¹?g333333ë?i   t   cmapt   interpolationt   bilineart   colorsg      à?t	   linewidthi   t   axs   %gt   gammat   weights   DEN optimization (R-free)(   g      à?g      à?g      à?(   t
   matplotlibR   t   lent   AssertionErrort   add_subplott   set_positiont   contourft   YlGnBut   contourt   colorbart
   set_xtickst   ranget
   set_ytickst   set_xticklabelst   set_yticklabelst
   set_xlabelt
   set_ylabelt	   set_title(   t   figureR   R	   t   zR   t   pt   csett   x(    (    sD   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/den/plot_weights.pyt   _plot_weights   s    :!$$c         C   s=   d d l  m } | j ƒ  } t | |  | | ƒ | j ƒ  d  S(   Niÿÿÿÿ(   t   pyplot(   R
   R!   R   R    t   show(   R   R	   R   t   pltR   (    (    sD   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/den/plot_weights.pyt   plot_weights_pyplot   s    c    	   8   C   s   dA dB dC dD dE dF dG dH dI dJ dK dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ d[ d\ d] d^ d_ d` da db dc dd de df dg dh di dj dk dl dm dn do dp dq dr ds dt du dv g6 }  t  t t g  |  D] } | d> ^ q¸ ƒ ƒ ƒ } t  t t g  |  D] } | d? ^ qç ƒ ƒ ƒ } g  | D] } g  ^ q} t | ƒ } x^ t t | ƒ ƒ D]J } xA t t | ƒ ƒ D]- } | | | } | | j |  | d@ ƒ qWWq>Wt | | | ƒ d  S(w   Ng        g      à?gmçû©ñÒÍ?g      ð?g^KÈ=›Í?g      @g˜LŒJêÌ?g      @gÇK7‰A`Í?g      $@gÞ“‡…ZÓÌ?g      9@gŽuqàÍ?g      I@gäƒžÍªÏÍ?g      Y@gÖVì/»'Ï?g      i@gÖÅm4€·Ð?gš™™™™™É?gZõ¹ÚŠýÍ?góÒo_Î?gõJY†8ÖÍ?g\ AñcÌÍ?gV-²Í?gÉv¾Ÿ/Í?gˆôÛ×�sÎ?g'1¬ZÐ?gš™™™™™Ù?g%ušÎ?g®Ø_vOÎ?g+ö—Ý“‡Í?g†ZÓ¼ãÍ?g´Yõ¹ÚŠÍ?gq=
×£pÍ?gà-� ø1Î?gž^)ËÏ?g333333ã?gÏfÕçj+Î?gGx$(Î?g{ƒ/L¦
Î?gDioð…ÉÌ?gâX·Ñ Î?g}®¶bÙÍ?gs×òAÏÎ?gš™™™™™é?gM„O¯”Í?gyX¨5Í;Î?g“©‚QIÍ?gÙÎ÷SãÍ?g' ‰°áéÍ?gb¡Ö4ï8Í?g§èH.ÿ!Í?gÅ °rh‘Í?g#J{ƒ/LÎ?g¶„|Ð³YÍ?g”ö_˜LÍ?gVŸ«­Ø_Î?g)ËÇº¸Í?gKY†8ÖÅÍ?gì/»'Í?i    i   i   (   g        g      à?gmçû©ñÒÍ?(   g        g      ð?g^KÈ=›Í?(   g        g      @g˜LŒJêÌ?(   g        g      @gÇK7‰A`Í?(   g        g      $@gÞ“‡…ZÓÌ?(   g        g      9@gŽuqàÍ?(   g        g      I@gäƒžÍªÏÍ?(   g        g      Y@gÖVì/»'Ï?(   g        g      i@gÖÅm4€·Ð?(   gš™™™™™É?g      à?gZõ¹ÚŠýÍ?(   gš™™™™™É?g      ð?góÒo_Î?(   gš™™™™™É?g      @gõJY†8ÖÍ?(   gš™™™™™É?g      @g\ AñcÌÍ?(   gš™™™™™É?g      $@gV-²Í?(   gš™™™™™É?g      9@gÉv¾Ÿ/Í?(   gš™™™™™É?g      I@g^KÈ=›Í?(   gš™™™™™É?g      Y@gˆôÛ×�sÎ?(   gš™™™™™É?g      i@g'1¬ZÐ?(   gš™™™™™Ù?g      à?g%ušÎ?(   gš™™™™™Ù?g      ð?g®Ø_vOÎ?(   gš™™™™™Ù?g      @g\ AñcÌÍ?(   gš™™™™™Ù?g      @g+ö—Ý“‡Í?(   gš™™™™™Ù?g      $@g†ZÓ¼ãÍ?(   gš™™™™™Ù?g      9@g´Yõ¹ÚŠÍ?(   gš™™™™™Ù?g      I@gq=
×£pÍ?(   gš™™™™™Ù?g      Y@gà-� ø1Î?(   gš™™™™™Ù?g      i@gž^)ËÏ?(   g333333ã?g      à?gÏfÕçj+Î?(   g333333ã?g      ð?gq=
×£pÍ?(   g333333ã?g      @g%ušÎ?(   g333333ã?g      @gGx$(Î?(   g333333ã?g      $@g{ƒ/L¦
Î?(   g333333ã?g      9@gDioð…ÉÌ?(   g333333ã?g      I@gâX·Ñ Î?(   g333333ã?g      Y@g}®¶bÙÍ?(   g333333ã?g      i@gs×òAÏÎ?(   gš™™™™™é?g      à?gM„O¯”Í?(   gš™™™™™é?g      ð?gyX¨5Í;Î?(   gš™™™™™é?g      @g“©‚QIÍ?(   gš™™™™™é?g      @gÙÎ÷SãÍ?(   gš™™™™™é?g      $@g' ‰°áéÍ?(   gš™™™™™é?g      9@gÇK7‰A`Í?(   gš™™™™™é?g      I@gb¡Ö4ï8Í?(   gš™™™™™é?g      Y@g§èH.ÿ!Í?(   gš™™™™™é?g      i@gÅ °rh‘Í?(   g      ð?g      à?g#J{ƒ/LÎ?(   g      ð?g      ð?g¶„|Ð³YÍ?(   g      ð?g      @g”ö_˜LÍ?(   g      ð?g      @gVŸ«­Ø_Î?(   g      ð?g      $@gZõ¹ÚŠýÍ?(   g      ð?g      9@g)ËÇº¸Í?(   g      ð?g      I@gKY†8ÖÅÍ?(   g      ð?g      Y@gì/»'Í?(   g      ð?g      i@g\ AñcÌÍ?(   t   sortedt   listt   setR   R   t   appendR$   (	   t   grid_resultsR   R   R	   t   valuest   nt   it   jt   k(    (    sD   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/den/plot_weights.pyt   exercise   s~    	//!t   __main__N(   t
   __future__R    R    R$   R/   t   __name__(    (    (    sD   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/den/plot_weights.pyt   <module>   s
   			C