From 525d2d6e835092ddd6a62d7b3a09589d66082e15 Mon Sep 17 00:00:00 2001
From: Hugo TRACHINO <hugo.trachino@pc127-05.pedago.ensiie.fr>
Date: Fri, 4 Oct 2019 14:12:13 +0200
Subject: [PATCH] ./follow lance

---
 sgi2018/src/follower/transform.h | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/sgi2018/src/follower/transform.h b/sgi2018/src/follower/transform.h
index d7a8f5a..f2409f2 100644
--- a/sgi2018/src/follower/transform.h
+++ b/sgi2018/src/follower/transform.h
@@ -17,7 +17,7 @@ namespace Transform
         // to generate a point from a parameter i, you can use  :
         ComplexMatrix res(xiVector.size(), startIndexes.size());
         for(unsigned int i = 0; i < xiVector.size(); i ++){
-          Point endPoint(startPoint.omega, startPoint.size, xiVector[startIndexes[i]]);
+          Point endPoint(startPoint.omega, startPoint.size, xiVector[i]);
           std::vector<unsigned int> indexesResult = Follow::adaptativStep(startPoint, endPoint, startIndexes, distanceThreshold, minstepsize);
           ComplexVector v(indexesResult.size());
           for(unsigned int j = 0; j < indexesResult.size(); j ++){
@@ -41,8 +41,9 @@ namespace Transform
         // NOTE: identical to the previous function but dumps all eigenvalues produced, not just targets, to outCsv
         //
         ComplexMatrix res(xiVector.size(), startIndexes.size());
+        #pragma omp parallel for
         for(unsigned int i = 0; i < xiVector.size(); i ++){
-          Point endPoint(startPoint.omega, startPoint.size, xiVector[startIndexes[i]]);
+          Point endPoint(startPoint.omega, startPoint.size, xiVector[i]);
           std::vector<unsigned int> indexesResult = Follow::adaptativStep(startPoint, endPoint, startIndexes, distanceThreshold, minstepsize);
           ComplexVector v(indexesResult.size());
           outCsv.write( (ComplexVector) phaseSpeed(startPoint.omega, startPoint.formatedEigenvalues) );
-- 
GitLab