matrixbuilder.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_BACKENDS_MATRIXBUILDER_HH
#define DUNE_FUFEM_BACKENDS_MATRIXBUILDER_HH
 
#ifndef DUNE_FUFEM_DISABLE_HEADER_DEPRECATION
#warning This header is deprecated and will be removed after 2.11. Use the ISTLMatrixBackend::PatternBuilder from dune-assembler instead.
#endif
 
#include <array>
#include <vector>
 
#include <dune/common/indices.hh>
#include <dune/common/hybridutilities.hh>
 
#include <dune/istl/matrix.hh>
#include <dune/istl/bcrsmatrix.hh>
#include <dune/istl/matrixindexset.hh>
#include <dune/istl/multitypeblockvector.hh>
#include <dune/istl/multitypeblockmatrix.hh>
 
namespace Dune::Fufem {
 
 
 
template<class Matrix>
class
MatrixBuilder;
 
 
 
template<class T, class A>
class
[[deprecated("This class is deprecated and will be removed after 2.11. Use the ISTLMatrixBackend::PatternBuilder from dune-assembler instead.")]]
MatrixBuilder<Dune::BCRSMatrix<T, A>>
{
public:
 
  using Matrix = Dune::BCRSMatrix<T, A>;
 
  MatrixBuilder(Matrix& matrix) :
    matrix_(matrix)
  {}
 
  template<class RowSizeInfo, class ColSizeInfo>
  void resize(const RowSizeInfo& rowSizeInfo, const ColSizeInfo& colSizeInfo)
  {
    indices_.resize(rowSizeInfo.size(), colSizeInfo.size());
  }
 
  template<class RowIndex, class ColIndex>
  void insertEntry(const RowIndex& rowIndex, const ColIndex& colIndex)
  {
    indices_.add(rowIndex[0], colIndex[0]);
  }
 
  void setupMatrix()
  {
    indices_.exportIdx(matrix_);
  }
 
private:
  Dune::MatrixIndexSet indices_;
  Matrix& matrix_;
};
 
 
 
template<class Row0, class... Rows>
class
[[deprecated("This class is deprecated and will be removed after 2.11. Use the ISTLMatrixBackend::PatternBuilder from dune-assembler instead.")]]
MatrixBuilder<MultiTypeBlockMatrix<Row0,Rows...>>
{
public:
 
  using Matrix = MultiTypeBlockMatrix<Row0,Rows...>;
 
  MatrixBuilder(Matrix& matrix) :
    matrix_(matrix)
  {}
 
  template<class RowSizeInfo, class ColSizeInfo>
  void resize(const RowSizeInfo& rowSizeInfo, const ColSizeInfo& colSizeInfo)
  {
    for(size_t i=0; i<rowBlocks_; i++)
      for(size_t j=0; j<colBlocks_; j++)
        indices_[i][j].resize(rowSizeInfo.size({i}), colSizeInfo.size({j}));
  }
 
  template<class RowIndex, class ColIndex>
  void insertEntry(const RowIndex& rowIndex, const ColIndex& colIndex)
  {
    indices_[rowIndex[0]][colIndex[0]].add(rowIndex[1], colIndex[1]);
  }
 
  void setupMatrix()
  {
    Hybrid::forEach(Hybrid::integralRange(Dune::index_constant<rowBlocks_>()), [&](auto&& i) {
      Hybrid::forEach(Hybrid::integralRange(Dune::index_constant<colBlocks_>()), [&](auto&& j) {
        indices_[i][j].exportIdx(matrix_[i][j]);
      });
    });
 
  }
 
private:
  static constexpr size_t rowBlocks_ = Matrix::N();
  static constexpr size_t colBlocks_ = Matrix::M();
  std::array<std::array<Dune::MatrixIndexSet,colBlocks_>,rowBlocks_> indices_;
  Matrix& matrix_;
};
 
 
template<class T, class A>
class
[[deprecated("This class is deprecated and will be removed after 2.11. Use the ISTLMatrixBackend::PatternBuilder from dune-assembler instead.")]]
MatrixBuilder<Dune::Matrix<Dune::BCRSMatrix<T, A>>>
{
public:
 
  using Matrix = Dune::Matrix<Dune::BCRSMatrix<T, A>>;
 
  MatrixBuilder(Matrix& matrix) :
    matrix_(matrix)
  {}
 
  template<class RowSizeInfo, class ColSizeInfo>
  void resize(const RowSizeInfo& rowSizeInfo, const ColSizeInfo& colSizeInfo)
  {
    auto rowBlocks = rowSizeInfo.size();
    auto colBlocks = colSizeInfo.size();
    indices_.resize(rowBlocks);
    for(size_t i=0; i<rowBlocks; i++)
    {
      indices_[i].resize(colBlocks);
      for(size_t j=0; j<colBlocks; j++)
        indices_[i][j].resize(rowSizeInfo.size({i}), colSizeInfo.size({j}));
    }
  }
 
  template<class RowIndex, class ColIndex>
  void insertEntry(const RowIndex& rowIndex, const ColIndex& colIndex)
  {
    indices_[rowIndex[0]][colIndex[0]].add(rowIndex[1], colIndex[1]);
  }
 
  void setupMatrix()
  {
    auto rowBlocks = indices_.size();
    auto colBlocks = indices_[0].size();
    matrix_.setSize(rowBlocks, colBlocks);
    for(size_t i=0; i<rowBlocks; i++)
      for(size_t j=0; j<colBlocks; j++)
        indices_[i][j].exportIdx(matrix_[i][j]);
  }
 
private:
  std::vector<std::vector<Dune::MatrixIndexSet>> indices_;
  Matrix& matrix_;
};
 
 
 
} // namespace Dune::Fufem
 
 
 
#endif // DUNE_FUFEM_BACKENDS_MATRIXBUILDER_HH