diff --git a/matlab/dynare_estimation_1.m b/matlab/dynare_estimation_1.m
index a0eb52924b210df5d63c48acbc458d8d9affa865..4ce3ff8aa502cdaae927c26b6ed36987c071de13 100644
--- a/matlab/dynare_estimation_1.m
+++ b/matlab/dynare_estimation_1.m
@@ -420,6 +420,9 @@ switch options_.MCMC_jumping_covariance
     end
 end
 
+bounds = prior_bounds(bayestopt_, options_.prior_trunc); %reset bounds as lb and ub must only be operational during mode-finding
+DSMH_sampler(objective_function,xparam1,bounds,dataset_,dataset_info,options_,M_,estim_params_,bayestopt_,oo_);
+
 if (any(bayestopt_.pshape  >0 ) && options_.mh_replic) || ...
         (any(bayestopt_.pshape >0 ) && options_.load_mh_file)  %% not ML estimation
     bounds = prior_bounds(bayestopt_, options_.prior_trunc); %reset bounds as lb and ub must only be operational during mode-finding
diff --git a/matlab/particles b/matlab/particles
index 84d213ea3e705f36871bc8f960d175a330fd71de..857168ddf8c2d386dac875da1ede31fd689be309 160000
--- a/matlab/particles
+++ b/matlab/particles
@@ -1 +1 @@
-Subproject commit 84d213ea3e705f36871bc8f960d175a330fd71de
+Subproject commit 857168ddf8c2d386dac875da1ede31fd689be309