coupling.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2018 - 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_non_matching_coupling
#define dealii_non_matching_coupling
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/function_lib.h>
#include <deal.II/base/quadrature.h>
 
#include <deal.II/dofs/dof_handler.h>
 
#include <deal.II/fe/component_mask.h>
#include <deal.II/fe/fe.h>
#include <deal.II/fe/fe_values.h>
#include <deal.II/fe/mapping_q1.h>
 
#include <deal.II/grid/grid_tools_cache.h>
 
#include <deal.II/lac/affine_constraints.h>
 
DEAL_II_NAMESPACE_OPEN
 
namespace NonMatching
{
  template <int dim0, int dim1, int spacedim, typename number = double>
  void
  create_coupling_sparsity_pattern(
    const DoFHandler<dim0, spacedim> &space_dh,
    const DoFHandler<dim1, spacedim> &immersed_dh,
    const Quadrature<dim1>           &quad,
    SparsityPatternBase              &sparsity,
    const AffineConstraints<number>  &constraints    = {},
    const ComponentMask              &space_comps    = {},
    const ComponentMask              &immersed_comps = {},
    const Mapping<dim0, spacedim>    &space_mapping =
      StaticMappingQ1<dim0, spacedim>::mapping,
    const Mapping<dim1, spacedim> &immersed_mapping =
      StaticMappingQ1<dim1, spacedim>::mapping,
    const AffineConstraints<number> &immersed_constraints =
      AffineConstraints<number>());
 
  template <int dim0, int dim1, int spacedim, typename number = double>
  void
  create_coupling_sparsity_pattern(
    const GridTools::Cache<dim0, spacedim> &cache,
    const DoFHandler<dim0, spacedim>       &space_dh,
    const DoFHandler<dim1, spacedim>       &immersed_dh,
    const Quadrature<dim1>                 &quad,
    SparsityPatternBase                    &sparsity,
    const AffineConstraints<number>        &constraints    = {},
    const ComponentMask                    &space_comps    = {},
    const ComponentMask                    &immersed_comps = {},
    const Mapping<dim1, spacedim>          &immersed_mapping =
      StaticMappingQ1<dim1, spacedim>::mapping,
    const AffineConstraints<number> &immersed_constraints =
      AffineConstraints<number>());
 
 
  template <int dim0, int dim1, int spacedim, typename Matrix>
  void
  create_coupling_mass_matrix(
    const DoFHandler<dim0, spacedim>                     &space_dh,
    const DoFHandler<dim1, spacedim>                     &immersed_dh,
    const Quadrature<dim1>                               &quad,
    Matrix                                               &matrix,
    const AffineConstraints<typename Matrix::value_type> &constraints =
      AffineConstraints<typename Matrix::value_type>(),
    const ComponentMask           &space_comps    = {},
    const ComponentMask           &immersed_comps = {},
    const Mapping<dim0, spacedim> &space_mapping =
      StaticMappingQ1<dim0, spacedim>::mapping,
    const Mapping<dim1, spacedim> &immersed_mapping =
      StaticMappingQ1<dim1, spacedim>::mapping,
    const AffineConstraints<typename Matrix::value_type> &immersed_constraints =
      AffineConstraints<typename Matrix::value_type>());
 
  template <int dim0, int dim1, int spacedim, typename Matrix>
  void
  create_coupling_mass_matrix(
    const GridTools::Cache<dim0, spacedim>               &cache,
    const DoFHandler<dim0, spacedim>                     &space_dh,
    const DoFHandler<dim1, spacedim>                     &immersed_dh,
    const Quadrature<dim1>                               &quad,
    Matrix                                               &matrix,
    const AffineConstraints<typename Matrix::value_type> &constraints =
      AffineConstraints<typename Matrix::value_type>(),
    const ComponentMask           &space_comps    = {},
    const ComponentMask           &immersed_comps = {},
    const Mapping<dim1, spacedim> &immersed_mapping =
      StaticMappingQ1<dim1, spacedim>::mapping,
    const AffineConstraints<typename Matrix::value_type> &immersed_constraints =
      AffineConstraints<typename Matrix::value_type>());
 
  template <int dim0, int dim1, int spacedim, typename Number = double>
  void
  create_coupling_sparsity_pattern(
    const double                           &epsilon,
    const GridTools::Cache<dim0, spacedim> &cache0,
    const GridTools::Cache<dim1, spacedim> &cache1,
    const DoFHandler<dim0, spacedim>       &dh0,
    const DoFHandler<dim1, spacedim>       &dh1,
    const Quadrature<dim1>                 &quad,
    SparsityPatternBase                    &sparsity,
    const AffineConstraints<Number> &constraints0 = AffineConstraints<Number>(),
    const ComponentMask             &comps0       = {},
    const ComponentMask             &comps1       = {});
 
  template <int dim0, int dim1, int spacedim, typename Matrix>
  void
  create_coupling_mass_matrix(
    Functions::CutOffFunctionBase<spacedim>              &kernel,
    const double                                         &epsilon,
    const GridTools::Cache<dim0, spacedim>               &cache0,
    const GridTools::Cache<dim1, spacedim>               &cache1,
    const DoFHandler<dim0, spacedim>                     &dh0,
    const DoFHandler<dim1, spacedim>                     &dh1,
    const Quadrature<dim0>                               &quadrature0,
    const Quadrature<dim1>                               &quadrature1,
    Matrix                                               &matrix,
    const AffineConstraints<typename Matrix::value_type> &constraints0 =
      AffineConstraints<typename Matrix::value_type>(),
    const ComponentMask &comps0 = {},
    const ComponentMask &comps1 = {});
} // namespace NonMatching
DEAL_II_NAMESPACE_CLOSE
 
#endif