[Asm] CompSparsPat now also computes NodeAdjTable

AssemblerLib/ComputeSparsityPattern.cpp

@@ -7,18 +7,22 @@
  *
  */
 
+#include "MeshLib/NodeAdjacencyTable.h"
+#include "LocalToGlobalIndexMap.h"
+
 #include "ComputeSparsityPattern.h"
 
 namespace AssemblerLib
 {
 
 SparsityPattern
-computeSparsityPattern(
-        MeshLib::NodeAdjacencyTable const& node_adjacency_table,
-        LocalToGlobalIndexMap const& dof_table,
+computeSparsityPattern(LocalToGlobalIndexMap const& dof_table,
         MeshLib::Mesh const& mesh
         )
 {
+    MeshLib::NodeAdjacencyTable node_adjacency_table;
+    node_adjacency_table.createTable(mesh.getNodes());
+
     std::vector<std::vector<GlobalIndexType> > global_idcs;
     global_idcs.reserve(mesh.getNNodes());
     for (std::size_t n=0; n<mesh.getNNodes(); ++n)

AssemblerLib/ComputeSparsityPattern.h

@@ -12,19 +12,18 @@
 
 #include <vector>
 
-#include "MeshLib/NodeAdjacencyTable.h"
-#include "LocalToGlobalIndexMap.h"
 #include "ProcessLib/NumericsConfig.h"
 
 namespace AssemblerLib
 {
 
+class LocalToGlobalIndexMap;
+
 using SparsityPattern = std::vector<GlobalIndexType>;
 
 /**
  * @brief Computes a sparsity pattern for the given inputs.
  *
- * @param node_adjacency_table describes the mesh topology
  * @param dof_table            maps mesh nodes to global indices
  * @param mesh                 mesh for which the two parameters above are defined
  *
@@ -32,7 +31,6 @@ using SparsityPattern = std::vector<GlobalIndexType>;
  */
 SparsityPattern
 computeSparsityPattern(
-        MeshLib::NodeAdjacencyTable const& node_adjacency_table,
         LocalToGlobalIndexMap const& dof_table,
         MeshLib::Mesh const& mesh
         );