stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixing bugs in previous commit

time-shift
Michel Juillard 2012-08-30 12:44:46 +02:00
parent 1bac2d34c3
commit e8f159f66b
2 changed files with 16 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/independent_metropolis_hastings_core.m
View File

@ -105,13 +105,19 @@ for b = fblck:nblck,
% this will not work if the master uses a random generator not
% available in the slave (different Matlab version or
% Matlab/Octave cluster). Therefor the trap.
set_dynare_seed(options_.DynareRandomStreams.algo);
% this set the random generator type (the seed is useless but
% needed by the function)
set_dynare_seed(options_.DynareRandomStreams.algo,...
options_.DynareRandomStreams.seed);
% this set the state
set_dynare_random_generator_state(record.Seeds(b).Unifor, ...
record.Seeds(b).Normal);
catch
% if the state set by master is incompatible with the slave, we
% only reseed
set_dynare_seed(options_.DynareRandomStreams.seed+b);
end
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'])

9
matlab/random_walk_metropolis_hastings_core.m
View File

@ -139,10 +139,17 @@ for b = fblck:nblck,
% this will not work if the master uses a random generator not
% available in the slave (different Matlab version or
% Matlab/Octave cluster). Therefor the trap.
set_dynare_seed(options_.DynareRandomStreams.algo);
% this set the random generator type (the seed is useless but
% needed by the function)
set_dynare_seed(options_.DynareRandomStreams.algo,...
options_.DynareRandomStreams.seed);
% this set the state
set_dynare_random_generator_state(record.Seeds(b).Unifor, ...
record.Seeds(b).Normal);
catch
% if the state set by master is incompatible with the slave, we
% only reseed
set_dynare_seed(options_.DynareRandomStreams.seed+b);
end
if (options_.load_mh_file~=0) && (fline(b)>1) && OpenOldFile(b)