PDSampler.jl icon indicating copy to clipboard operation
PDSampler.jl copied to clipboard

test for node with linear thinning

Open tlienart opened this issue 7 years ago • 0 comments

In the helper function ls_updatepq!, those lines:

bounce = fg.factors[fidx].nextevent(vcxf, vcvf)
acc    = bounce.dobounce(g, vcvf)
while !acc
    bounce = fg.factors[fidx].nextevent(vcxf, vcvf)
    acc    = bounce.dobounce(g, vcvf)
end
tauf = bounce.tau

the while !acc is for the case of thinning (corresponds to the accept/reject step in thinning). This is valid of course (same code as BPS) but so far the test case in LBPS has been a multivariate gaussian for which acc==true. For the sake of testing, one should

  • define a chain with simple factors
  • factors require thinning for sampling
  • test the helper function with that.

tlienart avatar Apr 10 '17 09:04 tlienart