vector_tools_integrate_difference.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2020 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_vector_tools_integrate_difference_h
#define dealii_vector_tools_integrate_difference_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/numerics/vector_tools_common.h>
 
DEAL_II_NAMESPACE_OPEN
 
template <int dim, int spacedim>
DEAL_II_CXX20_REQUIRES((concepts::is_valid_dim_spacedim<dim, spacedim>))
class DoFHandler;
 
template <int dim, typename Number>
class Function;
template <int dim, int spacedim>
class Mapping;
template <int dim>
class Quadrature;
 
template <int dim, int spacedim>
DEAL_II_CXX20_REQUIRES((concepts::is_valid_dim_spacedim<dim, spacedim>))
class Triangulation;
 
namespace hp
{
  template <int dim, int spacedim>
  class MappingCollection;
  template <int dim>
  class QCollection;
} // namespace hp
 
 
namespace VectorTools
{
  template <int dim, typename Number, class OutVector, int spacedim>
  DEAL_II_CXX20_REQUIRES(concepts::is_writable_dealii_vector_type<OutVector>)
  void integrate_difference(const Mapping<dim, spacedim>     &mapping,
                            const DoFHandler<dim, spacedim>  &dof,
                            const ReadVector<Number>         &fe_function,
                            const Function<spacedim, Number> &exact_solution,
                            OutVector                        &difference,
                            const Quadrature<dim>            &q,
                            const NormType                   &norm,
                            const Function<spacedim, double> *weight = nullptr,
                            const double                      exponent = 2.);
 
  template <int dim, typename Number, class OutVector, int spacedim>
  DEAL_II_CXX20_REQUIRES(concepts::is_writable_dealii_vector_type<OutVector>)
  void integrate_difference(const DoFHandler<dim, spacedim>  &dof,
                            const ReadVector<Number>         &fe_function,
                            const Function<spacedim, Number> &exact_solution,
                            OutVector                        &difference,
                            const Quadrature<dim>            &q,
                            const NormType                   &norm,
                            const Function<spacedim, double> *weight = nullptr,
                            const double                      exponent = 2.);
 
  template <int dim, typename Number, class OutVector, int spacedim>
  DEAL_II_CXX20_REQUIRES(concepts::is_writable_dealii_vector_type<OutVector>)
  void integrate_difference(const hp::MappingCollection<dim, spacedim> &mapping,
                            const DoFHandler<dim, spacedim>            &dof,
                            const ReadVector<Number>         &fe_function,
                            const Function<spacedim, Number> &exact_solution,
                            OutVector                        &difference,
                            const hp::QCollection<dim>       &q,
                            const NormType                   &norm,
                            const Function<spacedim, double> *weight = nullptr,
                            const double                      exponent = 2.);
 
  template <int dim, typename Number, class OutVector, int spacedim>
  DEAL_II_CXX20_REQUIRES(concepts::is_writable_dealii_vector_type<OutVector>)
  void integrate_difference(const DoFHandler<dim, spacedim>  &dof,
                            const ReadVector<Number>         &fe_function,
                            const Function<spacedim, Number> &exact_solution,
                            OutVector                        &difference,
                            const hp::QCollection<dim>       &q,
                            const NormType                   &norm,
                            const Function<spacedim, double> *weight = nullptr,
                            const double                      exponent = 2.);
 
  template <int dim, int spacedim, class InVector>
  DEAL_II_CXX20_REQUIRES(concepts::is_dealii_vector_type<InVector>)
  double compute_global_error(const Triangulation<dim, spacedim> &tria,
                              const InVector &cellwise_error,
                              const NormType &norm,
                              const double    exponent = 2.);
 
} // namespace VectorTools
 
DEAL_II_NAMESPACE_CLOSE
 
#endif // dealii_vector_tools_integrate_difference_h