psblas_sparsity_pattern.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2019 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_psblas_sparsity_h
#define dealii_psblas_sparsity_h
 
#include <deal.II/lac/sparsity_pattern_base.h>
 
#ifdef DEAL_II_WITH_PSBLAS
#  include <deal.II/lac/psblas_common.h>
#endif
 
DEAL_II_NAMESPACE_OPEN
 
#ifdef DEAL_II_WITH_PSBLAS
namespace PSCToolkitWrappers
{
 
  class SparsityPattern : public SparsityPatternBase
  {
  public:
    using size_type = ::types::global_dof_index;
 
    SparsityPattern() = default;
 
    SparsityPattern(const IndexSet &parallel_partitioning,
                    const MPI_Comm  communicator = MPI_COMM_WORLD);
 
    virtual ~SparsityPattern() override = default;
 
    template <typename ForwardIterator>
    void
    add_entries(const size_type row,
                ForwardIterator begin,
                ForwardIterator end,
                const bool      indices_are_sorted = false);
 
    virtual void
    add_row_entries(const size_type                  &row,
                    const ArrayView<const size_type> &columns,
                    const bool indices_are_sorted = false) override;
 
    void
    add(const size_type i, const size_type j);
 
    using SparsityPatternBase::add_entries;
 
    void
    compress();
 
  private:
    std::shared_ptr<psb_c_descriptor> psblas_descriptor;
 
    friend class SparseMatrix;
  };
 
 
} // namespace PSCToolkitWrappers
 
#endif // DEAL_II_WITH_PSBLAS
 
DEAL_II_NAMESPACE_CLOSE
 
#endif