From eaedc6c99abc98e83b105d946d3f458a01f61588 Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@jrc.ec.europa.eu>
Date: Tue, 25 Feb 2014 09:34:31 +0100
Subject: [PATCH] Fixed bug when the last saved metropolis block is complete
and a metropolis is continued with load_mh_file
(cherry picked from commit 5ae06a211c8d398833a5dd80666eef0c1b029f97)
---
matlab/metropolis_hastings_initialization.m | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/matlab/metropolis_hastings_initialization.m b/matlab/metropolis_hastings_initialization.m
index 7872e591e..8c6ed5a9f 100644
--- a/matlab/metropolis_hastings_initialization.m
+++ b/matlab/metropolis_hastings_initialization.m
@@ -243,13 +243,14 @@ elseif options_.load_mh_file && ~options_.mh_recover
LastLineNumber = record.LastLineNumber;
if LastLineNumber < MAX_nruns
NewFile = ones(nblck,1)*LastFileNumber;
+ fline = ones(nblck,1)*(LastLineNumber+1);
else
NewFile = ones(nblck,1)*(LastFileNumber+1);
+ fline = ones(nblck,1);
end
ilogpo2 = record.LastLogPost;
ix2 = record.LastParameters;
fblck = 1;
- fline = ones(nblck,1)*(LastLineNumber+1);
NumberOfPreviousSimulations = sum(record.MhDraws(:,1),1);
fprintf('Estimation::mcmc: I am writting a new mh-history file... ');
record.MhDraws = [record.MhDraws;zeros(1,3)];
--
GitLab