#include <tetmesh_faceted_surfaces.h>

Inheritance diagram for oomph::DiskWithTwoLayersTetMeshFacetedSurface:

Public Member Functions
	DiskWithTwoLayersTetMeshFacetedSurface (DiskLikeGeomObjectWithBoundaries *disk_parametrised_by_nonsingular_coordinates_pt, const unsigned &half_nsegment, const unsigned &first_one_based_boundary_id_for_disk, const unsigned &one_based_region_id_above_disk, const unsigned &one_based_region_id_below_disk, unsigned &last_one_based_boundary_id_for_disk, unsigned &last_one_based_boundary_id)

Public Member Functions inherited from oomph::DiskTetMeshFacetedSurface
	DiskTetMeshFacetedSurface (DiskLikeGeomObjectWithBoundaries *disk_parametrised_by_nonsingular_coordinates_pt, const unsigned &half_nsegment, const unsigned &first_one_based_boundary_id_for_disk, unsigned &last_one_based_boundary_id_for_disk)

	~DiskTetMeshFacetedSurface ()
	Destructor. More...

void	boundary_zeta01 (const unsigned &facet_id, const double &zeta_boundary, Vector< double > &zeta)

Public Member Functions inherited from oomph::TetMeshFacetedSurface
	TetMeshFacetedSurface ()
	Constructor: More...

virtual	~TetMeshFacetedSurface ()
	Empty destructor. More...

unsigned	nvertex () const
	Number of vertices. More...

unsigned	nfacet () const
	Number of facets. More...

unsigned	one_based_facet_boundary_id (const unsigned &j) const
	One-based boundary id of j-th facet. More...

unsigned	one_based_vertex_boundary_id (const unsigned &j) const
	First (of possibly multiple) one-based boundary id of j-th vertex. More...

double	vertex_coordinate (const unsigned &j, const unsigned &i) const
	i-th coordinate of j-th vertex More...

unsigned	nvertex_on_facet (const unsigned &j) const
	Number of vertices defining the j-th facet. More...

bool	boundaries_can_be_split_in_tetgen ()
	Test whether boundary can be split in tetgen. More...

void	enable_boundaries_can_be_split_in_tetgen ()
	Test whether boundaries can be split in tetgen. More...

void	disable_boundaries_can_be_split_in_tetgen ()
	Test whether boundaries can be split in tetgen. More...

TetMeshFacet *	facet_pt (const unsigned &j) const
	Pointer to j-th facet. More...

TetMeshVertex *	vertex_pt (const unsigned &j) const
	Pointer to j-th vertex. More...

DiskLikeGeomObjectWithBoundaries *	geom_object_with_boundaries_pt ()

void	output (std::ostream &outfile) const
	Output. More...

void	output (const std::string &filename) const
	Output. More...

virtual void	boundary_zeta12 (const unsigned &facet_id, const double &zeta_boundary, Vector< double > &zeta)

virtual void	boundary_zeta20 (const unsigned &facet_id, const double &zeta_boundary, Vector< double > &zeta)

Vector< unsigned >	vertex_index_in_tetgen (const unsigned &f)

Public Member Functions inherited from oomph::TetMeshFacetedClosedSurface
	TetMeshFacetedClosedSurface ()
	Constructor: More...

virtual	~TetMeshFacetedClosedSurface ()
	Empty destructor. More...

void	enable_faceted_volume_represents_hole_for_gmsh ()
	Declare closed surface to represent hole for gmsh. More...

void	disable_faceted_volume_represents_hole_for_gmsh ()
	Declare closed surface NOT to represent hole for gmsh. More...

bool	faceted_volume_represents_hole_for_gmsh () const
	Does closed surface represent hole for gmsh? More...

const double &	internal_point_for_tetgen (const unsigned &j, const unsigned &i) const
	i=th coordinate of the j-th internal point for tetgen More...

void	set_hole_for_tetgen (const Vector< double > &hole_point)
	Specify coordinate of hole for tetgen. More...

void	set_region_for_tetgen (const unsigned &region_id, const Vector< double > &region_point)

unsigned	ninternal_point_for_tetgen ()

const int &	region_id_for_tetgen (const unsigned &j) const

bool	internal_point_identifies_hole_for_tetgen (const unsigned &j)
	Is j-th internal point for tetgen associated with a hole? More...

bool	internal_point_identifies_region_for_tetgen (const unsigned &j)
	Is j-th internal point for tetgen associated with a region? More...

Additional Inherited Members
Protected Attributes inherited from oomph::DiskTetMeshFacetedSurface
unsigned	Nfacet
	Number of facets. More...

std::vector< bool >	Facet_is_on_boundary
	Is facet on boundary? More...

Vector< double >	Left_boundary_coordinate
	Left boundary coordinate of i-th facet. More...

Vector< double >	Right_boundary_coordinate
	Right boundary coordinate of i-th facet. More...

