songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

example added

This commit is contained in:
Bernd Gärtner 2006-03-10 12:34:32 +00:00
parent 0aa2469e6c
commit d9e740ff1a
2 changed files with 51 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
.gitattributes vendored
View File

@ -1011,6 +1011,7 @@ Mesh_2/doc_tex/Mesh_2/example-conform.eps -text svneol=unset#unset
Mesh_2/doc_tex/Mesh_2/example-conform.png -text svneol=unset#unset
Mesh_2/doc_tex/Mesh_2_ref/part_of_a_cluster.eps -text
Mesh_2/doc_tex/Mesh_2_ref/part_of_a_cluster.gif -text svneol=unset#unset
Min_annulus_d/examples/min_annulus_d.C -text
Min_circle_2/demo/Min_circle_2/min_circle_2.vcproj -text
Min_ellipse_2/demo/Min_ellipse_2/min_ellipse_2.vcproj -text
Min_quadrilateral_2/demo/Min_quadrilateral_2/min_quadrilateral_2.vcproj -text

50
Min_annulus_d/examples/min_annulus_d.C Normal file
View File

@ -0,0 +1,50 @@
// file: examples/Min_annulus_d/min_annulus_d.C
// computes the smallest enclosing annulus of two point
// sets on nested squares in R^2
#include <CGAL/Simple_cartesian.h>
#include <CGAL/Min_annulus_d.h>
#include <CGAL/Optimisation_d_traits_2.h>
typedef CGAL::Simple_cartesian<double> K;
typedef K::Point_2 Point;
typedef CGAL::Optimisation_d_traits_2<K> Traits;
typedef CGAL::Min_annulus_d<Traits> Min_annulus;
int main()
{
// points on the squares [-1,1]^2 and [-2,2]^2
Point P[8] = { Point(-1,-1), Point(-1,1), Point(1,-1), Point(1,1),
Point(-2,-2), Point(-2,2), Point(2,-2), Point(2,2)};
Min_annulus ma(P, P+8);
// get center of annulus
Min_annulus::Coordinate_iterator coord_it;
std::cout << "center:"; // homogeneous point, (0,0,1)
for (coord_it = ma.center_coordinates_begin();
coord_it != ma.center_coordinates_end();
++coord_it)
std::cout << " " << *coord_it;
std::cout << std::endl;
// get inner squared radius, 1^2+1^2 = 2
std::cout << "Inner squared radius: " <<
ma.squared_inner_radius_numerator() << " / " <<
ma.squared_radii_denominator() << " ~ " <<
ma.squared_inner_radius_numerator() /
ma.squared_radii_denominator() << std::endl;
// get outer squared radius, 2^2+2^2 = 8
std::cout << "Outer squared radius: " <<
ma.squared_outer_radius_numerator() << " / " <<
ma.squared_radii_denominator() << " ~ " <<
ma.squared_outer_radius_numerator() /
ma.squared_radii_denominator() << std::endl;
return 0;
}