MappedFEFaceValues Class Reference

#include <mappedfefacevalues.h>

Inheritance diagram for MappedFEFaceValues:

[legend]

Collaboration diagram for MappedFEFaceValues:

[legend]

List of all members.

Public Member Functions
	MappedFEFaceValues (FiniteElementInterface &fe, Index faceId)
	~MappedFEFaceValues ()
void	setPoints (const ArrayOfVecDouble &v, const VecDouble &weights)
void	setQuadrature (const dealii::Quadrature< 2 > &quad)
void	compute (MappingOrthoMesh &map, unsigned compute_flags)
virtual void	compute_values ()
virtual void	compute_grads (MappingOrthoMesh &map)
virtual void	compute_J (MappingOrthoMesh &map)
virtual void	compute_JxW (MappingOrthoMesh &map, const VecDouble &weights)
virtual void	compute_JxW (MappingOrthoMesh &map)
Protected Attributes
unsigned	_dir
VecIndex	_face_to_cell_map
VecDouble	e1
VecDouble	e2

---

Detailed Description

This class specializes the MappedFEValues to read values from the faces. It is used mainly to compute Neumman conditions at the boundaries and it is supposed to be used in the same way as the MappedFEValues.

There are some minor differences though. First in the setPts() method you pass actually the points in 2D and the face. These points are then mapped in the reference domain in 3D to the points in the exact face. You have a method to get the jacobian to make integrations on manifolds called J_Faces()

Definition at line 22 of file mappedfefacevalues.h.

---

Constructor & Destructor Documentation

MappedFEFaceValues::MappedFEFaceValues	(	FiniteElementInterface &	fe,
		Index	faceId
	)

Definition at line 7 of file mappedfefacevalues.cpp.

  :MappedFEValuesBase(fe,fe.face_to_cell_local_map(faceId).size())
{
  _dir = faceId;//Store the face of the face values

  //Store the mapping between the local face dofs
  //and the local cell dofs. This is obtained
  //fron the finite element.
  _face_to_cell_map= fe.face_to_cell_local_map(faceId);

  UnitCube ref;
  ref.face_base_vectors(faceId,&e1,&e2);

  
}

MappedFEFaceValues::~MappedFEFaceValues

(

)

[inline]

Definition at line 34 of file mappedfefacevalues.h.

{}

---

Member Function Documentation

void MappedFEFaceValues::compute	(	MappingOrthoMesh &	map,
		unsigned	compute_flags
	)

void MappedFEFaceValues::compute_grads

(

MappingOrthoMesh &

map

)

[virtual]

Compute the grads of the bases functions

Definition at line 69 of file mappedfefacevalues.cpp.

{
  throw new Exception("MappedFEFaceValues::compute_grads()  not implemented yet\n");
}

void MappedFEFaceValues::compute_J

(

MappingOrthoMesh &

map

)

[virtual]

Compute jacobians of the points list

Definition at line 76 of file mappedfefacevalues.cpp.

{
  static Matrix GradT(3,3); //A matrix 3x3
  VecDoubleRef pt(3);
  static VecDouble aux1(3);
  static VecDouble aux2(3);

  for(unsigned i=0;i<_pts.size();i++)
  {
    //get pt i
    _pts.getVecValues(i,&pt); //Get the point

    //get GradT
    map.T_grad(pt,GradT);

    GradT.vmult(aux1,e1);
    GradT.vmult(aux2,e2);
    
    _J(i) = NumericMethods::abs_vector_product(aux1,aux2);
  }

}

void MappedFEFaceValues::compute_JxW

(

MappingOrthoMesh &

map

)

[virtual]

Definition at line 115 of file mappedfefacevalues.cpp.

{
  compute_JxW(map,_weights);
}

void MappedFEFaceValues::compute_JxW	(	MappingOrthoMesh &	map,
		const VecDouble &	weights
	)			[virtual]

Definition at line 100 of file mappedfefacevalues.cpp.

{
  assert(weights.size() == _pts.size());
  compute_J(map);
  _JxW.term_mult(_J,weights);
}

void MappedFEFaceValues::compute_values

(

)

[virtual]

Compute the values of the bases function

Definition at line 49 of file mappedfefacevalues.cpp.

{
  unsigned iEnd = _pts.size();
  VecDoubleRef pt(3);
  unsigned n_local_dofs = _face_to_cell_map.size();

  for (unsigned i=0;i<iEnd;i++) //for all points in _pts list
  {
    _pts.getVecValues(i,&pt); //Get the point
    for (unsigned j=0;j<n_local_dofs;j++) //for all dofs
    {
      //evaluat the base function j on the point pt
      //and store in the matrix this->_Values;
      _Values(i,j)=_fe.shape_value(_face_to_cell_map[j],pt);
    }
  }

}

void MappedFEFaceValues::setPoints	(	const ArrayOfVecDouble &	v,
		const VecDouble &	weights
	)

Set the list of points

Reimplemented from MappedFEValuesBase.

Definition at line 23 of file mappedfefacevalues.cpp.

{
  
  //The array of points must be 2D points
  assert(v.vecs_size() ==2);

  ArrayOfVecDouble pts;
  pts.reinit(v.size(),3);

  VecDoubleRef pt3D(3);
  VecDouble    pt2D(2);
  UnitCube refDomain;

  for (unsigned i=0;i<v.size();i++)
  {
    pts.getVecValues(i,&pt3D);
    v.copyVecValues(i,&pt2D);
    refDomain.project_2d_points_into_faces(pt2D,&pt3D,static_cast<FaceDirection3D>(_dir));
  }
  pts.print();
  MappedFEValuesBase::setPoints(pts,weights);
}

void MappedFEFaceValues::setQuadrature

(

const dealii::Quadrature< 2 > &

quad

)

Definition at line 108 of file mappedfefacevalues.cpp.

{
  //printf("Caralho\n");
  MappedFEFaceValues::setPoints(quad.get_points(),quad.get_weights());
}

---

Member Data Documentation

unsigned MappedFEFaceValues::_dir [protected]

Definition at line 27 of file mappedfefacevalues.h.

VecIndex MappedFEFaceValues::_face_to_cell_map [protected]

Definition at line 28 of file mappedfefacevalues.h.

VecDouble MappedFEFaceValues::e1 [protected]

Definition at line 29 of file mappedfefacevalues.h.

VecDouble MappedFEFaceValues::e2 [protected]

Definition at line 29 of file mappedfefacevalues.h.

---

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

• mappedfefacevalues.h
• mappedfefacevalues.cpp