New samplers and mode compute methods

To test, I put a call in dynare_estimation1.m

Herbst_Schorfheide_sampler(TargetFun,xparam1,mh_bounds,dataset_,dataset_info,options_,M_,estim_params_,bayestopt_,oo_)
or
DSMH_sampler(TargetFun,xparam1,mh_bounds,dataset_,dataset_info,options_,M_,estim_params_,bayestopt_,oo_)

a priori everything is defined in Global_initialization.m with comprehensive notations (DSMH_ for DSMH) and HB_ for Herbst &Schorfheide