void
arrangement::compute_criticalCurves_type_II()
{
for (int i = 0; i < (int) convolutions_o.size(); ++i)
{
// for all convolution arrangements
Arrangement_2::Edge_iterator first = convolutions_o[i].edges_begin();
Arrangement_2::Edge_iterator previous = convolutions_o[i].edges_begin();
bool flag = true;
Arrangement_2 copy(convolutions_o[i]);
Walk_pl walk_pl(copy);
for (Arrangement_2::Edge_iterator edge = convolutions_o[i].edges_begin(); edge != convolutions_o[i].edges_end(); ++edge)
{
if (flag)
{
++edge;
flag = false;
}
if (previous->data().compare(edge->data()) != 0)
{
// add critical curve type II in the copy
double x = CGAL::to_double(edge->curve().source().x());
double y = CGAL::to_double(edge->curve().source().y());
double radius = r1r2;
Rat_point_2 center(x, y);
Rat_circle_2 circle(center, radius * radius);
Conic_curve_2 conic_arc(circle);
insert(nonCriticalRegions, conic_arc);
}
previous = edge;
}
if (previous->data().compare(first->data()) != 0)
{
// add critical curve type II in the copy
double x = CGAL::to_double(first->curve().source().x());
double y = CGAL::to_double(first->curve().source().y());
double radius = r1r2;
Rat_point_2 center(x, y);
Rat_circle_2 circle(center, radius * radius);
Conic_curve_2 conic_arc(circle);
insert(nonCriticalRegions, conic_arc);
}
// keep only wanted curves and put it in the arrangement
for (Arrangement_2::Edge_iterator e = nonCriticalRegions.edges_begin(); e != nonCriticalRegions.edges_end(); ++e)
if ((e->data().compare("") == 0))
keep_arc(e, copy, walk_pl);
for (Arrangement_2::Edge_iterator edge = convolutions_o[i].edges_begin(); edge != convolutions_o[i].edges_end(); ++edge)
insert(nonCriticalRegions, edge->curve());
}
}
void CriticalCurves::setParameters(double radius_1, double radius_2, Arrangements_2 insets_1, Arrangements_2 insets_2)
{
Arrangement_2_iterator inset_1 = insets_1.begin();
Arrangement_2_iterator inset_2 = insets_2.begin();
while (inset_1 != insets_1.end() && inset_2 != insets_2.end())
{
Arrangement_2 arrangement;
// Add the curves of the inset.
for (Edge_iterator edge = inset_1->edges_begin(); edge != inset_1->edges_end(); ++edge)
{
insert(arrangement, edge->curve());
}
// Add the critical curves of type I.
for (Edge_iterator edge = inset_2->edges_begin(); edge != inset_2->edges_end(); ++edge)
{
if (CGAL::COLLINEAR == edge->curve().orientation())
{
// Displaced a segment.
Nt_traits nt_traits;
Algebraic_ft factor = nt_traits.convert(Rational(radius_1) + Rational(radius_2));
Conic_point_2 source = edge->curve().source();
Conic_point_2 target = edge->curve().target();
Algebraic_ft delta_x = target.x() - source.x();
Algebraic_ft delta_y = target.y() - source.y();
Algebraic_ft length = nt_traits.sqrt(delta_x * delta_x + delta_y * delta_y);
Algebraic_ft translation_x = factor * delta_y / length;
Algebraic_ft translation_y = - factor * delta_x / length;
Conic_point_2 point_1(source.x() + translation_x, source.y() + translation_y);
Conic_point_2 point_2(target.x() + translation_x, target.y() + translation_y);
Algebraic_ft a = - delta_y;
Algebraic_ft b = delta_x;
Algebraic_ft c = factor * length - (source.y() * target.x() - source.x() * target.y());
X_monotone_curve_2 x_monotone_curve(a, b, c, point_1, point_2);
insert(arrangement, x_monotone_curve);
}
else
{
// Displaces an arc.
Rational two(2);
Rational four(4);
Rational r = edge->curve().r();
Rational s = edge->curve().s();
Rational t = edge->curve().t();
Rational u = edge->curve().u();
Rational v = edge->curve().v();
Rational w = edge->curve().w();
Nt_traits nt_traits;
Rational x_center = - u / (two * r);
Rational y_center = - v / (two * r);
Rat_point_2 rat_center(x_center, y_center);
Conic_point_2 center(nt_traits.convert(x_center), nt_traits.convert(y_center));
Rational radius = Rational(radius_1) + two * Rational(radius_2);
Algebraic_ft coefficient = nt_traits.convert(radius / Rational(radius_2));
Conic_point_2 source_1 = edge->curve().source();
Algebraic_ft x_source_2 = center.x() + coefficient * (source_1.x() - center.x());
Algebraic_ft y_source_2 = center.y() + coefficient * (source_1.y() - center.y());
Conic_point_2 source_2(x_source_2, y_source_2);
Conic_point_2 target_1 = edge->curve().target();
Algebraic_ft x_target_2 = center.x() + coefficient * (target_1.x() - center.x());
Algebraic_ft y_target_2 = center.y() + coefficient * (target_1.y() - center.y());
Conic_point_2 target_2(x_target_2, y_target_2);
Rat_circle_2 circle(rat_center, radius * radius);
Conic_arc_2 conic_arc(circle, CGAL::COUNTERCLOCKWISE, source_2, target_2);
insert(arrangement, conic_arc);
}
}
// Add the critical curves of type II.
for (Edge_iterator edge = inset_2->edges_begin(); edge != inset_2->edges_end(); ++edge)
{
double x = CGAL::to_double(edge->curve().source().x());
double y = CGAL::to_double(edge->curve().source().y());
double radius = radius_1 + radius_2;
Rat_point_2 center(x, y);
Rat_circle_2 circle(center, radius * radius);
Conic_arc_2 conic_arc(circle);
insert(arrangement, conic_arc);
}
// Remove the curves which are not include in the inset.
Objects objects;
Face_handle face;
for (Edge_iterator edge = arrangement.edges_begin(); edge != arrangement.edges_end(); ++edge)
{
CGAL::zone(*inset_1, edge->curve(), std::back_inserter(objects));
for (Object_iterator object = objects.begin(); object != objects.end(); ++object)
{
if (assign(face, *object))
{
if (face->is_unbounded())
{
remove_edge(arrangement, edge);
break;
}
}
}
objects.clear();
}
// Print essential information on the standard input.
std::cout << "Arrangement:" << std::endl;
std::cout << " Number of vertices: " << arrangement.number_of_vertices() << std::endl;
std::cout << " Number of edges : " << arrangement.number_of_edges() << std::endl;
std::cout << " Number of face : " << arrangement.number_of_faces() << std::endl;
this->critical_curves.push_back(arrangement);
++inset_1;
++inset_2;
}
// Commit changes.
emit(criticalCurvesChanged());
return;
}
void
arrangement::compute_ACScell()
{
/*
for (int i = 0; i < (int)convolutions.size(); ++i)
for (Arrangement_2::Edge_iterator edge = convolutions[i].edges_begin(); edge != convolutions[i].edges_end(); ++edge)
{
convolution_r_all = Arrangement_2(convolutions[i]);
insert(convolution_r_all, edge->curve());
Arrangement_2::Edge_iterator e = convolution_r_all.edges_end();
e->set_data(edge->data());
}
*/
/*
std::cout << "edge label convo : ";
for (Arrangement_2::Edge_iterator edge = convolution_r_all.edges_begin(); edge != convolution_r_all.edges_end(); ++edge)
std::cout << edge->data() << " ";
std::cout << std::endl;
*/
for (int i = 0; i < (int)point_in_faces.size(); ++i)
{
// do a copy
Arrangement_2 copy(convolutions[0]);
Observer observer(copy);
// add a circle
Rat_point_2 center(point_in_faces[i][0], point_in_faces[i][1]);
Rat_circle_2 circle(center, r1r2 * r1r2);
Conic_curve_2 conic_arc(circle);
insert(copy, conic_arc);
// then add rho label for arc
for (Arrangement_2::Edge_iterator edge = copy.edges_begin(); edge != copy.edges_end(); ++edge)
{
if (edge->data().compare("") == 0)
{
edge->set_data("rho_");
edge->twin()->set_data("rho_");
}
}
Walk_pl walk_pl(copy);
std::vector<std::vector<double> > points;
compute_pointInCells(copy, points);
for (int j = 0; j < (int)points.size(); ++j)
{
double _x = point_in_faces[i][0] - points[j][0];
double _y = point_in_faces[i][1] - points[j][1];
if (r1r2 < sqrt(_x * _x + _y * _y))
{
Rational x_(points[j][0]);
Rational y_(points[j][1]);
Conic_point_2 p(x_,y_);
Arrangement_2::Vertex_handle v = insert_point(copy, p, walk_pl);
try
{
ACScell cell(i);
Arrangement_2::Ccb_halfedge_circulator first_outer_ccb = v->face()->outer_ccb();
Arrangement_2::Ccb_halfedge_circulator outer_ccb = v->face()->outer_ccb();
do
{
// for retrieving ACScell Begin and End
insert(cell.arr, outer_ccb->curve());
// then continue
cell.addLabel(outer_ccb->data());
++outer_ccb;
} while (outer_ccb != first_outer_ccb);
ACScells.push_back(cell);
}
catch (const std::exception exn) {}
}
}
}
// clean ACScells
for (int i = 0; i < (int)ACScells.size(); ++i)
ACScells[i].cleanLabels();
// compute ACScells id
int cpt = 0;
int previous = 0;
for (int i = 0; i < (int)ACScells.size(); ++i)
{
if (previous != ACScells[i].NCR)
{
cpt = 0;
previous = ACScells[i].NCR;
}
ACScells[i].id = std::to_string(ACScells[i].NCR) + "." + std::to_string(cpt);
++cpt;
}
}
void
arrangement::compute_criticalCurves_type_I()
{
double radius_1 = ((double)manipulator_diametre)/2.0;
double radius_2 = ((double)target_diametre)/2.0;
for (int i = 0; i < (int) convolutions_o.size(); ++i)
{
Arrangement_2 copy(convolutions_o[i]);
Walk_pl walk_pl(copy);
// Add the critical curves of type I.
for (Arrangement_2::Edge_iterator edge = convolutions_o[i].edges_begin(); edge != convolutions_o[i].edges_end(); ++edge)
{
if (CGAL::COLLINEAR == edge->curve().orientation())
{
// Displaced a segment.
Nt_traits nt_traits;
Algebraic_ft factor = nt_traits.convert(Rational(radius_1) + Rational(radius_2));
Conic_point_2 source = edge->curve().source();
Conic_point_2 target = edge->curve().target();
Algebraic_ft delta_x = target.x() - source.x();
Algebraic_ft delta_y = target.y() - source.y();
Algebraic_ft length = nt_traits.sqrt(delta_x * delta_x + delta_y * delta_y);
Algebraic_ft translation_x = - factor * delta_y / length;
Algebraic_ft translation_y = factor * delta_x / length;
Conic_point_2 point_1(source.x() + translation_x, source.y() + translation_y);
Conic_point_2 point_2(target.x() + translation_x, target.y() + translation_y);
Algebraic_ft a = delta_y;
Algebraic_ft b = - delta_x;
Algebraic_ft c = factor * length + (source.y() * target.x() - source.x() * target.y());
X_monotone_curve_2 x_monotone_curve(a, b, c, point_1, point_2);
insert(nonCriticalRegions, x_monotone_curve);
}
else
{
// Displaces an arc.
Rational two(2);
Rational r = edge->curve().r();
Rational u = edge->curve().u();
Rational v = edge->curve().v();
Nt_traits nt_traits;
Rational x_center = - u / (two * r);
Rational y_center = - v / (two * r);
Rat_point_2 rat_center(x_center, y_center);
Conic_point_2 center(nt_traits.convert(x_center), nt_traits.convert(y_center));
Rational radius = Rational(radius_1) + two * Rational(radius_2);
Algebraic_ft coefficient = nt_traits.convert(radius / Rational(radius_2));
Conic_point_2 source_1 = edge->curve().source();
Algebraic_ft x_source_2 = center.x() + coefficient * (source_1.x() - center.x());
Algebraic_ft y_source_2 = center.y() + coefficient * (source_1.y() - center.y());
Conic_point_2 source_2(x_source_2, y_source_2);
Conic_point_2 target_1 = edge->curve().target();
Algebraic_ft x_target_2 = center.x() + coefficient * (target_1.x() - center.x());
Algebraic_ft y_target_2 = center.y() + coefficient * (target_1.y() - center.y());
Conic_point_2 target_2(x_target_2, y_target_2);
Rat_circle_2 circle(rat_center, radius * radius);
Conic_curve_2 conic_arc(circle, CGAL::CLOCKWISE, source_2, target_2);
insert(nonCriticalRegions, conic_arc);
}
}
// keep only wanted curves and put it in the arrangement
for (Arrangement_2::Edge_iterator e = nonCriticalRegions.edges_begin(); e != nonCriticalRegions.edges_end(); ++e)
if ((e->data().compare("") == 0))
keep_arc(e, copy, walk_pl);
for (Arrangement_2::Edge_iterator edge = convolutions_o[i].edges_begin(); edge != convolutions_o[i].edges_end(); ++edge)
insert(nonCriticalRegions, edge->curve());
}
}
