--- /dev/null
+## ---------------------------------------------------------------------
+##
+## Copyright (C) 2017 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.
+##
+## ---------------------------------------------------------------------
+
+#
+# Configuration for the CGAL library:
+#
+
+CONFIGURE_FEATURE(CGAL)
--- /dev/null
+## ---------------------------------------------------------------------
+##
+## Copyright (C) 2017 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.
+##
+## ---------------------------------------------------------------------
+
+#
+# Try to find the CGAL libraries
+#
+# This module exports
+#
+# CGAL_INCLUDE_DIRS
+#
+
+SET(CGAL_DIR "" CACHE PATH "An optional hint to a CGAL installation")
+SET_IF_EMPTY(CGAL_DIR "$ENV{CGAL_DIR}")
+
+IF(NOT "${CGAL_DIR}" STREQUAL "")
+ SET(CGAL_DIR ${CGAL_DIR})
+ENDIF()
+
+IF(DEAL_II_HAVE_CXX17)
+ # temporarily disable ${CMAKE_SOURCE_DIR}/cmake/modules for module lookup
+ LIST(REMOVE_ITEM CMAKE_MODULE_PATH ${CMAKE_SOURCE_DIR}/cmake/modules/)
+ FIND_PACKAGE(CGAL)
+ LIST(APPEND CMAKE_MODULE_PATH ${CMAKE_SOURCE_DIR}/cmake/modules/)
+ELSE(DEAL_II_HAVE_CXX17)
+ SET(CGAL_FOUND FALSE)
+ SET(CGAL_INCLUDE_DIRS "-NOTFOUND")
+ MESSAGE("-- WARNING: CGAL wrappers require cxx17. Disabling CGAL Support.")
+ENDIF(DEAL_II_HAVE_CXX17)
+
+DEAL_II_PACKAGE_HANDLE(CGAL
+ INCLUDE_DIRS REQUIRED CGAL_INCLUDE_DIRS
+ CLEAR CGAL_INCLUDE_DIRS
+ )
DEAL_II_ARPACK_WITH_PARPACK=1 \
DEAL_II_WITH_ASSIMP=1 \
DEAL_II_WITH_BOOST=1 \
+ DEAL_II_WITH_CGAL=1 \
DEAL_II_WITH_TASKFLOW=1 \
DEAL_II_WITH_COMPLEX_VALUES=1 \
DEAL_II_WITH_CUDA=1 \
--- /dev/null
+New: Added support for the CGAL library (www.cgal.org).
+<br>
+(Marco Feder, Luca Heltai, 2022/04/26)
#cmakedefine DEAL_II_WITH_ARPACK
#cmakedefine DEAL_II_WITH_ARBORX
#cmakedefine DEAL_II_WITH_ASSIMP
+#cmakedefine DEAL_II_WITH_CGAL
#cmakedefine DEAL_II_WITH_COMPLEX_VALUES
#cmakedefine DEAL_II_WITH_CUDA
#cmakedefine DEAL_II_WITH_GINKGO
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2020 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_utilities_h
+#define dealii_cgal_utilities_h
+
+#include <deal.II/base/config.h>
+
+#ifdef DEAL_II_WITH_CGAL
+# include <CGAL/Cartesian.h>
+# include <CGAL/Simple_cartesian.h>
+
+DEAL_II_NAMESPACE_OPEN
+/**
+ * Interface to the Computational Geometry Algorithm Library (CGAL).
+ *
+ * CGAL is a software project that provides easy access to efficient and
+ * reliable geometric algorithms. The library offers data structures and
+ * algorithms like triangulations, Voronoi diagrams, Boolean operations on
+ * polygons and polyhedra, point set processing, arrangements of curves, surface
+ * and volume mesh generation, geometry processing, alpha shapes, convex hull
+ * algorithms, shape reconstruction, AABB and KD trees...
+ *
+ * You can learn more about the CGAL library at https://www.cgal.org/
+ */
+namespace CGALWrappers
+{
+ /**
+ * Convert from deal.II Point to any compatible CGAL point.
+ *
+ * @tparam CGALPointType Any of the CGAL point types
+ * @tparam dim Dimension of the point
+ * @param [in] p An input deal.II Point<dim>
+ * @return CGALPointType A CGAL point
+ */
+ template <typename CGALPointType, int dim>
+ inline CGALPointType
+ to_cgal(const dealii::Point<dim> &p);
+
+ /**
+ * Convert from various CGAL point types to deal.II Point.
+ * @tparam dim Dimension of the point
+ * @tparam CGALPointType Any of the CGAL point types
+ * @param p
+ * @return dealii::Point<dim>
+ */
+ template <int dim, typename CGALPointType>
+ inline dealii::Point<dim>
+ to_dealii(const CGALPointType &p);
+} // namespace CGALWrappers
+
+# ifndef DOXYGEN
+// Template implementations
+namespace CGALWrappers
+{
+ template <typename CGALPointType, int dim>
+ inline CGALPointType
+ to_cgal(const dealii::Point<dim> &p)
+ {
+ constexpr int cdim = CGALPointType::Ambient_dimension::value;
+ static_assert(dim <= cdim, "Only dim <= cdim supported");
+ if constexpr (cdim == 1)
+ return CGALPointType(p[0]);
+ else if constexpr (cdim == 2)
+ return CGALPointType(p[0], dim > 1 ? p[1] : 0);
+ else if constexpr (cdim == 3)
+ return CGALPointType(p[0], dim > 1 ? p[1] : 0, dim > 2 ? p[2] : 0);
+ else
+ Assert(false, dealii::ExcNotImplemented());
+ return CGALPointType();
+ }
+
+
+
+ template <int dim, typename CGALPointType>
+ inline dealii::Point<dim>
+ to_dealii(const CGALPointType &p)
+ {
+ constexpr int cdim = CGALPointType::Ambient_dimension::value;
+ if constexpr (dim == 1)
+ return dealii::Point<dim>(CGAL::to_double(p.x()));
+ else if constexpr (dim == 2)
+ return dealii::Point<dim>(CGAL::to_double(p.x()),
+ cdim > 1 ? CGAL::to_double(p.y()) : 0);
+ else if constexpr (dim == 3)
+ return dealii::Point<dim>(CGAL::to_double(p.x()),
+ cdim > 1 ? CGAL::to_double(p.y()) : 0,
+ cdim > 2 ? CGAL::to_double(p.z()) : 0);
+ else
+ Assert(false, dealii::ExcNotImplemented());
+ }
+} // namespace CGALWrappers
+# endif
+
+DEAL_II_NAMESPACE_CLOSE
+
+#endif
+#endif
--- /dev/null
+CMAKE_MINIMUM_REQUIRED(VERSION 3.1.0)
+INCLUDE(../setup_testsubproject.cmake)
+PROJECT(testsuite CXX)
+IF(DEAL_II_WITH_CGAL)
+ DEAL_II_PICKUP_TESTS()
+ENDIF()
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2020 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.
+//
+// ---------------------------------------------------------------------
+
+// Test conversion from CGAL to deal.II Point and viceversa
+
+#include <deal.II/base/point.h>
+
+#include <CGAL/IO/io.h>
+#include <CGAL/Simple_cartesian.h>
+#include <deal.II/cgal/utilities.h>
+
+#include "../tests.h"
+
+using namespace CGALWrappers;
+
+int
+main()
+{
+ initlog();
+ using CGALPoint = CGAL::Point_3<CGAL::Simple_cartesian<double>>;
+ // Test conversion from deal.II Point to CGAL Point
+ {
+ const Point<3> p(1.0, 2.0, 3.0);
+ const auto cgal_point = to_cgal<CGALPoint>(p);
+ deallog << "CGAL Point: " << cgal_point << std::endl;
+ }
+
+ // Test conversion from CGAL Point to deal.II Point
+ {
+ const CGALPoint cgal_point(1.0, 2.0, 3.0);
+ const auto deal_ii_point = to_dealii<3>(cgal_point);
+ deallog << "deal.II Point: " << deal_ii_point << std::endl;
+ }
+}
--- /dev/null
+
+DEAL::CGAL Point: 1.00000 2.00000 3.00000
+DEAL::deal.II Point: 1.00000 2.00000 3.00000
projects / dealii.git / commitdiff