A 2D element map for a curved quadrilateral. More...

#include <elementMaps.hh>

Inheritance diagram for concepts::BlendingQuad2d:

Collaboration diagram for concepts::BlendingQuad2d:

Public Member Functions
	BlendingQuad2d (const MappingEdge2d edgemap0, const MappingEdge2d edgemap1, const MappingEdge2d edgemap2, const MappingEdge2d edgemap3=0, const Real2d vtx0=0, const Real2d vtx1=0, const Real2d vtx2=0, const Real2d vtx3=0)
	Constructor.
	BlendingQuad2d (const BlendingQuad2d &v)
	Copy constructor.
virtual BlendingQuad2d *	clone () const
	Returns a copy of the map.
virtual MappingEdge2d *	edge (const uint edge) const
	Returns a copy of the edge mapping of one edge.
virtual MapReal2d	hessian (uint i, const Real x, const Real y) const
	Returns the Hessian in a 2D linear map.
virtual MapReal2d	jacobian (const Real x, const Real y) const
	Returns the Jacobian in a 2D linear map.
virtual Real	jacobianDeterminant (const Real x, const Real y) const
	Returns determinant of the Jacobian.
virtual MapReal2d	jacobianInverse (const Real x, const Real y) const
	Returns the inverse of the Jacobian in a 2D linear map.
virtual Real	lineElement (const Real x, const uint edge) const
	Returns factor of differential element integrating over an edge.
virtual Real2d	operator() (Real x, Real y) const
	Returns a point in 2D mapped from the unit square [0,1]² onto the element in the original mesh.
virtual MappingQuad2d *	part (const Real2d x0, const Real2d y0) const
	Returns a new object of an element map for a rectangular part of the reference cell defined by the two points.
virtual bool	straight () const
	~BlendingQuad2d ()
Protected Member Functions
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
Private Member Functions
void	testJacobian_ () const
	Tests if the jacobian is positive.
Private Attributes
const MappingEdge2d *	edgemap_ [4]
	Map of the edges.
Real2d	vtx_ [4]
	Vector that contains the vertices.

Detailed Description

A 2D element map for a curved quadrilateral.

The mapping is defined by the Linear Blending Function Method (proposed by Gordon and Hall, 1971). Therefor the mapping of the edges has to be given.

Test:: test::TestBlending2D

Author:: , Kersten Schmidt, 2004

Definition at line 852 of file elementMaps.hh.

Constructor & Destructor Documentation

concepts::BlendingQuad2d::BlendingQuad2d	(	const MappingEdge2d *	edgemap0,
		const MappingEdge2d *	edgemap1,
		const MappingEdge2d *	edgemap2,
		const MappingEdge2d *	edgemap3 = `0`,
		const Real2d *	vtx0 = `0`,
		const Real2d *	vtx1 = `0`,
		const Real2d *	vtx2 = `0`,
		const Real2d *	vtx3 = `0`
	)

Constructor.

Pointers to edge mappings and vertex coordinates can be given. There has to be given at least

two opposite edges
two adjacent edges and a the other vertex
an edge and the two other vertices
four vertices

If more edges or vertices are given, the particular coordinates of the vertices has to match.

The edges or vertices have to be given counter clockwise. The first vertex of the edgemap0 is the first vertex of the quad, the second in the edge map the second ...

The direction of an edge has to be changed by its method inverse(), e.g.

std::auto_ptr<MappingEdge2d> invEdge = edge.inverse();
        BlendingQuad2d quad(invEdge, ...);

The class holds own copies of the given edge mappings and coordinates of the vertices.

First edge is $\xi_2=0$ , where $\xi_2$ is the second local coordinate.

concepts::BlendingQuad2d::BlendingQuad2d ( const BlendingQuad2d & v )

Copy constructor.

concepts::BlendingQuad2d::~BlendingQuad2d ( )

Member Function Documentation

virtual BlendingQuad2d* concepts::BlendingQuad2d::clone ( ) const [inline, virtual]

Returns a copy of the map.

Implements concepts::MappingQuad2d.

Definition at line 900 of file elementMaps.hh.

Here is the call graph for this function:

virtual MappingEdge2d* concepts::BlendingQuad2d::edge ( const uint edge ) const [inline, virtual]

Returns a copy of the edge mapping of one edge.

The edge mappings are directed counter-clockwise, i.e. for the lower edge from the left lower vertex to the right lower one.

Parameters:

edge	number of edge, i.e. 0 - lower, 1 - right, 2 - upper, 3 - left

Reimplemented from concepts::MappingQuad2d.

Definition at line 896 of file elementMaps.hh.

Here is the call graph for this function:

virtual MapReal2d concepts::BlendingQuad2d::hessian	(	uint	i,
		const Real	x,
		const Real	y
	)		const `[virtual]`

Returns the Hessian in a 2D linear map.

Parameters:

i	coordinate (0 = x, 1 = y)

Implements concepts::MappingQuad2d.

virtual std::ostream& concepts::BlendingQuad2d::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from concepts::MappingQuad2d.

virtual MapReal2d concepts::BlendingQuad2d::jacobian	(	const Real	x,
		const Real	y
	)		const `[virtual]`

Returns the Jacobian in a 2D linear map.

Implements concepts::MappingQuad2d.

virtual Real concepts::MappingQuad2d::jacobianDeterminant	(	const Real	x,
		const Real	y
	)		const `[virtual, inherited]`

Returns determinant of the Jacobian.

virtual MapReal2d concepts::MappingQuad2d::jacobianInverse	(	const Real	x,
		const Real	y
	)		const `[virtual, inherited]`

Returns the inverse of the Jacobian in a 2D linear map.

virtual Real concepts::BlendingQuad2d::lineElement	(	const Real	x,
		const uint	edge
	)		const `[virtual]`

Returns factor of differential element integrating over an edge.

The factor is either $\sqrt{J_{11}^2+J_{21}^2}$ or $\sqrt{J_{12}^2+J_{22}^2}$ .

The routine exists explicitly because it is more efficient to implement it than to use jacobian itself, i.e. for affine elements the line element is constant.

Parameters:

x	local variable on edge (in [0,1]), stands for $\xi_1$ or $\xi_2$
edge	number of edge, i.e. 0: $\xi_2 = 0$ , 1: $\xi_1 = 1$ , 2: $\xi_2 = 1$ , 3: $\xi_1 = 0$