3
\T� � $ @ sF d ddddddddd d
d d d
ddddddddddddddd d d d d d!d"d#g$Z eZd$Zd%Zd&Zd'd(lZd'd(lZd'd(l Z yd'd)l
m Z e dd*�Z
W n ek
r� d+d,� Z
Y nX dZdZdZdZdZdZdZdZd-Ze jd�kr�d1Zd1Zd�Zn d2Zd2Zd�Zeed0 Z G d3d� de �Z G d4d� de �Z G d5d� de �Z G d6d� de �Z!G d7d � d e e"�Z#G d8d� de �Z$G d9d� de e"�Z%G d:d
� d
e �Z&G d;d� de �Z'G dd
� d
e&e(�Z*G d?d� de&e(e)�Z+G d@d� de e,�Z-e e#e&e*e(e+e e)e-g Z.e!e e$e e%e e'e iZ/eeeeeeeefZ0y d'd(l1Z1W n. ek
�rd G dAdB� dBe2�Z3e3� Z1[3Y nX y
e1j4 W n> e5k
�r� e6e1j7� dC��r�e1j7� `8dDd � Z9dEd � Z:Y n6X e1j4� Z4e6e4dC��r�e4`8e4fdFd �Z:e4fdGd �Z9[1[4d�dHd �Z;G dId � d e2�Zj?e<� G dMdN� dNe2�Z@G dOd� de2�ZAG dPdQ� dQe2�ZBd�dRdS�ZCeDjEZFdTdU� ZGdVdW� ZHdXdY� ZIdZd[� ZJd�d]d^�ZKd_d`� ZLdadb� ZMG dcdd� dde2�ZNeN� jOZPd�dedf�ZQdgdh� ZRdidj� ZSdkdldmdndodpdqdrdsdt� fdudv�ZTd�dwdx�ZUd�dydz�ZVeAd{ee#e*e gg d|d�d0d'd}�ZWeAd~ee#e*e e e+gg d�ZXeAd~eg g d�ZYd'd(lZZZeZj[d�eZj\eZj]B �j^Z_eZj[d��j^Z`eZj[d��j^ZaeZj[d�eZj\eZjbB �Zc[Zy d'd(ldZeW n ek
�r� Y nX d�d�d��Zfd�d�� Zgd�d�� Zhd�d�d��Zid�d�� Zjd�d�� Zke� s r+ T� �? � l ��N�Zoi@�Tc @ s e Zd Zdd� ZdS )r c G s d S )Nr( )�self�contextr) r( r( r* �handle� s zDecimalException.handleN)�__name__�
__module__� __qualname__r1 r( r( r( r* r � s c @ s e Zd ZdS )r N)r2 r3 r4 r( r( r( r* r � s
c @ s e Zd Zdd� ZdS )r c G s, |r(t |d j|d jdd�}|j|�S tS )Nr& �nT)�_dec_from_triple�_sign�_int�_fix_nan�_NaN)r/ r0 r) �ansr( r( r* r1 � s
zInvalidOperation.handleN)r2 r3 r4 r1 r( r( r( r* r � s c @ s e Zd Zdd� ZdS )r c G s t S )N)r: )r/ r0 r) r( r( r* r1 s zConversionSyntax.handleN)r2 r3 r4 r1 r( r( r( r* r � s c @ s e Zd Zdd� ZdS )r
c G s t | S )N)�_SignedInfinity)r/ r0 �signr) r( r( r* r1 s zDivisionByZero.handleN)r2 r3 r4 r1 r( r( r( r* r
s
c @ s e Zd Zdd� ZdS )r c G s t S )N)r: )r/ r0 r) r( r( r* r1 s zDivisionImpossible.handleN)r2 r3 r4 r1 r( r( r( r* r s c @ s e Zd Zdd� ZdS )r c G s t S )N)r: )r/ r0 r) r( r( r* r1 ) s zDivisionUndefined.handleN)r2 r3 r4 r1 r( r( r( r* r ! s c @ s e Zd ZdS )r N)r2 r3 r4 r( r( r( r* r , s
c @ s e Zd Zdd� ZdS )r c G s t S )N)r: )r/ r0 r) r( r( r* r1 C s zInvalidContext.handleN)r2 r3 r4 r1 r( r( r( r* r 8 s c @ s e Zd ZdS )r N)r2 r3 r4 r( r( r( r* r F s
c @ s e Zd ZdS )r
N)r2 r3 r4 r( r( r( r* r
R s c @ s e Zd Zdd� ZdS )r c G s� |j ttttfkrt| S |dkrR|j tkr4t| S t|d|j |j |j d �S |dkr�|j t
krlt| S t|d|j |j |j d �S d S )Nr& �9r. ) �roundingr r r r r< r r6 �prec�Emaxr )r/ r0 r= r) r( r( r* r1 s s
zOverflow.handleN)r2 r3 r4 r1 r( r( r( r* r ] s c @ s e Zd ZdS )r N)r2 r3 r4 r( r( r( r* r � s
c @ s e Zd ZdS )r N)r2 r3 r4 r( r( r( r* r � s
c @ s e Zd Zefdd�ZdS )�
MockThreadingc C s
|j t S )N)�modules� __xname__)r/ �sysr( r( r* �local� s zMockThreading.localN)r2 r3 r4 rE rF r( r( r( r* rB � s rB �__decimal_context__c C s, | t ttfkr | j� } | j� | tj� _d S )N)r r r �copy� clear_flags� threading�current_threadrG )r0 r( r( r* r � s c C s4 y
t j� jS tk
r. t� } | t j� _| S X d S )N)rJ rK rG �AttributeErrorr )r0 r( r( r* r � s
c C s, y| j S tk
r& t� }|| _ |S X d S )N)rG rL r )�_localr0 r( r( r* r � s c C s( | t ttfkr | j� } | j� | |_d S )N)r r r rH rI rG )r0 rM r( r( r* r � s c C s | d krt � } t| �S )N)r �_ContextManager)Zctxr( r( r* r � s $ c
@ s� e Zd Zd�Zd�dd�Zed d
� �Zd d � Zd
d� Zd�dd�Z dd� Z
dd� Z dd� Z d�dd�Z
d�dd�Zd�dd �Zd�d d �Zd�d d �Zd�d!d"�Zd#d$� Zd%d&� Zd'd(� Zd)d*� Zd�d,d-�Zd�d.d/�Zd�d0d1�Z�d d2d3�Z�dd5d6�Z�dd7d8�Z e Z �dd9d:�Z �dd;d<�Z �dd=d>�Z e Z!�dd?d@�Z"dAdB� Z#�ddCdD�Z$�ddEdF�Z%�d dGdH�Z&�d
dIdJ�Z'�d dKdL�Z(�d dMdN�Z)�d
dOdP�Z*�ddQdR�Z+dSdT� Z,dUdV� Z-e-Z.dWdX� Z/e0e/�Z/dYdZ� Z1e0e1�Z1d[d\� Z2d]d^� Z3d_d`� Z4dadb� Z5dcdd� Z6dedf� Z7dgdh� Z8didj� Z9dkdl� Z:dmdn� Z;dodp� Ze6e7e8e9e:e;e
z Decimal._cmpc C s< t | |dd�\} }|tkr |S | j||�r.dS | j|�dkS )NT)� equality_opFr& )�_convert_for_comparison�NotImplementedr� r� )r/ r� r0 r( r( r* �__eq__� s zDecimal.__eq__c C s< t | |�\} }|tkr|S | j||�}|r.dS | j|�dk S )NFr& )r� r� r� r� )r/ r� r0 r; r( r( r* �__lt__� s zDecimal.__lt__c C s< t | |�\} }|tkr|S | j||�}|r.dS | j|�dkS )NFr& )r� r� r� r� )r/ r� r0 r; r( r( r* �__le__� s zDecimal.__le__c C s< t | |�\} }|tkr|S | j||�}|r.dS | j|�dkS )NFr& )r� r� r� r� )r/ r� r0 r; r( r( r* �__gt__� s zDecimal.__gt__c C s< t | |�\} }|tkr|S | j||�}|r.dS | j|�dkS )NFr& )r� r� r� r� )r/ r� r0 r; r( r( r* �__ge__� s zDecimal.__ge__c C s> t |dd�}| js |r0|jr0| j||�}|r0|S t| j|��S )NT)�raiseit)�_convert_otherrP r� r r� )r/ r� r0 r; r( r( r* �compare� s zDecimal.comparec C s� | j r4| j� rtd��n | j� r$tS | jr0t S tS | jdkrNtd| jt �}ntt
| j t �}t | j �| t }| dkr||n| }|dkr�dS |S )Nz"Cannot hash a signaling NaN value.r& �
r. r, r� �����)
rP r� rr �is_nan� _PyHASH_NANr7 � _PyHASH_INFrO �pow�_PyHASH_MODULUS�
_PyHASH_10INVrU r8 )r/ Zexp_hashZhash_r; r( r( r* �__hash__� s
zDecimal.__hash__c C s t | jttt| j��| j�S )N)r r7 rk ro rU r8 rO )r/ r( r( r* �as_tuple� s zDecimal.as_tuplec C s� | j r | j� rtd��ntd��| s(dS t| j�}| jdkrR|d| j d }}nr| j }x(|dkr�|d dkr�|d }|d8 }q\W | j }t|| @ j� d |�}|r�||L }||8 }d| |> }| j r�| }||fS )Nz#cannot convert NaN to integer ratioz(cannot convert Infinity to integer ratior& r. r� rz )r& r. )
rP r� rl �
OverflowErrorrU r8 rO �minr~ r7 )r/ r5 r� Zd5Zd2Zshift2r( r( r* r} � s,
zDecimal.as_integer_ratioc C s dt | � S )Nz
Decimal('%s'))r` )r/ r( r( r* �__repr__' s zDecimal.__repr__Fc C sd ddg| j }| jrL| jdkr&|d S | jdkr>|d | j S |d | j S | jt| j� }| jdkrt|dkrt|}n6|s~d
}n,| jd kr�|d
d d
}n|d
d d
}|dkr�d }d
d | | j }nN|t| j�k� r�| jd |t| j� }d}n | jd |� }d
| j|d � }||k�r*d}n*|d k�r:t� }ddg|j d|| }|| | | S )NrS rT rZ ZInfinityr5 �NaNr� r& � r. rQ r[ �.�e�Ez%+di����)r7 rP rO r8 rf r �capitals) r/ �engr0 r= �
leftdigits�dotplacerv rw rV r( r( r* �__str__, s:
zDecimal.__str__c C s | j d|d�S )NT)r� r0 )r� )r/ r0 r( r( r* �
to_eng_string` s zDecimal.to_eng_stringc C sT | j r| j|d�}|r|S |d kr(t� }| rB|jtkrB| j� }n| j� }|j|�S )N)r0 )rP r� r r? r �copy_abs� copy_negate�_fix)r/ r0 r; r( r( r* �__neg__i s
zDecimal.__neg__c C sT | j r| j|d�}|r|S |d kr(t� }| rB|jtkrB| j� }nt| �}|j|�S )N)r0 )rP r� r r? r r� r r� )r/ r0 r; r( r( r* �__pos__ s
zDecimal.__pos__Tc C sJ |s | j � S | jr&| j|d�}|r&|S | jr:| j|d�}n | j|d�}|S )N)r0 )r� rP r� r7 r� r� )r/ �roundr0 r; r( r( r* �__abs__� s zDecimal.__abs__c
C sl t |�}|tkr|S |d kr"t� }| js.|jr�| j||�}|rB|S | j� rr| j|jkrj|j� rj|jtd�S t | �S |j� r�t |�S t
| j |j �}d}|j t
kr�| j|jkr�d}| r�| r�t
| j|j�}|r�d}t|d|�}|j|�}|S | �s"t||j |j d �}|j||j �}|j|�}|S |�sZt|| j |j d �}| j||j �}|j|�}|S t| �}t|�}t|||j�\}}t� } |j|jk�r�|j|jk�r�t|d|�}|j|�}|S |j|jk �r�|| }}|jdk�r�d| _|j|j |_|_nd| _n&|jdk�r d| _d\|_|_nd| _|jdk�r@|j|j | _n|j|j | _|j| _t | �}|j|�}|S )Nz
-INF + INFr& r. rQ )r& r& )r� r� r rP r� r� r7 rd r r r� rO r? r r6 r� �maxr@ �_rescaleri �
_normalizer= rU rV )
r/ r� r0 r; rV Z negativezeror= �op1�op2r� r( r( r* �__add__� s|
zDecimal.__add__c C sH t |�}|tkr|S | js |jr6| j||d�}|r6|S | j|j� |d�S )N)r0 )r� r� rP r� r� r� )r/ r� r0 r; r( r( r* �__sub__ s zDecimal.__sub__c C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r� )r/ r� r0 r( r( r* �__rsub__ s zDecimal.__rsub__c C sD t |�}|tkr|S |d kr"t� }| j|jA }| js:|jr�| j||�}|rN|S | j� rn|sf|jtd�S t | S |j� r�| s�|jtd�S t | S | j
|j
}| s�| r�t |d|�}|j |�}|S | j
dkr�t ||j
|�}|j |�}|S |j
dk�rt || j
|�}|j |�}|S t| �}t|�}t |t|j|j �|�}|j |�}|S )Nz (+-)INF * 0z 0 * (+-)INFrQ �1)r� r� r r7 rP r� r� rd r r< rO r6 r� r8 ri r` rU )r/ r� r0 Z
resultsignr; Z resultexpr� r� r( r( r* �__mul__ sH
zDecimal.__mul__c C s� t |�}|tkrtS |d kr"t� }| j|jA }| js:|jr�| j||�}|rN|S | j� rj|j� rj|jtd�S | j� rzt | S |j� r�|jt
d� t |d|j � �S |s�| s�|jt
d�S |jtd|�S | s�| j|j }d}n�t|j�t| j� |j d }| j|j | }t| �}t|�} |dk�r:t|jd| | j�\}}
n t|j| jd| �\}}
|
�rt|d dk�r�|d7 }n<| j|j } x.|| k �r�|d dk�r�|d }|d7 }�q�W t |t|�|�}|j|�S )
Nz(+-)INF/(+-)INFzDivision by infinityrQ z0 / 0zx / 0r& r. r� rz )r� r� r r7 rP r� r� rd r r< r r6 �Etinyr r
rO rf r8 r@ ri �divmodrU r` r� ) r/ r� r0 r= r; rV �coeff�shiftr� r� � remainder� ideal_expr( r( r* � __truediv__Q sP
zDecimal.__truediv__c C s | j |j A }|j� r | j}nt| j|j�}| j� |j� }| sP|j� sP|dkrjt|dd�| j||j�fS ||jk�r t | �}t |�}|j
|j
kr�| j d|j
|j
9 _ n| j d|j
|j
9 _ t |j |j �\}} |d|j k �r t|t
|�d�t| j t
| �|�fS |jtd�}
|
|
fS )Nr, rQ r& r� z%quotient too large in //, % or divmodr� )r7 r� rO r� r� r6 r� r? r@ ri rV rU r� r` rd r ) r/ r� r0 r= r� �expdiffr� r� �q�rr; r( r( r* �_divide� s*
zDecimal._dividec C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r� )r/ r� r0 r( r( r* � __rtruediv__� s zDecimal.__rtruediv__c C s� t |�}|tkr|S |d kr"t� }| j||�}|r:||fS | j|jA }| j� r~|j� rj|jtd�}||fS t| |jtd�fS |s�| s�|jt d�}||fS |jt
d|�|jtd�fS | j ||�\}}|j |�}||fS )Nzdivmod(INF, INF)zINF % xz divmod(0, 0)zx // 0zx % 0)
r� r� r r� r7 r� rd r r< r r
r� r� )r/ r� r0 r; r= Zquotientr� r( r( r* �
__divmod__� s0
zDecimal.__divmod__c C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r� )r/ r� r0 r( r( r* � __rdivmod__� s zDecimal.__rdivmod__c C s� t |�}|tkr|S |d kr"t� }| j||�}|r6|S | j� rJ|jtd�S |sj| r^|jtd�S |jtd�S | j||�d }|j |�}|S )NzINF % xzx % 0z0 % 0r. )
r� r� r r� r� rd r r r� r� )r/ r� r0 r; r� r( r( r* �__mod__� s"
zDecimal.__mod__c C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r� )r/ r� r0 r( r( r* �__rmod__� s zDecimal.__rmod__c C s� |d krt � }t|dd�}| j||�}|r.|S | j� rB|jtd�S |sb| rV|jtd�S |jtd�S |j� r|t| �}|j|�S t | j
|j
�}| s�t | j d|�}|j|�S | j
� |j
� }||jd kr�|jt�S |d kr�| j||j�}|j|�S t| �}t|�}|j|jk�r(| jd |j|j 9 _n| jd |j|j 9 _t|j|j�\}} d| |d@ |jk�r~| |j8 } |d7 }|d |j k�r�|jt�S | j }
| d
k �r�d|
}
| } t |
t| �|�}|j|�S ) NT)r� zremainder_near(infinity, x)zremainder_near(x, 0)zremainder_near(0, 0)rQ r. r, r� r& r� )r r� r� r� rd r r r r� r� rO r6 r7 r� r@ r r� r? ri rV rU r� r` ) r/ r� r0 r; �ideal_exponentr� r� r� r� r� r= r( r( r* �remainder_near sZ
zDecimal.remainder_nearc C s� t |�}|tkr|S |d kr"t� }| j||�}|r6|S | j� rb|j� rR|jtd�S t| j|jA S |s�| r�|jt d| j|jA �S |jt
d�S | j ||�d S )Nz
INF // INFzx // 0z0 // 0r& ) r� r� r r� r� rd r r< r7 r
r r� )r/ r� r0 r; r( r( r* � __floordiv__L s$ zDecimal.__floordiv__c C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r� )r/ r� r0 r( r( r* �
__rfloordiv__h s zDecimal.__rfloordiv__c C s8 | j � r(| j� rtd��| jr"dnd}nt| �}t|�S )Nz%Cannot convert signaling NaN to floatz-nan�nan)r� r� rl r7 r` rp )r/ �sr( r( r* � __float__o s zDecimal.__float__c C st | j r(| j� rtd��n| j� r(td��d| j }| jdkrT|t| j� d| j S |t| jd | j� pjd� S d S )Nz Cannot convert NaN to integerz"Cannot convert infinity to integerr. r& r� rQ r� ) rP r� rl r� r� r7 rO rU r8 )r/ r� r( r( r* �__int__y s
zDecimal.__int__c C s | S )Nr( )r/ r( r( r* �real� s z Decimal.realc C s t d�S )Nr& )r )r/ r( r( r* �imag� s z Decimal.imagc C s | S )Nr( )r/ r( r( r* � conjugate� s zDecimal.conjugatec C s t t| ��S )N)�complexrp )r/ r( r( r* � __complex__� s zDecimal.__complex__c C sR | j }|j|j }t|�|krJ|t|�| d � jd�}t| j|| jd�S t| �S )NrQ T) r8 r@ �clamprf rg r6 r7 rO r )r/ r0 ZpayloadZmax_payload_lenr( r( r* r9 � s zDecimal._fix_nanc C sX | j r | j� r| j|�S t| �S |j� }|j� }| s�|j|g|j }tt | j
|�|�}|| j
krx|j t � t
| jd|�S t| �S t| j�| j
|j }||kr�|j td| j�}|j t� |j t� |S ||k }|r�|}| j
|k �r�t| j�| j
| } | dk �rt
| jd|d �} d} | j|j }
|
| | �} | jd | � �p>d} | dk�r~tt| �d �} t| �|jk�r~| d d� } |d7 }||k�r�|j td| j�}nt
| j| |�}| �r�|�r�|j t� |�r�|j t� | �r�|j t� |j t� |�s�|j t � |S |�r |j t� |jdk�rP| j
|k�rP|j t � | jd| j
| }
t
| j|
|�S t| �S )NrQ z
above Emaxr& r� r. r� )rP r� r9 r r� �EtoprA r� r� r� rO rd r r6 r7 rf r8 r@ r r r �_pick_rounding_functionr? r` rU r r
)r/ r0 r� r� �exp_maxZnew_expZexp_minr; Zself_is_subnormalrx Zrounding_method�changedr� r� r( r( r* r� � sn
z Decimal._fixc C s t | j|�rdS dS d S )Nr& r. r� )�
_all_zerosr8 )r/ r@ r( r( r* � _round_down s zDecimal._round_downc C s | j |� S )N)r� )r/ r@ r( r( r* � _round_up s zDecimal._round_upc C s* | j | dkrdS t| j |�r"dS dS d S )NZ56789r. r& r� )r8 r� )r/ r@ r( r( r* �_round_half_up s
zDecimal._round_half_upc C s t | j|�rdS | j|�S d S )Nr. r� )� _exact_halfr8 r� )r/ r@ r( r( r* �_round_half_down s zDecimal._round_half_downc C s8 t | j|�r*|dks&| j|d dkr*dS | j|�S d S )Nr& r. �02468r� )r� r8 r� )r/ r@ r( r( r* �_round_half_even# s zDecimal._round_half_evenc C s | j r| j|�S | j|� S d S )N)r7 r� )r/ r@ r( r( r* �_round_ceiling+ s
zDecimal._round_ceilingc C s | j s| j|�S | j|� S d S )N)r7 r� )r/ r@ r( r( r* � _round_floor2 s
zDecimal._round_floorc C s0 |r | j |d dkr | j|�S | j|� S d S )Nr. Z05)r8 r� )r/ r@ r( r( r* � _round_05up9 s
zDecimal._round_05up)r r r r r r r r c C sb |d k r2t |t�std��tdd| �}| j|�S | jrR| j� rJtd��ntd��t| j dt
��S )Nz+Second argument to round should be integralr& r� zcannot round a NaNzcannot round an infinity) r_ rU rr r6 �quantizerP r� rl r� r� r )r/ r5 rV r( r( r* � __round__K s /
zDecimal.__round__c C s0 | j r | j� rtd��ntd��t| jdt��S )Nzcannot round a NaNzcannot round an infinityr& )rP r� rl r� rU r� r )r/ r( r( r* � __floor__� s
zDecimal.__floor__c C s0 | j r | j� rtd��ntd��t| jdt��S )Nzcannot round a NaNzcannot round an infinityr& )rP r� rl r� rU r� r )r/ r( r( r* �__ceil__� s
zDecimal.__ceil__c C s t |dd�}t |dd�}| js$|jr�|d kr2t� }| jdkrJ|jtd| �S |jdkrb|jtd|�S | jdkrr| }nf|jdkr�|}nV| jdkr�|s�|jtd�S t| j|jA }n*|jdkr�| s�|jtd�S t| j|jA }n0t| j|jA t t
| j �t
|j � �| j|j �}|j ||�S ) NT)r� rY r� r5 rZ zINF * 0 in fmaz0 * INF in fma)
r� rP r rO rd r r< r7 r6 r` rU r8 r� )r/ r� Zthirdr0 �productr( r( r* �fma� s6
z Decimal.fmac C s� t |�}|tkr|S t |�}|tkr(|S |d kr6t� }| j� }|j� }|j� }|sZ|sZ|r�|dkrp|jtd| �S |dkr�|jtd|�S |dkr�|jtd|�S |r�| j|�S |r�|j|�S |j|�S | j� o�|j� o�|j� s�|jtd�S |dk r�|jtd�S |�s |jtd�S |j� |j k�r(|jtd�S | �rD| �rD|jtd�S |j
� �rTd}n| j }t t
|��}t| j� �}t|j� �} |j
| td |j|� | }x t| j�D ]}
t|d |�}�q�W t|| j
|�}t|t|�d�S )
Nr, r� z@pow() 3rd argument not allowed unless all arguments are integersr& zApow() 2nd argument cannot be negative when 3rd argument specifiedz pow() 3rd argument cannot be 0zSinsufficient precision: pow() 3rd argument must not have more than precision digitszXat least one of pow() 1st argument and 2nd argument must be nonzero; 0**0 is not definedr� )r� r� r r� rd r r9 �
_isintegerr� r@ �_isevenr7 rh rU ri �to_integral_valuer� rV �ranger6 r` ) r/ r� �modulor0 r� r� Z
modulo_is_nanr= �base�exponent�ir( r( r* �
_power_modulo� sl
zDecimal._power_moduloc C s t | �}|j|j }}x |d dkr6|d }|d7 }qW t |�}|j|j }}x |d dkrn|d }|d7 }qPW |dk�r ||9 }x |d dkr�|d }|d7 }q�W |dk r�d S |d| } |jdkr�| } |j� o�|jdk�r| jt|� }
t| |
|d �} nd} tddd| | | �S |jdk�r�|d } | dk�r�|| @ |k�rPd S t |�d }
|d
d }|t
t |��k�r~d S t |
| |�}
t || |�}|
d k�s�|d k�r�d S |
|k�r�d S d |
}n�| d k�r�t |�d
d }
t
d |
|�\}}|�r d S x$|d dk�r$|d }|
d8 }
�qW |d d }|t
t |��k�rHd S t |
| |�}
t || |�}|
d k�sx|d k�r|d S |
|k�r�d S d|
}nd S |d| k�r�d S |
| }tdt |�|�S |dk�r�|d| d }}n�|dk�r t
t t|| ���| k�r d S t |�}|dk�r>t
t t|�| ��| k�r>d S |d| }}x:|d |d k�oldkn �r�|d }|d }�qPW x:|d |d k�o�dkn �r�|d }|d }�q�W |dk�rt|dk�r�||k�r�d S t
||�\}}|dk�rd S dt |� | > }x>t
|||d �\}}||k�r>P n||d | | }�qW ||k�oh|dk�spd S |}|dk�r�||d t|� k�r�d S || }||9 }|d| k�r�d S t |�}|j� �r�|jdk�r�| jt|� }
t||
|t
|� �} nd} td|d| || �S )Nr� r& r. r� rQ r, � r� � �] �A rz � r[ �d )r, r r� r )ri rU rV r= r� r7 rO r� r6 �_nbitsrf r` �_decimal_lshift_exactr� rh � _log10_lb)r/ r� �p�x�xc�xe�y�yc�yer r� ZzerosZ
last_digitr� Zemaxr� ru r5 Zxc_bits�rem�ar� r� Zstr_xcr( r( r* � _power_exact( s� :
&&&&
zDecimal._power_exactc C s@ |d k r| j |||�S t|�}|tkr*|S |d kr8t� }| j||�}|rL|S |sd| s`|jtd�S tS d}| jdkr�|j � r�|j
� s�d}n| r�|jtd�S | j � } | s�|jdkr�t |dd�S t
| S | j� r�|jdkr�t
| S t |dd�S | tk�r�|j � �rZ|jdk�rd}n ||jk�r"|j}nt|�}| j| }|d|j k �rxd|j }|jt� n |jt� |jt� d|j }t |dd| |�S | j� }|j� �r�|jdk|dk k�r�t |dd�S t
| S d }d} | j� |j� }
|dk|jdkk�r |
tt|j��k�rHt |d|jd �}n,|j� } |
tt| ��k�rHt |d| d �}|d k�r�| j||jd �}|d k �r�|dk�r�t d|j|j�}d} |d k�r:|j} t | �}
|
j|
j }}t |�}|j|j }}|j dk�r�| }d }xJt ||||| | �\}}|d
d tt|��| d �r P |d 7 }�q�W t |t|�|�}| �r2|j � �r2t|j�|jk�r�|jd t|j� }t |j|jd| |j| �}|j � }|j!� xt"D ]}d|j#|< �q�W |j$|�}|jt� |j%t& �r�|jt'� |j%t( �r |jt(d |j� x:t't&ttt)fD ] }|j%| �r|j|� �qW n
|j$|�}|S )
Nz0 ** 0r& r. z+x ** y with x negative and y not an integerrQ r� FTr[ rz r� z
above Emax)*r r� r� r r� rd r �_Oner7 r� r� r� r6 r< r� r@ rU rO r r r� �_log10_exp_boundrf r` rA r� r r8 ri rV r= �_dpowerrH rI �_signals�trapsr� �flagsr
r r r )r/ r� r� r0 r; Z result_signZ
multiplierrV Zself_adj�exactZboundr� r r
r r r r r �extrar� r� Z
newcontextZ exceptionr( r( r* �__pow__ s�
"
zDecimal.__pow__c C s" t |�}|tkr|S |j| |d�S )N)r0 )r� r� r )r/ r� r0 r( r( r* �__rpow__� s zDecimal.__rpow__c C s� |d krt � }| jr(| j|d�}|r(|S | j|�}|j� r>|S |sPt|jdd�S |j|j� g|j }t
|j �}|j }x.|j |d dkr�||k r�|d7 }|d8 }qvW t|j|j d |� |�S )N)r0 rQ r& r. )
r rP r� r� r� r6 r7 rA r� r� rf r8 rO )r/ r0 r; �dupr� �endrV r( r( r* � normalize� s$
zDecimal.normalizec C s� t |dd�}|d krt� }|d kr(|j}| js4|jr|| j||�}|rH|S |j� sX| j� r||j� rp| j� rpt| �S |jtd�S |j � |j
ko�|j kn s�|jtd�S | s�t | j
d|j
�}|j|�S | j� }||j kr�|jtd�S ||j
d |jk�r|jtd�S | j|j
|�}|j� |j k�r0|jtd�S t|j�|jk�rN|jtd�S |�rn|j� |jk �rn|jt� |j
| j
k�r�|| k�r�|jt� |jt� |j|�}|S ) NT)r� zquantize with one INFz)target exponent out of bounds in quantizerQ z9exponent of quantize result too large for current contextr. z7quantize result has too many digits for current context)r� r r? rP r� r� r rd r r� rO rA r6 r7 r� r� r@ r� rf r8 �Eminr
r r )r/ rV r? r0 r; r� r( r( r* r�
sT
zDecimal.quantizec C sD t |dd�}| js|jr8| j� r(|j� p6| j� o6|j� S | j|jkS )NT)r� )r� rP r� � is_infiniterO )r/ r� r0 r( r( r* � same_quantumJ
s
zDecimal.same_quantumc C s� | j rt| �S | s t| jd|�S | j|krHt| j| jd| j| |�S t| j�| j | }|dk rzt| jd|d �} d}| j| }|| |�}| jd |� p�d}|dkr�tt |�d �}t| j||�S )NrQ r& r� r. )
rP r r6 r7 rO r8 rf r� r` rU )r/ rV r? rx Z
this_functionr� r� r( r( r* r� Y
s"
zDecimal._rescalec C sh |dkrt d��| js | r$t| �S | j| j� d | |�}|j� | j� krd|j|j� d | |�}|S )Nr& z'argument should be at least 1 in _roundr. )rl rP r r� r� )r/ �placesr? r; r( r( r* �_round{
s
zDecimal._roundc C s� | j r"| j|d�}|r|S t| �S | jdkr4t| �S | sFt| jdd�S |d krTt� }|d krb|j}| jd|�}|| kr�|j t
� |j t � |S )N)r0 r& rQ ) rP r� r rO r6 r7 r r? r� rd r r )r/ r? r0 r; r( r( r* �to_integral_exact�
s$
zDecimal.to_integral_exactc C s` |d krt � }|d kr |j}| jr>| j|d�}|r6|S t| �S | jdkrPt| �S | jd|�S d S )N)r0 r& )r r? rP r� r rO r� )r/ r? r0 r; r( r( r* r� �
s
zDecimal.to_integral_valuec C s� |d krt � }| jrB| j|d�}|r(|S | j� rB| jdkrBt| �S | sdt| jd| jd �}|j|�S | jdkrz|j t
d�S |j d }t | �}|j
d? }|j
d@ r�|jd }t| j�d? d }n|j}t| j�d d? }|| }|dkr�|d| 9 }d } nt|d| �\}}
|
} ||8 }d| } x(|| } | | k�r:P n | | d? } �q$W | �o\| | |k} | �r�|dk�r|| d| } n| d| 9 } ||7 }n| d
dk�r�| d7 } tdt| �|�}|j� }|jt�}
|j|�}|
|_|S ) N)r0 r& rQ r, r. zsqrt(-x), x > 0r� r Trz )r rP r� r� r7 r r6 rO r� rd r r@ ri rV rU rf r8 r� r` �
_shallow_copy�
_set_roundingr r? )r/ r0 r; r@ �opr� �c�lr� r r� r5 r� r? r( r( r* �sqrt�
s`
z Decimal.sqrtc C s� t |dd�}|d krt� }| js&|jr~| j� }|j� }|s>|r~|dkrX|dkrX| j|�S |dkrr|dkrr|j|�S | j||�S | j|�}|dkr�| j|�}|dkr�|}n| }|j|�S )NT)r� r. r& r� )r� r rP r� r� r� r� �
compare_total)r/ r� r0 �sn�onr, r; r( r( r* r� % s&
z Decimal.maxc C s� t |dd�}|d krt� }| js&|jr~| j� }|j� }|s>|r~|dkrX|dkrX| j|�S |dkrr|dkrr|j|�S | j||�S | j|�}|dkr�| j|�}|dkr�| }n|}|j|�S )NT)r� r. r& r� )r� r rP r� r� r� r� r/ )r/ r� r0 r0 r1 r, r; r( r( r* r� O s&
z Decimal.minc C s8 | j r
dS | jdkrdS | j| jd � }|dt|� kS )NFr& TrQ )rP rO r8 rf )r/ �restr( r( r* r� q s
zDecimal._isintegerc C s( | s| j dkrdS | jd| j dkS )Nr& Tr. r� r� )rO r8 )r/ r( r( r* r� z s zDecimal._isevenc C s. y| j t| j� d S tk
r( dS X d S )Nr. r& )rO rf r8 rr )r/ r( r( r* r� � s zDecimal.adjustedc C s | S )Nr( )r/ r( r( r* � canonical� s zDecimal.canonicalc C s. t |dd�}| j||�}|r |S | j||d�S )NT)r� )r0 )r� r� r� )r/ r� r0 r; r( r( r* �compare_signal� s
zDecimal.compare_signalc C sf t |dd�}| jr |j r tS | j r0|jr0tS | j}| j� }|j� }|sP|�r||kr�t| j�| jf}t|j�|jf}||k r�|r�tS tS ||kr�|r�tS tS tS |r�|dkr�tS |dkr�tS |dkr�tS |dkr�tS n4|dkr�tS |dkr�tS |dk�r tS |dk�rtS | |k �r tS | |k�r*tS | j|jk �rF|�rBtS tS | j|jk�rb|�r^tS tS tS )NT)r� r. r, ) r� r7 � _NegativeOner r� rf r8 �_ZerorO )r/ r� r0 r= Zself_nanZ other_nanZself_keyZ other_keyr( r( r* r/ � sf
zDecimal.compare_totalc C s& t |dd�}| j� }|j� }|j|�S )NT)r� )r� r� r/ )r/ r� r0 r� �or( r( r* �compare_total_mag� s zDecimal.compare_total_magc C s t d| j| j| j�S )Nr& )r6 r8 rO rP )r/ r( r( r* r� � s zDecimal.copy_absc C s2 | j rtd| j| j| j�S td| j| j| j�S d S )Nr& r. )r7 r6 r8 rO rP )r/ r( r( r* r� � s zDecimal.copy_negatec C s" t |dd�}t|j| j| j| j�S )NT)r� )r� r6 r7 r8 rO rP )r/ r� r0 r( r( r* � copy_sign� s
zDecimal.copy_signc C s� |d krt � }| j|d�}|r"|S | j� d
kr2tS | s:tS | j� dkrNt| �S |j}| j� }| jdkr�|t t
|j d d ��kr�t dd|j d �}�n0| jdkr�|t t
|j
� d d ��kr�t dd|j
� d �}n�| jdko�|| k �r
t ddd|d d | �}n�| jdk�rB|| d k �rBt dd|d | d �}n�t| �}|j|j }}|jdk�rj| }d}xFt|||| �\} }
| dd t t
| ��| d �r�P |d7 }�qpW t dt
| �|
�}|j� }|jt�} |j|�}| |_|S ) N)r0 r. r& r[ r� rQ r> rz r� r� )r r� r� r6 r r r@ r� r7 rf r` rA r6 r� ri rU rV r= �_dexpr) r* r r� r? ) r/ r0 r; r �adjr+ r, r� r r� rV r? r( r( r* rV sJ $( "
z Decimal.expc C s dS )NTr( )r/ r( r( r* � is_canonicalL s zDecimal.is_canonicalc C s | j S )N)rP )r/ r( r( r* � is_finiteT s zDecimal.is_finitec C s
| j dkS )NrZ )rO )r/ r( r( r* r$ \ s zDecimal.is_infinitec C s
| j dkS )Nr5 rY )r5 rY )rO )r/ r( r( r* r� ` s zDecimal.is_nanc C s, | j s | rdS |d kr t� }|j| j� kS )NF)rP r r# r� )r/ r0 r( r( r* � is_normald s
zDecimal.is_normalc C s
| j dkS )Nr5 )rO )r/ r( r( r* r� l s zDecimal.is_qnanc C s
| j dkS )Nr. )r7 )r/ r( r( r* � is_signedp s zDecimal.is_signedc C s
| j dkS )NrY )rO )r/ r( r( r* r� t s zDecimal.is_snanc C s, | j s | rdS |d kr t� }| j� |jk S )NF)rP r r� r# )r/ r0 r( r( r* � is_subnormalx s
zDecimal.is_subnormalc C s | j o| jdkS )NrQ )rP r8 )r/ r( r( r* �is_zero� s zDecimal.is_zeroc C s� | j t| j� d }|dkr4tt|d d ��d S |dkrXttd| d d ��d S t| �}|j|j }}|dkr�t|d| �}t|�}t|�t|� ||k S |ttd| | �� d S )Nr. � r� r, r& r� r� )rO rf r8 r` ri rU rV )r/ r; r+ r, r� �num�denr( r( r* �
_ln_exp_bound� s zDecimal._ln_exp_boundc
C s |d krt � }| j|d�}|r"|S | s*tS | j� dkr:tS | tkrFtS | jdkr\|jt d�S t
| �}|j |j }}|j
}|| j� d }x>t|||�}|ddttt|���| d r�P |d7 }q�W tt |dk �tt|��| �}|j� }|jt�} |j|�}| |_|S ) N)r0 r. zln of a negative valuer, rz r� r[ r& )r r� �_NegativeInfinityr� � _Infinityr r6 r7 rd r ri rU rV r@ rE �_dlogrf r` rh r6 r) r* r r� r? )
r/ r0 r; r+ r, r� r r& r� r? r( r( r* �ln� s:
$
z
Decimal.lnc C s� | j t| j� d }|dkr,tt|��d S |dkrHttd| ��d S t| �}|j|j }}|dkr�t|d| �}td| �}t|�t|� ||k d S td| | �}t|�| |dk d S ) Nr. r, r& r� �� Z231r� r� )rO rf r8 r` ri rU rV )r/ r; r+ r, r� rC rD r( r( r* r � s zDecimal._log10_exp_boundc
C sH |d krt � }| j|d�}|r"|S | s*tS | j� dkr:tS | jdkrP|jtd�S | jd dkr�| jdd � dt | j�d kr�t
| j t | j� d �}n�t | �}|j
|j }}|j}|| j� d }x>t|||�}|dd t tt|���| d r�P |d
7 }q�W tt
|dk �tt|��| �}|j� }|jt�} |j|�}| |_|S ) N)r0 r. zlog10 of a negative valuer& r� rQ r, rz r� r[ )r r� rF r� rG r7 rd r r8 rf r rO ri rU rV r@ r �_dlog10r` rh r6 r) r* r r� r? )
r/ r0 r; r+ r, r� r r& r� r? r( r( r* �log10� s:
. $
z
Decimal.log10c C sV | j |d�}|r|S |d kr"t� }| j� r.tS | s@|jtdd�S t| j� �}|j|�S )N)r0 zlogb(0)r. ) r� r r� rG rd r
r r� r� )r/ r0 r; r( r( r* �logb
s z Decimal.logbc C s8 | j dks| jdkrdS x| jD ]}|dkr dS q W dS )Nr& FZ01T)r7 rO r8 )r/ �digr( r( r* �
_islogical>
s zDecimal._islogicalc C s� |j t|� }|dkr$d| | }n|dk r<||j d � }|j t|� }|dkr`d| | }n|dk rx||j d � }||fS )Nr& rQ )r@ rf )r/ r0 �opa�opbZdifr( r( r* �
_fill_logicalL
s zDecimal._fill_logicalc C s~ |d krt � }t|dd�}| j� s.|j� r8|jt�S | j|| j|j�\}}djdd� t||�D ��}t d|j
d�pxdd�S )NT)r� rS c S s$ g | ] \}}t t|�t|�@ ��qS r( )r` rU )�.0r �br( r( r* �
g
s z'Decimal.logical_and..r& rQ ) r r� rO rd r rR r8 rn �zipr6 rg )r/ r� r0 rP rQ r� r( r( r* � logical_andY
s
zDecimal.logical_andc C s( |d krt � }| jtdd|j d�|�S )Nr& r� )r � logical_xorr6 r@ )r/ r0 r( r( r* �logical_invertj
s zDecimal.logical_invertc C s~ |d krt � }t|dd�}| j� s.|j� r8|jt�S | j|| j|j�\}}djdd� t||�D ��}t d|j
d�pxdd�S )NT)r� rS c S s$ g | ] \}}t t|�t|�B ��qS r( )r` rU )rS r rT r( r( r* rU
s z&Decimal.logical_or..r& rQ ) r r� rO rd r rR r8 rn rV r6 rg )r/ r� r0 rP rQ r� r( r( r* �
logical_orq
s
zDecimal.logical_orc C s~ |d krt � }t|dd�}| j� s.|j� r8|jt�S | j|| j|j�\}}djdd� t||�D ��}t d|j
d�pxdd�S )NT)r� rS c S s$ g | ] \}}t t|�t|�A ��qS r( )r` rU )rS r rT r( r( r* rU �
s z'Decimal.logical_xor..r& rQ ) r r� rO rd r rR r8 rn rV r6 rg )r/ r� r0 rP rQ r� r( r( r* rX �
s
zDecimal.logical_xorc C s� t |dd�}|d krt� }| js&|jr~| j� }|j� }|s>|r~|dkrX|dkrX| j|�S |dkrr|dkrr|j|�S | j||�S | j� j|j� �}|dkr�| j|�}|dkr�|}n| }|j|�S )NT)r� r. r& r� ) r� r rP r� r� r� r� r� r/ )r/ r� r0 r0 r1 r, r; r( r( r* �max_mag�
s&
zDecimal.max_magc C s� t |dd�}|d krt� }| js&|jr~| j� }|j� }|s>|r~|dkrX|dkrX| j|�S |dkrr|dkrr|j|�S | j||�S | j� j|j� �}|dkr�| j|�}|dkr�| }n|}|j|�S )NT)r� r. r& r� ) r� r rP r� r� r� r� r� r/ )r/ r� r0 r0 r1 r, r; r( r( r* �min_mag�
s&
zDecimal.min_magc C s� |d krt � }| j|d�}|r"|S | j� dkr2tS | j� dkrTtdd|j |j� �S |j� }|jt � |j
� | j |�}|| kr�|S | j tdd|j
� d �|�S )N)r0 r. r& r> r� r� )r r� r� rF r6 r@ r� rH r* r �_ignore_all_flagsr� r� r� )r/ r0 r; �new_selfr( r( r* �
next_minus�
s"
zDecimal.next_minusc C s� |d krt � }| j|d�}|r"|S | j� dkr2tS | j� dkrTtdd|j |j� �S |j� }|jt � |j
� | j |�}|| kr�|S | j tdd|j
� d �|�S )N)r0 r. r> r& r� r� )r r� r� rG r6 r@ r� rH r* r r] r� r� r� )r/ r0 r; r^ r( r( r* � next_plus�
s"
zDecimal.next_plusc C s� t |dd�}|d krt� }| j||�}|r.|S | j|�}|dkrJ| j|�S |dkr^| j|�}n
| j|�}|j� r�|jt d|j
� |jt � |jt � nD|j
� |jk r�|jt� |jt� |jt � |jt � |s�|jt� |S )NT)r� r& r. z Infinite result from next_towardr� )r� r r� r� r9 r` r_ r� rd r r7 r r r� r# r r
r )r/ r� r0 r; Z
comparisonr( r( r* � next_toward�
s4
zDecimal.next_towardc C s� | j � r dS | j� rdS | j� }|dkr,dS |d
kr8dS | j� rN| jrJdS dS |d kr\t� }| j|d�rv| jrrd S d
S | jr�d S d S d S )Nr� r� r. z +Infinityz -Infinityz-Zeroz+Zero)r0 z
-Subnormalz
+Subnormalz-Normalz+Normalr� )r� r� r� rA r7 r r@ )r/ r0 �infr( r( r* � number_class+ s, zDecimal.number_classc C s t d�S )Nr� )r )r/ r( r( r* �radixU s z
Decimal.radixc C s� |d krt � }t|dd�}| j||�}|r.|S |jdkrB|jt�S |j t|� ko^|jkn sn|jt�S | j� r~t | �S t|�}| j
}|jt |� }|dkr�d| | }n|dk r�|| d � }||d � |d |� }t | j
|jd�p�d| j�S )NT)r� r& rQ )r r� r� rO rd r r@ rU r� r r8 rf r6 r7 rg )r/ r� r0 r; �torot�rotdig�topadZrotatedr( r( r* �rotateY s,
"
zDecimal.rotatec C s� |d krt � }t|dd�}| j||�}|r.|S |jdkrB|jt�S d|j|j }d|j|j }|t|� kox|kn s�|jt�S | j � r�t
| �S t | j | j
| jt|� �}|j|�}|S )NT)r� r& r, r� )r r� r� rO rd r rA r@ rU r� r r6 r7 r8 r� )r/ r� r0 r; ZliminfZlimsupr� r( r( r* �scalebz s"
zDecimal.scalebc C s |d krt � }t|dd�}| j||�}|r.|S |jdkrB|jt�S |j t|� ko^|jkn sn|jt�S | j� r~t | �S t|�}| j
}|jt |� }|dkr�d| | }n|dk r�|| d � }|dk r�|d |� }n |d| }||j d � }t | j
|jd��p d| j�S )NT)r� r& rQ )r r� r� rO rd r r@ rU r� r r8 rf r6 r7 rg )r/ r� r0 r; re rf rg Zshiftedr( r( r* r� � s2
"
z
Decimal.shiftc C s | j t| �ffS )N)� __class__r` )r/ r( r( r* �
__reduce__� s zDecimal.__reduce__c C s t | �tkr| S | jt| ��S )N)�typer rj r` )r/ r( r( r* �__copy__� s zDecimal.__copy__c C s t | �tkr| S | jt| ��S )N)rl r rj r` )r/ �memor( r( r* � __deepcopy__� s zDecimal.__deepcopy__c C sN |d krt � }t||d�}| jrXt| j|�}t| j� �}|d dkrL|d7 }t|||�S |d d krvddg|j |d< |d dkr�t | j| j
| j d �} |j }|d }|d k �r|d dkr�| j
|d |�} nF|d d
kr�| j| |�} n*|d d k�rt| j
�|k�r| j
||�} | �rB| j d k�rB|d d
k�rB| jd |�} | j t| j
� } |d dk�r�| �r||d k �r|d | }
nd }
nB|d d
k�r�| }
n.|d d k�r�| j d k�r�| dk�r�| }
nd }
|
d k �r�d} d|
| j
} nP|
t| j
�k�r| j
d|
t| j
� } d} n"| j
d |
� �p$d} | j
|
d � } | |
}
t| j| | |
|�S )N)� _localeconvrl �%�g�Gr, � precision�eEr. zfF%ZgGr& r� rQ rS i����)r �_parse_format_specifierrP � _format_signr7 r` r� �
_format_alignr� r6 r8 rO r? r' r� rf �_format_number)r/ Z specifierr0 rp �specr= �bodyr? rt r� r� rv rw rV r( r( r* �
__format__� sZ
"
zDecimal.__format__)rO r8 r7 rP )rQ N)NN)N)N)N)N)N)N)FN)N)N)N)TN)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)NN)N)N)NN)N)NN)NN)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)NN)�r2 r3 r4 � __slots__r^ � classmethodrq r� r� r� r� r� r� r� r� r� r� r� r� r� r� r} r� r� r� r� r� r� r� �__radd__r� r� r� �__rmul__r� r� r� r� r� r� r� r� r� r� r� r� � __trunc__r� �propertyr� r� r� r9 r� r� r� r� r� r� r� r� r� �dictr� r� r� r� r� r r r r r" r� r% r� r' r( r� � to_integralr. r� r� r� r� r� r3 r4 r/ r8 r� r� r9 rV r< r= r$ r� r> r� r? r� r@ rA rE rI r rL rM rO rR rW rY rZ rX r[ r\ r_ r` ra rc rd rh ri r� rk rm ro r| r( r( r( r* r + s
(
!@
2
4
V 7 ;! $ K
f > , U n Y = " c * " I K 2 3
. * ! 'Fc C s& t jt�}| |_||_||_||_|S )N)r] r^ r r7 r8 rO rP )r= Z coefficientr Zspecialr/ r( r( r* r6 s
r6 c @ s$ e Zd Zdd� Zdd� Zdd� ZdS )rN c C s |j � | _d S )N)rH � new_context)r/ r� r( r( r* �__init__9 s z_ContextManager.__init__c C s t � | _t| j� | jS )N)r �
saved_contextr r� )r/ r( r( r* � __enter__; s
z_ContextManager.__enter__c C s t | j� d S )N)r r� )r/ �t�v�tbr( r( r* �__exit__? s z_ContextManager.__exit__N)r2 r3 r4 r� r� r� r( r( r( r* rN 3 s rN c @ s� e Zd Zd�dd�Zdd� Zdd� Zdd � Zd
d � Zd d
� Zdd� Z dd� Z
dd� Z dd� Z dd� Z
e
Zd�dd�Zdd� Zd d � Zd d � ZdZd d!� Zd"d#� Zd$d%� Zd�d'd(�Zd)d*� Zd+d,� Zd-d.� Zd/d0� Zd1d2� Z d3d4� Z d5d6� Z d7d8� Z d9d:� Z d;d<� Z!d=d>� Z"d?d@� Z#dAdB� Z$dCdD� Z%dEdF� Z&dGdH� Z'dIdJ� Z(dKdL� Z)dMdN� Z*dOdP� Z+dQdR� Z,dSdT� Z-dUdV� Z.dWdX� Z/dYdZ� Z0d[d\� Z1d]d^� Z2d_d`� Z3dadb� Z4dcdd� Z5dedf� Z6dgdh� Z7didj� Z8dkdl� Z9dmdn� Z:dodp� Z;dqdr� Zdwdx� Z?dydz� Z@d{d|� ZAd}d~� ZBdd�� ZCd�d�� ZDd�d�� ZEd�d�� ZFd�d�d��ZGd�d�� ZHd�d�� ZId�d�� ZJd�d�� ZKd�d�� ZLd�d�� ZMd�d�� ZNd�d�� ZOd�d�� ZPd�d�� ZQd�d�� ZRd�d�� ZSd�d�� ZTd�d�� ZUeUZVdS )�r Nc
s> yt }
W n tk
r Y nX |d k r*|n|
j| _|d k r>|n|
j| _|d k rR|n|
j| _|d k rf|n|
j| _|d k rz|n|
j| _|d k r�|n|
j| _| d kr�g | _n| | _�d kr�|
j j
� | _ n.t �t �s�t �fdd�t
� D ��| _ n�| _ � d k�r
t jt
d�| _n0t � t ��s4t � fdd�t
� D ��| _n� | _d S )Nc 3 s | ]}|t |� k�fV qd S )N)rU )rS r� )r r( r* � n s z#Context.__init__..r& c 3 s | ]}|t |� k�fV qd S )N)rU )rS r� )r r( r* r� u s )r � NameErrorr@ r? r# rA r� r� �_ignored_flagsr rH r_ r� r �fromkeysr ) r/ r@ r? r# rA r� r� r r r� Zdcr( )r r r* r� U s.
zContext.__init__c C s� t |t�std| ��|dkr<||kr�td||||f ��nJ|dkrb||k r�td||||f ��n$||k sr||kr�td||||f ��tj| ||�S )Nz%s must be an integerz-infz %s must be in [%s, %d]. got: %srb z %s must be in [%d, %s]. got: %sz %s must be in [%d, %d]. got %s)r_ rU rr rl r] � __setattr__)r/ �namert ZvminZvmaxr( r( r* �_set_integer_checky s
zContext._set_integer_checkc C sh t |t�std| ��x |D ]}|tkr td| ��q W x tD ]}||kr>td| ��q>W tj| ||�S )Nz%s must be a signal dictz %s is not a valid signal dict)r_ r� rr r �KeyErrorr] r� )r/ r� r� �keyr( r( r* �_set_signal_dict� s
zContext._set_signal_dictc C s� |dkr| j ||dd�S |dkr0| j ||dd�S |dkrH| j ||dd�S |dkr`| j ||dd�S |d krx| j ||dd�S |d
kr�|tkr�td | ��tj| ||�S |d ks�|d
kr�| j||�S |dkr�tj| ||�S td| ��d S )Nr@ r. rb r# z-infr& rA r� r� r? z%s: invalid rounding moder r r� z.'decimal.Context' object has no attribute '%s')r� �_rounding_modesrr r] r� r� rL )r/ r� rt r( r( r* r� � s( zContext.__setattr__c C s t d| ��d S )Nz%s cannot be deleted)rL )r/ r� r( r( r* � __delattr__� s zContext.__delattr__c C sN dd� | j j� D �}dd� | jj� D �}| j| j| j| j| j| j| j ||ffS )Nc S s g | ]\}}|r|�qS r( r( )rS �sigr� r( r( r* rU � s z&Context.__reduce__..c S s g | ]\}}|r|�qS r( r( )rS r� r� r( r( r* rU � s )
r �itemsr rj r@ r? r# rA r� r� )r/ r r r( r( r* rk � s
zContext.__reduce__c C s| g }|j dt| � � dd� | jj� D �}|j ddj|� d � dd� | jj� D �}|j ddj|� d � dj|�d S )
NzrContext(prec=%(prec)d, rounding=%(rounding)s, Emin=%(Emin)d, Emax=%(Emax)d, capitals=%(capitals)d, clamp=%(clamp)dc S s g | ]\}}|r|j �qS r( )r2 )rS r r� r( r( r* rU � s z$Context.__repr__..zflags=[z, �]c S s g | ]\}}|r|j �qS r( )r2 )rS r� r� r( r( r* rU � s ztraps=[�))rm �varsr r� rn r )r/ r� �namesr( r( r* r� � s zContext.__repr__c C s x| j D ]}d| j |< qW d S )Nr& )r )r/ �flagr( r( r* rI � s zContext.clear_flagsc C s x| j D ]}d| j |< qW d S )Nr& )r )r/ r� r( r( r* � clear_traps� s zContext.clear_trapsc
C s. t | j| j| j| j| j| j| j| j| j � }|S )N)
r r@ r? r# rA r� r� r r r� )r/ �ncr( r( r* r) � s zContext._shallow_copyc
C s6 t | j| j| j| j| j| j| jj� | j j� | j
� }|S )N) r r@ r? r# rA r� r� r rH r r� )r/ r� r( r( r* rH � s
z Context.copyc G sZ t j||�}|| jkr(|� j| f|�� S d| j|< | j| sN|� j| f|�� S ||��d S )Nr. )�_condition_map�getr� r1 r r )r/ Z conditionZ explanationr) �errorr( r( r* rd � s
zContext._raise_errorc C s
| j t� S )N)�
_ignore_flagsr )r/ r( r( r* r] � s zContext._ignore_all_flagsc G s | j t|� | _ t|�S )N)r� rj )r/ r r( r( r* r� � s zContext._ignore_flagsc G s< |r t |d ttf�r |d }x|D ]}| jj|� q$W d S )Nr& )r_ rk rj r� �remove)r/ r r� r( r( r* �
_regard_flags� s
zContext._regard_flagsc C s t | j| j d �S )Nr. )rU r# r@ )r/ r( r( r* r� s z
Context.Etinyc C s t | j| j d �S )Nr. )rU rA r@ )r/ r( r( r* r� s z Context.Etopc C s | j }|| _ |S )N)r? )r/ rl r? r( r( r* r* s zContext._set_roundingrQ c C sj t |t�r*||j� ks d|kr*| jtd�S t|| d�}|j� r`t|j�| j | j
kr`| jtd�S |j | �S )NrR zAtrailing or leading whitespace and underscores are not permitted.)r0 z diagnostic info too long in NaN) r_ r` rb rd r r r� rf r8 r@ r� r� )r/ rC r� r( r( r* �create_decimal# s zContext.create_decimalc C s t j|�}|j| �S )N)r rq r� )r/ r r� r( r( r* �create_decimal_from_float4 s
z!Context.create_decimal_from_floatc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* rh F s z Context.absc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* �add[ s
z Context.addc C s t |j| ��S )N)r` r� )r/ r r( r( r* �_applyp s zContext._applyc C s t |t�std��|j� S )Nz,canonical requires a Decimal as an argument.)r_ r rr r3 )r/ r r( r( r* r3 s s
zContext.canonicalc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� � s ! zContext.comparec C s t |dd�}|j|| d�S )NT)r� )r0 )r� r4 )r/ r rT r( r( r* r4 � s zContext.compare_signalc C s t |dd�}|j|�S )NT)r� )r� r/ )r/ r rT r( r( r* r/ � s zContext.compare_totalc C s t |dd�}|j|�S )NT)r� )r� r8 )r/ r rT r( r( r* r8 � s zContext.compare_total_magc C s t |dd�}|j� S )NT)r� )r� r� )r/ r r( r( r* r� � s
zContext.copy_absc C s t |dd�}t|�S )NT)r� )r� r )r/ r r( r( r* � copy_decimal� s
zContext.copy_decimalc C s t |dd�}|j� S )NT)r� )r� r� )r/ r r( r( r* r� s
zContext.copy_negatec C s t |dd�}|j|�S )NT)r� )r� r9 )r/ r rT r( r( r* r9 s zContext.copy_signc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* �divide+ s
zContext.dividec C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* �
divide_intP s
zContext.divide_intc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* r� g s
zContext.divmodc C s t |dd�}|j| d�S )NT)r� )r0 )r� rV )r/ r r( r( r* rV | s z Context.expc C s t |dd�}|j||| d�S )NT)r� )r0 )r� r� )r/ r rT r, r( r( r* r� � s z Context.fmac C s t |t�std��|j� S )Nz/is_canonical requires a Decimal as an argument.)r_ r rr r< )r/ r r( r( r* r< � s
zContext.is_canonicalc C s t |dd�}|j� S )NT)r� )r� r= )r/ r r( r( r* r= � s zContext.is_finitec C s t |dd�}|j� S )NT)r� )r� r$ )r/ r r( r( r* r$ � s zContext.is_infinitec C s t |dd�}|j� S )NT)r� )r� r� )r/ r r( r( r* r� � s
zContext.is_nanc C s t |dd�}|j| d�S )NT)r� )r0 )r� r> )r/ r r( r( r* r> � s zContext.is_normalc C s t |dd�}|j� S )NT)r� )r� r� )r/ r r( r( r* r� s zContext.is_qnanc C s t |dd�}|j� S )NT)r� )r� r? )r/ r r( r( r* r? s zContext.is_signedc C s t |dd�}|j� S )NT)r� )r� r� )r/ r r( r( r* r� $ s
zContext.is_snanc C s t |dd�}|j| d�S )NT)r� )r0 )r� r@ )r/ r r( r( r* r@ 4 s zContext.is_subnormalc C s t |dd�}|j� S )NT)r� )r� rA )r/ r r( r( r* rA J s zContext.is_zeroc C s t |dd�}|j| d�S )NT)r� )r0 )r� rI )r/ r r( r( r* rI [ s z
Context.lnc C s t |dd�}|j| d�S )NT)r� )r0 )r� rL )r/ r r( r( r* rL q s z
Context.log10c C s t |dd�}|j| d�S )NT)r� )r0 )r� rM )r/ r r( r( r* rM � s z Context.logbc C s t |dd�}|j|| d�S )NT)r� )r0 )r� rW )r/ r rT r( r( r* rW � s zContext.logical_andc C s t |dd�}|j| d�S )NT)r� )r0 )r� rY )r/ r r( r( r* rY � s zContext.logical_invertc C s t |dd�}|j|| d�S )NT)r� )r0 )r� rZ )r/ r rT r( r( r* rZ � s zContext.logical_orc C s t |dd�}|j|| d�S )NT)r� )r0 )r� rX )r/ r rT r( r( r* rX � s zContext.logical_xorc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� s z Context.maxc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r[ )r/ r rT r( r( r* r[ & s zContext.max_magc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� 7 s z Context.minc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r\ )r/ r rT r( r( r* r\ R s zContext.min_magc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* �minusc s z
Context.minusc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* �multiplyt s
zContext.multiplyc C s t |dd�}|j| d�S )NT)r� )r0 )r� r_ )r/ r r( r( r* r_ � s zContext.next_minusc C s t |dd�}|j| d�S )NT)r� )r0 )r� r` )r/ r r( r( r* r` � s zContext.next_plusc C s t |dd�}|j|| d�S )NT)r� )r0 )r� ra )r/ r rT r( r( r* ra � s zContext.next_towardc C s t |dd�}|j| d�S )NT)r� )r0 )r� r" )r/ r r( r( r* r" � s zContext.normalizec C s t |dd�}|j| d�S )NT)r� )r0 )r� rc )r/ r r( r( r* rc � s / zContext.number_classc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* �plus) s z Context.plusc C s: t |dd�}|j||| d�}|tkr2td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r r� rr )r/ r rT r� r� r( r( r* �power: s
I z
Context.powerc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� � s 7 zContext.quantizec C s t d�S )Nr� )r )r/ r( r( r* rd � s z
Context.radixc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* r� � s
zContext.remainderc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� � s zContext.remainder_nearc C s t |dd�}|j|| d�S )NT)r� )r0 )r� rh )r/ r rT r( r( r* rh s zContext.rotatec C s t |dd�}|j|�S )NT)r� )r� r% )r/ r rT r( r( r* r% 1 s zContext.same_quantumc C s t |dd�}|j|| d�S )NT)r� )r0 )r� ri )r/ r rT r( r( r* ri I s zContext.scalebc C s t |dd�}|j|| d�S )NT)r� )r0 )r� r� )r/ r rT r( r( r* r� \ s z
Context.shiftc C s t |dd�}|j| d�S )NT)r� )r0 )r� r. )r/ r r( r( r* r. z s z Context.sqrtc C s8 t |dd�}|j|| d�}|tkr0td| ��n|S d S )NT)r� )r0 z Unable to convert %s to Decimal)r� r� r� rr )r/ r rT r� r( r( r* �subtract� s
zContext.subtractc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* r� � s zContext.to_eng_stringc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* �
to_sci_string� s zContext.to_sci_stringc C s t |dd�}|j| d�S )NT)r� )r0 )r� r( )r/ r r( r( r* r( � s zContext.to_integral_exactc C s t |dd�}|j| d�S )NT)r� )r0 )r� r� )r/ r r( r( r* r� � s zContext.to_integral_value) NNNNNNNNN)N)rQ )N)Wr2 r3 r4 r� r� r� r� r� rk r� rI r� r) rH rm rd r] r� r� r� r� r� r* r� r� rh r� r� r3 r� r4 r/ r8 r� r� r� r9 r� r� r� rV r� r< r= r$ r� r> r� r? r� r@ rA rI rL rM rW rY rZ rX r� r[ r� r\ r� r� r_ r` ra r" rc r� r� r� rd r� r� rh r% ri r� r. r� r� r� r( r� r� r( r( r( r* r B s�
"
$#
%
#2
P:&" c @ s& e Zd Zd Zd
dd�Zdd� ZeZdS ) ri r= rU rV Nc C sf |d kr d | _ d| _d | _nFt|t�rD|j| _ t|j�| _|j| _n |d | _ |d | _|d | _d S )Nr& r. r, )r= rU rV r_ r r7 r8 rO )r/ rt r( r( r* r� s
z_WorkRep.__init__c C s d| j | j| jf S )Nz (%r, %r, %r))r= rU rV )r/ r( r( r* r� ( s z_WorkRep.__repr__)r= rU rV )N)r2 r3 r4 r} r� r� r� r( r( r( r* ri s
ri c C s� | j |j k r|}| }n| }|}tt|j��}tt|j��}|j td|| d � }||j d |k rpd|_||_ | jd|j |j 9 _|j |_ | |fS )Nr. r, r� r� )rV rf r` rU r� )r� r� r@ Ztmpr� Ztmp_lenZ other_lenrV r( r( r* r� / s r� c C sb | dkr dS |dkr | d| S t t| ��}t|�t|jd�� }|| k rPd S | d| S d S )Nr& r� rQ )r` rh rf �rstrip)r5 r� Zstr_nZval_nr( r( r* r
O s r
c C sF | dks|dkrt d��d}x$||kr@||| | d? }}q W |S )Nr& z3Both arguments to _sqrt_nearest should be positive.r. )rl )r5 r rT r( r( r* �
_sqrt_nearestd s
r� c C s2 d|> | |? }}|d| |d @ |d@ |k S )Nr. r, r( )r
r� rT r� r( r( r* �_rshift_nearests s r� c C s& t | |�\}}|d| |d@ |k S )Nr, r. )r� )r rT r� r� r( r( r* � _div_nearest{ s r� r c C s� | | }d}xn||kr*t |�|| > |ksF||krzt |�|| ? |krzt|| d> |t||t||� |� �}|d7 }qW tdtt|�� d| � }t||�}t||�}x0t|d dd�D ] }t||�t|| |� }q�W t|| |�S )Nr& r. r� r[ i����r� )rh r� r� r� rU rf r` r� ) r
�M�Lr �R�TZyshift�wr� r( r( r* �_ilog� s
r� c
C s� |d7 }t t| ��}|| || dk }|dkr�d| }|| | }|dkrZ| d| 9 } nt| d| �} t| |�}t|�}t|| |�}|| } nd}t|d| �} t| | d�S )Nr, r. r& r� r )rf r` r� r� �
_log10_digits)
r, r� r r- r r� r� �log_dZlog_10Z log_tenpowerr( r( r* rK � s
rK c C s� |d7 }t t| ��}|| || dk }|dkrr|| | }|dkrR| d| 9 } nt| d| �} t| d| �}nd}|r�t tt|���d }|| dkr�t|t|| � d| �}q�d}nd}t|| d�S )Nr, r. r& r� r )rf r` r� r� rh r� ) r, r� r r- r r� r� r Z f_log_tenr( r( r* rH � s" rH c @ s e Zd Zdd� Zdd� ZdS )�
_Log10Memoizec C s
d| _ d S )NZ/23025850929940456840179914546843642076011014886)rx )r/ r( r( r* r� s z_Log10Memoize.__init__c C s� |dk rt d��|t| j�kr�d}xLd|| d }tttd| |�d��}|| d � d| krdP |d7 }q$W |jd�d d � | _t| jd |d � �S )
Nr& zp should be nonnegativer[ r� r, r rQ r. r� )rl rf rx r` r� r� r� rU )r/ r r r� rx r( r( r* � getdigits s z_Log10Memoize.getdigitsN)r2 r3 r4 r� r� r( r( r( r* r� s r� c C s� t | |> | �}tdtt|�� d| � }t| |�}||> }x.t|d dd�D ]}t| || || �}qRW x6t|d dd �D ]"}||d > }t||| |�}q�W || S )
Nr� r[ r. r& r, i����r� r� r� )r rU rf r` r� r� ) r
r� r� r� r� r ZMshiftr r� r( r( r* �_iexp& s
r� c C s� |d7 }t d|tt| �� d �}|| }|| }|dkrH| d| }n| d| }t|t|��\}}t|d| �}tt|d| �d�|| d fS )Nr, r& r. r� i� r[ )r� rf r` r� r� r� r� ) r, r� r r r� r� ZcshiftZquotr r( r( r* r: K s r: c C s� t tt|���| }t| ||| d �}|| }|dkrJ|| d| }nt|| d| �}|dkr�t t| ��| dk|dkkr�d|d d d| } }
q�d| d | } }
n,t||d |d �\} }
t| d�} |
d7 }
| |
fS )Nr. r& r� )rf r` rh rH r� r: ) r r r r r rT Zlxcr� Zpcr� rV r( r( r* r o s
r r �F �5 �( � rB � r� rz ) r� �2�3�4�5�6�7�8r> c C s0 | dkrt d��t| �}dt|� ||d S )Nr& z0The argument to _log10_lb should be nonnegative.r )rl r` rf )r, Z
correctionZstr_cr( r( r* r � s r c C sL t | t�r| S t | t�r t| �S |r8t | t�r8tj| �S |rHtd| ��tS )Nz Unable to convert %s to Decimal)r_ r rU rp rq rr r� )r� r� Z allow_floatr( r( r* r� � s
r� c C s� t |t�r| |fS t |tj�rR| jsDt| jtt| j �|j
�| j �} | t|j �fS |rrt |tj
�rr|jdkrr|j}t |t�r�t� }|r�d|jt< n |jtd� | tj|�fS ttfS )Nr& r. z;strict semantics for mixing floats and Decimals are enabled)r_ r �_numbersZRationalrP r6 r7 r` rU r8 � denominatorrO � numeratorZComplexr� r� rp r r r rd rq r� )r/ r� r� r0 r( r( r* r� � s$
r� r i?B )r@ r? r r rA r# r� r� r\ )r@ r? r r a� # A numeric string consists of:
# \s*
(?P[-+])? # an optional sign, followed by either...
(
(?=\d|\.\d) # ...a number (with at least one digit)
(?P\d*) # having a (possibly empty) integer part
(\.(?P\d*))? # followed by an optional fractional part
(E(?P[-+]?\d+))? # followed by an optional exponent, or...
|
Inf(inity)? # ...an infinity, or...
|
(?Ps)? # ...an (optionally signaling)
NaN # NaN
(?P\d*) # with (possibly empty) diagnostic info.
)
# \s*
\Z
z0*$z50*$z�\A
(?:
(?P.)?
(?P[<>=^])
)?
(?P[-+ ])?
(?P\#)?
(?P0)?
(?P(?!0)\d+)?
(?P,)?
(?:\.(?P0|(?!0)\d+))?
(?P[eEfFgGn%])?
\Z
c C s� t j| �}|d kr td| ��|j� }|d }|d }|d d k |d< |d rv|d k rbtd| ��|d k rvtd| ��|p|d|d< |p�d|d< |d d kr�d
|d < t|d p�d �|d < |d
d k r�t|d
�|d
< |d
dkr�|d d ks�|d dkr�d|d
< |d dk�rfd|d< |d k�r&tj� }|d d k �r@td| ��|d |d< |d |d< |d |d< n*|d d k�r|d|d< ddg|d< d|d< |S )NzInvalid format specifier: �fill�align�zeropadz7Fill character conflicts with '0' in format specifier: z2Alignment conflicts with '0' in format specifier: � �>r= rT � minimumwidthrQ rt r& rl ZgGnr. r5 rr �
thousands_sepzJExplicit thousands separator conflicts with 'n' type in format specifier: �grouping�
decimal_pointrS r[ r� )� _parse_format_specifier_regex�matchrl � groupdictrU �_locale�
localeconv)� format_specrp ru Z format_dictr� r� r( r( r* rv D sN
rv c C s� |d }|d }||t | � t |� }|d }|dkrF| | | }nj|dkr\|| | }nT|dkrr| | | }n>|dkr�t |�d }|d |� | | ||d � }ntd ��|S )
Nr� r� r� �u s~
&
.
^
0",#
%$+
*
P
%
)