added calls to re-initialize persistent variable 'penalty' is

dsge_likelihood() and dsge_likelihood_hh() before calling the
optimizer

matlab/dynare_estimation_1.m

@@ -140,6 +140,8 @@ end
 
 %% Estimation of the posterior mode or likelihood mode
 if ~isequal(options_.mode_compute,0) && ~options_.mh_posterior_mode_estimation
+    % initializing persistent variable 'penalty' in dsge_likelihood()
+    dsge_likelihood();
     switch options_.mode_compute
       case 1
         optim_options = optimset('display','iter','LargeScale','off', ...

matlab/newrat.m

@@ -55,6 +55,8 @@ htol0=htol_base;
 gibbstol=length(BayesInfo.pshape)/50; %25;
 
 func_hh = str2func([func2str(func0),'_hh']);
+% initialize persistent variable 'penalty' in func_hh()
+func_hh();
 fval0=feval(func0,x,DynareDataset,DynareOptions,Model,EstimatedParameters,BayesInfo,DynareResults);
 fval=fval0;