³ò
O?ÇGc           @   sÄ   d  d k  Td  d k Td  d k Td „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z	 d	 „  Z
 d
 „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d S(   iÿÿÿÿ(   t   *c         C   sa   t  |  ƒ } t | | ƒ } x? t d | ƒ D]. } |  | } | t | ƒ } | | | | <q+ W| S(   Ni    (   t   lent   make_zero_matrixt   ranget   conj(   t   psit   nt   rhot   jt   cjt   cj2(    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt
   densop_std   s     
c   
   	   C   sÂ   t  |  ƒ } t | | ƒ } t |  | ƒ } x‘ t d | ƒ D]€ } | | } | t | ƒ } x] t d | ƒ D]L } xC t d | ƒ D]2 }	 | | |	 c | | | | | | |	 7<q€ Wqj Wq: W| S(   Ni    (   R   R   t   basis_coeffs_onR   R   (
   R   t   basisR   R   t   coeffst   kt   ckt   ck2t   iR   (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt	   densop_on#   s     
  8c      	   C   sÑ   t  |  ƒ } t | | ƒ } t |  | ƒ } t | | ƒ } x‘ t d | ƒ D]€ } | | } | t | ƒ } x] t d | ƒ D]L }	 xC t d | ƒ D]2 }
 | |	 |
 c | | | |	 | | |
 7<q Wqy WqI W| S(   Ni    (   R   R   t   basis_coeffst   linear_combinationR   R   (   R   R   R   R   R   t   wR   R   R   R   R   (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   densop0   s     
  8c      	   C   s¡   t  | d ƒ } t | | ƒ } x{ t d | ƒ D]j } xa t d | ƒ D]P } xG t d | ƒ D]6 } | | | c |  | | | | | | | 7<q[ WqE Wq/ W| S(   Ni    (   R   R   R   (   t   pst   psisR   R   R   R   R   (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   densopp>   s       <c           C   s
   t  d ƒ S(   Ni   (   t   make_identity_matrix(    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_IJ   s    c          C   s)   t  d ƒ }  t |  |  g |  |  g g ƒ S(   Ng      à?(   t   sqrtt   sackmat(   t   s(    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_HO   s    	c           C   s   t  d d g d d g g ƒ S(   Ni    i   (   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_nW   s    	c           C   s   t  d d g d d g g ƒ S(   Ni   i    (   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_nbar]   s    	c           C   s   t  d d g d d g g ƒ S(   Ni    i   (   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Xd   s    	c           C   s   t  d d g d d g g ƒ S(   Ni    i   iÿÿÿÿ(   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Yi   s    	c           C   s   t  d d g d d g g ƒ S(   Ni   i    iÿÿÿÿ(   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Zn   s    	c           C   s   t  d d g d d g g ƒ S(   Ni    i   (   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_sXv   s    	c           C   s   t  d d g d d g g ƒ S(   Ni    y              ð¿y              ð?(   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_sY{   s    	c           C   s   t  d d g d d g g ƒ S(   Ni   i    iÿÿÿÿ(   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_sZ€   s    	c           C   sF   t  d d d d g d d d d g d d d d g d d d d g g ƒ S(   Ni   i    (   R   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Sˆ   s
    c         C   s   t  t ƒ  |  | ƒ S(   N(   t
   multikroniR%   (   R   t   N(    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Zi   s    c         C   sµ   t  ƒ  } t ƒ  } t | |  | ƒ } t | | | ƒ } t | |  | ƒ } t | | | ƒ } t ƒ  }	 t |	 |  | ƒ }
 t |	 | | ƒ } | | | | |
 | | | | | } | S(   N(   R!   R"   R*   R#   (   R   R   R+   R   t   nbart   nit   njt   nbarit   nbarjt   Xt   Xit   Xjt   Sij(    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Sij•   s    			*c   
      C   sc   t  ƒ  } t ƒ  } t | |  | ƒ } t | |  | ƒ } t ƒ  } t | | | ƒ } | | | }	 |	 S(   N(   R!   R"   R*   R#   (
   R   R   R+   R   R-   R.   R0   R2   R4   t   Cij(    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pyt   get_Cij´   s    			N(   t
   sackmatc_mt
   cplxreal_mt   mathR   R   R   R   R   R    R!   R"   R#   R$   R%   R&   R'   R(   R)   R,   R6   R8   (    (    (    s/   /u5/kerl/pub_http_internet/doc/qnt/code/qc_m.pys   <module>	   s(   


																		