MappedFEValues Class Reference

#include <mappedfevalues.h>

Inheritance diagram for MappedFEValues:

[legend]

Collaboration diagram for MappedFEValues:

[legend]

List of all members.

Public Member Functions
	MappedFEValues (FiniteElementInterface &fe)
void	setPoints (const dealii::Quadrature< 3 > &quad)
void	setPoints (const ArrayOfVecDouble &pts, const VecDouble &weights)
void	compute (MappingOrthoMesh &map, unsigned compute_flags)
void	compute_values ()
void	compute_grads (MappingOrthoMesh &map)
void	compute_J (MappingOrthoMesh &map)
void	compute_JxW (MappingOrthoMesh &map, const VecDouble &weights)
void	compute_JxW (MappingOrthoMesh &map)
	~MappedFEValues ()

---

Detailed Description

The sole purpose of this class is to get values of the FE in the deformed cell.

The user first set the points (on the reference domain) for which the FE values will be computed and then call the many compute methods of this class like computeShapeValues() and computeShapeGrad(), etc.

After that, the user can get the values computed through methods like shape_value(Index pt,Index i) or shape_grad(Index pt,Index i), etc..... where pt is the index in the vector of points passed to this class and i the base function associated with shape function i of the element.

Definition at line 30 of file mappedfevalues.h.

---

Constructor & Destructor Documentation

MappedFEValues::MappedFEValues

(

FiniteElementInterface &

fe

)

Definition at line 95 of file mappedfevalues.cpp.

  :MappedFEValuesBase(fe,fe.n_local_dofs())
 {

 }

MappedFEValues::~MappedFEValues

(

)

[inline]

Definition at line 52 of file mappedfevalues.h.

{}

---

Member Function Documentation

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

Definition at line 102 of file mappedfevalues.cpp.

{
  m_compute_flags=compute_flags;
  if (compute_flags | COMPUTE_VALUES)
    compute_values();
  if  (compute_flags | COMPUTE_GRADS)
    compute_grads(map);
  if  (compute_flags | COMPUTE_PTS)
    compute_points(map);
  if  (compute_flags | COMPUTE_JxW)
    compute_JxW(map,_weights);
  else if (compute_flags | COMPUTE_J)
    compute_J(map);

}

void MappedFEValues::compute_grads

(

MappingOrthoMesh &

map

)

Compute the gras of the bases functions

Definition at line 34 of file mappedfevalues.cpp.

{
  unsigned iEnd = _pts.size();
  unsigned n_local_dofs = _fe.n_local_dofs();
  VecDoubleRef pt(3);
  static VecDouble vGradRef(3);
  Matrix invTGrad(3,3);
  for (unsigned iPt=0;iPt<iEnd;iPt++) //for all points in _pts list
  {
    _pts.getVecValues(iPt,&pt); //Get the point
    map.Inv_T_Grad(pt,invTGrad);

    for (unsigned j=0;j<n_local_dofs;j++) //for all dofs
    {
      _fe.shape_grad(j,pt,vGradRef);
      invTGrad.Tvmult(_Grad[iPt*m_n_dofs + j],vGradRef);
    }
  }
}

void MappedFEValues::compute_J

(

MappingOrthoMesh &

map

)

Compute jacobians of the points list

Definition at line 57 of file mappedfevalues.cpp.

{
  static Matrix GradT(3,3); //A matrix 3x3
  VecDoubleRef pt(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);
    _J(i)=NumericMethods::det3x3(GradT);
  }
}

void MappedFEValues::compute_JxW

(

MappingOrthoMesh &

map

)

Definition at line 87 of file mappedfevalues.cpp.

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

void MappedFEValues::compute_JxW	(	MappingOrthoMesh &	map,
		const VecDouble &	weights
	)

Compute the jacobians calling the compute_J() method and the calculate the jacobians multiplied by weights.

Parameters:

weights

Returns:

Definition at line 80 of file mappedfevalues.cpp.

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

void MappedFEValues::compute_values

(

)

Compute the values of the basis function

Compute the values of the bases functions in the deformed domain

Definition at line 15 of file mappedfevalues.cpp.

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

  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
    {
      //evaluates the base function j on the point pt
      //and store in the matrix this->_Values;
      _Values(i,j)=_fe.shape_value(j,pt);
    }
  }
}

void MappedFEValues::setPoints	(	const ArrayOfVecDouble &	pts,
		const VecDouble &	weights
	)

Reimplemented from MappedFEValuesBase.

Definition at line 119 of file mappedfevalues.cpp.

{
  MappedFEValuesBase::setPoints(pts,weights);
}

void MappedFEValues::setPoints

(

const dealii::Quadrature< 3 > &

quad

)

Definition at line 5 of file mappedfevalues.cpp.

{
  MappedFEValuesBase::setPoints(quad.get_points(),quad.get_weights());
}

---

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

• mappedfevalues.h
• mappedfevalues.cpp