sparsity_pattern_base.cc

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------
 
#include <deal.II/lac/sparsity_pattern_base.h>
 
#include <boost/container/small_vector.hpp>
 
#include <algorithm>
#include <utility>
 
DEAL_II_NAMESPACE_OPEN
 
void
SparsityPatternBase::add_entries(
  const ArrayView<const std::pair<size_type, size_type>> &inputs)
{
  // We always want to sort so that we can use the optimized add_row_entries()
  // function
  boost::container::small_vector<std::pair<size_type, size_type>, 128> entries(
    inputs.begin(), inputs.end());
  std::sort(entries.begin(), entries.end());
  boost::container::small_vector<size_type, 128> columns;
 
  auto entry = entries.begin();
  while (entry != entries.end())
    {
      const auto row     = entry->first;
      auto       row_end = entry;
      while (row_end != entries.end() && row_end->first == row)
        ++row_end;
 
      columns.resize(0);
      columns.reserve(row_end - entry);
      // Simultaneously transform and uniquify
      columns.push_back(entry->second);
      ++entry;
      while (entry != row_end)
        {
          if (columns.back() != entry->second)
            columns.push_back(entry->second);
          ++entry;
        }
 
      add_row_entries(row,
                      make_array_view(columns.begin(), columns.end()),
                      true);
    }
}
 
DEAL_II_NAMESPACE_CLOSE