Vector< unsigned >	Boundary_id
	ID of boundary the i-th facet is located on. More...

unsigned	Nelement_on_disk_boundary
	Number of elements on disk boundary. More...

TriangleMesh< TPoissonElement< 2, 2 > > *	Tri_mesh_pt
	Mesh used to facet-ise (discretise) disk. More...

std::map< Node , TetMeshVertex >	Equivalent_vertex_pt
	Mapping between nodes and vertices. More...

Protected Attributes inherited from oomph::TetMeshFacetedSurface
Vector< TetMeshVertex * >	Vertex_pt
	Vector pointers to vertices. More...

Vector< TetMeshFacet * >	Facet_pt
	Vector of pointers to facets. More...

bool	Boundaries_can_be_split_in_tetgen

Vector< Vector< unsigned > >	Facet_vertex_index_in_tetgen

DiskLikeGeomObjectWithBoundaries *	Geom_object_with_boundaries_pt
	GeomObject with boundaries associated with this surface. More...

Detailed Description

//////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// TetMeshFacetedClosedSurface that defines disk with two layers on top and bottom

Constructor & Destructor Documentation

◆ DiskWithTwoLayersTetMeshFacetedSurface()

oomph::DiskWithTwoLayersTetMeshFacetedSurface::DiskWithTwoLayersTetMeshFacetedSurface	(	DiskLikeGeomObjectWithBoundaries *	disk_parametrised_by_nonsingular_coordinates_pt,
		const unsigned &	half_nsegment,
		const unsigned &	first_one_based_boundary_id_for_disk,
		const unsigned &	one_based_region_id_above_disk,
		const unsigned &	one_based_region_id_below_disk,
		unsigned &	last_one_based_boundary_id_for_disk,
		unsigned &	last_one_based_boundary_id
	)

inline

Constructor: Pass pointer to GeomObject (with boundaries, and parametrised by coordinates without coordinate singularities, i.e. not polars, say) that defines the shape of the disk. Other args specify:

half the number of segments on perimeter of disk
first one-based boundary ID to be used to enumerate the boundaries on the disk
one based region ID for region above the disk
one based region ID for region below the disk

Computes:

last one-based boundary id used to enumerate boundaries on disk
last one-based boundary id used to enumerate any boundaries

