DealBase Class Reference

#include <dealbase.h>

Inheritance diagram for DealBase:

[legend]

List of all members.

Public Member Functions
	DealBase ()
	DealBase (const Triangulation< DIM > &tria)
	~DealBase ()
Triangulation< DIM > &	getTriangulation () const
Cell	getCell (int index)
Cell	endCell () const
Cell	beginCell () const
Face	beginFace ()
Face	endFace ()
void	getFaceBarycenter (Face &face, Point3D &p)
Point3D	getFaceBarycenter (Face &face)
double	getCentralValue (const Cell &cell, const VecDouble &values) const
double	getFaceValue (const Face &face, const VecDouble &values)
double	getFaceValue (const Cell &cell, FaceDirection3D dir, const VecDouble &values)
unsigned	getFaceValueIndex (const Cell &cell, FaceDirection3D dir)
bool	getFaceFlag (const Face &face, const VecBool &values)
Cell	getCellAtPoint (const Point< DIM > &p)
double	getNeighborValue (const Cell &cell, CellDirection3D dir, const VecDouble &cValues)
double	getVerticeValue (const Cell &cell, VertexDirection3D vertex, const VecDouble &values)
double	getVerticeValue (unsigned vIndex, const VecDouble &values)
unsigned	getCentralPointIndex (Cell &cell)
void	getVerticePoint (unsigned vIndex, Point< DIM > &X) const
Cell	getEndCell ()
Cell	getAdjCell (Cell &cell, CellDirection3D dir)
Cell	getAdjCell (Cell &cell, CellDirection3D dir1, CellDirection3D dir2)
bool	getCellFlag (const Cell &cell, const VecBool &cVFlags) const
bool	getVerticeFlag (const Cell &cell, VertexDirection3D dir, const VecBool &vVFlags) const
void	setCentralValue (const Cell &cell, double dd, VecDouble &values)
void	setCentralValue (DoFCell &cell, double dd, VecDouble &values)
void	setFaceValue (const Face &face, double dd, VecDouble &values)
void	setFaceValue (const Cell &cell, FaceDirection3D dir, double dd, VecDouble &values)
void	setFaceFlag (const Face &face, bool b, VecBool &values)
void	setCellFlag (const Cell &cell, const bool flag, VecBool &cVFlags) const
void	setVerticeFlag (const Cell &cell, VertexDirection3D dir, const bool flag, VecBool &vVFlags) const
double	getL2NormAtCells (const VecDouble &cValues) const
double	getIntegralAtCells (const VecDouble &cValues) const
unsigned	numCells () const
unsigned	numFaces () const
unsigned	numVertices () const
double	faceArea (Face face)
int	getAdjCellIndex (unsigned iCell, CellDirection3D dir)
double	getAdjCellValue (unsigned iCell, CellDirection3D dir, const VecDouble &v)
int	getFaceIndex (unsigned iCell, FaceDirection3D dir)
double	getFaceValue (unsigned iCell, FaceDirection3D dir, const VecDouble &v)
bool	isCellAtBoundary (unsigned iCell)
bool	hasAdjCell (unsigned iCell, CellDirection3D dir)
void	setCentralValuesFromFunction (Function3D &f, VecDouble &cValues, unsigned component=0)
void	setCentralValuesInFunctionDomain (Function3D &f, MapIntDouble &values, unsigned component=0)
void	setVerticesValuesFromFunction (Function< DIM > &f, VecDouble &vValues, unsigned deg=0)
void	setFacesValuesFromFunction (Function3D &f, VecDouble &fValues, unsigned deg=0, bool bOnlyPrescribedBoundary=false)
void	setFacesValuesFromFunction (Function3D &f, Matrix &MValues, bool bOnlyPrescribedBoundary=false)
void	tagFacesInDomain (Function3D &f, VecTag &tags, char value)
void	tagFacesInDomain (Function3D &f, VecBool &tags)
void	setVerticeValue (Cell &cell, VertexDirection3D vertex, double dd, VecDouble &values)
void	setMaterialIdOfPrescribedBoundaryCells (GeneralFunctionInterface &f, unsigned material_id)
void	setMaterialIdFromFunctionDomain (GeneralFunctionInterface &f, unsigned material_id)
void	getCellsInFunctionDomain (Function3D &f, VecIndex &vec)
void	printMatrixMaple (FullMatrix< double > &M)
void	printMatrixMaple (SparseMatrix< double > &M)
void	printVertices ()
void	printCells ()
void	printFaces ()
void	printPoint (const Point< DIM > &P, std::ostream &out)
void	printVector (const VecDouble &vSol)
void	printPoint (const Point< DIM > &p)
void	printCentralValues (const VecDouble &cValues)
void	printCentralValues (const VecDouble &cValues1, const VecDouble &cValues2)
void	printCentralValues (const VecDouble &cValues1, const VecDouble &cValues2, const VecDouble &cValues3, const VecDouble &cValues4)
void	printFacesValuesAtBoundary (const VecDouble &fValues, const VecDouble &fValues2)
void	printFacesValuesAtBoundary (const VecDouble &fValues1, const VecDouble &fValues2, const VecDouble &fValues3, const VecDouble &fValues4)
void	printFacesValuesAtBoundary (const VecDouble &fValues)
void	printCell (Cell cell, std::ostream &out=std::cout)
void	printCells (VecIndex &vec, std::ostream &out=std::cout)
bool	isFirstCell (Cell &cell)
bool	isValid (Cell &cell)
bool	hasNeighbor (Cell &cell, CellDirection3D dir)
void	flagPrescribedFaces (Function3D &f, VecBool &vFlags)
void	setPrescribedBoundaryId (Function3D &f, unsigned boundary_id)
double	getDX (Cell &cell)
double	getDZ (Cell &cell)
double	getDY (Cell &cell)
Static Public Member Functions
static void	printTriangulation (Triangulation< 3 > &tria, std::ostream &out=std::cout)
Static Public Attributes
static const int	INVALID_INDEX = -1
Protected Member Functions
void	setTriangulation (Triangulation< DIM > &trig)
Protected Attributes
int	n_cells
Private Attributes
Triangulation< DIM > *	m_pTriangulation
Cell	END_CELL

