Allow to perform neighborhood_width sampling around calibrated values;

matlab/dynare_sensitivity.m

@@ -48,6 +48,8 @@ if isfield(options_gsa,'graph_format'),
 end
 if isfield(options_gsa,'mode_file'),
     options_.mode_file=options_gsa.mode_file;
+elseif isfield(options_gsa,'neighborhood_width') && options_gsa_.neighborhood_width>0,
+    options_.mode_file='';    
 end
 
 options_.order = 1;

matlab/gsa/stab_map_.m

@@ -241,7 +241,11 @@ if fload==0,
         %             end
         %         end
         %load([fname_,'_mode'])
+        if neighborhood_width>0 && isempty(options_.mode_file),
+            xparam1 = get_all_parameters(estim_params_,M_);
+        else
         eval(['load ' options_.mode_file '.mat;']);
+        end
         if neighborhood_width>0,
             for j=1:nshock,
                 lpmat0(:,j) = randperm(Nsam)'./(Nsam+1); %latin hypercube