gibbs_sample and update_gibbs_sample in algorithm 3.10 mutate the argument a. So maybe it would be clearer to call them gibbs_sample! and update_gibbs_sample!, according to julia convention of functions mutating their arguments.
Possibly also removing the explicit assignment a=gibbs_sample(...).
mykelk
commented
3 years ago
gibbs_sampleandupdate_gibbs_samplein algorithm 3.10 mutate the argumenta. So maybe it would be clearer to call themgibbs_sample!andupdate_gibbs_sample!, according to julia convention of functions mutating their arguments. Possibly also removing the explicit assignmenta=gibbs_sample(...).