additional_data.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2022 - 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_cgal_additional_data_h
#define dealii_cgal_additional_data_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/exceptions.h>
 
#ifdef DEAL_II_WITH_CGAL
 
#  include <CGAL/version.h>
#  if CGAL_VERSION_MAJOR >= 6
#    include <CGAL/Installation/internal/disable_deprecation_warnings_and_errors.h>
#  endif
#  include <CGAL/Mesh_facet_topology.h>
 
#  include <limits>
 
DEAL_II_NAMESPACE_OPEN
namespace CGALWrappers
{
  enum class FacetTopology
  {
    facet_vertices_on_surface = CGAL::FACET_VERTICES_ON_SURFACE,

    facet_vertices_on_same_surface_patch =
      CGAL::FACET_VERTICES_ON_SAME_SURFACE_PATCH,

    facet_vertices_on_same_surface_patch_with_adjacency_check =
      CGAL::FACET_VERTICES_ON_SAME_SURFACE_PATCH_WITH_ADJACENCY_CHECK,
  };

  template <int dim>
  struct AdditionalData
  {
    double edge_size;

    double facet_angle;

    double facet_size;
    double facet_distance;

    FacetTopology facet_topology;

    double cell_radius_edge_ratio;
    double cell_size;

    AdditionalData(
      double        edge_s  = std::numeric_limits<double>::max(),
      double        facet_a = 0.,
      double        facet_s = 0.,
      double        facet_d = 0.,
      FacetTopology facet_t =
        ::CGALWrappers::FacetTopology::facet_vertices_on_surface,
      double cell_radius_edge_r = 0.,
      double cell_s             = 0.)
    {
      AssertThrow(
        dim == 3,
        ExcMessage(
          "These struct can be instantiated with 3d Triangulations only."));
      edge_size              = edge_s;
      facet_angle            = facet_a;
      facet_size             = facet_s;
      facet_distance         = facet_d;
      facet_topology         = facet_t;
      cell_radius_edge_ratio = cell_radius_edge_r;
      cell_size              = cell_s;
    }
  };

  template <>
  struct AdditionalData<2>
  {
    double angular_bound;

    double radius_bound;

    double distance_bound;

    AdditionalData(double angular_b  = 0.,
                   double radius_b   = 0.,
                   double distance_b = 0.)
    {
      angular_bound  = angular_b;
      radius_bound   = radius_b;
      distance_bound = distance_b;
    }
  };
 
} // namespace CGALWrappers
DEAL_II_NAMESPACE_CLOSE
 
#else
 
DEAL_II_NAMESPACE_OPEN
namespace CGALWrappers
{
  template <int dim>
  struct AdditionalData
  {};
} // namespace CGALWrappers
 
DEAL_II_NAMESPACE_CLOSE
#endif
 
#endif