Commit d6993cec authored by Marco Ratto's avatar Marco Ratto
Browse files

1) trap the case where the model does not solve for any parameter values in the MC sample;

2) fixed initialization of alpha2_stab;
parent 22e0b9cc
...@@ -128,7 +128,7 @@ options_gsa = set_default_option(options_gsa,'Nsam',2048); ...@@ -128,7 +128,7 @@ options_gsa = set_default_option(options_gsa,'Nsam',2048);
options_gsa = set_default_option(options_gsa,'load_redform',0); options_gsa = set_default_option(options_gsa,'load_redform',0);
options_gsa = set_default_option(options_gsa,'load_rmse',0); options_gsa = set_default_option(options_gsa,'load_rmse',0);
options_gsa = set_default_option(options_gsa,'load_stab',0); options_gsa = set_default_option(options_gsa,'load_stab',0);
options_gsa = set_default_option(options_gsa,'alpha2_stab',0.3); options_gsa = set_default_option(options_gsa,'alpha2_stab',0);
options_gsa = set_default_option(options_gsa,'ksstat',0.1); options_gsa = set_default_option(options_gsa,'ksstat',0.1);
options_gsa = set_default_option(options_gsa,'pvalue_ks',0.001); options_gsa = set_default_option(options_gsa,'pvalue_ks',0.001);
options_gsa = set_default_option(options_gsa,'pvalue_corr',0.001); options_gsa = set_default_option(options_gsa,'pvalue_corr',0.001);
...@@ -218,6 +218,11 @@ end ...@@ -218,6 +218,11 @@ end
if options_gsa.stab && ~options_gsa.ppost, if options_gsa.stab && ~options_gsa.ppost,
x0 = stab_map_(OutputDirectoryName,options_gsa); x0 = stab_map_(OutputDirectoryName,options_gsa);
if isempty(x0),
disp(' ')
disp('Sensitivity computations stopped: no parameter set provided a unique solution')
return
end
end end
% reduced form % reduced form
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment