coupling.h

// ---------------------------------------------------------------------
//
// Copyright (C) 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 at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#ifndef dealii_non_matching_coupling
#define dealii_non_matching_coupling

#include <deal.II/base/config.h>
#include <deal.II/base/quadrature.h>

#include <deal.II/grid/grid_tools_cache.h>

#include <deal.II/fe/component_mask.h>
#include <deal.II/fe/mapping_q1.h>
#include <deal.II/fe/fe.h>

#include <deal.II/dofs/dof_handler.h>

#include <deal.II/lac/constraint_matrix.h>

DEAL_II_NAMESPACE_OPEN

namespace NonMatching
{
  template<int dim0, int dim1, int spacedim, typename Sparsity>
  void create_coupling_sparsity_pattern(const DoFHandler<dim0, spacedim> &space_dh,
                                        const DoFHandler<dim1, spacedim> &immersed_dh,
                                        const Quadrature<dim1>           &quad,
                                        Sparsity                         &sparsity,
                                        const ConstraintMatrix           &constraints = ConstraintMatrix(),
                                        const ComponentMask              &space_comps = ComponentMask(),
                                        const ComponentMask              &immersed_comps = ComponentMask(),
                                        const Mapping<dim0, spacedim>    &space_mapping = StaticMappingQ1<dim0,spacedim>::mapping,
                                        const Mapping<dim1, spacedim>    &immersed_mapping = StaticMappingQ1<dim1, spacedim>::mapping);

  template<int dim0, int dim1, int spacedim, typename Sparsity>
  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,
                                        Sparsity                         &sparsity,
                                        const ConstraintMatrix           &constraints = ConstraintMatrix(),
                                        const ComponentMask              &space_comps = ComponentMask(),
                                        const ComponentMask              &immersed_comps = ComponentMask(),
                                        const Mapping<dim1, spacedim>    &immersed_mapping = StaticMappingQ1<dim1, spacedim>::mapping);


  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 ConstraintMatrix           &constraints = ConstraintMatrix(),
                                   const ComponentMask              &space_comps = ComponentMask(),
                                   const ComponentMask              &immersed_comps = ComponentMask(),
                                   const Mapping<dim0, spacedim>    &space_mapping = StaticMappingQ1<dim0,spacedim>::mapping,
                                   const Mapping<dim1, spacedim>    &immersed_mapping = StaticMappingQ1<dim1, spacedim>::mapping);

  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 ConstraintMatrix                 &constraints = ConstraintMatrix(),
                                   const ComponentMask                    &space_comps = ComponentMask(),
                                   const ComponentMask                    &immersed_comps = ComponentMask(),
                                   const Mapping<dim1, spacedim>          &immersed_mapping = StaticMappingQ1<dim1, spacedim>::mapping);
}
DEAL_II_NAMESPACE_CLOSE

#endif