diff --git a/sgi2018/src/follower/transform.h b/sgi2018/src/follower/transform.h index d7a8f5a81568391efd000022ed57bdfcde00aad2..f2409f24d8c82d6f9029f5dd803b358785142bba 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) );