Home

Research

Divulgazione

Programs

Books

OEIS

Teaching2425

Links

Papers: List Papers; (with Abstracts); Curriculum (in Italian): long version ; short version; in English: long version. Google Scholar profile. ResearchGate page. Orcid ID. Scopus Author ID. Web of Science Researcher ID, Mathematical Reviews page, Zentralblatt page, IRIS-CINECA bibliometric parameters (italian ASN) [2023]. English C1 badge.

A) specialfunct_package-final.gp: Pari/GP script. It can be used via gp2c. Please read the description at the top of the gp-file. Toward the bottom of the gp-file you will find some correctness and performances tests.
It contains several functions:

1) logΓ(x), for x>0;
2) ψ(x), for x>0;
3) Bateman G-function, for x>0;
4) Hurwitz ζ(s,x), for x>0, s>1, (in another file named hurwitz_forfft.gp there's version with x in (0,1) for FFT applications);
4*) ONLY in another file: Hurwitz ζ'(s,x),x in (0,1), s>1 (ONLY in another file hurwitz_forfft.gp there's version with x in (0,1) for FFT applications);
5) Dirichlet β(s), s>1;
6) Dirichlet β'(s), s>1;
7) Dirichlet β values (s), s>1 (it computes β'/ β(s), β(s) and β'(s) at the same time);
8) Catalan constant computation.

b) hurwitz_forfft.gp: Pari/GP script. It can be used via gp2c. Please read the description at the top of the gp-file. Toward the bottom of the gp-file you will find some correctness and performances tests. All these programs require in input a file named primroot.res containing:
on the first line: a prime number q; on the second line: a primitive root g of q.
It contains several functions:

1) zetah(s, k₁, k₂), s>1 (for FFT applications; results saved on file); it computes: ζ(s,(g^k mod q)/q) for k₁≤k≤k₂;
2) zetahprime(s, k₁, k₂), s>1 (for FFT applications; results saved on file) the first derivative of the Hurwitz ζ function: ζ'(s,(g^k mod q)/q) for k₁≤k≤k₂;
3) zetahtotal(s, k₁, k₂), s>1 (for FFT applications; results saved on file); ζ(s,(g^k mod q)/q) and ζ'(s,(g^k mod q)/q) for k₁≤k≤k₂;
4) hurwitz_DIF_even (compute both ζ-Hurwitz and ζ'-Hurwitz with reflection fornulae; even case; Propositions 9 and 12 of [1]; results saved on file);
it computes ζ(s,(g^k mod q)/q) + ζ(s,1-(g^k mod q)/q) for k₁≤k≤k₂ and ζ'(s,(g^k mod q)/q) + ζ'(s,1-(g^k mod q)/q) for k₁≤k≤k₂;
5) hurwitz_DIF_odd (compute both ζ-Hurwitz and ζ'-Hurwitz with reflection fornulae; odd case; Propositions 9 and 12 of [1]; results saved on file);
it computes ζ(s,(g^k mod q)/q) - ζ(s,1-(g^k mod q)/q) for k₁≤k≤k₂ and ζ'(s,(g^k mod q)/q) - ζ'(s,1-(g^k mod q)/q) for k₁≤k≤k₂.

C) eta.gp: Pari/GP script. It can be used via gp2c. Please read the description at the top of the gp-file. Toward the bottom of the gp-file you will find some correctness and performances tests. It contains the scripts to compute the alternating Hurwitz zeta-function η(s,x); for x>0, s>1.

[1] A. Languasco - A unified strategy to compute some special functions of number-theoretic interest - J. Number Theory 247 (2023), 118--161.