impl_reparam_mesh.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_REPARAM_MESH_HPP
#define QUGAR_IMPL_REPARAM_MESH_HPP
 
 
#include <qugar/bbox.hpp>
#include <qugar/cart_grid_tp.hpp>
#include <qugar/domain_function.hpp>
#include <qugar/point.hpp>
#include <qugar/reparam_mesh.hpp>
#include <qugar/tensor_index_tp.hpp>
#include <qugar/tolerance.hpp>
 
#include <cstddef>
#include <functional>
#include <span>
#include <string>
#include <vector>
 
 
namespace qugar::impl {
 
template<int dim, int range> class ImplReparamMesh : public ReparamMesh<dim, range>
{
public:
  explicit ImplReparamMesh(int order);
 
 
  void generate_wirebasket(const std::vector<std::reference_wrapper<const ImplicitFunc<range>>> &impl_funcs,
    const BoundBox<range> &domain,
    const Tolerance &tol);
 
  void generate_wirebasket(const std::vector<std::reference_wrapper<const ImplicitFunc<range>>> &impl_funcs,
    const std::vector<int> &cell_ids,
    const BoundBox<range> &domain,
    const Tolerance &tol);
 
private:
  void generate_wirebasket(const std::vector<int> &cell_ids,
    const BoundBox<range> &domain,
    const Tolerance &tol,
    const std::function<bool(const Point<range> &)> &check_face,
    const std::function<bool(const Point<range> &)> &check_internal);
 
public:
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == dim)
  void orient_cells_positively();
 
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == dim)
  void orient_cells_positively(const std::vector<int> &cell_ids);
 
private:
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == (dim + 1))
  void orient_levelset_cells_positively(const std::vector<int> &cell_ids,
    const std::function<Point<range>(const Point<range> &)> &outer_normal_computer);
 
public:
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == (dim + 1))
  void orient_levelset_cells_positively(
    const std::vector<std::reference_wrapper<const ImplicitFunc<range>>> &impl_funcs);
 
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == (dim + 1))
  void orient_levelset_cells_positively(
    const std::vector<std::reference_wrapper<const ImplicitFunc<range>>> &impl_funcs,
    const std::vector<int> &cell_ids);
 
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == (dim + 1))
  void orient_facet_cells_positively(int local_facet_id);
 
  template<int dim_aux = dim>
    requires(dim_aux == dim && range == (dim + 1))
  void orient_facet_cells_positively(int local_facet_id, const std::vector<int> &cell_ids);
};
 
 
}// namespace qugar::impl
 
#endif// QUGAR_IMPL_REPARAM_MESH_HPP