From d2de20ef215b386c7a4dbedc437a884595a909a3 Mon Sep 17 00:00:00 2001
From: astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>
Date: Thu, 6 Jan 2011 10:35:24 +0000
Subject: Parallelized for the powell method

git-svn-id: file:///home/svn/opt_utilities@159 ed2142bd-67ad-457f-ba7c-d818d4011675
---
 methods/linmin/linmin.hpp        | 2 ++
 methods/powell/powell_method.hpp | 3 +++
 2 files changed, 5 insertions(+)

(limited to 'methods')

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];
-- 
cgit v1.2.2