From c12942d8aa11af3242fb0e9c4581cce7c0829f44 Mon Sep 17 00:00:00 2001 From: Matt Graham Date: Fri, 1 Dec 2023 18:31:23 +0000 Subject: [PATCH] Fix Lorenz 63 FillArrays bug --- test/models/lorenz63.jl | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/test/models/lorenz63.jl b/test/models/lorenz63.jl index 338c112..d960a17 100644 --- a/test/models/lorenz63.jl +++ b/test/models/lorenz63.jl @@ -2,7 +2,6 @@ module Lorenz63 using Base.Threads using Distributions -using FillArrays using HDF5 using Random using PDMats @@ -57,7 +56,7 @@ function init( Tsit5(); save_everystep=false ) - for u in eachcol(Matrix{S}(undef, 3, n_tasks)) + for u in eachcol(zeros(S, 3, n_tasks)) ] state_dimension = 3 observation_dimension = length(parameters.observed_indices) @@ -67,9 +66,9 @@ function init( ( MvNormal(m, isa(s, Vector) ? PDiagMat(s.^2) : ScalMat(length(m), s.^2)) for (m, s) in ( - (Ones{S}(state_dimension), parameters.initial_state_std), - (Zeros{S}(state_dimension), parameters.state_noise_std), - (Zeros{T}(observation_dimension), parameters.observation_noise_std), + (ones(S, state_dimension), parameters.initial_state_std), + (zeros(S, state_dimension), parameters.state_noise_std), + (zeros(T, observation_dimension), parameters.observation_noise_std), ) )... )