Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TracyWidom beta=4 seems off by a constant #83

Open
alanedelman opened this issue Nov 7, 2024 · 2 comments
Open

TracyWidom beta=4 seems off by a constant #83

alanedelman opened this issue Nov 7, 2024 · 2 comments

Comments

@alanedelman
Copy link
Contributor

alanedelman commented Nov 7, 2024

I compared with TracyWidomBeta.jl and also the following random airy operator at beta=4
and those two matched. The code in this package is off by 1.14 or something, if it were the sqrt(2)
(1.41 not 1.14) I might have understood why. Perhaps something in the normalization?
Happy to understand better.

function stochastic(β,n)
    h=n^-(1/3)
    x=h:h:10
    N=length(x)
    d=(-2/h^2 .- x) + sqrt(4/(h*β))*randn(N) # diagonal
    e=fill(1/h^2,N-1) # subdiagonal
    eigmax(SymTridiagonal(d,e))
end
@dlfivefifty
Copy link
Member

I suspect the code is still what @jiahao wrote a decade ago.... so I don't know if anyone will remember why

@alanedelman
Copy link
Contributor Author

alanedelman commented Jan 16, 2025

so there is an oddball convention when beta=4 that seems to be out there
Bornemann mentions this in his footnote on page 4 of
https://arxiv.org/pdf/0904.1581.

That number 1.14 mentioned earlier is in fact 2^(1/6) ≈ 1.12.

This convention is in wikipedia and mathematica.

IMHO, it is a thoroughly ridiculous idea because it doesn't generalize to other betas though there are some arguments in favor of the convention.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants