diff options
| -rw-r--r-- | methods/linmin/linmin.hpp | 2 | ||||
| -rw-r--r-- | methods/powell/powell_method.hpp | 3 |
2 files changed, 5 insertions, 0 deletions
diff --git a/methods/linmin/linmin.hpp b/methods/linmin/linmin.hpp index ddf2a8e..9ec1dd8 100644 --- a/methods/linmin/linmin.hpp +++ b/methods/linmin/linmin.hpp @@ -53,6 +53,7 @@ namespace opt_utilities pT xt; opt_eq(xt,p1); +#pragma omp parallel for for(size_t i=0;i<get_size(xt);++i) { //get_element(xt,i)+=x*get_element((pT)xi1,i); @@ -88,6 +89,7 @@ namespace opt_utilities //cout<<xx<<endl; fret=brent(ax,xx,bx,fadpt,TOL,xmin); //cout<<xmin<<endl; +#pragma omp parallel for for(j=0;j<n;++j) { //get_element(xi,j)*=xmin; diff --git a/methods/powell/powell_method.hpp b/methods/powell/powell_method.hpp index 041fa93..294e101 100644 --- a/methods/powell/powell_method.hpp +++ b/methods/powell/powell_method.hpp @@ -116,6 +116,7 @@ namespace opt_utilities del=0.0; for(i=0;i<n;++i) { +#pragma omp parallel for for(j=0;j<n;++j) { //get_element(xit,j)=xi[j][i]; @@ -138,6 +139,7 @@ namespace opt_utilities std::cerr<<"powell exceeding maximun iterations."<<std::endl; return; } +#pragma omp parallel for for(j=0;j<n;++j) { //get_element(ptt,j)=T(2.)*get_element(p,j)-get_element(pt,j); @@ -155,6 +157,7 @@ namespace opt_utilities if(t<T(0.)) { linmin(p,xit,fret,*p_fo); +#pragma omp parallel for for(j=0;j<n;++j) { xi[j][ibig-1]=xi[j][n-1]; |