Adding the ranks topology arg to SynGlobalDofs

src/field/DistributedField.cpp

@@ -391,7 +391,7 @@ namespace gmshfem::field
   }
 
   template< class T_Scalar, field::Form T_Form >
-  void DistributedField< T_Scalar, T_Form >::syncGlobalDofs(std::vector< unsigned long long > &localToGlobal, std::vector< unsigned long long > &readIDs)
+  void DistributedField< T_Scalar, T_Form >::syncGlobalDofs(std::vector< unsigned long long > &localToGlobal, std::vector< unsigned long long > &readIDs, const std::optional<std::vector<int>>&  neighboringRanks)
   {
   #ifdef HAVE_MPI
 

src/field/DistributedField.h

@@ -61,7 +61,7 @@ namespace gmshfem::field
     virtual void preProMPI(const std::optional<std::vector<int>>&  neighboringRanks = std::nullopt) override;
     virtual std::unordered_set<dofs::Dof*> getAllOwnedDofs() override;
     virtual std::unordered_set<dofs::Dof*> getNonOwnedDofs() override;
-    virtual void syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs) override;
+    virtual void syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs, const std::optional<std::vector<int>>&  neighboringRanks = std::nullopt) override;
     void dumpOwnedKeys();
 
    public:

src/field/FieldInterface.cpp

@@ -169,7 +169,7 @@ namespace gmshfem::field
   }
 
   template< class T_Scalar >
-  void FieldInterface< T_Scalar >::syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs)
+  void FieldInterface< T_Scalar >::syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs, const std::optional<std::vector<int>>&  neighboringRanks)
   {
     // Do nothing
   }

src/field/FieldInterface.h

@@ -244,7 +244,7 @@ namespace gmshfem::field
     virtual void preProMPI(const std::optional<std::vector<int>>&  neighboringRanks = std::nullopt);
     virtual std::unordered_set<dofs::Dof*> getAllOwnedDofs();
     virtual std::unordered_set<dofs::Dof*> getNonOwnedDofs();
-    virtual void syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs);
+    virtual void syncGlobalDofs(std::vector<unsigned long long>& localToGlobal, std::vector<unsigned long long>& readIDs, const std::optional<std::vector<int>>&  neighboringRanks = std::nullopt);
 
 
    public:

src/problem/Formulation.cpp

@@ -1242,7 +1242,7 @@ namespace gmshfem::problem
     timeSyncGlobal.tick();
     /* A dof knows its gmsh ID and its global index -> let the field handle the sharing */
     for (auto [tag, fieldInterfacePtr] : _unknownFields) {
-      fieldInterfacePtr->syncGlobalDofs(localToGlobalMap, localIdOfNonOwned);
+      fieldInterfacePtr->syncGlobalDofs(localToGlobalMap, localIdOfNonOwned, neighboringRanks);
     }
     common::GmshFem::BarrierMPI();