impl_unfitted_domain.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_IMPL_UNFITTED_DOMAIN_HPP
#define QUGAR_IMPL_UNFITTED_DOMAIN_HPP
 
 
 
#include <qugar/bbox.hpp>
#include <qugar/cart_grid_tp.hpp>
#include <qugar/domain_function.hpp>
#include <qugar/unfitted_domain.hpp>
#include <qugar/unfitted_domain_binary_part.hpp>
 
#include <array>
#include <cassert>
#include <cstddef>
#include <functional>
#include <memory>
#include <optional>
#include <vector>
 
namespace qugar::impl {
 
 
template<int dim> class UnfittedImplDomain : public UnfittedDomain<dim>
{
public:
  using FuncSign = qugar::impl::FuncSign;
  using FacetsStatus = std::array<ImmersedFacetStatus, static_cast<std::size_t>(dim) * 2>;
  using GridPtr = std::shared_ptr<const CartGridTP<dim>>;
  using FuncPtr = std::shared_ptr<const ImplicitFunc<dim>>;
  using KDTree = UnfittedBinarySpacePart<dim>;
  using KDTreePtr = std::shared_ptr<UnfittedBinarySpacePart<dim>>;
 
  explicit UnfittedImplDomain(const FuncPtr phi, GridPtr grid);
 
  explicit UnfittedImplDomain(const FuncPtr phi, GridPtr grid, const std::vector<std::int64_t> &cells);
 
  [[nodiscard]] FuncPtr get_impl_func() const;
 
private:
  FuncPtr phi_;
 
  void create_decomposition(KDTree &tree,
    const std::function<FuncSign(const BoundBox<dim> &)> &func_sign,
    const std::optional<std::vector<std::int64_t>> &target_cells);
 
  void classify_undetermined_sign_cell(KDTree &tree);
};
 
 
}// namespace qugar::impl
 
 
#endif// QUGAR_IMPL_UNFITTED_DOMAIN_HPP