cut_quadrature.hpp

Go to the documentation of this file.

// --------------------------------------------------------------------------
//
// Copyright (C) 2025-present by Pablo Antolin
//
// This file is part of the QUGaR library.
//
// SPDX-License-Identifier:    MIT
//
// --------------------------------------------------------------------------
 
#ifndef QUGAR_LIBRARY_CUT_QUADRATURE_HPP
#define QUGAR_LIBRARY_CUT_QUADRATURE_HPP
 
 
 
#include <qugar/point.hpp>
#include <qugar/types.hpp>
#include <qugar/unfitted_domain.hpp>
 
#include <cstddef>
#include <cstdint>
#include <memory>
#include <vector>
 
 
namespace qugar {
 
enum ImmersedStatus : std::uint8_t { cut, full, empty };
 
template<int dim> struct CutCellsQuad
{
  std::vector<std::int64_t> cells;
  std::vector<int> n_pts_per_cell;
 
  std::vector<Point<dim>> points;
  std::vector<real> weights;
 
  void reserve(const int n_cells, const int n_tot_pts);
};
 
template<int dim> struct CutIsoBoundsQuad
{
  std::vector<std::int64_t> cells;
  std::vector<int> local_facet_ids;
  std::vector<int> n_pts_per_facet;
 
  std::vector<Point<dim>> points;
  std::vector<real> weights;
 
  void reserve(const int n_cells, const int n_tot_pts);
};
 
template<int dim> struct CutUnfBoundsQuad
{
  std::vector<std::int64_t> cells;
  std::vector<int> n_pts_per_cell;
 
  std::vector<Point<dim>> points;
  std::vector<real> weights;
  std::vector<Point<dim>> normals;
 
  void reserve(const int n_cells, const int n_tot_pts);
};
 
template<int dim>
std::shared_ptr<const CutCellsQuad<dim>>
  create_quadrature(const UnfittedDomain<dim> &unf_domain, const std::vector<std::int64_t> &cells, int n_pts_dir);
 
template<int dim>
std::shared_ptr<const CutUnfBoundsQuad<dim>> create_unfitted_bound_quadrature(const UnfittedDomain<dim> &unf_domain,
  const std::vector<std::int64_t> &cells,
  int n_pts_dir,
  bool include_facet_unf_bdry,
  bool exclude_ext_bdry);
 
template<int dim>
std::shared_ptr<const CutIsoBoundsQuad<dim - 1>> create_facets_quadrature_interior_integral(
  const UnfittedDomain<dim> &unf_domain,
  const std::vector<std::int64_t> &cells,
  const std::vector<int> &facets,
  int n_pts_dir);
 
template<int dim>
std::shared_ptr<const CutIsoBoundsQuad<dim - 1>> create_facets_quadrature_exterior_integral(
  const UnfittedDomain<dim> &unf_domain,
  const std::vector<std::int64_t> &cells,
  const std::vector<int> &facets,
  int n_pts_dir);
 
 
}// namespace qugar
 
#endif// QUGAR_LIBRARY_CUT_QUADRATURE_HPP