GRINS::ComputeDistanceResidual Class Reference

#include <distance_function.h>

Public Member Functions
	ComputeDistanceResidual (const libMesh::Elem *belem, const libMesh::Point *point)
	~ComputeDistanceResidual ()
void	operator() (const libMesh::DenseVector< libMesh::Real > &U, libMesh::DenseVector< libMesh::Real > &F)

Private Attributes
const libMesh::Elem &	_belem
const unsigned int	_dim
const libMesh::Point &	_p
libMesh::AutoPtr< libMesh::FEBase >	_fe
libMesh::FEBase *	fe
libMesh::DenseVector< libMesh::Real >	Up
libMesh::DenseVector< libMesh::Real >	Fp
libMesh::DenseVector< libMesh::Real >	Um
libMesh::DenseVector< libMesh::Real >	Fm

Detailed Description

This class provides the functionality to compute the residual required when computing the distance function. The residual here corresponds to the derivative of the distance between the field point and the boundary point wrt the boundary parameterization. When this residual is zero, the distance is minimized, and we have found the point we're looking for.

Definition at line 188 of file distance_function.h.

Constructor & Destructor Documentation

GRINS::ComputeDistanceResidual::ComputeDistanceResidual ( const libMesh::Elem *  belem, const libMesh::Point *  point  )

inline

Definition at line 193 of file distance_function.h.

References _belem.

                                                                                 :
      _belem(*belem),
      _dim(_belem.dim()+1),
      _p(*point),
      _fe (libMesh::FEBase::build(_belem.dim(), libMesh::FEType(libMesh::FIRST, libMesh::LAGRANGE))),
      fe  (_fe.get())
    {
      // only supporting QUAD4 this way for now
      libmesh_assert(_belem.type()==libMesh::QUAD4);
    }

GRINS::ComputeDistanceResidual::~ComputeDistanceResidual ( )

inline

Definition at line 205 of file distance_function.h.

{}

Member Function Documentation

void GRINS::ComputeDistanceResidual::operator() ( const libMesh::DenseVector< libMesh::Real > &  U, libMesh::DenseVector< libMesh::Real > &  F  )

inline

Definition at line 208 of file distance_function.h.

References _belem, _dim, _p, and fe.

    {

      libmesh_assert(U.size()==_dim-1);
      libmesh_assert(F.size()==_dim-1);

      libMesh::DenseVector<libMesh::Real> xx(_dim);
      xx.zero();

      libMesh::DenseMatrix<libMesh::Real> xx_U(_dim, _dim-1);
      xx_U.zero();

      std::vector<libMesh::Point> xi(1);

      xi[0](0) = U(0);
      xi[0](1) = U(1);

      // grab basis functions (evaluated at qpts)
      const std::vector<std::vector<libMesh::Real> > &basis = fe->get_phi();
      const std::vector<std::vector<libMesh::RealGradient> > &dbasis = fe->get_dphi();

      // reinitialize finite element data at xi
      fe->reinit(&_belem, &xi);

      // interpolate location and derivatives
      for (unsigned int inode=0; inode<_belem.n_nodes(); ++inode)
        for (unsigned int idim=0; idim<_dim; ++idim)
          {
            xx(idim) += _belem.point(inode)(idim) * basis[inode][0];

            for (unsigned int jdim=0; jdim<_dim-1; ++jdim)
              xx_U(idim, jdim) += _belem.point(inode)(idim) * dbasis[inode][0](jdim);
          }


      // form the residual
      F.zero();
      for (unsigned int ires=0; ires<_dim-1; ++ires)
        for (unsigned int idim=0; idim<_dim; ++idim)
          F(ires) += ( xx(idim) - _p(idim) ) * xx_U(idim, ires);


    }

Member Data Documentation

const libMesh::Elem& GRINS::ComputeDistanceResidual::_belem

private

Definition at line 256 of file distance_function.h.

Referenced by ComputeDistanceResidual(), and operator()().

const unsigned int GRINS::ComputeDistanceResidual::_dim

private

Definition at line 257 of file distance_function.h.

Referenced by operator()().

libMesh::AutoPtr<libMesh::FEBase> GRINS::ComputeDistanceResidual::_fe

private

Definition at line 260 of file distance_function.h.

const libMesh::Point& GRINS::ComputeDistanceResidual::_p

private

Definition at line 258 of file distance_function.h.

Referenced by operator()().

libMesh::FEBase* GRINS::ComputeDistanceResidual::fe

private

Definition at line 261 of file distance_function.h.

Referenced by operator()().

libMesh::DenseVector<libMesh::Real> GRINS::ComputeDistanceResidual::Fm

private

Definition at line 264 of file distance_function.h.

libMesh::DenseVector<libMesh::Real> GRINS::ComputeDistanceResidual::Fp

private

Definition at line 264 of file distance_function.h.

libMesh::DenseVector<libMesh::Real> GRINS::ComputeDistanceResidual::Um

private

Definition at line 264 of file distance_function.h.

libMesh::DenseVector<libMesh::Real> GRINS::ComputeDistanceResidual::Up

private

Definition at line 264 of file distance_function.h.

---

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

• src/utilities/include/grins/distance_function.h