Commit 39432dca authored by Plaszczynski Stephane's avatar Plaszczynski Stephane

use 2.38^2/dim as default for scale

parent 32c7dfca
...@@ -50,7 +50,7 @@ int main(int argc, char *argv[]) ...@@ -50,7 +50,7 @@ int main(int argc, char *argv[])
} }
string ar_file=string("ar_vs_length_")+fileName; string ar_file=string("ar_vs_length_")+fileName;
string scale_file=string("scale_")+fileName; string scale_file=string("scale_")+fileName;
string cor_file=string("prop_cov_")+fileName; string cor_file=string("last_step_")+fileName;
cout << "output files:" <<endl; cout << "output files:" <<endl;
cout << ar_file << " " << scale_file << " " << cor_file << " " << endl; cout << ar_file << " " << scale_file << " " << cor_file << " " << endl;
...@@ -63,7 +63,7 @@ int main(int argc, char *argv[]) ...@@ -63,7 +63,7 @@ int main(int argc, char *argv[])
const string proposal_cov=parser.params.find<string>("proposal_cov",""); const string proposal_cov=parser.params.find<string>("proposal_cov","");
const int t0=parser.params.find<int>("t0",2000); const int t0=parser.params.find<int>("t0",2000);
const int ts=parser.params.find<int>("ts",10000); const int ts=parser.params.find<int>("ts",10000);
double scale=parser.params.find<double>("scale",0.001); double scale=parser.params.find<double>("scale",-1);
const bool do_move=parser.params.find<bool>("do_move",false); const bool do_move=parser.params.find<bool>("do_move",false);
bool verbose=parser.params.find<bool>("verbose",false); bool verbose=parser.params.find<bool>("verbose",false);
...@@ -83,6 +83,11 @@ int main(int argc, char *argv[]) ...@@ -83,6 +83,11 @@ int main(int argc, char *argv[])
cout <<"dim=" <<dim <<endl; cout <<"dim=" <<dim <<endl;
//SP 25/06/19: compute scale if undefined depending on dim
if (scale<0) {
scale=2.38*2.38/dim;
cout << "refining scale to " << scale << endl;
}
//--------defining seeds----------------------------------------------- //--------defining seeds-----------------------------------------------
unsigned int seed(0); unsigned int seed(0);
if (parser.params.param_present("seed")){ if (parser.params.param_present("seed")){
......
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