numcxx::SimpleGrid Class Reference

Description

Class containing data for simple grid data structure.

Examples:

40-stationary-heat-fe.cxx, 41-stationary-heat-fv.cxx, 42-convtest-fem-sin.cxx, 43-convtest-fvm-sin.cxx, 44-transient-heat-fe.cxx, and 46-nonlin-fvm.cxx.

Definition at line 19 of file simplegrid.hxx.

Public Member Functions
	SimpleGrid ()
	Trivial construtor. More...
	SimpleGrid (const SimpleGrid &g)=delete
	SimpleGrid (const Geometry &geometry, const char *triangle_flags)
	Construct simple grid from geometry via library call to Triangle mesh generator. More...
const TArray2< double > &	get_points () const
	Get array of point coordinates. More...
const TArray2< int > &	get_cells () const
	Get array of point indices describing cells. More...
const TArray1< int > &	get_cellregions () const
	Get array region markers. More...
const TArray2< int > &	get_bfaces () const
	Get array of point indices describing boundary faces. More...
const TArray1< int > &	get_bfaceregions () const
	Get array of boundary markers. More...
const int	spacedim () const
	Return dimension of space. More...
const int	griddim () const
	Return dimension of grid. More...
const int	ncells () const
	Return number of cells. More...
const int	npoints () const
	Return number of points. More...
const int	nbfaces () const
	Return number of boundary faces. More...
void	calc_hminmax (double &hmin, double &hmax) const
	Calculate some grid data. More...

Static Public Member Functions
static std::shared_ptr< SimpleGrid >	create (const Geometry &geometry, const std::string triangle_flags)
	Static constructor from reference to geometry. More...
static std::shared_ptr< SimpleGrid >	create (std::shared_ptr< Geometry > geometry, const char *triangle_flags)
	Static constructor from shared pointer to geometry. More...

Private Attributes
std::shared_ptr< TArray2< double > >	points =nullptr
std::shared_ptr< TArray2< int > >	cells =nullptr
std::shared_ptr< TArray1< int > >	cellregions =nullptr
std::shared_ptr< TArray2< int > >	bfaces =nullptr
std::shared_ptr< TArray1< int > >	bfaceregions =nullptr

Constructor & Destructor Documentation

numcxx::SimpleGrid::SimpleGrid ( )

inline

Trivial construtor.

Definition at line 24 of file simplegrid.hxx.

{};

numcxx::SimpleGrid::SimpleGrid ( const SimpleGrid &  g )

delete

numcxx::SimpleGrid::SimpleGrid	(	const Geometry &	geometry,
		const char *	triangle_flags
	)

Construct simple grid from geometry via library call to Triangle mesh generator.

The triangle_flags must contain the z flag

