CGAL users discussion list

[tang <>]

Dear Dr. Rineau,

Thank you very much for your kindly reply. I still want to simplify the
surface mesh.

Since the CGAL can only cope with short edge, I add a short edge to the flat
triangle (2*area(triangle)/max(length)<minlength), and then, I will let the
CGAL to simplify the mesh. 

my idea is as follows:

1) add short edge to flag triangles (the function 'splitsmalltriangle');

2) simplify the mesh by removing short edges.

The steps of funciton splitsmalltriangle is as as follows:

for every facet (triangle):
  calculate the area of the triangle;
  calculate lengths of 3 edges and find the maximal length;
  if(any length &lt; minlength) continue;
  if(2*area/max(length) > minlength) continue;
  get the projection of A (say it to B) on maximal edge, where A is the
opposite vertex of maximal edge;
  split the triangle by edge AB;
end

Now I don't know how to split the facet (triangle) by this edge, could you
please give me some hints?


Thanks,
Zhanghong Tang


BTW: the code posted last time has a mistake, it should be like this:
                dmax=max(max(d12, d23), d31);
                if(area/dmax>edgelengthmin)continue;
                /// get the projection of pi to longest edge
                if(dmax==d12)
                ...




--
View this message in context: 
http://cgal-discuss.949826.n4.nabble.com/Error-when-building-CGAL4-0-2-demo-examples-64-bit-version-tp4655937p4655954.html
Sent from the cgal-discuss mailing list archive at Nabble.com.