JuliaMath / JuliaMath/RandomMatrices.jl
TracyWidom beta=4 seems off by a constant
Đang mở
- Ngôn ngữ chính
- Julia
- Star
- 94
- Fork
- 41
- Chỉ số merge pull request
- Không có pull request nào được merge trong 30 ngày
Mô tả
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
```
Hướng dẫn đóng góp
Chưa lập chỉ mục được hướng dẫn đóng góp cho kho mã nguồn này
Đánh giá
Issue này chưa được đánh giá.