From a7d640be84a8b25730504824e4429c211f39188c Mon Sep 17 00:00:00 2001 From: odunbar Date: Sun, 17 May 2026 15:42:36 -0700 Subject: [PATCH 1/5] fixed bug with 1d --- src/Utilities.jl | 2 +- test/MarkovChainMonteCarlo/runtests.jl | 14 ++++++++++++++ 2 files changed, 15 insertions(+), 1 deletion(-) diff --git a/src/Utilities.jl b/src/Utilities.jl index 01b8d02bf..f9a7a7f6f 100644 --- a/src/Utilities.jl +++ b/src/Utilities.jl @@ -1163,7 +1163,7 @@ function create_noise_injector( end C = cov(prior) - m = reshape(mean(prior), :, 1) + m = (ndims(prior) > 1) ? reshape(mean(prior), :, 1) : [mean(prior)] E = Matrix(E) enc_m = E * m + b diff --git a/test/MarkovChainMonteCarlo/runtests.jl b/test/MarkovChainMonteCarlo/runtests.jl index 5bf6d05ee..46751d9e5 100644 --- a/test/MarkovChainMonteCarlo/runtests.jl +++ b/test/MarkovChainMonteCarlo/runtests.jl @@ -624,6 +624,20 @@ end @test noise_injector.use_noise # check for noise_injector_threshold @test noise_injector.scaling == 0.5 + # check 1D + input_dim = 1 + n_samples = 10 + prior_1d = constrained_gaussian("1d-check", 0, 1, -Inf, 5) + in_data = sample(prior_1d, n_samples) + out_data = sample(prior_1d, n_samples) + io_pairs_1d = PairedDataContainer(in_data, out_data, data_are_columns = true) + + # lossless encoding + lossless_sch = create_encoder_schedule((minmax_scale(), "in")) + initialize_and_encode_with_schedule!(lossless_sch, io_pairs_1d; prior_cov = cov(prior_1d)) + + noise_injector = create_noise_injector(lossless_sch, prior_1d, 0.0, 0.5) + end From ec1325ff6ca0e2656b6f1aa7b42b8ac61a000e41 Mon Sep 17 00:00:00 2001 From: odunbar Date: Sun, 17 May 2026 16:32:54 -0700 Subject: [PATCH 2/5] typo --- test/MarkovChainMonteCarlo/runtests.jl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/test/MarkovChainMonteCarlo/runtests.jl b/test/MarkovChainMonteCarlo/runtests.jl index 46751d9e5..925ef1f27 100644 --- a/test/MarkovChainMonteCarlo/runtests.jl +++ b/test/MarkovChainMonteCarlo/runtests.jl @@ -8,6 +8,7 @@ using CalibrateEmulateSample.EnsembleKalmanProcesses using CalibrateEmulateSample.MarkovChainMonteCarlo const MCMC = MarkovChainMonteCarlo using CalibrateEmulateSample.ParameterDistributions +const PD = CalibrateEmulateSample.ParameterDistributions using CalibrateEmulateSample.Emulators using CalibrateEmulateSample.DataContainers using CalibrateEmulateSample.Utilities @@ -628,8 +629,8 @@ end input_dim = 1 n_samples = 10 prior_1d = constrained_gaussian("1d-check", 0, 1, -Inf, 5) - in_data = sample(prior_1d, n_samples) - out_data = sample(prior_1d, n_samples) + in_data = PD.sample(prior_1d, n_samples) + out_data = PD.sample(prior_1d, n_samples) io_pairs_1d = PairedDataContainer(in_data, out_data, data_are_columns = true) # lossless encoding From 2da150ba5e4dfdffd186a9e9339b7136deee2758 Mon Sep 17 00:00:00 2001 From: odunbar Date: Sun, 17 May 2026 16:36:02 -0700 Subject: [PATCH 3/5] typo --- test/MarkovChainMonteCarlo/runtests.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/MarkovChainMonteCarlo/runtests.jl b/test/MarkovChainMonteCarlo/runtests.jl index 925ef1f27..ab954c002 100644 --- a/test/MarkovChainMonteCarlo/runtests.jl +++ b/test/MarkovChainMonteCarlo/runtests.jl @@ -8,7 +8,7 @@ using CalibrateEmulateSample.EnsembleKalmanProcesses using CalibrateEmulateSample.MarkovChainMonteCarlo const MCMC = MarkovChainMonteCarlo using CalibrateEmulateSample.ParameterDistributions -const PD = CalibrateEmulateSample.ParameterDistributions +const PD = ParameterDistributions using CalibrateEmulateSample.Emulators using CalibrateEmulateSample.DataContainers using CalibrateEmulateSample.Utilities From 892d2c29408bc8e0f91d7b48751ce8de6dd701a0 Mon Sep 17 00:00:00 2001 From: odunbar Date: Sun, 17 May 2026 17:20:55 -0700 Subject: [PATCH 4/5] julia formatting removes [x;;] to [x] notation... --- src/Utilities.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/Utilities.jl b/src/Utilities.jl index f9a7a7f6f..21001a828 100644 --- a/src/Utilities.jl +++ b/src/Utilities.jl @@ -1163,7 +1163,7 @@ function create_noise_injector( end C = cov(prior) - m = (ndims(prior) > 1) ? reshape(mean(prior), :, 1) : [mean(prior)] + m = (ndims(prior) > 1) ? reshape(mean(prior), :, 1) : fill(mean(prior),1,1) E = Matrix(E) enc_m = E * m + b From e0a37fa100961d304590cffd748e8a4b7012d3d2 Mon Sep 17 00:00:00 2001 From: odunbar Date: Sun, 17 May 2026 17:27:50 -0700 Subject: [PATCH 5/5] run format --- src/Utilities.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/Utilities.jl b/src/Utilities.jl index 21001a828..2ebce7011 100644 --- a/src/Utilities.jl +++ b/src/Utilities.jl @@ -1163,7 +1163,7 @@ function create_noise_injector( end C = cov(prior) - m = (ndims(prior) > 1) ? reshape(mean(prior), :, 1) : fill(mean(prior),1,1) + m = (ndims(prior) > 1) ? reshape(mean(prior), :, 1) : fill(mean(prior), 1, 1) E = Matrix(E) enc_m = E * m + b