| Server IP : 35.80.110.71 / Your IP : 216.73.216.221 Web Server : Apache/2.4.58 (Ubuntu) System : Linux ip-172-31-21-44 6.17.0-1019-aws #19~24.04.1-Ubuntu SMP Tue Jun 23 18:53:06 UTC 2026 x86_64 User : ubuntu ( 1000) PHP Version : 8.3.31 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : ON | Perl : ON | Python : OFF | Sudo : ON | Pkexec : OFF Directory : /lib/python3/dist-packages/mpmath/functions/__pycache__/ |
Upload File : |
�
�"`h� � � � d Z ddlZ G d� de� ZddlmZ d� Zd� Zd � Zd
� Z dd�Z
dgfd�Zdd
�Zdd�Z
dd�Zdd�Zedd�� Zedd�� Zy)aI
---------------------------------------------------------------------
.. sectionauthor:: Juan Arias de Reyna <arias@us.es>
This module implements zeta-related functions using the Riemann-Siegel
expansion: zeta_offline(s,k=0)
* coef(J, eps): Need in the computation of Rzeta(s,k)
* Rzeta_simul(s, der=0) computes Rzeta^(k)(s) and Rzeta^(k)(1-s) simultaneously
for 0 <= k <= der. Used by zeta_offline and z_offline
* Rzeta_set(s, derivatives) computes Rzeta^(k)(s) for given derivatives, used by
z_half(t,k) and zeta_half
* z_offline(w,k): Z(w) and its derivatives of order k <= 4
* z_half(t,k): Z(t) (Riemann Siegel function) and its derivatives of order k <= 4
* zeta_offline(s): zeta(s) and its derivatives of order k<= 4
* zeta_half(1/2+it,k): zeta(s) and its derivatives of order k<= 4
* rs_zeta(s,k=0) Computes zeta^(k)(s) Unifies zeta_half and zeta_offline
* rs_z(w,k=0) Computes Z^(k)(w) Unifies z_offline and z_half
----------------------------------------------------------------------
This program uses Riemann-Siegel expansion even to compute
zeta(s) on points s = sigma + i t with sigma arbitrary not
necessarily equal to 1/2.
It is founded on a new deduction of the formula, with rigorous
and sharp bounds for the terms and rest of this expansion.
More information on the papers:
J. Arias de Reyna, High Precision Computation of Riemann's
Zeta Function by the Riemann-Siegel Formula I, II
We refer to them as I, II.
In them we shall find detailed explanation of all the
procedure.
The program uses Riemann-Siegel expansion.
This is useful when t is big, ( say t > 10000 ).
The precision is limited, roughly it can compute zeta(sigma+it)
with an error less than exp(-c t) for some constant c depending
on sigma. The program gives an error when the Riemann-Siegel
formula can not compute to the wanted precision.
� Nc � � e Zd Zd� Zy)�RSCachec � � ddi i g| _ y )Nr �
)� _rs_cache)�ctxs �9/usr/lib/python3/dist-packages/mpmath/functions/rszeta.py�__init__zRSCache.__init__6 s � ��B��B���
� N)�__name__�
__module__�__qualname__r
� r r r r 5 s � �(r r � )�defunc �X � |dz }|dz }t | j d|dz z � d|z dz | j |� z
� }| j d|z � }t | j d|z � | j |� dz |z � }|| _ i }| j |d<