1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
// Copyright (c) 2006 INRIA Sophia-Antipolis (France).
// All rights reserved.
//
// This file is part of CGAL (www.cgal.org).
// You can redistribute it and/or modify it under the terms of the GNU
// General Public License as published by the Free Software Foundation,
// either version 3 of the License, or (at your option) any later version.
//
// Licensees holding a valid commercial license may use this file in
// accordance with the commercial license agreement provided with the software.
//
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
//
// $URL: svn+ssh://scm.gforge.inria.fr/svn/cgal/branches/next/Mesh_3/examples/Mesh_3/implicit_functions.h $
// $Id: implicit_functions.h 67610 2012-02-06 11:39:56Z lrineau $
//
//
// Author(s) : Steve OUDOT, Laurent RINEAU
///////////////// Definitions of several famous surfaces /////////////////
double sphere_function (double, double, double); // (c=(0,0,0), r=1)
double ellipsoid_function (double, double, double); // (c=(0,0,0), r=1)
double torus_function (double, double, double); // (c=(0,0,0), r=2)
double chair_function (double, double, double); // (c=(0,0,0), r=6)
double tanglecube_function (double, double, double); // (c=(0,0,0), r=4)
double octic_function (double, double, double); // (c=(0,0,0), r=2)
double heart_function (double, double, double); // (c=(0,0,0), r=2)
double klein_function (double, double, double); // (c=(0,0,0), r=4)
double ring_function (double, double, double); // (c=(0,0,0), r=?)
double false_knot_function (double, double, double); // (c=(0,0,0), r=1)
double knot1_function (double, double, double); // (c=(0,0,0), r=4)
double knot2_function (double, double, double); // (c=(0,0,0), r=4)
double knot3_function (double, double, double); // (c=(0,0,0), r=4)
double cube_function (double, double, double); // (c=(0,0,0), r=2)
template <int Sq_radius>
double sphere_function (double x, double y, double z) // (c=(0,0,0), r=Sq_radius)
{
double x2=x*x, y2=y*y, z2=z*z;
return (x2+y2+z2)/Sq_radius - 1;
}
template <typename FT, typename Point>
class FT_to_point_function_wrapper : public std::unary_function<Point, FT>
{
typedef FT (*Implicit_function)(FT, FT, FT);
Implicit_function function;
public:
FT_to_point_function_wrapper(Implicit_function f) : function(f) {}
FT operator()(Point p) const { return function(p.x(), p.y(), p.z()); }
};
|
