Skip to content
Snippets Groups Projects
Commit fe7db2f9 authored by FROGE Ewen's avatar FROGE Ewen
Browse files

Update file Compute_Innovation.jl

parent ad1c464d
No related branches found
No related tags found
No related merge requests found
......@@ -83,8 +83,8 @@ y = (y .- mean(y, dims=2)) ./ std(y, dims=2) #Normalization over each realisatio
Tmin = p + 1 # First point for computing innovation
Analogs = Float64.(y[1, :])
KNNpreds = zeros(Nreal, T - Tmin + 1)
Threads.@threads for i in 1:Nreal
real = 2 * i + 2 # Pick every two realisations, starting at the third
Threads.@threads for i in 1:Nreal
real = 2 * i + 2 #Pick one every two realisation to prevent correlation, and start at the third realisation to prevent correlation to analog dictionary
println(Threads.threadid())
Realisation = Float64.(y[real, :])
KNNpreds[i, :] = Preds(Realisation, Analogs, Tmin, T, p, k)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment