Skip to content
Snippets Groups Projects
Commit 3ca9485e authored by Marco Ratto's avatar Marco Ratto
Browse files

bug fix when calling gsa

set default order=1
set analytic derivation =1
parent 83afb256
No related branches found
No related tags found
No related merge requests found
...@@ -61,7 +61,7 @@ options_ident = set_default_option(options_ident,'advanced',0); ...@@ -61,7 +61,7 @@ options_ident = set_default_option(options_ident,'advanced',0);
options_ident = set_default_option(options_ident,'normalize_jacobians',1); options_ident = set_default_option(options_ident,'normalize_jacobians',1);
options_ident = set_default_option(options_ident,'lik_init',1); options_ident = set_default_option(options_ident,'lik_init',1);
if options_ident.gsa_sample_file, if options_ident.gsa_sample_file,
GSAFolder = checkpath('gsa'); GSAFolder = checkpath('gsa',M_.dname);
if options_ident.gsa_sample_file==1, if options_ident.gsa_sample_file==1,
load([GSAFolder,filesep,fname_,'_prior'],'lpmat','lpmat0','istable'); load([GSAFolder,filesep,fname_,'_prior'],'lpmat','lpmat0','istable');
elseif options_ident.gsa_sample_file==2, elseif options_ident.gsa_sample_file==2,
...@@ -105,11 +105,13 @@ nlags = options_ident.ar; ...@@ -105,11 +105,13 @@ nlags = options_ident.ar;
periods = options_ident.periods; periods = options_ident.periods;
replic = options_ident.replic; replic = options_ident.replic;
useautocorr = options_ident.useautocorr; useautocorr = options_ident.useautocorr;
options_.order=1;
options_.ar=nlags; options_.ar=nlags;
options_.prior_mc = options_ident.prior_mc; options_.prior_mc = options_ident.prior_mc;
options_.options_ident = options_ident; options_.options_ident = options_ident;
options_.Schur_vec_tol = 1.e-8; options_.Schur_vec_tol = 1.e-8;
options_.nomoments=0; options_.nomoments=0;
options_.analytic_derivation=1;
options_ = set_default_option(options_,'datafile',[]); options_ = set_default_option(options_,'datafile',[]);
options_.mode_compute = 0; options_.mode_compute = 0;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment