From 2ea112446f6e4e67587c2e24b532112328b7a7a5 Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@jrc.ec.europa.eu>
Date: Fri, 25 Sep 2015 17:09:07 +0200
Subject: [PATCH] for slice, always initialize empty invhess unless a new chain
 is initialized from a mode file

---
 matlab/check_posterior_sampler_options.m | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/matlab/check_posterior_sampler_options.m b/matlab/check_posterior_sampler_options.m
index 5f0cd7b7d3..acaff8da94 100644
--- a/matlab/check_posterior_sampler_options.m
+++ b/matlab/check_posterior_sampler_options.m
@@ -73,6 +73,10 @@ if strcmp(options_.posterior_sampling_method,'slice')
    
     if options_.load_mh_file,
         posterior_sampler_options.slice_initialize_with_mode = 0;
+    else
+        if ~posterior_sampler_options.slice_initialize_with_mode,
+            posterior_sampler_options.invhess=[];
+        end
     end    
 
     if posterior_sampler_options.rotated,
-- 
GitLab