triangulate(tri_options.c_str(), &in, &out, &vorout);
// TEMPORARY
- write_out_data(&out);
+ // write_out_data(&out);
// now copy the results back into the corresponding FGTriangle
// structures
if ( fabs(p0.x() - p1.x()) > 3 * FG_EPSILON ) {
// use y = mx + b
- // sort these in a sensable order
+ // sort these in a sensible order
if ( p0.x() > p1.x() ) {
Point3D tmp = p0;
p0 = p1;
y_err = fabs(current->y() - (m * current->x() + b));
- if ( y_err < 10 * FG_EPSILON ) {
+ if ( y_err < 20 * FG_EPSILON ) {
// cout << "FOUND EXTRA SEGMENT NODE (Y)" << endl;
// cout << p0 << " < " << *current << " < "
// << p1 << endl;
x_err = fabs(current->x() - p0.x());
// cout << " found a potential point, x err = "
// << x_err << endl;
- if ( x_err < 10*FG_EPSILON ) {
+ if ( x_err < 20 * FG_EPSILON ) {
// cout << "FOUND EXTRA SEGMENT NODE (X)" << endl;
// cout << p0 << " < " << *current << " < "
// << p1 << endl;