---

Detailed Description

MeshDeal3D This class is a Wrapper providing a class Mesh used by the finite volume methods. In background this class call methods from the class Triangulation of the library deal.II

Definition at line 25 of file dealbase.h.

---

Constructor & Destructor Documentation

DealBase::DealBase

(

)

[inline]

Definition at line 38 of file dealbase.h.

{m_pTriangulation = NULL;n_cells=-1;}

DealBase::DealBase

(

const Triangulation< DIM > &

tria

)

[inline]

Definition at line 39 of file dealbase.h.

{setTriangulation((Triangulation<DIM>&) tria);}

DealBase::~DealBase

(

)

[inline]

Definition at line 40 of file dealbase.h.

{}

---

Member Function Documentation

double DealBase::faceArea

(

Face

face

)

Definition at line 902 of file dealbase.cpp.

{
  return NumericMethods::squareMeasure(face->vertex(BL_VERTEX),
                                       face->vertex(BR_VERTEX),
                                       face->vertex(UR_VERTEX),
                                       face->vertex(UL_VERTEX));
}

void DealBase::flagPrescribedFaces	(	Function3D &	f,
		VecBool &	vFlags
	)

Definition at line 427 of file dealbase.cpp.

{
  assert(vFlags.size() == numFaces());
  for (Face face=beginFace();face!=endFace();face++)
  {
    Point3D p = face->center();
    setFaceFlag(face,face->at_boundary() && (f.isInDomain(p)),vFlags);
  }
}

double DealBase::getDX

(

Cell &

cell

)

[inline]

Definition at line 157 of file dealbase.h.

{return  cell->vertex(VERTEX_100)(0)-cell->vertex(VERTEX_000)(0);}

double DealBase::getDY

(

Cell &

cell

)

[inline]

Definition at line 159 of file dealbase.h.

{return  cell->vertex(VERTEX_010)(1)-cell->vertex(VERTEX_000)(1);}

double DealBase::getDZ

(

Cell &

cell

)

[inline]

Definition at line 158 of file dealbase.h.

{return  cell->vertex(VERTEX_001)(2)-cell->vertex(VERTEX_000)(2);}

double DealBase::getIntegralAtCells

(

const VecDouble &

cValues

)

const

Reimplemented in DealOrthoMesh.

double DealBase::getL2NormAtCells

(

const VecDouble &

cValues

)

const

Reimplemented in DealOrthoMesh.

Triangulation<DIM>& DealBase::getTriangulation

(

)

const [inline]

Definition at line 41 of file dealbase.h.

{assert(m_pTriangulation);return *m_pTriangulation;}

bool DealBase::hasNeighbor	(	Cell &	cell,
		CellDirection3D	dir
	)			[inline]

Definition at line 154 of file dealbase.h.

{return (cell->neighbor_index(dir) != DealBase::INVALID_INDEX);}

bool DealBase::isFirstCell

(

Cell &

cell

)

[inline]

Definition at line 152 of file dealbase.h.

{return cell == getTriangulation().begin_active();}

bool DealBase::isValid

(

Cell &

cell

)

[inline]

Definition at line 153 of file dealbase.h.

{return (cell.state() == IteratorState::valid);} 

unsigned DealBase::numCells

(

)

const [inline]

Definition at line 87 of file dealbase.h.

{return n_cells;}

unsigned DealBase::numFaces

(

)

const [inline]

Definition at line 88 of file dealbase.h.

{return getTriangulation().n_active_faces();}

unsigned DealBase::numVertices

(

)

const [inline]

Definition at line 89 of file dealbase.h.

{return getTriangulation().n_vertices();}

