immersed_surface_quadrature.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2017 - 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_immersed_surface_quadrature
#define dealii_non_matching_immersed_surface_quadrature
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/point.h>
#include <deal.II/base/quadrature.h>
#include <deal.II/base/subscriptor.h>
#include <deal.II/base/tensor.h>
 
#include <vector>
 
DEAL_II_NAMESPACE_OPEN
namespace NonMatching
{
  template <int dim, int spacedim = dim>
  class ImmersedSurfaceQuadrature : public Quadrature<dim>
  {
  public:
    ImmersedSurfaceQuadrature() = default;
 
    ImmersedSurfaceQuadrature(const std::vector<Point<dim>>          &points,
                              const std::vector<double>              &weights,
                              const std::vector<Tensor<1, spacedim>> &normals);
 
    void
    clear();
 
    void
    push_back(const Point<dim>          &point,
              const double               weight,
              const Tensor<1, spacedim> &normal);
 
    const Tensor<1, spacedim> &
    normal_vector(const unsigned int i) const;
 
    const std::vector<Tensor<1, spacedim>> &
    get_normal_vectors() const;
 
  protected:
    std::vector<Tensor<1, spacedim>> normals;
  };
 
} // namespace NonMatching
DEAL_II_NAMESPACE_CLOSE
 
#endif