songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #3567 from goran-w/goran-w-patch-3566 

Modified output_to_medit() for consistent facet orientation
This commit is contained in:
Laurent Rineau 2019-01-18 16:33:32 +01:00
parent 45066e7d91 509594f21d
commit 73ea2ee4ea
1 changed files with 17 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
Mesh_3/include/CGAL/IO/File_medit.h
View File

@ -815,27 +815,28 @@ output_to_medit(std::ostream& os,
fit != c3t3.facets_in_complex_end();
++fit)
{
for (int i=0; i<4; i++)
{
if (i != fit->second)
{
const Vertex_handle& vh = (*fit).first->vertex(i);
os << V[vh] << ' ';
}
}
typename C3T3::Facet f = (*fit);
// Apply priority among subdomains, to get consistent facet orientation per subdomain-pair interface.
if (f.first->subdomain_index() < f.first->neighbor(f.second)->subdomain_index())
f = tr.mirror_facet(f);
// Get facet vertices in CCW order.
Vertex_handle vh1 = f.first->vertex((f.second + 1) % 4);
Vertex_handle vh2 = f.first->vertex((f.second + 2) % 4);
Vertex_handle vh3 = f.first->vertex((f.second + 3) % 4);
// Facet orientation also depends on parity.
if (f.second % 2 != 0)
std::swap(vh2, vh3);
os << V[vh1] << ' ' << V[vh2] << ' ' << V[vh3] << ' ';
os << get(facet_pmap, *fit) << '\n';
// Print triangle again if needed
if ( print_each_facet_twice )
{
for (int i=0; i<4; i++)
{
if (i != fit->second)
{
const Vertex_handle& vh = (*fit).first->vertex(i);
os << V[vh] << ' ';
}
}
os << V[vh1] << ' ' << V[vh2] << ' ' << V[vh3] << ' ';
os << get(facet_twice_pmap, *fit) << '\n';
}
}