private static void RenderStrip(Tesselator tess, HalfEdge halfEdge, int size)
{
/* Render as many CCW triangles as possible in a strip starting from
* edge "e". The strip *should* contain exactly "size" triangles
* (otherwise we've goofed up somewhere).
*/
tess.CallBegin(TriangleListType.TriangleStrip);
tess.CallVertex(halfEdge.originVertex.ClientIndex);
tess.CallVertex(halfEdge.directionVertex.ClientIndex);
while (!halfEdge.leftFace.Marked())
{
halfEdge.leftFace.marked = true;
--size;
halfEdge = halfEdge.Dprev;
tess.CallVertex(halfEdge.originVertex.ClientIndex);
if (halfEdge.leftFace.Marked())
{
break;
}
halfEdge.leftFace.marked = true;
--size;
halfEdge = halfEdge.nextEdgeCCWAroundOrigin;
tess.CallVertex(halfEdge.directionVertex.ClientIndex);
}
if (size != 0)
{
throw new Exception();
}
tess.CallEnd();
}
public Dictionary(Tesselator tesselator)
{
Node nodeHead;
nodeHead = this.head;
nodeHead.Key = null;
nodeHead.next = nodeHead;
nodeHead.prev = nodeHead;
this.tesseator = tesselator;
}
public Dictionary(Tesselator tesselator)
{
Node nodeHead = _head;
nodeHead.Key = null;
nodeHead.next = nodeHead;
nodeHead.prev = nodeHead;
_tess = tesselator;
}
public Dictionary(Tesselator tesselator)
{
Node nodeHead;
nodeHead = this.head;
nodeHead.Key = null;
nodeHead.next = nodeHead;
nodeHead.prev = nodeHead;
this.tesseator = tesselator;
}
private void RenderTriangle(Tesselator tess, HalfEdge e, int size)
{
/* Just add the triangle to a triangle list, so we can render all
* the separate triangles at once.
*/
if (size != 1)
{
throw new Exception();
}
Face.AddToTrail(ref e.leftFace, ref this.lonelyTriList);
}
public void Connect(Tesselate.Tesselator tesselator, bool setEdgeFlag)
{
tesselator.callBegin = BeginCallBack;
tesselator.callEnd = EndCallBack;
tesselator.callVertex = VertexCallBack;
tesselator.callCombine = CombineCallBack;
if (setEdgeFlag)
{
tesselator.callEdgeFlag = EdgeFlagCallBack;
}
}
public void Connect(List <Vertex> vertextList,
Tesselate.Tesselator tesselator,
Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
this._inputVertextList = vertextList;
tesselator.WindingRule = windingRule;
NeedEdgeFlag = setEdgeFlag;
tesselator.SetListener(this);
}
public void Connect(Tesselate.Tesselator tesselator, Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
tesselator.callBegin = new Tesselate.Tesselator.CallBeginDelegate(BeginCallBack);
tesselator.callEnd = new Tesselate.Tesselator.CallEndDelegate(EndCallBack);
tesselator.callVertex = new Tesselate.Tesselator.CallVertexDelegate(VertexCallBack);
tesselator.callCombine = new Tesselate.Tesselator.CallCombineDelegate(CombineCallBack);
tesselator.windingRule = windingRule;
if (setEdgeFlag)
{
tesselator.callEdgeFlag += new Tesselate.Tesselator.CallEdgeFlagDelegate(EdgeFlagCallBack);
}
}
/// <summary>
/// connect to actual Tesselator
/// </summary>
/// <param name="tesselator"></param>
/// <param name="setEdgeFlag"></param>
public void Connect(Tesselator tesselator, bool setEdgeFlag)
{
NeedEdgeFlag = setEdgeFlag;
tesselator.SetListener(this);
//tesselator.callBegin = OnBegin;
//tesselator.callEnd = OnEnd;
//tesselator.callVertex = OnVertex;
//tesselator.callCombine = OnCombine;
//if (setEdgeFlag)
//{
// tesselator.callEdgeFlag = OnEdgeFlag;
//}
}
public void Connect(List <Vertex> vertextList, Tesselate.Tesselator tesselator, Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
this.inputVertextList = vertextList;
tesselator.callBegin += new Tesselate.Tesselator.CallBeginDelegate(BeginCallBack);
tesselator.callEnd += new Tesselate.Tesselator.CallEndDelegate(EndCallBack);
tesselator.callVertex += new Tesselate.Tesselator.CallVertexDelegate(VertexCallBack);
tesselator.callCombine += new Tesselate.Tesselator.CallCombineDelegate(CombineCallBack);
tesselator.windingRule = windingRule;
if (setEdgeFlag)
{
tesselator.callEdgeFlag += new Tesselate.Tesselator.CallEdgeFlagDelegate(EdgeFlagCallBack);
}
}
FrameBuffer _currentFrameBuffer;//default = null, system provide frame buffer
internal CanvasGL2d(int canvasW, int canvasH)
{
this.canvasW = canvasW;
this.canvasH = canvasH;
////setup viewport size
int max = Math.Max(canvasW, canvasH);
////square viewport
orthoView = MyMat4.ortho(0, max, 0, max, 0, 1);
//-----------------------------------------------------------------------
shaderRes = new CanvasToShaderSharedResource();
shaderRes.OrthoView = orthoView;
//-----------------------------------------------------------------------
basicFillShader = new BasicFillShader(shaderRes);
smoothLineShader = new SmoothLineShader(shaderRes);
rectFillShader = new RectFillShader(shaderRes);
gdiImgTextureShader = new GdiImageTextureShader(shaderRes);
gdiImgTextureWithWhiteTransparentShader = new GdiImageTextureWithWhiteTransparentShader(shaderRes);
blurShader = new BlurShader(shaderRes);
glesTextureShader = new OpenGLESTextureShader(shaderRes);
invertAlphaFragmentShader = new InvertAlphaLineSmoothShader(shaderRes); //used with stencil ***
// tessListener.Connect(tess,
//Tesselate.Tesselator.WindingRuleType.Odd, true);
conv3x3TextureShader = new Conv3x3TextureShader(shaderRes);
msdfShader = new DrawingGL.MultiChannelSdf(shaderRes);
msdfSubPixelRenderingShader = new DrawingGL.MultiChannelSubPixelRenderingSdf(shaderRes);
sdfShader = new DrawingGL.SingleChannelSdf(shaderRes);
//----
Tesselator tess = new Tesselator();
tess.WindingRule = Tesselator.WindingRuleType.Odd;
tessTool = new TessTool(tess);
//-----------------------------------------------------------------------
//GL.Enable(EnableCap.CullFace);
//GL.FrontFace(FrontFaceDirection.Cw);
//GL.CullFace(CullFaceMode.Back);
GL.Enable(EnableCap.Blend);
GL.BlendFunc(BlendingFactorSrc.SrcAlpha, BlendingFactorDest.OneMinusSrcAlpha);
GL.ClearColor(1, 1, 1, 1);
//-------------------------------------------------------------------------------
GL.Viewport(0, 0, canvasW, canvasH);
}
void RunTest(int instructionStreamIndex, Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
Tesselate.Tesselator tesselator = new Tesselate.Tesselator();
tesselator.callBegin += new Tesselate.Tesselator.CallBeginDelegate(BeginCallBack);
tesselator.callEnd += new Tesselate.Tesselator.CallEndDelegate(EndCallBack);
tesselator.callVertex += new Tesselate.Tesselator.CallVertexDelegate(VertexCallBack);
tesselator.callCombine += new Tesselate.Tesselator.CallCombineDelegate(CombineCallBack);
tesselator.windingRule = windingRule;
if (setEdgeFlag)
{
tesselator.callEdgeFlag += new Tesselate.Tesselator.CallEdgeFlagDelegate(EdgeFlagCallBack);
}
ParseStreamForTesselator(tesselator, instructionStreamIndex);
}
private void RunTest(int instructionStreamIndex, Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
Tesselate.Tesselator tesselator = new Tesselate.Tesselator();
tesselator.callBegin += BeginCallBack;
tesselator.callEnd += EndCallBack;
tesselator.callVertex += VertexCallBack;
tesselator.callCombine += CombineCallBack;
tesselator.windingRule = windingRule;
if (setEdgeFlag)
{
tesselator.callEdgeFlag = EdgeFlagCallBack;
}
ParseStreamForTesselator(tesselator, instructionStreamIndex);
}
public void BeginCallBack(Tesselator.TriangleListType type)
{
switch (type)
{
case Tesselator.TriangleListType.Triangles:
Gl.glBegin(Gl.GL_TRIANGLES);
break;
case Tesselator.TriangleListType.TriangleFan:
Gl.glBegin(Gl.GL_TRIANGLE_FAN);
break;
case Tesselator.TriangleListType.TriangleStrip:
Gl.glBegin(Gl.GL_TRIANGLE_STRIP);
break;
}
}
public void BeginCallBack(Tesselator.TriangleListType type)
{
switch (type)
{
case Tesselator.TriangleListType.Triangles:
GL.Begin(BeginMode.Triangles);
break;
case Tesselator.TriangleListType.TriangleFan:
GL.Begin(BeginMode.TriangleFan);
break;
case Tesselator.TriangleListType.TriangleStrip:
GL.Begin(BeginMode.TriangleStrip);
break;
}
}
public void ParseStreamForTesselator(Tesselate.Tesselator tesselator, int instructionStreamIndex)
{
m_VertexList.Clear();
m_CurrentOutput = 0;
string[] instructionStream = m_InsructionStream[instructionStreamIndex];
for (int curInstruction = 0; curInstruction < instructionStream.Length; curInstruction++)
{
switch (instructionStream[curInstruction])
{
case "BP":
tesselator.BeginPolygon();
break;
case "BC":
tesselator.BeginContour();
break;
case "V":
double x = Convert.ToDouble(instructionStream[curInstruction + 1]);
double y = Convert.ToDouble(instructionStream[curInstruction + 2]);
curInstruction += 2;
double[] coords = new double[3];
coords[0] = x;
coords[1] = y;
tesselator.AddVertex(coords, m_VertexList.Count);
m_VertexList.Add(new Vertex(x, y));
break;
case "EC":
tesselator.EndContour();
break;
case "EP":
tesselator.EndPolygon();
break;
default:
throw new Exception();
}
}
}
void RunTest(int instructionStreamIndex, Tesselator.WindingRuleType windingRule, bool setEdgeFlag)
{
Tesselate.Tesselator tesselator = new Tesselate.Tesselator();
TessTest.TessListener listner = new TessTest.TessListener();
tesselator.WindingRule = windingRule;
listner.Connect(tesselator, setEdgeFlag);
//tesselator.callBegin += new Tesselate.Tesselator.CallBeginDelegate(BeginCallBack);
//tesselator.callEnd += new Tesselate.Tesselator.CallEndDelegate(EndCallBack);
//tesselator.callVertex += new Tesselate.Tesselator.CallVertexDelegate(VertexCallBack);
//tesselator.callCombine += new Tesselate.Tesselator.CallCombineDelegate(CombineCallBack);
//tesselator.windingRule = windingRule;
//if (setEdgeFlag)
//{
// tesselator.callEdgeFlag += new Tesselate.Tesselator.CallEdgeFlagDelegate(EdgeFlagCallBack);
//}
//ParseStreamForTesselator(tesselator, instructionStreamIndex);
}
public void BeginCallBack(Tesselator.TriangleListType type)
{
//Console.WriteLine("begin: " + type.ToString());
//Assert.IsTrue(GetNextOutputAsString() == "B");
//switch (type)
//{
// case Tesselator.TriangleListType.Triangles:
// Assert.IsTrue(GetNextOutputAsString() == "TRI");
// break;
// case Tesselator.TriangleListType.TriangleFan:
// Assert.IsTrue(GetNextOutputAsString() == "FAN");
// break;
// case Tesselator.TriangleListType.TriangleStrip:
// Assert.IsTrue(GetNextOutputAsString() == "STRIP");
// break;
// default:
// throw new Exception("unknown TriangleListType '" + type.ToString() + "'.");
//}
}
static void RenderFan(Tesselator tess, HalfEdge e, int size)
{
/* Render as many CCW triangles as possible in a fan starting from
* edge "e". The fan *should* contain exactly "size" triangles
* (otherwise we've goofed up somewhere).
*/
tess.CallBegin(Tesselator.TriangleListType.TriangleFan);
tess.CallVertex(e.originVertex.clientIndex);
tess.CallVertex(e.directionVertex.clientIndex);
while (!e.leftFace.Marked())
{
e.leftFace.marked = true;
--size;
e = e.nextEdgeCCWAroundOrigin;
tess.CallVertex(e.directionVertex.clientIndex);
}
if (size != 0)
{
throw new Exception();
}
tess.CallEnd();
}
static bool CheckForLeftSplice(Tesselator tess, ActiveRegion regUp)
/*
* Check the upper and lower edge of "regUp", to make sure that the
* eUp.Dst is above eLo, or eLo.Dst is below eUp (depending on which
* destination is rightmost).
*
* Theoretically, this should always be true. However, splitting an edge
* into two pieces can change the results of previous tests. For example,
* suppose at one point we checked eUp and eLo, and decided that eUp.Dst
* is barely above eLo. Then later, we split eLo into two edges (eg. from
* a splice operation like this one). This can change the result of
* the test so that now eUp.Dst is incident to eLo, or barely below it.
* We must correct this condition to maintain the dictionary invariants
* (otherwise new edges might get inserted in the wrong place in the
* dictionary, and bad stuff will happen).
*
* We fix the problem by just splicing the offending vertex into the
* other edge.
*/
{
ActiveRegion regLo = RegionBelow(regUp);
HalfEdge eUp = regUp.upperHalfEdge;
HalfEdge eLo = regLo.upperHalfEdge;
HalfEdge e;
if (eUp.directionVertex.VertEq(eLo.directionVertex))
{
throw new Exception();
}
if (eUp.directionVertex.VertLeq(eLo.directionVertex))
{
if (ContourVertex.EdgeSign(eUp.directionVertex, eLo.directionVertex, eUp.originVertex) < 0)
{
return false;
}
/* eLo.Dst is above eUp, so splice eLo.Dst into eUp */
regUp.RegionAbove().dirty = regUp.dirty = true;
e = Mesh.meshSplitEdge(eUp);
Mesh.meshSplice(eLo.otherHalfOfThisEdge, e);
e.leftFace.isInterior = regUp.inside;
}
else
{
if (ContourVertex.EdgeSign(eLo.directionVertex, eUp.directionVertex, eLo.originVertex) > 0) return false;
/* eUp.Dst is below eLo, so splice eUp.Dst into eLo */
regUp.dirty = regLo.dirty = true;
e = Mesh.meshSplitEdge(eLo);
Mesh.meshSplice(eUp.nextEdgeCCWAroundLeftFace, eLo.otherHalfOfThisEdge);
e.rightFace.isInterior = regUp.inside;
}
return true;
}
static void GetIntersectData(Tesselator tess, ContourVertex isect,
ContourVertex orgUp, ContourVertex dstUp,
ContourVertex orgLo, ContourVertex dstLo)
/*
* We've computed a new intersection point, now we need a "data" pointer
* from the user so that we can refer to this new vertex in the
* rendering callbacks.
*/
{
int[] data4 = new int[4];
double[] weights4 = new double[4];
data4[0] = orgUp.clientIndex;
data4[1] = dstUp.clientIndex;
data4[2] = orgLo.clientIndex;
data4[3] = dstLo.clientIndex;
isect.C_0 = isect.C_1 = isect.C_2 = 0;
VertexWeights(isect, orgUp, dstUp, out weights4[0], out weights4[1]);
VertexWeights(isect, orgLo, dstLo, out weights4[2], out weights4[3]);
CallCombine(tess, isect, data4, weights4, true);
}
static bool CheckForRightSplice(Tesselator tess, ActiveRegion regUp)
/*
* Check the upper and lower edge of "regUp", to make sure that the
* eUp.Org is above eLo, or eLo.Org is below eUp (depending on which
* origin is leftmost).
*
* The main purpose is to splice right-going edges with the same
* dest vertex and nearly identical slopes (ie. we can't distinguish
* the slopes numerically). However the splicing can also help us
* to recover from numerical errors. For example, suppose at one
* point we checked eUp and eLo, and decided that eUp.Org is barely
* above eLo. Then later, we split eLo into two edges (eg. from
* a splice operation like this one). This can change the result of
* our test so that now eUp.Org is incident to eLo, or barely below it.
* We must correct this condition to maintain the dictionary invariants.
*
* One possibility is to check these edges for intersection again
* (ie. CheckForIntersect). This is what we do if possible. However
* CheckForIntersect requires that tess.currentSweepVertex lies between eUp and eLo,
* so that it has something to fall back on when the intersection
* calculation gives us an unusable answer. So, for those cases where
* we can't check for intersection, this routine fixes the problem
* by just splicing the offending vertex into the other edge.
* This is a guaranteed solution, no matter how degenerate things get.
* Basically this is a combinatorial solution to a numerical problem.
*/
{
ActiveRegion regLo = RegionBelow(regUp);
HalfEdge eUp = regUp.upperHalfEdge;
HalfEdge eLo = regLo.upperHalfEdge;
if (eUp.originVertex.VertLeq(eLo.originVertex))
{
if (ContourVertex.EdgeSign(eLo.directionVertex, eUp.originVertex, eLo.originVertex) > 0)
{
return false;
}
/* eUp.Org appears to be below eLo */
if (!eUp.originVertex.VertEq(eLo.originVertex))
{
/* Splice eUp.Org into eLo */
Mesh.meshSplitEdge(eLo.otherHalfOfThisEdge);
Mesh.meshSplice(eUp, eLo.Oprev);
regUp.dirty = regLo.dirty = true;
}
else if (eUp.originVertex != eLo.originVertex)
{
/* merge the two vertices, discarding eUp.Org */
tess.vertexPriorityQue.Delete(eUp.originVertex.priorityQueueHandle);
//pqDelete(tess.pq, eUp.Org.pqHandle); /* __gl_pqSortDelete */
SpliceMergeVertices(tess, eLo.Oprev, eUp);
}
}
else
{
if (ContourVertex.EdgeSign(eUp.directionVertex, eLo.originVertex, eUp.originVertex) < 0)
{
return false;
}
/* eLo.Org appears to be above eUp, so splice eLo.Org into eUp */
regUp.RegionAbove().dirty = regUp.dirty = true;
Mesh.meshSplitEdge(eUp.otherHalfOfThisEdge);
Mesh.meshSplice(eLo.Oprev, eUp);
}
return true;
}
static void ConnectLeftDegenerate(Tesselator tess, ActiveRegion regUp, ContourVertex vEvent)
/*
* The currentSweepVertex vertex lies exactly on an already-processed edge or vertex.
* Adding the new vertex involves splicing it into the already-processed
* part of the mesh.
*/
{
HalfEdge e, eTopLeft, eTopRight, eLast;
ActiveRegion reg;
e = regUp.upperHalfEdge;
if (e.originVertex.VertEq(vEvent))
{
/* e.Org is an unprocessed vertex - just combine them, and wait
* for e.Org to be pulled from the queue
*/
SpliceMergeVertices(tess, e, vEvent.edgeThisIsOriginOf);
return;
}
if (!e.directionVertex.VertEq(vEvent))
{
/* General case -- splice vEvent into edge e which passes through it */
Mesh.meshSplitEdge(e.otherHalfOfThisEdge);
if (regUp.fixUpperEdge)
{
/* This edge was fixable -- delete unused portion of original edge */
Mesh.DeleteHalfEdge(e.nextEdgeCCWAroundOrigin);
regUp.fixUpperEdge = false;
}
Mesh.meshSplice(vEvent.edgeThisIsOriginOf, e);
SweepEvent(tess, vEvent); /* recurse */
return;
}
/* vEvent coincides with e.Dst, which has already been processed.
* Splice in the additional right-going edges.
*/
regUp = TopRightRegion(regUp);
reg = RegionBelow(regUp);
eTopRight = reg.upperHalfEdge.otherHalfOfThisEdge;
eTopLeft = eLast = eTopRight.nextEdgeCCWAroundOrigin;
if (reg.fixUpperEdge)
{
/* Here e.Dst has only a single fixable edge going right.
* We can delete it since now we have some real right-going edges.
*/
if (eTopLeft == eTopRight)
{
throw new Exception(); /* there are some left edges too */
}
DeleteRegion(reg);
Mesh.DeleteHalfEdge(eTopRight);
eTopRight = eTopLeft.Oprev;
}
Mesh.meshSplice(vEvent.edgeThisIsOriginOf, eTopRight);
if (!eTopLeft.EdgeGoesLeft())
{
/* e.Dst had no left-going edges -- indicate this to AddRightEdges() */
eTopLeft = null;
}
AddRightEdges(tess, regUp, eTopRight.nextEdgeCCWAroundOrigin, eLast, eTopLeft, true);
}
static void ConnectRightVertex(Tesselator tess, ActiveRegion regUp, HalfEdge eBottomLeft)
/*
* Purpose: connect a "right" vertex vEvent (one where all edges go left)
* to the unprocessed portion of the mesh. Since there are no right-going
* edges, two regions (one above vEvent and one below) are being merged
* into one. "regUp" is the upper of these two regions.
*
* There are two reasons for doing this (adding a right-going edge):
* - if the two regions being merged are "inside", we must add an edge
* to keep them separated (the combined region would not be monotone).
* - in any case, we must leave some record of vEvent in the dictionary,
* so that we can merge vEvent with features that we have not seen yet.
* For example, maybe there is a vertical edge which passes just to
* the right of vEvent; we would like to splice vEvent into this edge.
*
* However, we don't want to connect vEvent to just any vertex. We don''t
* want the new edge to cross any other edges; otherwise we will create
* intersection vertices even when the input data had no self-intersections.
* (This is a bad thing; if the user's input data has no intersections,
* we don't want to generate any false intersections ourselves.)
*
* Our eventual goal is to connect vEvent to the leftmost unprocessed
* vertex of the combined region (the union of regUp and regLo).
* But because of unseen vertices with all right-going edges, and also
* new vertices which may be created by edge intersections, we don''t
* know where that leftmost unprocessed vertex is. In the meantime, we
* connect vEvent to the closest vertex of either chain, and mark the region
* as "fixUpperEdge". This flag says to delete and reconnect this edge
* to the next processed vertex on the boundary of the combined region.
* Quite possibly the vertex we connected to will turn out to be the
* closest one, in which case we won''t need to make any changes.
*/
{
HalfEdge eNew;
HalfEdge eTopLeft = eBottomLeft.nextEdgeCCWAroundOrigin;
ActiveRegion regLo = RegionBelow(regUp);
HalfEdge eUp = regUp.upperHalfEdge;
HalfEdge eLo = regLo.upperHalfEdge;
bool degenerate = false;
if (eUp.directionVertex != eLo.directionVertex)
{
CheckForIntersect(tess, regUp);
}
/* Possible new degeneracies: upper or lower edge of regUp may pass
* through vEvent, or may coincide with new intersection vertex
*/
if (eUp.originVertex.VertEq(tess.currentSweepVertex))
{
Mesh.meshSplice(eTopLeft.Oprev, eUp);
regUp = TopLeftRegion(regUp);
eTopLeft = RegionBelow(regUp).upperHalfEdge;
FinishLeftRegions(tess, RegionBelow(regUp), regLo);
degenerate = true;
}
if (eLo.originVertex.VertEq(tess.currentSweepVertex))
{
Mesh.meshSplice(eBottomLeft, eLo.Oprev);
eBottomLeft = FinishLeftRegions(tess, regLo, null);
degenerate = true;
}
if (degenerate)
{
AddRightEdges(tess, regUp, eBottomLeft.nextEdgeCCWAroundOrigin, eTopLeft, eTopLeft, true);
return;
}
/* Non-degenerate situation -- need to add a temporary, fixable edge.
* Connect to the closer of eLo.Org, eUp.Org.
*/
if (eLo.originVertex.VertLeq(eUp.originVertex))
{
eNew = eLo.Oprev;
}
else
{
eNew = eUp;
}
eNew = Mesh.meshConnect(eBottomLeft.Lprev, eNew);
/* Prevent cleanup, otherwise eNew might disappear before we've even
* had a chance to mark it as a temporary edge.
*/
AddRightEdges(tess, regUp, eNew, eNew.nextEdgeCCWAroundOrigin, eNew.nextEdgeCCWAroundOrigin, false);
eNew.otherHalfOfThisEdge.regionThisIsUpperEdgeOf.fixUpperEdge = true;
WalkDirtyRegions(tess, regUp);
}
public PolygonTesselator()
{
m_Tesselator = new Tesselate.Tesselator();
m_Tesselator.callBegin += new Tesselator.CallBeginDelegate(BeginCallBack);
m_Tesselator.callEnd += new Tesselator.CallEndDelegate(EndCallBack);
m_Tesselator.callVertex += new Tesselator.CallVertexDelegate(VertexCallBack);
//m_Tesselator.callEdgeFlag += new Tesselator.CallEdgeFlagDelegate(EdgeFlagCallBack);
m_Tesselator.callCombine += new Tesselator.CallCombineDelegate(CombineCallBack);
}
static void DoneEdgeDict(Tesselator tess)
{
ActiveRegion reg;
int fixedEdges = 0;
/* __GL_DICTLISTKEY */
/* __GL_DICTLISTMIN */
while ((reg = tess.edgeDictionary.GetMinNode().Key) != null)
{
/*
* At the end of all processing, the dictionary should contain
* only the two sentinel edges, plus at most one "fixable" edge
* created by ConnectRightVertex().
*/
if (!reg.sentinel)
{
if (!reg.fixUpperEdge)
{
throw new System.Exception();
}
if (++fixedEdges != 1)
{
throw new System.Exception();
}
}
if (reg.windingNumber != 0)
{
throw new Exception();
}
DeleteRegion(reg);
}
tess.edgeDictionary = null;
}
static void AddSentinel(Tesselator tess, double t)
/*
* We add two sentinel edges above and below all other edges,
* to avoid special cases at the top and bottom.
*/
{
HalfEdge halfEdge;
ActiveRegion activeRedion = new ActiveRegion();
halfEdge = tess.mesh.MakeEdge();
halfEdge.originVertex.x = SENTINEL_COORD;
halfEdge.originVertex.y = t;
halfEdge.directionVertex.x = -SENTINEL_COORD;
halfEdge.directionVertex.y = t;
tess.currentSweepVertex = halfEdge.directionVertex; /* initialize it */
activeRedion.upperHalfEdge = halfEdge;
activeRedion.windingNumber = 0;
activeRedion.inside = false;
activeRedion.fixUpperEdge = false;
activeRedion.sentinel = true;
activeRedion.dirty = false;
activeRedion.upperHalfEdgeDictNode = tess.edgeDictionary.Insert(activeRedion); /* __gl_dictListInsertBefore */
}
static void SweepEvent(Tesselator tess, ContourVertex vEvent)
/*
* Does everything necessary when the sweep line crosses a vertex.
* Updates the mesh and the edge dictionary.
*/
{
ActiveRegion regUp, reg;
HalfEdge e, eTopLeft, eBottomLeft;
tess.currentSweepVertex = vEvent; /* for access in EdgeLeq() */
/* Check if this vertex is the right endpoint of an edge that is
* already in the dictionary. In this case we don't need to waste
* time searching for the location to insert new edges.
*/
e = vEvent.edgeThisIsOriginOf;
while (e.regionThisIsUpperEdgeOf == null)
{
e = e.nextEdgeCCWAroundOrigin;
if (e == vEvent.edgeThisIsOriginOf)
{
/* All edges go right -- not incident to any processed edges */
ConnectLeftVertex(tess, vEvent);
return;
}
}
/* Processing consists of two phases: first we "finish" all the
* active regions where both the upper and lower edges terminate
* at vEvent (ie. vEvent is closing off these regions).
* We mark these faces "inside" or "outside" the polygon according
* to their winding number, and delete the edges from the dictionary.
* This takes care of all the left-going edges from vEvent.
*/
regUp = TopLeftRegion(e.regionThisIsUpperEdgeOf);
reg = RegionBelow(regUp);
eTopLeft = reg.upperHalfEdge;
eBottomLeft = FinishLeftRegions(tess, reg, null);
/* Next we process all the right-going edges from vEvent. This
* involves adding the edges to the dictionary, and creating the
* associated "active regions" which record information about the
* regions between adjacent dictionary edges.
*/
if (eBottomLeft.nextEdgeCCWAroundOrigin == eTopLeft)
{
/* No right-going edges -- add a temporary "fixable" edge */
ConnectRightVertex(tess, regUp, eBottomLeft);
}
else
{
AddRightEdges(tess, regUp, eBottomLeft.nextEdgeCCWAroundOrigin, eTopLeft, eTopLeft, true);
}
}
public void Connect(Tesselate.Tesselator tesselator, bool setEdgeFlag)
{
NeedEdgeFlag = setEdgeFlag;
tesselator.SetListener(this);
}
static ActiveRegion AddRegionBelow(Tesselator tess,
ActiveRegion regAbove,
HalfEdge eNewUp)
/*
* Add a new active region to the sweep line, *somewhere* below "regAbove"
* (according to where the new edge belongs in the sweep-line dictionary).
* The upper edge of the new region will be "eNewUp".
* Winding number and "inside" flag are not updated.
*/
{
ActiveRegion regNew = new ActiveRegion();
regNew.upperHalfEdge = eNewUp;
/* __gl_dictListInsertBefore */
regNew.upperHalfEdgeDictNode = tess.edgeDictionary.InsertBefore(regAbove.upperHalfEdgeDictNode, regNew);
regNew.fixUpperEdge = false;
regNew.sentinel = false;
regNew.dirty = false;
eNewUp.regionThisIsUpperEdgeOf = regNew;
return regNew;
}
static void DonePriorityQ(Tesselator tess)
{
tess.vertexPriorityQue = null; /* __gl_pqSortDeletePriorityQ */
}
static void InitPriorityQue(Tesselator tess)
/*
* Insert all vertices into the priority queue which determines the
* order in which vertices cross the sweep line.
*/
{
MiniCollection.MaxFirstList<ContourVertex> priorityQue = tess.vertexPriorityQue = new MiniCollection.MaxFirstList<ContourVertex>();
ContourVertex vertexHead = tess.mesh.vertexHead;
for (ContourVertex curVertex = vertexHead.nextVertex; curVertex != vertexHead; curVertex = curVertex.nextVertex)
{
priorityQue.Add(out curVertex.priorityQueueHandle, curVertex);
}
}
static void RemoveDegenerateEdges(Tesselator tess)
{
// Remove zero-length edges, and contours with fewer than 3 vertices.
HalfEdge edgeHead = tess.mesh.halfEdgeHead;
HalfEdge nextHalfEdge;
for (HalfEdge currentEdge = edgeHead.nextHalfEdge; currentEdge != edgeHead; currentEdge = nextHalfEdge)
{
nextHalfEdge = currentEdge.nextHalfEdge;
HalfEdge nextEdgeCCWAroundLeftFace = currentEdge.nextEdgeCCWAroundLeftFace;
if (currentEdge.originVertex.VertEq(currentEdge.directionVertex)
&& currentEdge.nextEdgeCCWAroundLeftFace.nextEdgeCCWAroundLeftFace != currentEdge)
{
// Zero-length edge, contour has at least 3 edges
SpliceMergeVertices(tess, nextEdgeCCWAroundLeftFace, currentEdge); /* deletes e.Org */
Mesh.DeleteHalfEdge(currentEdge); /* e is a self-loop */
currentEdge = nextEdgeCCWAroundLeftFace;
nextEdgeCCWAroundLeftFace = currentEdge.nextEdgeCCWAroundLeftFace;
}
if (nextEdgeCCWAroundLeftFace.nextEdgeCCWAroundLeftFace == currentEdge)
{
// Degenerate contour (one or two edges)
if (nextEdgeCCWAroundLeftFace != currentEdge)
{
if (nextEdgeCCWAroundLeftFace == nextHalfEdge || nextEdgeCCWAroundLeftFace == nextHalfEdge.otherHalfOfThisEdge)
{
nextHalfEdge = nextHalfEdge.nextHalfEdge;
}
Mesh.DeleteHalfEdge(nextEdgeCCWAroundLeftFace);
}
if (currentEdge == nextHalfEdge || currentEdge == nextHalfEdge.otherHalfOfThisEdge) { nextHalfEdge = nextHalfEdge.nextHalfEdge; }
Mesh.DeleteHalfEdge(currentEdge);
}
}
}
/* __gl_computeInterior( tess ) computes the planar arrangement specified
* by the given contours, and further subdivides this arrangement
* into regions. Each region is marked "inside" if it belongs
* to the polygon, according to the rule given by tess.windingRule.
* Each interior region is guaranteed be monotone.
*/
public static int ComputeInterior(Tesselator tess)
/*
* __gl_computeInterior( tess ) computes the planar arrangement specified
* by the given contours, and further subdivides this arrangement
* into regions. Each region is marked "inside" if it belongs
* to the polygon, according to the rule given by tess.windingRule.
* Each interior region is guaranteed be monotone.
*/
{
ContourVertex vertex, vertexNext;
/* Each vertex defines an currentSweepVertex for our sweep line. Start by inserting
* all the vertices in a priority queue. Events are processed in
* lexicographic order, ie.
*
* e1 < e2 iff e1.x < e2.x || (e1.x == e2.x && e1.y < e2.y)
*/
RemoveDegenerateEdges(tess);
InitPriorityQue(tess);
InitEdgeDict(tess);
while (!tess.vertexPriorityQue.IsEmpty)
{
vertex = tess.vertexPriorityQue.DeleteMin();
for (;;)
{
if (!tess.vertexPriorityQue.IsEmpty)
{
vertexNext = tess.vertexPriorityQue.FindMin(); /* __gl_pqSortMinimum */
}
else
{
vertexNext = null;
}
if (vertexNext == null || !vertexNext.Equal2D(vertex)) break;
/* Merge together all vertices at exactly the same location.
* This is more efficient than processing them one at a time,
* simplifies the code (see ConnectLeftDegenerate), and is also
* important for correct handling of certain degenerate cases.
* For example, suppose there are two identical edges A and B
* that belong to different contours (so without this code they would
* be processed by separate sweep events). Suppose another edge C
* crosses A and B from above. When A is processed, we split it
* at its intersection point with C. However this also splits C,
* so when we insert B we may compute a slightly different
* intersection point. This might leave two edges with a small
* gap between them. This kind of error is especially obvious
* when using boundary extraction (GLU_TESS_BOUNDARY_ONLY).
*/
vertexNext = tess.vertexPriorityQue.DeleteMin(); /* __gl_pqSortExtractMin*/
SpliceMergeVertices(tess, vertex.edgeThisIsOriginOf, vertexNext.edgeThisIsOriginOf);
}
SweepEvent(tess, vertex);
}
/* Set tess.currentSweepVertex for debugging purposes */
/* __GL_DICTLISTKEY */
/* __GL_DICTLISTMIN */
tess.currentSweepVertex = tess.edgeDictionary.GetMinNode().Key.upperHalfEdge.originVertex;
DoneEdgeDict(tess);
DonePriorityQ(tess);
if (!RemoveDegenerateFaces(tess.mesh))
{
return 0;
}
tess.mesh.CheckMesh();
return 1;
}
public TessTool(Tesselate.Tesselator tess)
{
this.tess = tess;
this.tessListener = new TessListener2();
tessListener.Connect(tess, true);
}
static bool CheckForIntersect(Tesselator tess, ActiveRegion regUp)
/*
* Check the upper and lower edges of the given region to see if
* they intersect. If so, create the intersection and add it
* to the data structures.
*
* Returns true if adding the new intersection resulted in a recursive
* call to AddRightEdges(); in this case all "dirty" regions have been
* checked for intersections, and possibly regUp has been deleted.
*/
{
ActiveRegion regLo = RegionBelow(regUp);
HalfEdge eUp = regUp.upperHalfEdge;
HalfEdge eLo = regLo.upperHalfEdge;
ContourVertex orgUp = eUp.originVertex;
ContourVertex orgLo = eLo.originVertex;
ContourVertex dstUp = eUp.directionVertex;
ContourVertex dstLo = eLo.directionVertex;
double tMinUp, tMaxLo;
ContourVertex isect = new ContourVertex();
ContourVertex orgMin;
HalfEdge e;
if (dstLo.VertEq(dstUp))
{
throw new Exception();
}
if (ContourVertex.EdgeSign(dstUp, tess.currentSweepVertex, orgUp) > 0)
{
throw new Exception();
}
if (ContourVertex.EdgeSign(dstLo, tess.currentSweepVertex, orgLo) < 0)
{
throw new Exception();
}
if (orgUp == tess.currentSweepVertex || orgLo == tess.currentSweepVertex)
{
throw new Exception();
}
if (regUp.fixUpperEdge || regLo.fixUpperEdge)
{
throw new Exception();
}
if (orgUp == orgLo)
{
return false; /* right endpoints are the same */
}
tMinUp = Math.Min(orgUp.y, dstUp.y);
tMaxLo = Math.Max(orgLo.y, dstLo.y);
if (tMinUp > tMaxLo)
{
return false; /* t ranges do not overlap */
}
if (orgUp.VertLeq(orgLo))
{
if (ContourVertex.EdgeSign(dstLo, orgUp, orgLo) > 0)
{
return false;
}
}
else
{
if (ContourVertex.EdgeSign(dstUp, orgLo, orgUp) < 0)
{
return false;
}
}
EdgeIntersect(dstUp, orgUp, dstLo, orgLo, ref isect);
// The following properties are guaranteed:
if (!(Math.Min(orgUp.y, dstUp.y) <= isect.y))
{
throw new System.Exception();
}
if (!(isect.y <= Math.Max(orgLo.y, dstLo.y)))
{
throw new System.Exception();
}
if (!(Math.Min(dstLo.x, dstUp.x) <= isect.x))
{
throw new System.Exception();
}
if (!(isect.x <= Math.Max(orgLo.x, orgUp.x)))
{
throw new System.Exception();
}
if (isect.VertLeq(tess.currentSweepVertex))
{
/* The intersection point lies slightly to the left of the sweep line,
* so move it until it''s slightly to the right of the sweep line.
* (If we had perfect numerical precision, this would never happen
* in the first place). The easiest and safest thing to do is
* replace the intersection by tess.currentSweepVertex.
*/
isect.x = tess.currentSweepVertex.x;
isect.y = tess.currentSweepVertex.y;
}
/* Similarly, if the computed intersection lies to the right of the
* rightmost origin (which should rarely happen), it can cause
* unbelievable inefficiency on sufficiently degenerate inputs.
* (If you have the test program, try running test54.d with the
* "X zoom" option turned on).
*/
orgMin = orgUp.VertLeq(orgLo) ? orgUp : orgLo;
if (orgMin.VertLeq(isect))
{
isect.x = orgMin.x;
isect.y = orgMin.y;
}
if (isect.VertEq(orgUp) || isect.VertEq(orgLo))
{
/* Easy case -- intersection at one of the right endpoints */
CheckForRightSplice(tess, regUp);
return false;
}
if ((!dstUp.VertEq(tess.currentSweepVertex)
&& ContourVertex.EdgeSign(dstUp, tess.currentSweepVertex, isect) >= 0)
|| (!dstLo.VertEq(tess.currentSweepVertex)
&& ContourVertex.EdgeSign(dstLo, tess.currentSweepVertex, isect) <= 0))
{
/* Very unusual -- the new upper or lower edge would pass on the
* wrong side of the sweep currentSweepVertex, or through it. This can happen
* due to very small numerical errors in the intersection calculation.
*/
if (dstLo == tess.currentSweepVertex)
{
/* Splice dstLo into eUp, and process the new region(s) */
Mesh.meshSplitEdge(eUp.otherHalfOfThisEdge);
Mesh.meshSplice(eLo.otherHalfOfThisEdge, eUp);
regUp = TopLeftRegion(regUp);
eUp = RegionBelow(regUp).upperHalfEdge;
FinishLeftRegions(tess, RegionBelow(regUp), regLo);
AddRightEdges(tess, regUp, eUp.Oprev, eUp, eUp, true);
return true;
}
if (dstUp == tess.currentSweepVertex)
{
/* Splice dstUp into eLo, and process the new region(s) */
Mesh.meshSplitEdge(eLo.otherHalfOfThisEdge);
Mesh.meshSplice(eUp.nextEdgeCCWAroundLeftFace, eLo.Oprev);
regLo = regUp;
regUp = TopRightRegion(regUp);
e = RegionBelow(regUp).upperHalfEdge.Rprev;
regLo.upperHalfEdge = eLo.Oprev;
eLo = FinishLeftRegions(tess, regLo, null);
AddRightEdges(tess, regUp, eLo.nextEdgeCCWAroundOrigin, eUp.Rprev, e, true);
return true;
}
/* Special case: called from ConnectRightVertex. If either
* edge passes on the wrong side of tess.currentSweepVertex, split it
* (and wait for ConnectRightVertex to splice it appropriately).
*/
if (ContourVertex.EdgeSign(dstUp, tess.currentSweepVertex, isect) >= 0)
{
regUp.RegionAbove().dirty = regUp.dirty = true;
Mesh.meshSplitEdge(eUp.otherHalfOfThisEdge);
eUp.originVertex.x = tess.currentSweepVertex.x;
eUp.originVertex.y = tess.currentSweepVertex.y;
}
if (ContourVertex.EdgeSign(dstLo, tess.currentSweepVertex, isect) <= 0)
{
regUp.dirty = regLo.dirty = true;
Mesh.meshSplitEdge(eLo.otherHalfOfThisEdge);
eLo.originVertex.x = tess.currentSweepVertex.x;
eLo.originVertex.y = tess.currentSweepVertex.y;
}
/* leave the rest for ConnectRightVertex */
return false;
}
/* General case -- split both edges, splice into new vertex.
* When we do the splice operation, the order of the arguments is
* arbitrary as far as correctness goes. However, when the operation
* creates a new face, the work done is proportional to the size of
* the new face. We expect the faces in the processed part of
* the mesh (ie. eUp.Lface) to be smaller than the faces in the
* unprocessed original contours (which will be eLo.Oprev.Lface).
*/
Mesh.meshSplitEdge(eUp.otherHalfOfThisEdge);
Mesh.meshSplitEdge(eLo.otherHalfOfThisEdge);
Mesh.meshSplice(eLo.Oprev, eUp);
eUp.originVertex.x = isect.x;
eUp.originVertex.y = isect.y;
tess.vertexPriorityQue.Add(out eUp.originVertex.priorityQueueHandle, eUp.originVertex); /* __gl_pqSortInsert */
GetIntersectData(tess, eUp.originVertex, orgUp, dstUp, orgLo, dstLo);
regUp.RegionAbove().dirty = regUp.dirty = regLo.dirty = true;
return false;
}
static void ConnectLeftVertex(Tesselator tess, ContourVertex vEvent)
/*
* Purpose: connect a "left" vertex (one where both edges go right)
* to the processed portion of the mesh. Let R be the active region
* containing vEvent, and let U and L be the upper and lower edge
* chains of R. There are two possibilities:
*
* - the normal case: split R into two regions, by connecting vEvent to
* the rightmost vertex of U or L lying to the left of the sweep line
*
* - the degenerate case: if vEvent is close enough to U or L, we
* merge vEvent into that edge chain. The sub-cases are:
* - merging with the rightmost vertex of U or L
* - merging with the active edge of U or L
* - merging with an already-processed portion of U or L
*/
{
ActiveRegion regUp, regLo, reg;
HalfEdge eUp, eLo, eNew;
ActiveRegion tmp = new ActiveRegion();
/* assert( vEvent.anEdge.Onext.Onext == vEvent.anEdge ); */
/* Get a pointer to the active region containing vEvent */
tmp.upperHalfEdge = vEvent.edgeThisIsOriginOf.otherHalfOfThisEdge;
/* __GL_DICTLISTKEY */
/* __gl_dictListSearch */
regUp = Dictionary.dictSearch(tess.edgeDictionary, tmp).Key;
regLo = RegionBelow(regUp);
eUp = regUp.upperHalfEdge;
eLo = regLo.upperHalfEdge;
/* Try merging with U or L first */
if (ContourVertex.EdgeSign(eUp.directionVertex, vEvent, eUp.originVertex) == 0)
{
ConnectLeftDegenerate(tess, regUp, vEvent);
return;
}
/* Connect vEvent to rightmost processed vertex of either chain.
* e.Dst is the vertex that we will connect to vEvent.
*/
reg = eLo.directionVertex.VertLeq(eUp.directionVertex) ? regUp : regLo;
if (regUp.inside || reg.fixUpperEdge)
{
if (reg == regUp)
{
eNew = Mesh.meshConnect(vEvent.edgeThisIsOriginOf.otherHalfOfThisEdge, eUp.nextEdgeCCWAroundLeftFace);
}
else
{
HalfEdge tempHalfEdge = Mesh.meshConnect(eLo.Dnext, vEvent.edgeThisIsOriginOf);
eNew = tempHalfEdge.otherHalfOfThisEdge;
}
if (reg.fixUpperEdge)
{
FixUpperEdge(reg, eNew);
}
else
{
ComputeWinding(tess, AddRegionBelow(tess, regUp, eNew));
}
SweepEvent(tess, vEvent);
}
else
{
/* The new vertex is in a region which does not belong to the polygon.
* We don''t need to connect this vertex to the rest of the mesh.
*/
AddRightEdges(tess, regUp, vEvent.edgeThisIsOriginOf, vEvent.edgeThisIsOriginOf, null, true);
}
}
static void ComputeWinding(Tesselator tess, ActiveRegion reg)
{
reg.windingNumber = reg.RegionAbove().windingNumber + reg.upperHalfEdge.winding;
reg.inside = tess.IsWindingInside(reg.windingNumber);
}
static void InitEdgeDict(Tesselator tess)
/*
* We maintain an ordering of edge intersections with the sweep line.
* This order is maintained in a dynamic dictionary.
*/
{
/* __gl_dictListNewDict */
tess.edgeDictionary = new Dictionary(tess);
AddSentinel(tess, -SENTINEL_COORD);
AddSentinel(tess, SENTINEL_COORD);
}
static void FinishRegion(Tesselator tess, ActiveRegion reg)
/*
* Delete a region from the sweep line. This happens when the upper
* and lower chains of a region meet (at a vertex on the sweep line).
* The "inside" flag is copied to the appropriate mesh face (we could
* not do this before -- since the structure of the mesh is always
* changing, this face may not have even existed until now).
*/
{
HalfEdge e = reg.upperHalfEdge;
Face f = e.leftFace;
f.isInterior = reg.inside;
f.halfEdgeThisIsLeftFaceOf = e; // optimization for mesh.TessellateMonoRegion()
DeleteRegion(reg);
}
static HalfEdge FinishLeftRegions(Tesselator tess,
ActiveRegion regFirst, ActiveRegion regLast)
/*
* We are given a vertex with one or more left-going edges. All affected
* edges should be in the edge dictionary. Starting at regFirst.eUp,
* we walk down deleting all regions where both edges have the same
* origin vOrg. At the same time we copy the "inside" flag from the
* active region to the face, since at this point each face will belong
* to at most one region (this was not necessarily true until this point
* in the sweep). The walk stops at the region above regLast; if regLast
* is null we walk as far as possible. At the same time we relink the
* mesh if necessary, so that the ordering of edges around vOrg is the
* same as in the dictionary.
*/
{
ActiveRegion reg, regPrev;
HalfEdge e, ePrev;
regPrev = regFirst;
ePrev = regFirst.upperHalfEdge;
while (regPrev != regLast)
{
regPrev.fixUpperEdge = false; /* placement was OK */
reg = RegionBelow(regPrev);
e = reg.upperHalfEdge;
if (e.originVertex != ePrev.originVertex)
{
if (!reg.fixUpperEdge)
{
/* Remove the last left-going edge. Even though there are no further
* edges in the dictionary with this origin, there may be further
* such edges in the mesh (if we are adding left edges to a vertex
* that has already been processed). Thus it is important to call
* FinishRegion rather than just DeleteRegion.
*/
FinishRegion(tess, regPrev);
break;
}
/* If the edge below was a temporary edge introduced by
* ConnectRightVertex, now is the time to fix it.
*/
e = Mesh.meshConnect(ePrev.Lprev, e.otherHalfOfThisEdge);
FixUpperEdge(reg, e);
}
/* Relink edges so that ePrev.Onext == e */
if (ePrev.nextEdgeCCWAroundOrigin != e)
{
Mesh.meshSplice(e.Oprev, e);
Mesh.meshSplice(ePrev, e);
}
FinishRegion(tess, regPrev); /* may change reg.eUp */
ePrev = reg.upperHalfEdge;
regPrev = reg;
}
return ePrev;
}
public TessTool(Tesselator tess)
{
_tess = tess;
_tessListener = new TessListener();
_tessListener.Connect(tess, true);
}
static void AddRightEdges(Tesselator tess, ActiveRegion regUp,
HalfEdge eFirst, HalfEdge eLast, HalfEdge eTopLeft,
bool cleanUp)
/*
* Purpose: insert right-going edges into the edge dictionary, and update
* winding numbers and mesh connectivity appropriately. All right-going
* edges share a common origin vOrg. Edges are inserted CCW starting at
* eFirst; the last edge inserted is eLast.Oprev. If vOrg has any
* left-going edges already processed, then eTopLeft must be the edge
* such that an imaginary upward vertical segment from vOrg would be
* contained between eTopLeft.Oprev and eTopLeft; otherwise eTopLeft
* should be null.
*/
{
ActiveRegion reg, regPrev;
HalfEdge e, ePrev;
bool firstTime = true;
/* Insert the new right-going edges in the dictionary */
e = eFirst;
do
{
if (!e.originVertex.VertLeq(e.directionVertex))
{
throw new Exception();
}
AddRegionBelow(tess, regUp, e.otherHalfOfThisEdge);
e = e.nextEdgeCCWAroundOrigin;
} while (e != eLast);
/* Walk *all* right-going edges from e.Org, in the dictionary order,
* updating the winding numbers of each region, and re-linking the mesh
* edges to match the dictionary ordering (if necessary).
*/
if (eTopLeft == null)
{
eTopLeft = RegionBelow(regUp).upperHalfEdge.Rprev;
}
regPrev = regUp;
ePrev = eTopLeft;
for (;;)
{
reg = RegionBelow(regPrev);
e = reg.upperHalfEdge.otherHalfOfThisEdge;
if (e.originVertex != ePrev.originVertex) break;
if (e.nextEdgeCCWAroundOrigin != ePrev)
{
/* Unlink e from its current position, and relink below ePrev */
Mesh.meshSplice(e.Oprev, e);
Mesh.meshSplice(ePrev.Oprev, e);
}
/* Compute the winding number and "inside" flag for the new regions */
reg.windingNumber = regPrev.windingNumber - e.winding;
reg.inside = tess.IsWindingInside(reg.windingNumber);
/* Check for two outgoing edges with same slope -- process these
* before any intersection tests (see example in __gl_computeInterior).
*/
regPrev.dirty = true;
if (!firstTime && CheckForRightSplice(tess, regPrev))
{
AddWinding(e, ePrev);
DeleteRegion(regPrev);
Mesh.DeleteHalfEdge(ePrev);
}
firstTime = false;
regPrev = reg;
ePrev = e;
}
regPrev.dirty = true;
if (regPrev.windingNumber - e.winding != reg.windingNumber)
{
throw new Exception();
}
if (cleanUp)
{
/* Check for intersections between newly adjacent edges. */
WalkDirtyRegions(tess, regPrev);
}
}
// routine to render this primitive
public void CallRender(Tesselator tess, HalfEdge edge, int data)
{
render(tess, edge, data);
}
static void CallCombine(Tesselator tess, ContourVertex intersectionVertex, int[] vertexIndexArray, double[] vertexWeights, bool needed)
{
/* Copy coord data in case the callback changes it. */
double c0 = intersectionVertex.C_0;
double c1 = intersectionVertex.C_1;
double c2 = intersectionVertex.C_2;
intersectionVertex.clientIndex = 0;
tess.CallCombine(c0, c1, c2, vertexIndexArray, vertexWeights, out intersectionVertex.clientIndex);
if (intersectionVertex.clientIndex == 0)
{
if (!needed)
{
intersectionVertex.clientIndex = vertexIndexArray[0];
}
else
{
/* The only fatal error is when two edges are found to intersect,
* but the user has not provided the callback necessary to handle
* generated intersection points.
*/
throw new Exception("You need to provided a callback to handle generated intersection points.");
}
}
}
public void ParseStreamForTesselator(Tesselate.Tesselator tesselator, int instructionStreamIndex)
{
m_VertexList.Clear();
m_CurrentOutput = 0;
string[] instructionStream = m_InsructionStream[instructionStreamIndex];
for (int curInstruction = 0; curInstruction < instructionStream.Length; curInstruction++)
{
switch (instructionStream[curInstruction])
{
case "BP":
tesselator.BeginPolygon();
break;
case "BC":
tesselator.BeginContour();
break;
case "V":
double x = Convert.ToDouble(instructionStream[curInstruction + 1]);
double y = Convert.ToDouble(instructionStream[curInstruction + 2]);
curInstruction += 2;
tesselator.AddVertex(x, y, 0, m_VertexList.Count);
m_VertexList.Add(new Vertex(x, y));
break;
case "EC":
tesselator.EndContour();
break;
case "EP":
tesselator.EndPolygon();
break;
default:
throw new Exception();
}
}
//---------------
//test
//------------------
//int a = m_VertexList.Count;
////test draw
//Bitmap bmp = new Bitmap(300, 300);
//Graphics g = Graphics.FromImage(bmp);
//g.Clear(Color.White);
//int j = m_VertexList.Count;
//int lim = j - 1;
//for (int i = 0; i < lim; i++)
//{
// var v1 = m_VertexList[i];
// var v2 = m_VertexList[i + 1];
// g.DrawLine(Pens.Black,
// new PointF((float)v1.m_X, (float)v1.m_Y),
// new PointF((float)v2.m_X, (float)v2.m_Y));
//}
//g.Dispose();
//bmp.Save("tess" + (totalSampleCount++) + ".bmp");
}
static void SpliceMergeVertices(Tesselator tess, HalfEdge e1, HalfEdge e2)
/*
* Two vertices with idential coordinates are combined into one.
* e1.Org is kept, while e2.Org is discarded.
*/
{
int[] data4 = new int[4];
double[] weights4 = new double[] { 0.5f, 0.5f, 0, 0 };
data4[0] = e1.originVertex.clientIndex;
data4[1] = e2.originVertex.clientIndex;
CallCombine(tess, e1.originVertex, data4, weights4, false);
Mesh.meshSplice(e1, e2);
}
public static bool EdgeLeq(Tesselator tess, ActiveRegion reg1, ActiveRegion reg2)
/*
* Both edges must be directed from right to left (this is the canonical
* direction for the upper edge of each region).
*
* The strategy is to evaluate a "t" value for each edge at the
* current sweep line position, given by tess.currentSweepVertex. The calculations
* are designed to be very stable, but of course they are not perfect.
*
* Special case: if both edge destinations are at the sweep event,
* we sort the edges by slope (they would otherwise compare equally).
*/
{
ContourVertex currentSweepVertex = tess.currentSweepVertex;
HalfEdge e1, e2;
double t1, t2;
e1 = reg1.upperHalfEdge;
e2 = reg2.upperHalfEdge;
if (e1.directionVertex == currentSweepVertex)
{
if (e2.directionVertex == currentSweepVertex)
{
/* Two edges right of the sweep line which meet at the sweep currentSweepVertex.
* Sort them by slope.
*/
if (e1.originVertex.VertLeq(e2.originVertex))
{
return ContourVertex.EdgeSign(e2.directionVertex, e1.originVertex, e2.originVertex) <= 0;
}
return ContourVertex.EdgeSign(e1.directionVertex, e2.originVertex, e1.originVertex) >= 0;
}
return ContourVertex.EdgeSign(e2.directionVertex, currentSweepVertex, e2.originVertex) <= 0;
}
if (e2.directionVertex == currentSweepVertex)
{
return ContourVertex.EdgeSign(e1.directionVertex, currentSweepVertex, e1.originVertex) >= 0;
}
/* General case - compute signed distance *from* e1, e2 to currentSweepVertex */
t1 = ContourVertex.EdgeEval(e1.directionVertex, currentSweepVertex, e1.originVertex);
t2 = ContourVertex.EdgeEval(e2.directionVertex, currentSweepVertex, e2.originVertex);
return (t1 >= t2);
}
static void WalkDirtyRegions(Tesselator tess, ActiveRegion regUp)
/*
* When the upper or lower edge of any region changes, the region is
* marked "dirty". This routine walks through all the dirty regions
* and makes sure that the dictionary invariants are satisfied
* (see the comments at the beginning of this file). Of course
* new dirty regions can be created as we make changes to restore
* the invariants.
*/
{
ActiveRegion regLo = RegionBelow(regUp);
HalfEdge eUp, eLo;
for (;;)
{
/* Find the lowest dirty region (we walk from the bottom up). */
while (regLo.dirty)
{
regUp = regLo;
regLo = RegionBelow(regLo);
}
if (!regUp.dirty)
{
regLo = regUp;
regUp = regUp.RegionAbove();
if (regUp == null || !regUp.dirty)
{
/* We've walked all the dirty regions */
return;
}
}
regUp.dirty = false;
eUp = regUp.upperHalfEdge;
eLo = regLo.upperHalfEdge;
if (eUp.directionVertex != eLo.directionVertex)
{
/* Check that the edge ordering is obeyed at the Dst vertices. */
if (CheckForLeftSplice(tess, regUp))
{
/* If the upper or lower edge was marked fixUpperEdge, then
* we no longer need it (since these edges are needed only for
* vertices which otherwise have no right-going edges).
*/
if (regLo.fixUpperEdge)
{
DeleteRegion(regLo);
Mesh.DeleteHalfEdge(eLo);
regLo = RegionBelow(regUp);
eLo = regLo.upperHalfEdge;
}
else if (regUp.fixUpperEdge)
{
DeleteRegion(regUp);
Mesh.DeleteHalfEdge(eUp);
regUp = regLo.RegionAbove();
eUp = regUp.upperHalfEdge;
}
}
}
if (eUp.originVertex != eLo.originVertex)
{
if (eUp.directionVertex != eLo.directionVertex
&& !regUp.fixUpperEdge && !regLo.fixUpperEdge
&& (eUp.directionVertex == tess.currentSweepVertex || eLo.directionVertex == tess.currentSweepVertex))
{
/* When all else fails in CheckForIntersect(), it uses tess.currentSweepVertex
* as the intersection location. To make this possible, it requires
* that tess.currentSweepVertex lie between the upper and lower edges, and also
* that neither of these is marked fixUpperEdge (since in the worst
* case it might splice one of these edges into tess.currentSweepVertex, and
* violate the invariant that fixable edges are the only right-going
* edge from their associated vertex).
*/
if (CheckForIntersect(tess, regUp))
{
/* WalkDirtyRegions() was called recursively; we're done */
return;
}
}
else
{
/* Even though we can't use CheckForIntersect(), the Org vertices
* may violate the dictionary edge ordering. Check and correct this.
*/
CheckForRightSplice(tess, regUp);
}
}
if (eUp.originVertex == eLo.originVertex && eUp.directionVertex == eLo.directionVertex)
{
/* A degenerate loop consisting of only two edges -- delete it. */
AddWinding(eLo, eUp);
DeleteRegion(regUp);
Mesh.DeleteHalfEdge(eUp);
regUp = regLo.RegionAbove();
}
}
}