Move delete_seeds to ext

ext/AdvancedPSLibtaskExt.jl

@@ -72,7 +72,7 @@ function AdvancedPS.fork(trace::LibtaskTrace, isref::Bool=false)
  newtrace = copy(trace)
  update_rng!(newtrace)
  isref && AdvancedPS.delete_retained!(newtrace.model.f)
- isref && AdvancedPS.delete_seeds!(newtrace)
+ isref && delete_seeds!(newtrace)
 
  # add backward reference
  addreference!(newtrace.model.ctask.task, newtrace)
@@ -204,4 +204,15 @@ function replay(particle::AdvancedPS.Particle)
  return trace
 end
 
+"""
+ delete_seeds!(particle::Particle)
+
+Truncate the seed history from the `particle` rng. When forking the reference Particle
+we need to keep the seeds up to the current model iteration but generate new seeds
+and random values afterward.
+"""
+function delete_seeds!(particle::AdvancedPS.Particle)
+ return particle.rng.keys = particle.rng.keys[1:(particle.rng.count - 1)]
+end
+
 end

src/model.jl

@@ -19,17 +19,6 @@ Base.copy(trace::Trace) = Trace(copy(trace.model), deepcopy(trace.rng))
 
 function observe end
 
-"""
- delete_seeds!(particle::Particle)
-
-Truncate the seed history from the `particle` rng. When forking the reference Particle
-we need to keep the seeds up to the current model iteration but generate new seeds
-and random values afterward.
-"""
-function delete_seeds!(particle::Particle)
- return particle.rng.keys = particle.rng.keys[1:(particle.rng.count - 1)]
-end
-
 """
  gen_refseed!(particle::Particle)
 

test/pgas.jl

@@ -106,4 +106,15 @@
 
  @test vals1 ≈ vals2
  end
+
+ @testset "smc sampler" begin
+ model = BaseModel(Params(0.9, 0.32, 1))
+ npart = 10
+
+ sampler = AdvancedPS.SMC(npart)
+ chains = sample(model, sampler)
+
+ @test length(chains.trajectories) == npart
+ @test length(chains.trajectories[1].model.X) == 3
+ end
 end