Merge pull request #438 from ReactiveBayes/rxinfer-update-with-refactor-cvi-projection

test(fix): no marginalsamples option anymore

Author: bvdmitri

Project.toml

@@ -1,7 +1,7 @@
 name = "RxInfer"
 uuid = "86711068-29c9-4ff7-b620-ae75d7495b3d"
 authors = ["Bagaev Dmitry <d.v.bagaev@tue.nl> and contributors"]
-version = "4.1.0"
+version = "4.2.0"
 
 [deps]
 BayesBase = "b4ee3484-f114-42fe-b91c-797d54a0c67e"
@@ -56,7 +56,7 @@ Preferences = "1.4.3"
 PrettyTables = "2"
 ProgressMeter = "1.0.0"
 Random = "1.9"
-ReactiveMP = "~5.1.0"
+ReactiveMP = "~5.2.0"
 Reexport = "1.2.0"
 Rocket = "1.8.0"
 Static = "0.8.10, 1"

codemeta.json

@@ -9,12 +9,12 @@
     "downloadUrl": "https://github.com/reactivebayes/RxInfer.jl/releases",
     "issueTracker": "https://github.com/reactivebayes/RxInfer.jl/issues",
     "name": "RxInfer.jl",
-    "version": "4.1.0",
+    "version": "4.2.0",
     "description": "Julia package for automated, scalable and efficient Bayesian inference on factor graphs with reactive message passing. ",
     "applicationCategory": "Statistics",
     "developmentStatus": "active",
     "readme": "https://rxinfer.ml",
-    "softwareVersion": "4.1.0",
+    "softwareVersion": "4.2.0",
     "keywords": [
         "Bayesian inference",
         "message passing",

test/ext/ProjectionExt/inference_with_projection_tests.jl

@@ -411,9 +411,15 @@ end
         model = mymodel(C = C), data = (y = y,), meta = mymeta(), constraints = myconstraints(), initialization = myinitialization(), free_energy = true, iterations = 40
     )
 
-    @test mean(result.posteriors[:a][end]) ≈ a atol = 3e-2
-    @test foo(mean(result.posteriors[:a][end]), mean(result.posteriors[:b][end])) ≈ foo(a, b) atol = 3e-2
-    @test mean(result.posteriors[:μ][end]) ≈ foo(a, b) atol = 3e-2
+    conf_bound_a = 3 * std(result.posteriors[:a][end])
+    @test mean(result.posteriors[:a][end]) - conf_bound_a < a < mean(result.posteriors[:a][end]) + conf_bound_a
+
+    conf_bound_b = 3 * std(result.posteriors[:b][end])
+    @test mean(result.posteriors[:b][end]) - conf_bound_b < b < mean(result.posteriors[:b][end]) + conf_bound_b
+
+    @test mean(result.posteriors[:a][end]) ≈ a atol = 1e-2
+    @test foo(mean(result.posteriors[:a][end]), mean(result.posteriors[:b][end])) ≈ foo(a, b) atol = 1e-2
+    @test mean(result.posteriors[:μ][end]) ≈ foo(a, b) atol = 2e-2
     @test first(result.free_energy) > last(result.free_energy)
     @test count(<(0), diff(result.free_energy)) > 0.95
 
@@ -469,7 +475,7 @@ end
     end
 
     @meta function mymeta()
-        foo() -> CVIProjection(rng = StableRNG(42), marginalsamples = 20, outsamples = 5)
+        foo() -> CVIProjection(rng = StableRNG(42), sampling_strategy = FullSampling(10), outsamples = 5)
     end
 
     result = infer(