Skip to Content.
Sympa Menu

cgal-discuss - [cgal-discuss] 3D meshing - Subdomain index of a cell and it's mirror

Subject: CGAL users discussion list

List archive

[cgal-discuss] 3D meshing - Subdomain index of a cell and it's mirror


Chronological Thread 
  • From: Ramin H <>
  • To:
  • Subject: [cgal-discuss] 3D meshing - Subdomain index of a cell and it's mirror
  • Date: Mon, 3 May 2010 12:43:42 -0400
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:from:date:message-id:subject:to:content-type; b=WiGkxd/26j4FVfEzZpdkolIj+8cWmwh4dx4UmVruRSX9Wx6GaPNHmdiree9pbH3O5k ObLwbmcapbALl1PRhoNiVbvAf8mmC6DA/9ASSGxt98SJvigzRVvz7axOBkb1f7pK+4+B 5Hn0Uu4oSDm+Bdetmptv1uu1jGznQ0x2HSFKs=

Hi,
I still don't know everything about C3T3 even though I have read the documentation. In a 3D triangulation, let's say the cell representing a facet does not belong to domain (subdomain=0). If we find the mirror cell of this facet, should it belong to the domain (subdomain > 0) ?

To test this idea, I wrote the following code:

  for(Facet_iterator facet_it = c3t3_.facets_begin() ;
      facet_it != c3t3_.facets_end() ;
      ++facet_it)
  {
    Facet facet = *facet_it;
    const Facet f = c3t3_.triangulation().mirror_facet(facet);
    if( (int)(facet.first->subdomain_index()) == (int)(f.first->subdomain_index()) )
    {
      std::cerr << "Cell and it's mirror have equal subindex: " << (int)(f.first->subdomain_index()) << std::endl;
    }
  }

Here is some sample output:
Cell and it's mirror have equal subindex: 1
Cell and it's mirror have equal subindex: 1
Cell and it's mirror have equal subindex: 0
Cell and it's mirror have equal subindex: 1
Cell and it's mirror have equal subindex: 1
Cell and it's mirror have equal subindex: 1
Cell and it's mirror have equal subindex: 0
Cell and it's mirror have equal subindex: 1

Is this correct behavior? If Facet_iterator is iterating over the facets on the surface, why a cell representing that facet and it's mirror both can have subdomain index equal to 1 or 0?

Thanks,
-Ramin

PS. I have attached the input polyhedral domain as an .off file.


Attachment: _domain_to_mesh.off
Description: Binary data




Archive powered by MHonArc 2.6.16.

Top of Page