void DealBase::printCell	(	Cell	cell,
		std::ostream &	out = std::cout
	)

Definition at line 597 of file dealbase.cpp.

{
  char str[500];
  sprintf(str,"==== Cell %d ====\n",cell->index());
  out << str;
  out << "Barycenter: " << cell->barycenter() << endl;
  //     cout << "0 = " <<  cell->vertex(0) << endl;
  //     cout << "1 = " <<  cell->vertex(1) << endl;
  //     cout << "2 = " <<  cell->vertex(2) << endl;
  //     cout << "3 = " <<  cell->vertex(3) << endl;
  //     cout << "4 = " <<  cell->vertex(4) << endl;
  //     cout << "5 = " <<  cell->vertex(5) << endl;
  //     cout << "6 = " <<  cell->vertex(6) << endl;
  //     cout << "7 = " <<  cell->vertex(7) << endl;

  out << "V000 = "  << cell->vertex_index(VERTEX_000) <<  ") " << cell->vertex(VERTEX_000)  << endl;
  out << "V100 = " << cell->vertex_index(VERTEX_100) <<  ") " << cell->vertex(VERTEX_100) << endl;
  out << "V110 = " << cell->vertex_index(VERTEX_110) <<  ") " << cell->vertex(VERTEX_110) << endl;
  out << "V010 = " << cell->vertex_index(VERTEX_010) <<  ") " << cell->vertex(VERTEX_010) << endl;
  out << "V001 = " << cell->vertex_index(VERTEX_001) <<  ") " << cell->vertex(VERTEX_001) << endl;
  out << "V101 = " << cell->vertex_index(VERTEX_101) <<  ") " << cell->vertex(VERTEX_101) << endl;
  out << "V111 = " << cell->vertex_index(VERTEX_111) <<  ") " << cell->vertex(VERTEX_111) << endl;
  out << "V011 = " << cell->vertex_index(VERTEX_011) <<  ") " << cell->vertex(VERTEX_011) << endl;

  out << "RIGHT_FACE  " << cell->face_index (RIGHT_FACE) << ")" << endl;
  out << "LEFT_FACE   " << cell->face_index  (LEFT_FACE) << ")"  << endl;
  out << "FRONT_FACE  " << cell->face_index (FRONT_FACE) << ")" << endl;
  out << "BACK_FACE   " << cell->face_index  (BACK_FACE) << ")"  << endl;
  out << "BOTTOM_FACE " << cell->face_index(BOTTOM_FACE) << ")"<< endl;
  out << "UP_FACE     " << cell->face_index    (UP_FACE) << ")"    << endl;


}

void DealBase::printCells	(	VecIndex &	vec,
		std::ostream &	out = std::cout
	)

Definition at line 966 of file dealbase.cpp.

{
  for (unsigned i = 0;i<vec.size();i++)
  {
    out <<  i << ")";
    printCell(getCell(vec[i]),out);
  }
}

void DealBase::setPrescribedBoundaryId	(	Function3D &	f,
		unsigned	boundary_id
	)

Objetivo: Configurar o ID das faces das celulas localizadas na fronteira e que tem seu valor prescrito pela funcao f. Em outras palavras quando f.isInDomain(face) == PRESCRIBED_BOUNDARY;

Parametros: f = Funcao que designa fronteira. boundary_id = Valor a ser adicionado na face na fronteira. Retorno:

Definition at line 266 of file dealbase.cpp.

{
  Cell cell = getTriangulation().begin_active();
  Cell endc = getTriangulation().end();

  for (;cell!=endc;cell++)
  {
    for (unsigned face_no=0;face_no < GeometryInfo<DIM>::faces_per_cell;face_no++)
    {
      Face face = cell->face(face_no);
      if (face->at_boundary())
      {
        bool markFace = true;
        for (unsigned ivertex=0;ivertex< GeometryInfo<DIM>::vertices_per_face;ivertex++)
        {
          if (!f.isInDomain(face->vertex(ivertex)))
          {
            markFace = false;
            break;
          }
        }
        if (markFace)
        {
          face->set_boundary_indicator(boundary_id);
        }
      }
    }
  }
}

void DealBase::setTriangulation

(

Triangulation< DIM > &

trig

)

[protected]

Objetivo: Configurar o membro m_pTriangulation.

Parametros:

Retorno:

Definition at line 14 of file dealbase.cpp.

{
  m_pTriangulation = &trig;
  END_CELL = getTriangulation().end();
  n_cells=getTriangulation().n_active_cells();
}

---

Member Data Documentation

Cell DealBase::END_CELL [private]

Definition at line 29 of file dealbase.h.

const int DealBase::INVALID_INDEX = -1 [static]

Definition at line 35 of file dealbase.h.

Triangulation<DIM>* DealBase::m_pTriangulation [private]

Definition at line 28 of file dealbase.h.

int DealBase::n_cells [protected]

Definition at line 32 of file dealbase.h.

---

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

• dealbase.h
• dealbase.cpp