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

Update file Compute_Innovation.jl

parent b8fb10aa
No related branches found
No related tags found
No related merge requests found
......@@ -83,15 +83,13 @@ 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)
@time begin
Threads.@threads for i in 1:Nreal
real = 2 * i + 2 # Pick every two realisations, starting at the third
println(Threads.threadid())
Realisation = Float64.(y[real, :])
KNNpreds[i, :] = Preds(Realisation, Analogs, Tmin, T, p, k)
end
end
# Save associated velocity and time
V = zeros(Nreal, T - Tmin + 1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment