Skip to content
Snippets Groups Projects
Commit 46e715ba authored by Mark Jenkinson's avatar Mark Jenkinson
Browse files

Put in std::sqrt

parent 57c2bf36
No related branches found
No related tags found
No related merge requests found
...@@ -373,7 +373,7 @@ void scg(ColumnVector& x,const gEvalFunction& func){ ...@@ -373,7 +373,7 @@ void scg(ColumnVector& x,const gEvalFunction& func){
if(kappa<eps){ if(kappa<eps){
break; break;
} }
sigma=sigma0/sqrt(kappa); sigma=sigma0/std::sqrt(kappa);
xplus = x + sigma*d; xplus = x + sigma*d;
gplus=func.g_evaluate(xplus);gevals++; gplus=func.g_evaluate(xplus);gevals++;
gamma = (d.t()*(gplus - gradnew)).AsScalar()/sigma; gamma = (d.t()*(gplus - gradnew)).AsScalar()/sigma;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment