matrix_lib.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2002 - 2018 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.
//
// ---------------------------------------------------------------------
 
#ifndef dealii_matrix_lib_h
#define dealii_matrix_lib_h
 
#include <deal.II/base/subscriptor.h>
 
#include <deal.II/lac/solver_richardson.h>
#include <deal.II/lac/vector_memory.h>
 
DEAL_II_NAMESPACE_OPEN
 
template <typename number>
class Vector;
template <typename number>
class BlockVector;
template <typename number>
class SparseMatrix;
 
class DEAL_II_DEPRECATED MeanValueFilter : public Subscriptor
{
public:
  using size_type = types::global_dof_index;
 
  MeanValueFilter(const size_type component = numbers::invalid_size_type);
 
  template <typename number>
  void
  filter(Vector<number> &v) const;
 
  template <typename number>
  void
  filter(BlockVector<number> &v) const;
 
  template <typename number>
  void
  vmult(Vector<number> &dst, const Vector<number> &src) const;
 
  template <typename number>
  void
  vmult_add(Vector<number> &dst, const Vector<number> &src) const;
 
  template <typename number>
  void
  vmult(BlockVector<number> &dst, const BlockVector<number> &src) const;
 
  template <typename number>
  void
  vmult_add(BlockVector<number> &dst, const BlockVector<number> &src) const;
 
 
  template <typename VectorType>
  void
  Tvmult(VectorType &, const VectorType &) const;
 
  template <typename VectorType>
  void
  Tvmult_add(VectorType &, const VectorType &) const;
 
private:
  const size_type component;
};
 
 
 
//---------------------------------------------------------------------------
 
 
template <typename VectorType>
inline void
MeanValueFilter::Tvmult(VectorType &, const VectorType &) const
{
  Assert(false, ExcNotImplemented());
}
 
 
template <typename VectorType>
inline void
MeanValueFilter::Tvmult_add(VectorType &, const VectorType &) const
{
  Assert(false, ExcNotImplemented());
}
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif