JuliaMath / JuliaMath/SpecialFunctions.jl

Generalized zeta function cutoff causes discontinuity

Open
#488 0 comments 1 reaction 0 assignees View on GitHub
Dominant language
Julia
Stars
381
Forks
113
Avg merge
1d 23h
Merged PRs (30d)
1

Description

Version Information:
Julia 1.11.4
SpecialFunctions v2.5.0

`zeta(s, z)` gives incorrect values in the region $\Re (z) \geq 7 + \Re (s - 1) + |\Im (s) |$. For example:

```
julia> zeta(-13+7im, 0.5) # should be ~350.218774132674 + 259.188383123717im
-3371.4397140865763 + 2764.0891272966132im
```

The function jumps sharply at the boundary of the affected region and no longer behaves as expected. You can see the effects in these plots (using z=0.5):

![Image](https://github.com/user-attachments/assets/62669f54-1723-45bc-bc5d-758cfae9ad60)

If my understanding is correct, the Euler-Maclaurin expansion used should be well-defined for $\Re(s) \gt -17$ or so. However, I think the initial sum is prematurely truncated to zero terms when the cutoff is checked.

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.