From f77b101d7b828172d04732ad1ba74adc7926b806 Mon Sep 17 00:00:00 2001 From: Michel Juillard Date: Wed, 29 Aug 2012 22:12:50 +0200 Subject: [PATCH] fixing bugs in random generator handling code --- matlab/random_walk_metropolis_hastings_core.m | 3 ++- matlab/set_dynare_random_generator_state.m | 2 +- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/matlab/random_walk_metropolis_hastings_core.m b/matlab/random_walk_metropolis_hastings_core.m index 7ff9afb0c..217328a93 100644 --- a/matlab/random_walk_metropolis_hastings_core.m +++ b/matlab/random_walk_metropolis_hastings_core.m @@ -135,7 +135,8 @@ jloop=0; JSUM = 0; for b = fblck:nblck, jloop=jloop+1; - set_dynare_seed(options_.DynareRandomStreams.seed+b); + set_dynare_random_generator_state(record.Seeds(b).Unifor, ... + record.Seeds(b).Normal); if (options_.load_mh_file~=0) && (fline(b)>1) && OpenOldFile(b) load(['./' MhDirectoryName '/' ModelName '_mh' int2str(NewFile(b)) ... '_blck' int2str(b) '.mat']) diff --git a/matlab/set_dynare_random_generator_state.m b/matlab/set_dynare_random_generator_state.m index 3b9d3c3a4..cbe3c4bb4 100644 --- a/matlab/set_dynare_random_generator_state.m +++ b/matlab/set_dynare_random_generator_state.m @@ -46,7 +46,7 @@ function [state_u,state_n] = get_dynare_random_generator_state(state_u,state_n) 'don''t understand the origin of the problem, ' ... 'please, contact Dynare''s development team.']) end - s_new = RandStream(s.type,'State',state_u); + s.State = state_u; if matlab_ver_less_than('7.12') RandStream.setDefaultStream(s); else