Definition at line 79 of file simplegrid.cxx.

  {
    if (strchr(flags,'z')==nullptr)
    {
      throw std::runtime_error("numcxx: triangulate: Missing 'z' flag"); 
    }
    
    struct triangle::triangulateio in,out;
    int dbg=0;
        
    out.pointlist=0;
    out.pointattributelist=0;
    out.pointmarkerlist=0;
    out.numberofpoints=0;
    out.numberofpointattributes=0;
    out.trianglelist=0;
    out.triangleattributelist=0;
    out.trianglearealist=0;
    out.neighborlist=0;
    out.numberoftriangles=0;
    out.numberofcorners=0;
    out.numberoftriangleattributes=0;
    out.segmentlist=0;
    out.segmentmarkerlist=0;
    out.numberofsegments=0;
    out.holelist=0;
    out.numberofholes=0;
    out.regionlist=0;
    out.numberofregions=0;
    out.edgelist=0;
    out.edgemarkerlist=0;
    out.normlist=0;
    out.numberofedges=0;
        
    in.pointlist=0;
    in.pointattributelist=0;
    in.pointmarkerlist=0;
    in.numberofpoints=0;
    in.numberofpointattributes=0;
    in.trianglelist=0;
    in.triangleattributelist=0;
    in.trianglearealist=0;
    in.neighborlist=0;
    in.numberoftriangles=0;
    in.numberofcorners=0;
    in.numberoftriangleattributes=0;
    in.segmentlist=0;
    in.segmentmarkerlist=0;
    in.numberofsegments=0;
    in.holelist=0;
    in.numberofholes=0;
    in.regionlist=0;
    in.numberofregions=0;
    in.edgelist=0;
    in.edgemarkerlist=0;
    in.normlist=0;
    in.numberofedges=0;
        
        
    // take  point list
    if (geometry.points)
    {
      in.pointlist=geometry.points->data();
      in.numberofpoints=geometry.points->shape(0);
      in.numberofpointattributes=0;
    }
    else
    {
      throw std::runtime_error("numcxx: triangulate: Missing point list"); 
    }
    // take bfaces if defined
    if (geometry.bfaces)
    {
      if (!geometry.bfaceregions)
        throw std::runtime_error("numcxx: triangulate: Missing bfaceregions"); 

      in.segmentlist=geometry.bfaces->data();
      in.segmentmarkerlist=geometry.bfaceregions->data();
      in.numberofsegments=geometry.bfaces->shape(0);
    }
    else
    {
      throw std::runtime_error("numcxx: triangulate: Missing bface list"); 
    }
        
    // create region and hole information for triangle
    // Holes are marked with region number <0
    if (geometry.regionpoints)
    {
      int nregs,nholes,ireg,ihole;
      int i;
      nregs=nholes=0;
      if (!geometry.regionnumbers)
      {
        throw std::runtime_error("numcxx: triangulate: Missing region numbers"); 
      }
      if (!geometry.regionvolumes)
      {
        throw std::runtime_error("numcxx: triangulate: Missing region volumes"); 
      }
      for (i=0;i<geometry.regionnumbers->shape(0);i++)
      {
        if ((*geometry.regionnumbers)(i)>=0)
          nregs++; 
        else 
          nholes++;
      }
      in.numberofholes=nholes;
      in.numberofregions=nregs;
      if(nholes)
        in.holelist=(double*)malloc(2*nholes*sizeof(double));
      in.regionlist=(double*)malloc(4*nregs*sizeof(double));
      ireg=ihole=0;
            
      for (i=0;i<geometry.regionnumbers->shape(0);i++)
        if ((*geometry.regionnumbers)(i)>=0)
        {
          in.regionlist[ireg+0]=(*geometry.regionpoints)(i,0);
          in.regionlist[ireg+1]=(*geometry.regionpoints)(i,1);
          in.regionlist[ireg+2]=(*geometry.regionnumbers)(i);
          in.regionlist[ireg+3]=(*geometry.regionvolumes)(i);
          ireg+=4;
        }                
        else
        {
          in.holelist[ihole+0]=(*geometry.regionpoints)(i,0); 
          in.holelist[ihole+1]=(*geometry.regionpoints)(i,1); 
          ihole+=2;
        }
            
    }    
        
    
    if (dbg)
    {
      printf("in:\n");
      print_triangulateio(&in);
    }
        
        
    //  Do the damn thing!
    triangle::triangulate((char*)flags,&in,&out,NULL);
        
    if (dbg)
    {
      printf("out (flags=%s):\n",flags);
      print_triangulateio(&out);
    }
        
    // This triangle specific input information is not needed anymore.
    if(in.holelist) triangle::trifree(in.holelist);
    if(in.regionlist) triangle::trifree(in.regionlist);
        
    
    // Take point and cell list as is. 
    points=std::make_shared<numcxx::TArray2<double>>(out.numberofpoints,2,out.pointlist,[](double*p){triangle::trifree(p);});
    cells=std::make_shared<numcxx::TArray2<int>>(out.numberoftriangles,3,out.trianglelist,[](int*p){triangle::trifree(p);});
        

    // Our cellregions are integer, not double as in triangle. 
    // Furthermore, they are mandatory in our process, so we
    // set them to 0 if they are not available
    {
      int *triattr,i;
            
      cellregions=std::make_shared<numcxx::TArray1<int>>(out.numberoftriangles);
      if (out.triangleattributelist)
        for (i=0;i<out.numberoftriangles;i++) 
          (*cellregions)(i)=(int)out.triangleattributelist[i];
      else
        for (i=0;i<out.numberoftriangles;i++) 
          (*cellregions)(i)=0;
    
      if (in.triangleattributelist) triangle::trifree(in.triangleattributelist);
      if (out.triangleattributelist) triangle::trifree(out.triangleattributelist);
    }
        
    if (out.segmentlist)
    {
      bfaces=std::make_shared<numcxx::TArray2<int>>(out.numberofsegments,2,out.segmentlist,[](int*p){triangle::trifree(p);});
      bfaceregions=std::make_shared<numcxx::TArray1<int>>(out.numberofsegments,out.segmentmarkerlist,[](int*p){triangle::trifree(p);});
    }
    else
    {
      throw std::runtime_error("numcxx: triangulate: Missing segment list"); 
    }

  }

Here is the call graph for this function:

Member Function Documentation

static std::shared_ptr<SimpleGrid> numcxx::SimpleGrid::create ( const Geometry &  geometry, const std::string  triangle_flags  )

inlinestatic

Static constructor from reference to geometry.

Examples:

30-triangle.cxx, and 31-cut.cxx.

