./follow lance

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) );