vector_tools_constraints.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2020 - 2024 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_constraints_h
#define dealii_vector_tools_constraints_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/dofs/dof_handler.h>
 
#include <deal.II/multigrid/mg_constrained_dofs.h>
 
#include <map>
#include <set>
 
DEAL_II_NAMESPACE_OPEN
 
template <typename number>
class AffineConstraints;
template <int dim, int spacedim>
struct StaticMappingQ1;
template <int dim, typename Number>
class Function;
template <int dim, int spacedim>
class Mapping;
 
 
namespace VectorTools
{
  template <int dim, int spacedim>
  void
  compute_nonzero_normal_flux_constraints(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    const std::map<types::boundary_id, const Function<spacedim, double> *>
                                 &function_map,
    AffineConstraints<double>    &constraints,
    const Mapping<dim, spacedim> &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const bool use_manifold_for_normal = true);
 
  template <int dim, int spacedim>
  void
  compute_nonzero_normal_flux_constraints_on_level(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    const std::map<types::boundary_id, const Function<spacedim, double> *>
                                 &function_map,
    AffineConstraints<double>    &constraints,
    const Mapping<dim, spacedim> &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const IndexSet    &refinement_edge_indices = IndexSet(),
    const unsigned int level                   = numbers::invalid_unsigned_int,
    const bool         use_manifold_for_normal = true);
 
  template <int dim, int spacedim, typename number>
  void
  compute_no_normal_flux_constraints(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    AffineConstraints<number>          &constraints,
    const Mapping<dim, spacedim>       &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const bool use_manifold_for_normal = true);
 
  template <int dim, int spacedim>
  void
  compute_no_normal_flux_constraints_on_level(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    AffineConstraints<double>          &constraints,
    const Mapping<dim, spacedim>       &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const IndexSet    &refinement_edge_indices = IndexSet(),
    const unsigned int level                   = numbers::invalid_unsigned_int,
    const bool         use_manifold_for_normal = true);
 
  template <int dim, int spacedim>
  void
  compute_nonzero_tangential_flux_constraints(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    const std::map<types::boundary_id, const Function<spacedim, double> *>
                                 &function_map,
    AffineConstraints<double>    &constraints,
    const Mapping<dim, spacedim> &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const bool use_manifold_for_normal = true);
 
  template <int dim, int spacedim>
  void
  compute_normal_flux_constraints(
    const DoFHandler<dim, spacedim>    &dof_handler,
    const unsigned int                  first_vector_component,
    const std::set<types::boundary_id> &boundary_ids,
    AffineConstraints<double>          &constraints,
    const Mapping<dim, spacedim>       &mapping =
      (ReferenceCells::get_hypercube<dim>()
#ifndef _MSC_VER
         .template get_default_linear_mapping<dim, spacedim>()
#else
         .ReferenceCell::get_default_linear_mapping<dim, spacedim>()
#endif
         ),
    const bool use_manifold_for_normal = true);
 
} // namespace VectorTools
 
DEAL_II_NAMESPACE_CLOSE
 
#endif // dealii_vector_tools_constraints_h