Identify region for tetgen

                                         :
  DiskTetMeshFacetedSurface(disk_parametrised_by_nonsingular_coordinates_pt,
                            half_nsegment,
                            first_one_based_boundary_id_for_disk,
                            last_one_based_boundary_id_for_disk)
   {
  
    // Currently: Number of facets = number of elements on disk.
    // Set bounding region numbers for existing facets:
    unsigned nel=Nfacet;
    for (unsigned e=0;e<Nfacet;e++)
     {
      Facet_pt[e]->set_one_based_adjacent_region_id(
       one_based_region_id_above_disk);
      Facet_pt[e]->set_one_based_adjacent_region_id(
       one_based_region_id_below_disk);
     }
  
    // Additional facets: one square facet moving up/down from each
    // side and two lids.
    Nfacet+=2*Nelement_on_disk_boundary+2;
    Facet_pt.resize(Nfacet);
    Facet_is_on_boundary.resize(Nfacet,false);
  
    // Find sorted sequence of vertices going around the
    // perimeter of the disk
    Vector<TetMeshVertex*> sorted_perimeter_vertex_pt;
    std::set<TetMeshFacet*> boundary_facet_pt;
    bool first=true;
    for (unsigned e=0;e<nel;e++)
     {
      if (Facet_is_on_boundary[e])
       {
        boundary_facet_pt.insert(Facet_pt[e]);
        if (first)
         {
          TetMeshVertex* vertex0_pt=Facet_pt[e]->vertex_pt(0);
          sorted_perimeter_vertex_pt.push_back(vertex0_pt);
          first=false;
         }
       }
     }
    while (boundary_facet_pt.size()!=1)
     {
      for (std::set<TetMeshFacet*>::iterator it=boundary_facet_pt.begin();
           it!=boundary_facet_pt.end();it++)
       {
        unsigned n_last=sorted_perimeter_vertex_pt.size();
        if (sorted_perimeter_vertex_pt[n_last-1]==
            (*it)->vertex_pt(0))
         {
          sorted_perimeter_vertex_pt.push_back((*it)->vertex_pt(1));
          boundary_facet_pt.erase(*it);
          break;
         }
       }
     }
  
    //Find mean position
    double z_mean=0.0;
    double z_min= DBL_MAX;
    double z_max=-DBL_MAX;
    unsigned n=sorted_perimeter_vertex_pt.size();
    for (unsigned j=0;j<n;j++)
     {
      double z=sorted_perimeter_vertex_pt[j]->x(2);
      z_mean+=z;
      if (z>z_max) z_max=z;
      if (z<z_min) z_min=z;
     }
    z_mean/=double(n);
    double layer_thickness=0.6*(z_max-z_min);
    
    // Loop over upper/lower layer
    unsigned facet_count=nel;
    int offset_sign=-1;
    unsigned one_based_boundary_id=last_one_based_boundary_id_for_disk+1;
    unsigned n_vertex_on_disk=Vertex_pt.size();
    unsigned n_vertex_on_disk_perimeter=sorted_perimeter_vertex_pt.size();
    for (unsigned i_offset=0;i_offset<2;i_offset++)
     {
      unsigned one_based_region_id=one_based_region_id_below_disk;
      if (i_offset==1) one_based_region_id=one_based_region_id_above_disk;
  
      for (unsigned j=0;j<n_vertex_on_disk_perimeter;j++)
       {
        Vector<double> posn(3);
        posn[0]=sorted_perimeter_vertex_pt[j]->x(0);
        posn[1]=sorted_perimeter_vertex_pt[j]->x(1);
        posn[2]=z_mean+double(offset_sign)*layer_thickness;
        Vertex_pt.push_back(new TetMeshVertex(posn));
       }
      unsigned count=0;
      for (unsigned j=0;j<n_vertex_on_disk_perimeter-1;j++)
       {
        Facet_pt[facet_count]=new TetMeshFacet(4);
        Facet_pt[facet_count]->set_one_based_boundary_id(one_based_boundary_id);
        Facet_pt[facet_count]->set_vertex_pt(0,sorted_perimeter_vertex_pt[j]);
        Facet_pt[facet_count]->set_vertex_pt(1,sorted_perimeter_vertex_pt[j+1]);
        Facet_pt[facet_count]->
         set_vertex_pt(2,Vertex_pt[n_vertex_on_disk+
                                   i_offset*n_vertex_on_disk_perimeter+j+1]);
        Facet_pt[facet_count]->
         set_vertex_pt(3,Vertex_pt[n_vertex_on_disk+
                                   i_offset*n_vertex_on_disk_perimeter+j  ]);
        Facet_pt[facet_count]->
         set_one_based_adjacent_region_id(one_based_region_id);
        count++;
        facet_count++;
        one_based_boundary_id++;
       }
      Facet_pt[facet_count]=new TetMeshFacet(4);
      Facet_pt[facet_count]->set_one_based_boundary_id(one_based_boundary_id);
      Facet_pt[facet_count]->set_vertex_pt(0,sorted_perimeter_vertex_pt[count]);
      Facet_pt[facet_count]->set_vertex_pt(1,sorted_perimeter_vertex_pt[0]);
      Facet_pt[facet_count]->
       set_vertex_pt(2,Vertex_pt[n_vertex_on_disk+
                                 i_offset*n_vertex_on_disk_perimeter+0]);
      Facet_pt[facet_count]->
       set_vertex_pt(3,Vertex_pt[n_vertex_on_disk+
                                 i_offset*n_vertex_on_disk_perimeter+count]);
      Facet_pt[facet_count]->
       set_one_based_adjacent_region_id(one_based_region_id);
  
  
      // Closing lid
      facet_count++;
      one_based_boundary_id++;
      Facet_pt[facet_count]=new TetMeshFacet(n_vertex_on_disk_perimeter);
      Facet_pt[facet_count]->set_one_based_boundary_id(one_based_boundary_id);
      for (unsigned j=0;j<n_vertex_on_disk_perimeter;j++)
       {
        Facet_pt[facet_count]->
         set_vertex_pt(j,Vertex_pt[n_vertex_on_disk+
                                   i_offset*n_vertex_on_disk_perimeter+j]);
       }
      Facet_pt[facet_count]->
       set_one_based_adjacent_region_id(one_based_region_id);
  
      // Last one used (gets assigned twice but last one sticks!)
      last_one_based_boundary_id=one_based_boundary_id;
  
      // Bump
      facet_count++;
      one_based_boundary_id++;
  
      Vector<double> region_point(3,0.0);
      region_point[2]=z_mean+0.5*double(offset_sign)*layer_thickness;
      set_region_for_tetgen(one_based_region_id-1,region_point);
  
      offset_sign*=-1;
     }
   }

References e(), oomph::DiskTetMeshFacetedSurface::Facet_is_on_boundary, oomph::TetMeshFacetedSurface::Facet_pt, j, n, oomph::DiskTetMeshFacetedSurface::Nelement_on_disk_boundary, oomph::DiskTetMeshFacetedSurface::Nfacet, oomph::TetMeshFacetedClosedSurface::set_region_for_tetgen(), and oomph::TetMeshFacetedSurface::Vertex_pt.

The documentation for this class was generated from the following file:

tetmesh_faceted_surfaces.h

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DiskWithTwoLayersTetMeshFacetedSurface()