o
    à8Va  ã                   @   sÊ   d Z ddlmZ ddlmZ ddlmZmZmZm	Z	 ddl
mZmZmZmZmZmZmZmZmZmZmZmZmZmZ ddlmZ eG dd„ deeƒƒZeG d	d
„ d
eeƒƒZejZeG dd„ deƒƒZ dS )zReal and complex elements. é    )ÚDomainElement)Úpublic)ÚPythonMPContextÚ_mpfÚ_mpcÚ	_constant)ÚMPZ_ONEÚfzeroÚfoneÚfinfÚfninfÚfnanÚround_nearestÚmpf_mulÚrepr_dpsÚ	int_typesÚfrom_intÚ
from_floatÚfrom_strÚto_rational)Úmpqc                   @   ó2   e Zd ZdZdZdd„ Zedd„ eƒZdd„ Zd	S )
ÚRealElementzAn element of a real domain. ©Z__mpf__c                 C   ó
   || _ d S ©Nr   ©ÚselfÚval© r   ú@/usr/lib/python3/dist-packages/sympy/polys/domains/mpelements.pyÚ_set_mpf   ó   
zRealElement._set_mpfc                 C   ó   | j S r   r   ©r   r   r   r    Ú<lambda>   ó    zRealElement.<lambda>c                 C   ó   | j jS r   ©ÚcontextZ_parentr$   r   r   r    Úparent   ó   zRealElement.parentN)	Ú__name__Ú
__module__Ú__qualname__Ú__doc__Ú	__slots__r!   ÚpropertyÚ_mpf_r*   r   r   r   r    r      ó    r   c                   @   r   )
ÚComplexElementz An element of a complex domain. ©Z__mpc__c                 C   r   r   r5   r   r   r   r    Ú_set_mpc"   r"   zComplexElement._set_mpcc                 C   r#   r   r5   r$   r   r   r    r%   %   r&   zComplexElement.<lambda>c                 C   r'   r   r(   r$   r   r   r    r*   '   r+   zComplexElement.parentN)	r,   r-   r.   r/   r0   r6   r1   Z_mpc_r*   r   r   r   r    r4      r3   r4   c                   @   sb   e Zd Zddd„Zdd„ Zdd	„ Zd
d„ Zdd„ Zedd„ ƒZ	edd„ ƒZ
ddd„Zddd„ZdS )Ú	MPContexté5   NFc                 C   sX  |t g| _|d u r|  |¡ n|  |¡ t| _t| _| jt| jg| j_	| jt| jg| j_	|r3| | j_
n| | j_
t| _| jt| jg| j_	| | j_
| j| j| jg| _d| _d| _|d u r`|  ¡ | _n|du rht| _n|  |¡| _|  | j¡| _| js|d| _ntd| j ƒ| _|  t¡| _|  t¡| _|  ttf¡| _|  t¡| _|  t¡| _ |  t!¡| _"d S )NTFi@B é   )#r   Ú_prec_roundingZ	_set_precZ_set_dpsr   Zmpfr4   ZmpcÚnewZ_ctxdatar)   r   ZconstantÚtypesZtrap_complexZprettyÚ	_make_tolÚtolr	   Ú_convert_tolÚmake_mpfÚ	toleranceÚ	max_denomÚintZzeror
   ZoneZmake_mpcÚjr   Úinfr   Zninfr   Únan)ÚctxÚprecZdpsr>   Úrealr   r   r    Ú__init__/   s@   


zMPContext.__init__c                 C   s    d}dt d| j df}t||ƒS )N)r   é   é   é   r   r9   )r   rH   r   )rG   ZhundredZepsr   r   r    r=   ^   s   
zMPContext._make_tolc                 C   s   |   |  ¡ ¡S r   )r@   r=   ©rG   r   r   r    Úmake_tolc   s   zMPContext.make_tolc                 C   s`   t |tƒr	t|ƒS t |tƒrt|ƒS t|dƒr|jS | j\}}t |tƒr*t	|||ƒS t
d| ƒ‚)Nr2   zexpected a real number, got %s)Ú
isinstancer   r   Úfloatr   Úhasattrr2   r:   Ústrr   Ú
ValueError)rG   r>   rH   Zroundingr   r   r    r?   f   s   




zMPContext._convert_tolc                 C   s   t dt|ƒ ƒ‚)Nzcannot create mpf from )Ú	TypeErrorÚrepr)rG   ÚxÚstringsr   r   r    Ú_convert_fallbackr   s   zMPContext._convert_fallbackc                 C   s
   t | jƒS r   )r   Z_precrN   r   r   r    Ú_repr_digitsu   s   
zMPContext._repr_digitsc                 C   r#   r   )Z_dpsrN   r   r   r    Ú_str_digitsy   s   zMPContext._str_digitsTc                 C   sü   t |jƒ\}}|r|| jkr||fS d\}}}}||}	}
	 |	|
 }|||  }|| jkr.n|||||  |f\}}}}|
|	||
  }	}
q| j| | }t||ƒ}t|||  |||  ƒ}t||ƒ}|rh|sl||fS t|| ƒt|| ƒkr{|jS |jS )N)r   r9   r9   r   )r   r2   rB   r   ÚabsZ_mpq_)rG   ÚsÚlimitÚpÚqZp0Zq0Zp1Zq1ÚnÚdÚaZq2ÚkZnumberZbound1Zbound2r   r   r    r   }   s,   

ú

zMPContext.to_rationalc           	      C   s¢   |   |¡}|d u r|d u r| jp|  ¡  }}|d u r |   |¡}n	|d u r)|   |¡}t|| ƒ}||kr5dS t|ƒ}t|ƒ}||k rI|| }||kS || }||kS )NT)ZconvertrA   rO   r\   )	rG   r]   ÚtZrel_epsZabs_epsZdiffZabssZabstÚerrr   r   r    Úalmosteq›   s"   

ÿzMPContext.almosteq)r8   NNF)T)NN)r,   r-   r.   rJ   r=   rO   r?   rY   r1   rZ   r[   r   rg   r   r   r   r    r7   ,   s    
/


r7   N)!r/   Z!sympy.polys.domains.domainelementr   Zsympy.utilitiesr   Zmpmath.ctx_mp_pythonr   r   r   r   Zmpmath.libmpr   r	   r
   r   r   r   r   r   r   r   r   r   r   r   Zmpmath.rationalr   r   r4   ÚobjectÚ__new__r;   r7   r   r   r   r    Ú<module>   s    @