Ńņ
AüKc           @   s2   d  d k  l Z d  d k Z d   Z d   Z d S(   i’’’’(   t   divisionNc      
   C   s°  | |  j  o | |  }  } n t  |  } | |  |  } | | | } | | |  } d t i d | | |  } t i |  }	 t i |  }
 |	 d | |
 d | d |
 |	 | | |  |  <|
 d | |	 d | d |
 |	 | | | | <d | | |  <xd t |   D]V } | |  | } | | | } |	 | |
 | | |  | <|	 | |
 | | | | <qWxk t |  d |  D]V } | | |  } | | | } |	 | |
 | | | |  <|	 | |
 | | | | <q}Wxk t | d |  D]V } | | |  } | | | } |	 | |
 | | | |  <|	 | |
 | | | | <qėWxd t |  D]V } | |  | } | | | } |	 | |
 | | |  | <|	 | |
 | | | | <qRWd  S(   Ng      ą?i   i    i   (   t   lent   matht   atan2t   cost   sint   range(   t   pt   qt   At   Vt   nt   appt   aqqt   apqt   phit   ct   st   it   aipt   aiqt   apit   aqit   vipt   viq(    (    sG   /usr/users/fedorov/public_html/Numeric/Programs/Python/jacobi/jacobi.pyt   rotate   sH    22    c      
   C   s  t  |   } t  |  d  } | | j o d Sg  } t |  D]B } | g  } t |  D]! } | | | j o d n d qY ~ q? ~ } |  } d }	 d }
 d } x³ |
 d j o„ d }
 x t |  D] } x t | d |  D]l } t | | |  |	 t | | |  t | | |  j o' d }
 t | | | |  | d 7} qß qß WqÅ Wq„ Wg  } t |  D] } | | | | qi~ } | | | f S(   Ni    s   not square matrixi   gź-q=(   R   R   t   absR   (   t   MR   t   mt   _[1]t   jt   _[2]R   R
   R	   t   epst   rotatedt   sweepst   rR   t   _[3]t   E(    (    sG   /usr/users/fedorov/public_html/Numeric/Programs/Python/jacobi/jacobi.pyt   jacobi+   s,    \   ?/(   t
   __future__R    R   R   R&   (    (    (    sG   /usr/users/fedorov/public_html/Numeric/Programs/Python/jacobi/jacobi.pyt   <module>   s   	&