mirror of https://github.com/CGAL/cgal
55 lines
2.4 KiB
C++
55 lines
2.4 KiB
C++
#include <CGAL/Cartesian.h>
|
|
#include <CGAL/point_generators_2.h>
|
|
#include <CGAL/MP_Float.h>
|
|
#include <CGAL/Algebraic_kernel_for_circles_2_2.h>
|
|
#include <CGAL/Circular_kernel_2.h>
|
|
|
|
typedef CGAL::Quotient<CGAL::MP_Float> NT;
|
|
typedef CGAL::Cartesian<NT> Linear_k;
|
|
typedef CGAL::Algebraic_kernel_for_circles_2_2<NT> Algebraic_k;
|
|
typedef CGAL::Circular_kernel_2<Linear_k,Algebraic_k> Circular_k;
|
|
typedef Circular_k::Point_2 Point_2;
|
|
typedef Circular_k::Circle_2 Circle_2;
|
|
typedef Circular_k::Circular_arc_2 Circular_arc_2;
|
|
typedef Circular_k::Line_arc_2 Line_arc_2;
|
|
|
|
template <typename T>
|
|
double prob_2() {
|
|
CGAL::Random_points_in_square_2<Point_2> g(1.0);
|
|
double prob = 0.0;
|
|
for (int i = 0; i < 10000; i++) {
|
|
Point_2 p1, p2, p3, p4, p5, p6;
|
|
p1 = *g++; p2 = *g++; p3 = *g++;
|
|
p4 = *g++; p5 = *g++; p6 = *g++;
|
|
if(Linear_k().orientation_2_object()(Linear_k::Point_2(p1.x(), p1.y()),
|
|
Linear_k::Point_2(p2.x(), p2.y()),
|
|
Linear_k::Point_2(p3.x(), p3.y())) !=
|
|
CGAL::COUNTERCLOCKWISE) std::swap(p1, p3);
|
|
T o1 = T(p1, p2, p3);
|
|
if(Linear_k().orientation_2_object()(Linear_k::Point_2(p4.x(), p4.y()),
|
|
Linear_k::Point_2(p5.x(), p5.y()),
|
|
Linear_k::Point_2(p6.x(), p6.y())) !=
|
|
CGAL::COUNTERCLOCKWISE) std::swap(p4, p6);
|
|
T o2 = T(p4, p5, p6);
|
|
std::vector< CGAL::Object > res;
|
|
Circular_k().intersect_2_object()(o1, o2, std::back_inserter(res));
|
|
prob += (res.size() != 0) ? 1.0 : 0.0;
|
|
}
|
|
return prob/10000.0;
|
|
}
|
|
|
|
int main()
|
|
{
|
|
std::cout << "What is the probability that two arcs formed by" << std::endl;
|
|
std::cout << "three random counterclockwise-oriented points on" << std::endl;
|
|
std::cout << "an unity square intersect?" << std::endl;
|
|
std::cout << "The probability is: " << prob_2<Circular_arc_2>() <<
|
|
std::endl << std::endl;
|
|
std::cout << "And what about the probability that two circles formed by"
|
|
<< std::endl;
|
|
std::cout << "three random counterclockwise-oriented points on" << std::endl;
|
|
std::cout << "an unity square intersect?" << std::endl;
|
|
std::cout << "The probability is: " << prob_2<Circle_2>() << std::endl;
|
|
return 0;
|
|
};
|