void Edges::clip_line(Edge *e) {
Site *s1, *s2;
double x1,x2,y1,y2;
if(e -> a == 1.0 && e ->b >= 0.0) {
s1 = e -> ep[1];
s2 = e -> ep[0];
}
else {
s1 = e -> ep[0];
s2 = e -> ep[1];
}
#ifdef THROW_LESS_AWAY
if(e -> a == 1.0) {
if (s1!=(Site *)0) {
y1 = s1->coord.y;
if (y1 > bounds.t) {
y1 = bounds.t;
x1 = e -> c - e -> b * y1;
}
else if (y1 >= bounds.b) x1 = s1->coord.x;
else {
y1 = bounds.b;
x1 = e -> c - e -> b * y1;
}
}
else {
y1 = bounds.b;
x1 = e -> c - e -> b * y1;
}
if (s2!=(Site *)0) {
y2 = s2->coord.y;
if(y2<bounds.b) {
y2 = bounds.b;
x2 = e -> c - e -> b * y2;
}
else if (y2 <= bounds.t) x2 = s2->coord.x;
else {
y2 = bounds.t;
x2 = e -> c - e -> b * y2;
}
}
else {
y2 = bounds.t;
x2 = e -> c - e -> b * y2;
}
if ((x1> bounds.r && x2>bounds.r) || (x1<bounds.l&&x2<bounds.l))
return;
if(x1> bounds.r)
{ x1 = bounds.r; y1 = (e -> c - x1)/e -> b;};
if(x1<bounds.l)
{ x1 = bounds.l; y1 = (e -> c - x1)/e -> b;};
if(x2>bounds.r)
{ x2 = bounds.r; y2 = (e -> c - x2)/e -> b;};
if(x2<bounds.l)
{ x2 = bounds.l; y2 = (e -> c - x2)/e -> b;};
}
else {
if (s1!=(Site *)0) {
x1 = s1->coord.x;
if(x1>bounds.r) {
x1 = bounds.r;
y1 = e -> c - e -> a * x1;
}
else if (x1 >= bounds.l) y1 = s1->coord.y;
else {
x1 = bounds.l;
y1 = e -> c - e -> a * x1;
}
}
else {
x1 = bounds.l;
y1 = e -> c - e -> a * x1;
}
if (s2!=(Site *)0) {
x2 = s2->coord.x;
if(x2<bounds.l) {
x2 = bounds.l;
y2 = e -> c - e -> a * x2;
}
else if (x2 <= bounds.r) y2 = s2->coord.y;
else {
x2 = bounds.r;
y2 = e -> c - e -> a * x2;
}
}
else {
x2 = bounds.r;
y2 = e -> c - e -> a * x2;
}
if ((y1> bounds.t && y2>bounds.t) || (y1<bounds.b && y2<bounds.b))
return;
if(y1> bounds.t)
{ y1 = bounds.t; x1 = (e -> c - y1)/e -> a;};
if(y1<bounds.b)
{ y1 = bounds.b; x1 = (e -> c - y1)/e -> a;};
if(y2>bounds.t)
{ y2 = bounds.t; x2 = (e -> c - y2)/e -> a;};
if(y2<bounds.b)
{ y2 = bounds.b; x2 = (e -> c - y2)/e -> a;};
}
#else
if(e -> a == 1.0) {
if (s1!=(Site *)0) {
y1 = s1->coord.y;
if (y1 > bounds.t)
return;
else if (y1 >= bounds.b) x1 = s1->coord.x;
else {
y1 = bounds.b;
x1 = e -> c - e -> b * y1;
}
}
else {
y1 = bounds.b;
x1 = e -> c - e -> b * y1;
}
if (s2!=(Site *)0) {
y2 = s2->coord.y;
if(y2<bounds.b)
return;
else if (y2 <= bounds.t) x2 = s2->coord.x;
else {
y2 = bounds.t;
x2 = e -> c - e -> b * y2;
}
}
else {
y2 = bounds.t;
x2 = e -> c - e -> b * y2;
}
if ((x1> bounds.r && x2>bounds.r) || (x1<bounds.l&&x2<bounds.l))
return;
if(x1> bounds.r)
{ x1 = bounds.r; y1 = (e -> c - x1)/e -> b;};
if(x1<bounds.l)
{ x1 = bounds.l; y1 = (e -> c - x1)/e -> b;};
if(x2>bounds.r)
{ x2 = bounds.r; y2 = (e -> c - x2)/e -> b;};
if(x2<bounds.l)
{ x2 = bounds.l; y2 = (e -> c - x2)/e -> b;};
}
else {
if (s1!=(Site *)0) {
x1 = s1->coord.x;
if(x1>bounds.r)
return;
else if (x1 >= bounds.l) y1 = s1->coord.y;
else {
x1 = bounds.l;
y1 = e -> c - e -> a * x1;
}
}
else {
x1 = bounds.l;
y1 = e -> c - e -> a * x1;
}
if (s2!=(Site *)0) {
x2 = s2->coord.x;
if(x2<bounds.l)
return;
else if (x2 <= bounds.r) y2 = s2->coord.y;
else {
x2 = bounds.r;
y2 = e -> c - e -> a * x2;
}
}
else {
x2 = bounds.r;
y2 = e -> c - e -> a * x2;
}
if ((y1> bounds.t && y2>bounds.t) || (y1<bounds.b && y2<bounds.b))
return;
if(y1> bounds.t)
{ y1 = bounds.t; x1 = (e -> c - y1)/e -> a;};
if(y1<bounds.b)
{ y1 = bounds.b; x1 = (e -> c - y1)/e -> a;};
if(y2>bounds.t)
{ y2 = bounds.t; x2 = (e -> c - y2)/e -> a;};
if(y2<bounds.b)
{ y2 = bounds.b; x2 = (e -> c - y2)/e -> a;};
}
#endif
doSeg(e,x1,y1,x2,y2);
#ifdef VORLINES_CLIPEDGES
Layout *l = infos.nodes.front().layoutN->g;
Line seg;
seg.degree = 1;
seg.push_back(Coord(x1,y1));
seg.push_back(Coord(x2,y2));
gd<Drawn>(l).push_back(seg);
#endif
#ifdef STANDALONE
if (doPS) line (x1,y1,x2,y2);
#endif
}
void clip_line(Edge * e)
{
Site *s1, *s2;
double x1, x2, y1, y2;
if (e->a == 1.0 && e->b >= 0.0) {
s1 = e->ep[1];
s2 = e->ep[0];
} else {
s1 = e->ep[0];
s2 = e->ep[1];
}
if (e->a == 1.0) {
if (s1 != (Site *) NULL) {
y1 = s1->coord.y;
if (y1 > pymax)
return;
else if (y1 >= pymin)
x1 = s1->coord.x;
else {
y1 = pymin;
x1 = e->c - e->b * y1;
}
} else {
y1 = pymin;
x1 = e->c - e->b * y1;
}
if (s2 != (Site *) NULL) {
y2 = s2->coord.y;
if (y2 < pymin)
return;
else if (y2 <= pymax)
x2 = s2->coord.x;
else {
y2 = pymax;
x2 = e->c - e->b * y2;
}
} else {
y2 = pymax;
x2 = e->c - e->b * y2;
}
if (((x1 > pxmax) & (x2 > pxmax)) | ((x1 < pxmin) & (x2 < pxmin)))
return;
if (x1 > pxmax) {
x1 = pxmax;
y1 = (e->c - x1) / e->b;
};
if (x1 < pxmin) {
x1 = pxmin;
y1 = (e->c - x1) / e->b;
};
if (x2 > pxmax) {
x2 = pxmax;
y2 = (e->c - x2) / e->b;
};
if (x2 < pxmin) {
x2 = pxmin;
y2 = (e->c - x2) / e->b;
};
} else {
if (s1 != (Site *) NULL) {
x1 = s1->coord.x;
if (x1 > pxmax)
return;
else if (x1 >= pxmin)
y1 = s1->coord.y;
else {
x1 = pxmin;
y1 = e->c - e->a * x1;
}
} else {
x1 = pxmin;
y1 = e->c - e->a * x1;
}
if (s2 != (Site *) NULL) {
x2 = s2->coord.x;
if (x2 < pxmin)
return;
else if (x2 <= pxmax)
y2 = s2->coord.y;
else {
x2 = pxmax;
y2 = e->c - e->a * x2;
}
} else {
x2 = pxmax;
y2 = e->c - e->a * x2;
}
if (((y1 > pymax) & (y2 > pymax)) | ((y1 < pymin) & (y2 < pymin)))
return;
if (y1 > pymax) {
y1 = pymax;
x1 = (e->c - y1) / e->a;
};
if (y1 < pymin) {
y1 = pymin;
x1 = (e->c - y1) / e->a;
};
if (y2 > pymax) {
y2 = pymax;
x2 = (e->c - y2) / e->a;
};
if (y2 < pymin) {
y2 = pymin;
x2 = (e->c - y2) / e->a;
};
}
doSeg(e, x1, y1, x2, y2);
#ifdef STANDALONE
if (doPS)
line(x1, y1, x2, y2);
#endif
}
