From a82b89966ac518ad15f0f5a74de7aa034add53c9 Mon Sep 17 00:00:00 2001
From: kradchen <chenjuijun@realcan.cn>
Date: Tue, 30 May 2023 17:15:14 +0800
Subject: [PATCH] Improve filterTransmission performance.

---
 src/transmissionReconstruction/dataFilter/dataFilter.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/transmissionReconstruction/dataFilter/dataFilter.cpp b/src/transmissionReconstruction/dataFilter/dataFilter.cpp
index 6cefb4d..a05b982 100644
--- a/src/transmissionReconstruction/dataFilter/dataFilter.cpp
+++ b/src/transmissionReconstruction/dataFilter/dataFilter.cpp
@@ -36,6 +36,7 @@ namespace Recon {
     size.getData()[3] = aMSensData.size();
     auto idx = sub2ind(size, {aVrnBlock, aVrlBlock, aVsnBlock, aVslBlock});
     auto transData = zeros(idx.getDataSize(), 1, 1);
+    #pragma omp parallel for
     for (size_t i = 0; i < idx.getDataSize(); i++) {
         auto index = (size_t)(idx.getData()[i] - 1);
         auto sliceIndex = index / aMSensData0.getDataSize();
@@ -57,6 +58,7 @@ namespace Recon {
         auto inbetweenAngle = acosd(dot(aMSenderNormalBlock, aMReceiverNormalBlock) /
                                     (vecnorm(aMSenderNormalBlock, Norm2, 1) *
                                     vecnorm(aMReceiverNormalBlock, Norm2, 1)));
+        #pragma omp parallel for
         for (size_t i = 0; i < transData.getDataSize(); i++)
         {
             transData.getData()[i] = (inbetweenAngle.getData()[i]<aAngleLowerLimit || inbetweenAngle.getData()[i] > aAngleUpperLimit)?0.0:1.0;