swh:1:snp:78e145aa8174e576786284475a76cf6f187b3475
Tip revision: ff74e7f6e89013c6b08d86369d90f7c35c330a2b authored by Laurent Rineau on 31 July 2020, 13:52:12 UTC
Merge branch 'releases/CGAL-4.14-branch' into releases/CGAL-5.0-branch
Merge branch 'releases/CGAL-4.14-branch' into releases/CGAL-5.0-branch
Tip revision: ff74e7f
Skin_surface_traits_3.h
// Copyright (c) 2005 Rijksuniversiteit Groningen (Netherlands)
// All rights reserved.
//
// This file is part of CGAL (www.cgal.org).
//
// $URL$
// $Id$
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
//
//
// Author(s) : Nico Kruithof <Nico@cs.rug.nl>
#ifndef CGAL_SKIN_SURFACE_TRAITS_3_H
#define CGAL_SKIN_SURFACE_TRAITS_3_H
#include <CGAL/license/Skin_surface_3.h>
#include <CGAL/predicates/predicates_for_mixed_complex_3.h>
namespace CGAL {
/** Input: a list of n weighted points p_1...p_n and a query point x.
There is a plane separating the mixed cell defined by p_1...p_n-1
and the mixed cell defined by p_1...p_n. The predicate tests
whether x lies on the same side of this plane as the mixed cell
defined by p_1...p_n-1 (NEGATIVE), on the plane (ZERO) or on the
opposite side (POSITIVE).
**/
template <class K>
class Side_of_mixed_cell_3
{
public:
typedef typename K::FT FT;
typedef typename K::Point_3 Point_3;
typedef typename K::Weighted_point_3 Weighted_point_3;
typedef CGAL::Sign result_type;
Side_of_mixed_cell_3(const FT &shrink) : s(shrink) { }
result_type operator()(const Weighted_point_3 &p1,
const Weighted_point_3 &p2,
const Point_3 &x) const
{
return side_of_mixed_cellC3(p1.x(),p1.y(),p1.z(),p1.weight(),
p2.x(),p2.y(),p2.z(),p2.weight(),
x.x(),x.y(),x.z(),
s);
}
result_type operator()(const Weighted_point_3 &p1,
const Weighted_point_3 &p2,
const Weighted_point_3 &p3,
const Point_3 &x) const
{
return side_of_mixed_cellC3(p1.x(),p1.y(),p1.z(),p1.weight(),
p2.x(),p2.y(),p2.z(),p2.weight(),
p3.x(),p3.y(),p3.z(),p3.weight(),
x.x(),x.y(),x.z(),
s);
}
result_type operator()(const Weighted_point_3 &p1,
const Weighted_point_3 &p2,
const Weighted_point_3 &p3,
const Weighted_point_3 &p4,
const Point_3 &x) const
{
return side_of_mixed_cellC3(p1.x(),p1.y(),p1.z(),p1.weight(),
p2.x(),p2.y(),p2.z(),p2.weight(),
p3.x(),p3.y(),p3.z(),p3.weight(),
p4.x(),p4.y(),p4.z(),p4.weight(),
x.x(),x.y(),x.z(),
s);
}
private:
FT s;
};
/** Input: Two weighted points
Computes the anchor point of a Delaunay center and a Voronoi center
**/
template <class K>
class Construct_anchor_point_3
{
public:
typedef typename K::FT FT;
typedef typename K::Point_3 Point_3;
typedef Point_3 result_type;
Construct_anchor_point_3(const FT &shrink) : s(shrink) {}
result_type operator()(const Point_3 &p_del,
const Point_3 &p_vor) const
{
return Point_3((1-s)*p_del.x() + s*p_vor.x(),
(1-s)*p_del.y() + s*p_vor.y(),
(1-s)*p_del.z() + s*p_vor.z());
}
private:
FT s;
};
template <class K_>
class Skin_surface_traits_base_3
: public K_
{
public:
typedef K_ Kernel;
typedef Skin_surface_traits_base_3<Kernel> Self;
typedef typename Kernel::FT FT;
typedef CGAL::Side_of_mixed_cell_3<Self> Side_of_mixed_cell_3;
typedef CGAL::Construct_anchor_point_3<Self> Construct_anchor_point_3;
Skin_surface_traits_base_3() : s(-1) { }
Skin_surface_traits_base_3(FT s) : s(s) { }
void set_shrink(FT s_) {
s = s_;
}
FT get_shrink() const {
return s;
}
Side_of_mixed_cell_3 side_of_mixed_cell_3_object() const
{
CGAL_assertion((s>0) && (s<1));
return Side_of_mixed_cell_3(get_shrink());
}
Construct_anchor_point_3 construct_anchor_point_3_object() const {
return Construct_anchor_point_3(get_shrink());
}
private:
FT s;
};
// We need to introduce a "traits_base_3" class in order to get the
// specialization for Exact_predicates_inexact_constructions_kernel to work,
// otherwise there is a cycle in the derivation.
template < class K, bool UseFilteredPredicates=K::Has_filtered_predicates >
class Skin_surface_traits_3
: public Skin_surface_traits_base_3<K>
{
typedef Skin_surface_traits_base_3<K> Base;
public:
Skin_surface_traits_3() {}
Skin_surface_traits_3(typename Base::FT s) : Base(s) {}
enum { Has_filtered_predicates=false };
};
} //namespace CGAL
// Now specialize for Filtered_kernel<CK>, to get
// the filtered traits automatically.
#include <CGAL/Skin_surface_filtered_traits_3.h>
#include <CGAL/Filtered_kernel.h>
namespace CGAL {
// Just FK would be nicer, but VC 2005 messes it up with an "FK" in a base class when compiling degenerate_test.cpp
template < typename Sst3FK >
class Skin_surface_traits_3 < Sst3FK, true >
: public Skin_surface_filtered_traits_3 < Sst3FK >
{
typedef Skin_surface_filtered_traits_3 < Sst3FK > Base;
public:
typedef Sst3FK Kernel;
Skin_surface_traits_3() {}
Skin_surface_traits_3(typename Base::FT s) : Base(s) { }
};
} //namespace CGAL
#endif // CGAL_SKIN_SURFACE_TRAITS_3_H