From b4bfdd83e5fdf310e54af8e8127dc30ace50a804 Mon Sep 17 00:00:00 2001
From: StefanoPetrilli <stefanop_1999@hotmail.it>
Date: Sat, 21 May 2022 10:22:58 -0500
Subject: [PATCH] computeCandidateMate and isAlreadyMatched

---
 amgprec/impl/aggregator/MatchBoxPC.h          | 19 +++++
 ...mEdgesLinearSearchMesgBndlSmallMateCMP.cpp | 73 ++++++++++++++-----
 2 files changed, 74 insertions(+), 18 deletions(-)

diff --git a/amgprec/impl/aggregator/MatchBoxPC.h b/amgprec/impl/aggregator/MatchBoxPC.h
index 94ea7ea8..73908b9b 100644
--- a/amgprec/impl/aggregator/MatchBoxPC.h
+++ b/amgprec/impl/aggregator/MatchBoxPC.h
@@ -157,6 +157,25 @@ inline MilanLongInt firstComputeCandidateMate(MilanLongInt adj1,
                                               MilanLongInt* verLocInd,
                                               MilanReal* edgeLocWeight);
 
+inline bool isAlreadyMatched(MilanLongInt k,
+                                MilanLongInt* verLocInd,
+                                MilanLongInt StartIndex,
+                                MilanLongInt EndIndex,
+                                vector <MilanLongInt> &GMate,
+                                MilanLongInt* Mate,
+                                map <MilanLongInt, MilanLongInt> &Ghost2LocalMap);
+
+inline MilanLongInt computeCandidateMate(MilanLongInt adj1,
+                                         MilanLongInt adj2,
+                                         MilanReal* edgeLocWeight,
+                                         MilanLongInt k,
+                                         MilanLongInt* verLocInd,
+                                         MilanLongInt StartIndex,
+                                         MilanLongInt EndIndex,
+                                         vector <MilanLongInt> &GMate,
+                                         MilanLongInt* Mate,
+                                         map <MilanLongInt, MilanLongInt> &Ghost2LocalMap);
+
 void dalgoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP
         (
                 MilanLongInt NLVer, MilanLongInt NLEdge,
diff --git a/amgprec/impl/aggregator/algoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP.cpp b/amgprec/impl/aggregator/algoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP.cpp
index 84d05e08..e73c7939 100644
--- a/amgprec/impl/aggregator/algoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP.cpp
+++ b/amgprec/impl/aggregator/algoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP.cpp
@@ -443,13 +443,7 @@ void dalgoDistEdgeApproxDomEdgesLinearSearchMesgBndlSmallMateCMP(
                 w = -1;
                 heaviestEdgeWt = MilanRealMin; //Assign the smallest Value possible first LDBL_MIN
                 for (k = adj1; k < adj2; k++) {
-                    if ((verLocInd[k] < StartIndex) || (verLocInd[k] > EndIndex)) { //Is it a ghost vertex?
-                        if (GMate[Ghost2LocalMap[verLocInd[k]]] >= 0)// Already matched
-                            continue;
-                    } else { //A local vertex
-                        if (Mate[verLocInd[k] - StartIndex] >= 0) // Already matched
-                            continue;
-                    }
+                    if (isAlreadyMatched(k, verLocInd, StartIndex, EndIndex, GMate, Mate, Ghost2LocalMap)) continue;
 
                     if ((edgeLocWeight[k] > heaviestEdgeWt) ||
                         ((edgeLocWeight[k] == heaviestEdgeWt) && (w < verLocInd[k]))) {
@@ -1533,22 +1527,66 @@ inline MilanLongInt firstComputeCandidateMate(MilanLongInt adj1,
     return w;
 }
 
-/*
-inline MilanLongInt firstComputeCandidateMate(MilanLongInt adj1,
+/**
+ * //TODO documentation
+ * @param k
+ * @param verLocInd
+ * @param StartIndex
+ * @param EndIndex
+ * @param GMate
+ * @param Mate
+ * @param Ghost2LocalMap
+ * @return
+ */
+inline bool isAlreadyMatched(MilanLongInt k,
+                                MilanLongInt* verLocInd,
+                                MilanLongInt StartIndex,
+                                MilanLongInt EndIndex,
+                                vector <MilanLongInt> &GMate,
+                                MilanLongInt* Mate,
+                                map <MilanLongInt, MilanLongInt> &Ghost2LocalMap
+                                ) {
+
+    if ((verLocInd[k] < StartIndex) || (verLocInd[k] > EndIndex)) { //Is it a ghost vertex?
+        if (GMate[Ghost2LocalMap[verLocInd[k]]] >= 0)// Already matched
+            return true;
+    } else { //A local vertex
+        if (Mate[verLocInd[k] - StartIndex] >= 0) // Already matched
+            return true;
+    }
+
+    return false;
+}
+
+/**
+ * //TODO documentation
+ * @param adj1
+ * @param adj2
+ * @param edgeLocWeight
+ * @param k
+ * @param verLocInd
+ * @param StartIndex
+ * @param EndIndex
+ * @param GMate
+ * @param Mate
+ * @param Ghost2LocalMap
+ * @return
+ */
+inline MilanLongInt computeCandidateMate(MilanLongInt adj1,
                                               MilanLongInt adj2,
+                                              MilanReal* edgeLocWeight,
+                                              MilanLongInt k,
                                               MilanLongInt* verLocInd,
-                                              MilanReal* edgeLocWeight)
+                                              MilanLongInt StartIndex,
+                                              MilanLongInt EndIndex,
+                                              vector <MilanLongInt> &GMate,
+                                              MilanLongInt* Mate,
+                                              map <MilanLongInt, MilanLongInt> &Ghost2LocalMap)
 {
     MilanInt w = -1;
     MilanReal heaviestEdgeWt = MilanRealMin; //Assign the smallest Value possible first LDBL_MIN
     for (k = adj1; k < adj2; k++) {
-        if ((verLocInd[k] < StartIndex) || (verLocInd[k] > EndIndex)) { //Is it a ghost vertex?
-            if (GMate[Ghost2LocalMap[verLocInd[k]]] >= 0)// Already matched
-                continue;
-        } else { //A local vertex
-            if (Mate[verLocInd[k] - StartIndex] >= 0) // Already matched
-                continue;
-        }
+        if (isAlreadyMatched(k, verLocInd, StartIndex, EndIndex, GMate, Mate, Ghost2LocalMap)) continue;
 
         if ((edgeLocWeight[k] > heaviestEdgeWt) ||
             ((edgeLocWeight[k] == heaviestEdgeWt) && (w < verLocInd[k]))) {
@@ -1558,7 +1596,6 @@ inline MilanLongInt firstComputeCandidateMate(MilanLongInt adj1,
     } //End of for loop
     return w;
 }
- */
 #endif
 
 #endif
\ No newline at end of file