Definition at line 35 of file simplegrid.hxx.

{ return std::make_shared<SimpleGrid>(geometry,triangle_flags.c_str());}

static std::shared_ptr<SimpleGrid> numcxx::SimpleGrid::create ( std::shared_ptr< Geometry >  geometry, const char *  triangle_flags  )

inlinestatic

Static constructor from shared pointer to geometry.

Definition at line 38 of file simplegrid.hxx.

{ return std::make_shared<SimpleGrid>(*geometry,triangle_flags);}

const TArray2<double>& numcxx::SimpleGrid::get_points ( ) const

inline

Get array of point coordinates.

Examples:

42-convtest-fem-sin.cxx, and 43-convtest-fvm-sin.cxx.

Definition at line 42 of file simplegrid.hxx.

{ return *points;}; 

const TArray2<int>& numcxx::SimpleGrid::get_cells ( ) const

inline

Get array of point indices describing cells.

Definition at line 45 of file simplegrid.hxx.

{ return *cells;};

const TArray1<int>& numcxx::SimpleGrid::get_cellregions ( ) const

inline

Get array region markers.

Definition at line 48 of file simplegrid.hxx.

{return  *cellregions;};

const TArray2<int>& numcxx::SimpleGrid::get_bfaces ( ) const

inline

Get array of point indices describing boundary faces.

Definition at line 51 of file simplegrid.hxx.

{ return *bfaces;};

const TArray1<int>& numcxx::SimpleGrid::get_bfaceregions ( ) const

inline

Get array of boundary markers.

Definition at line 54 of file simplegrid.hxx.

{ return *bfaceregions;};

const int numcxx::SimpleGrid::spacedim ( ) const

inline

Return dimension of space.

Definition at line 57 of file simplegrid.hxx.

{ return points->shape(1);}

const int numcxx::SimpleGrid::griddim ( ) const

inline

Return dimension of grid.

Definition at line 60 of file simplegrid.hxx.

{ return bfaces->shape(1);}

const int numcxx::SimpleGrid::ncells ( ) const

inline

Return number of cells.

Definition at line 63 of file simplegrid.hxx.

{return cells->shape(0);}

const int numcxx::SimpleGrid::npoints ( ) const

inline

Return number of points.

Examples:

40-stationary-heat-fe.cxx, 41-stationary-heat-fv.cxx, 42-convtest-fem-sin.cxx, 43-convtest-fvm-sin.cxx, 44-transient-heat-fe.cxx, and 46-nonlin-fvm.cxx.

Definition at line 66 of file simplegrid.hxx.

{return points->shape(0);}

const int numcxx::SimpleGrid::nbfaces ( ) const

inline

Return number of boundary faces.

Definition at line 69 of file simplegrid.hxx.

{return bfaces->shape(0);}

Here is the call graph for this function:

void numcxx::SimpleGrid::calc_hminmax	(	double &	hmin,
		double &	hmax
	)		const

Calculate some grid data.

Examples:

42-convtest-fem-sin.cxx, and 43-convtest-fvm-sin.cxx.

Definition at line 55 of file simplegrid.cxx.

  {
    hmax=0.0;
    hmin=1.0e100;
    auto cells=*SimpleGrid::cells;
    auto points=*SimpleGrid::points;

    for (int icell=0; icell<ncells();icell++)
    {
      double h=dist(points, cells(icell,0), cells(icell,1));
      hmin=std::min(h,hmin);
      hmax=std::max(h,hmax);
   
      h=dist(points, cells(icell,1), cells(icell,2));
      hmin=std::min(h,hmin);
      hmax=std::max(h,hmax);

      h=dist(points, cells(icell,0), cells(icell,2));
      hmin=std::min(h,hmin);
      hmax=std::max(h,hmax);
    }

  }

Here is the call graph for this function:

Member Data Documentation

std::shared_ptr<TArray2<double> > numcxx::SimpleGrid::points =nullptr

private

Definition at line 75 of file simplegrid.hxx.

std::shared_ptr<TArray2<int> > numcxx::SimpleGrid::cells =nullptr

private

Definition at line 76 of file simplegrid.hxx.

std::shared_ptr<TArray1<int> > numcxx::SimpleGrid::cellregions =nullptr

private

Definition at line 77 of file simplegrid.hxx.

std::shared_ptr<TArray2<int> > numcxx::SimpleGrid::bfaces =nullptr

private

Definition at line 78 of file simplegrid.hxx.

std::shared_ptr<TArray1<int> > numcxx::SimpleGrid::bfaceregions =nullptr

private

Definition at line 79 of file simplegrid.hxx.

---

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

• simplegrid.hxx
• simplegrid.cxx