additional_data.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.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/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