Closed mmulvahill closed 8 years ago
Results weren't repeatable on any system. The global counters for the mh algos were retaining there last value from one run to the next. Changed from global to vars local to mcmc() and passes as pointers to each mh algo. Working now. See commit https://github.com/mmulvahill/pulsatile/commit/1b00457786917aab1ebda7409015ea5a234ee9f2 and branch: fix-repeatability
Check if need to use GetRNGstate() and PutRNGstate() after each RNG call -- currently doing only Get prior to first RNG call and Put after last.