#1443 Add some comments while looking a bit how to make FromVertexMatching...

#1443 Add some comments while looking a bit how to make FromVertexMatching work for a run from prepartitioned data (this will encompass storing directly the matrix behin the operator).

Sources/FiniteElement/FiniteElementSpace/Internal/DofProgramWiseIndexListPerVertexCoordIndexList.cpp

@@ -33,9 +33,9 @@ namespace MoReFEM
         
             DofProgramWiseIndexListPerVertexCoordIndexList
             ::DofProgramWiseIndexListPerVertexCoordIndexList(const unsigned int unique_id,
-                                                         const FEltSpace& felt_space,
-                                                         const NumberingSubset& numbering_subset)
-            :  unique_id_parent(unique_id),
+                                                             const FEltSpace& felt_space,
+                                                             const NumberingSubset& numbering_subset)
+            : unique_id_parent(unique_id),
             felt_space_(felt_space),
             numbering_subset_(numbering_subset)
             {
@@ -105,7 +105,6 @@ namespace MoReFEM
                                 {
                                     assert(!(!dof_ptr));
                                     const auto& dof = *dof_ptr;
-                                    std::cout << "DofProgramWiseIndexListPerVertexCoordIndexList" << std::endl;
                                     index_list.push_back(dof.GetProgramWiseIndex(numbering_subset));
                                     
                                     ++Ndof;

Sources/FiniteElement/FiniteElementSpace/Internal/DofProgramWiseIndexListPerVertexCoordIndexList.hpp

@@ -68,10 +68,6 @@ namespace MoReFEM
             // ============================
 
 
-
-
-
-
             /*!
              * \brief An internal class used to init correctly a VertexMatching interpolator.
              *
@@ -88,6 +84,10 @@ namespace MoReFEM
              *
              * \attention This class is intended at a very specific point of GodOfDof initialization (namely in GodOfDof::Init2(),
              * after some data have been initialized but before any processor-wise reduction occurred).
+             *
+             * \attention Due to the remarks above, this type of object is not adapted at all when operating a run from prepartitioned data (in which case
+             * by design we never get access to some of the data required here). The idea is to store directly the interpolation matrix in \a VertexMatching -
+             * which is a much more economical way to store the required data.
              */
             class DofProgramWiseIndexListPerVertexCoordIndexList
             : public ::MoReFEM::Crtp::UniqueId<DofProgramWiseIndexListPerVertexCoordIndexList, UniqueIdNS::AssignationMode::manual>

Sources/FiniteElement/FiniteElementSpace/Internal/DofProgramWiseIndexListPerVertexCoordIndexListManager.cpp

@@ -45,7 +45,7 @@ namespace MoReFEM
                      const unsigned int felt_space_index,
                      const unsigned int numbering_subset_index)
             {
-                // This method will be called from all GodOfDofs, but only the once that encompass felt_space
+                // This method will be called from all GodOfDofs, but only the one that encompass felt_space
                 // is allowed to proceed.
                 if (!god_of_dof.IsFEltSpace(felt_space_index))
                     return;

Sources/Geometry/Interpolator/VertexMatching.cpp

@@ -39,9 +39,6 @@ namespace MoReFEM
                 FilesystemNS::File::Read(in, filename, __FILE__, __LINE__);
                 
                 std::string line;
-
-                std::cout << "VM FILE = " << filename << std::endl;
-                
                 
                 while (getline(in ,line))
                 {

Sources/OperatorInstances/NonConformInterpolator/FromVertexMatching/FromVertexMatching.hpp

@@ -187,7 +187,7 @@ namespace MoReFEM
             /*!
              * \brief Pattern of the matrix.
              *
-             * \note It is stored only
+             * \note It is stored only if constructor requires it explicitly - if not it is reset to nullptr once used.
              */
             Wrappers::Petsc::MatrixPattern::const_unique_ptr matrix_pattern_ = nullptr;