Flow123d  release_3.0.0-506-g34af125
Public Types | Public Member Functions | Static Public Member Functions | Static Public Attributes | Private Member Functions | Static Private Member Functions | Static Private Attributes | List of all members
RefElement< dim > Class Template Reference

#include <ref_element.hh>

Collaboration diagram for RefElement< dim >:
Collaboration graph
[legend]

Public Types

typedef arma::vec::fixed< dim > LocalPoint
 
typedef arma::vec::fixed< dim+1 > BaryPoint
 
typedef arma::vec::fixed< dim > FaceBaryPoint
 
typedef std::vector< BaryPointBarycentricUnitVec
 
typedef const std::vector< LocalPoint > & CentersList
 

Public Member Functions

template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<>
const std::vector< std::vector< std::vector< unsigned int > > > nodes_of_subelements
 
template<>
const std::vector< std::vector< std::vector< unsigned int > > > nodes_of_subelements
 
template<>
const std::vector< std::vector< std::vector< unsigned int > > > nodes_of_subelements
 
template<>
vec::fixed< 1 > normal_vector (unsigned int sid)
 
template<>
vec::fixed< 2 > normal_vector (unsigned int sid)
 
template<>
vec::fixed< 3 > normal_vector (unsigned int sid)
 
template<>
auto clip (const BaryPoint &barycentric) -> BaryPoint
 
template<>
double side_measure (unsigned int sid)
 
template<>
double side_measure (unsigned int sid)
 
template<>
double side_measure (unsigned int sid)
 
template<>
unsigned int line_between_faces (unsigned int f1, unsigned int f2)
 
template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<>
const unsigned int side_permutations [][n_nodes_per_side]
 
template<unsigned int subdim>
auto interpolate (arma::vec::fixed< subdim+1 > coord, int sub_simplex_idx) -> BaryPoint
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<>
unsigned int count ()
 
template<unsigned int OutDim, unsigned int InDim>
const IdxVector< (InDim >OutDim?InDim+1:dim-InDim) > interact_ (unsigned int i)
 
template<template< unsigned int OutDim, unsigned int InDim > class TInteraction, unsigned int OutDim, unsigned int InDim>
const IdxVector< (InDim >OutDim?InDim+1:dim-InDim) > interact (TInteraction< OutDim, InDim > interaction)
 

Static Public Member Functions

static LocalPoint node_coords (unsigned int nid)
 
static LocalPoint normal_vector (unsigned int sid)
 
static BaryPoint clip (const BaryPoint &barycentric)
 
static unsigned int normal_orientation (unsigned int sid)
 
static double side_measure (unsigned int sid)
 
static unsigned int oposite_node (unsigned int sid)
 
static unsigned int line_between_faces (unsigned int f1, unsigned int f2)
 
static unsigned int permutation_index (unsigned int p[n_nodes_per_side])
 
static BaryPoint local_to_bary (const LocalPoint &lp)
 Converts from local to barycentric coordinates. More...
 
static LocalPoint bary_to_local (const BaryPoint &bp)
 Converts from barycentric to local coordinates. More...
 
static BarycentricUnitVec make_bary_unit_vec ()
 
static FaceBaryPoint barycentric_on_face (const BaryPoint &barycentric, unsigned int i_face)
 
static CentersList centers_of_subelements (unsigned int sub_dim)
 
static std::pair< unsigned int, unsigned int > zeros_positions (const BaryPoint &barycentric, double tolerance=std::numeric_limits< double >::epsilon()*2)
 
template<unsigned int subdim>
static unsigned int topology_idx (unsigned int zeros_positions)
 
template<unsigned int subdim>
static unsigned int count ()
 
template<unsigned int subdim>
static arma::mat::fixed< dim+1, subdim+1 > bary_coords (unsigned int sid)
 
template<unsigned int subdim>
static BaryPoint interpolate (arma::vec::fixed< subdim+1 > coord, int sub_simplex_idx)
 
static BaryPoint line_barycentric_interpolation (BaryPoint first_coords, BaryPoint second_coords, double first_theta, double second_theta, double theta)
 
template<template< unsigned int OutDim, unsigned int InDim > class TInteraction, unsigned int OutDim, unsigned int InDim>
static const IdxVector< (InDim >OutDim?InDim+1:dim-InDim) > interact (TInteraction< OutDim, InDim > interaction)
 

Static Public Attributes

static const unsigned int n_sides = dim + 1
 Number of sides. More...
 
static const unsigned int n_nodes = dim + 1
 Number of nodes. More...
 
static const unsigned int n_nodes_per_side = dim
 Number of nodes on one side. More...
 
static const unsigned int n_lines_per_node = dim
 Number of lines with one common node. More...
 
static const unsigned int n_nodes_per_line = 2
 Number of nodes in one line. More...
 
static const unsigned int n_sides_per_line = 2
 Number of sides with one common line. dim == 3. More...
 
static const unsigned int n_sides_per_node = dim
 Number of sides with one common line. More...
 
static const unsigned int n_lines_per_side = (unsigned int)((dim * (dim - 1)) / 2)
 Number of lines on boundary of one side. More...
 
static const unsigned int n_lines = (unsigned int)((dim * (dim + 1)) / 2)
 Number of lines, i.e. object of dimension dim-2 on the boundary of the reference element. More...
 
static const std::vector< std::vector< std::vector< unsigned int > > > nodes_of_subelements
 
static const unsigned int n_side_permutations = (dim+1)*(2*dim*dim-5*dim+6)/6
 
static const unsigned int side_permutations [n_side_permutations][n_nodes_per_side]
 

Private Member Functions

template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 1 > node_lines_ []
 
template<>
const IdxVector< 2 > node_lines_ []
 
template<>
const IdxVector< 3 > node_lines_ []
 
template<>
const IdxVector< 3 > side_nodes_ []
 
template<>
const IdxVector< 3 > node_sides_ []
 
template<>
const IdxVector< 2 > line_sides_ []
 
template<>
const IdxVector< 3 > side_lines_ []
 
template<>
const IdxVector< 2 > topology_zeros_ []
 
template<>
const IdxVector< 3 > topology_zeros_ []
 
template<>
const IdxVector< 6 > topology_zeros_ []
 
template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 2 > line_nodes_ []
 
template<>
const IdxVector< 1 > node_lines_ []
 
template<>
const IdxVector< 2 > node_lines_ []
 
template<>
const IdxVector< 3 > node_lines_ []
 
template<>
const IdxVector< 3 > side_nodes_ []
 
template<>
const IdxVector< 3 > node_sides_ []
 
template<>
const IdxVector< 2 > line_sides_ []
 
template<>
const IdxVector< 3 > side_lines_ []
 
template<>
const IdxVector< 2 > topology_zeros_ []
 
template<>
const IdxVector< 3 > topology_zeros_ []
 
template<>
const IdxVector< 6 > topology_zeros_ []
 
template<>
const IdxVector< 1 > interact_ (unsigned int i)
 This function is for "side_nodes" - for given side, give me nodes (0->0, 1->1). More...
 
template<>
const IdxVector< 2 > interact_ (unsigned int i)
 For line i {0}, give me indices of its nodes. More...
 
template<>
const IdxVector< 2 > interact_ (unsigned int i)
 For line i {0,1,2}, give me indices of its nodes. More...
 
template<>
const IdxVector< 2 > interact_ (unsigned int i)
 For line i {0,1,2,3,4,5}, give me indices of its nodes. More...
 
template<>
const IdxVector< 1 > interact_ (unsigned int i)
 For node i {0,1}, give me indices of lines. More...
 
template<>
const IdxVector< 2 > interact_ (unsigned int i)
 For node i {0,1,2}, give me indices of lines. More...
 
template<>
const IdxVector< 3 > interact_ (unsigned int i)
 For node i {0,1,2,3}, give me indices of lines. More...
 
template<>
const IdxVector< 3 > interact_ (unsigned int i)
 For side i {0,1,2}, give me indices of its nodes. More...
 
template<>
const IdxVector< 3 > interact_ (unsigned int i)
 For node i {0,1,2,3}, give me indices of sides. More...
 
template<>
const IdxVector< 2 > interact_ (unsigned int i)
 For line i {0,1,2,3}, give me indices of sides. More...
 
template<>
const IdxVector< 3 > interact_ (unsigned int i)
 For side i {0,1,2}, give me indices of its lines. More...
 

Static Private Member Functions

template<unsigned int OutDim, unsigned int InDim>
static const IdxVector< (InDim >OutDim?InDim+1:dim-InDim) > interact_ (unsigned int index)
 Internal part of the interact function. More...
 

Static Private Attributes

static const IdxVector< n_nodes_per_lineline_nodes_ [n_lines]
 For given line, returns its nodes indices. More...
 
static const IdxVector< n_lines_per_nodenode_lines_ [n_nodes]
 For given node, returns lines indices. More...
 
static const IdxVector< n_nodes_per_sideside_nodes_ [n_sides]
 For given side, returns nodes indices. For dim == 3. More...
 
static const IdxVector< n_sides_per_nodenode_sides_ [n_nodes]
 For given node, returns sides indices. For dim == 3. More...
 
static const IdxVector< n_sides_per_lineline_sides_ [n_lines]
 For given line, returns sides indices. For dim == 3. More...
 
static const IdxVector< n_lines_per_sideside_lines_ [n_sides]
 For given side, returns lines indices. For dim == 3. More...
 
static const IdxVector<(n_lines > n_nodes)?n_lines:n_nodes > topology_zeros_ [dim+1]
 

Detailed Description

template<unsigned int dim>
class RefElement< dim >

Definition at line 162 of file ref_element.hh.

Member Typedef Documentation

template<unsigned int dim>
typedef std::vector<BaryPoint> RefElement< dim >::BarycentricUnitVec

Definition at line 290 of file ref_element.hh.

template<unsigned int dim>
typedef arma::vec::fixed<dim+1> RefElement< dim >::BaryPoint

Barycentric coordinates.

e.g. coordinates (a,b,c) on triangle with vertices X, Y, Z represents a point: a*X+b*Y+c*Z

Definition at line 172 of file ref_element.hh.

template<unsigned int dim>
typedef const std::vector<LocalPoint>& RefElement< dim >::CentersList

Definition at line 305 of file ref_element.hh.

template<unsigned int dim>
typedef arma::vec::fixed<dim> RefElement< dim >::FaceBaryPoint

Definition at line 173 of file ref_element.hh.

template<unsigned int dim>
typedef arma::vec::fixed<dim> RefElement< dim >::LocalPoint

Definition at line 165 of file ref_element.hh.

Member Function Documentation

template<unsigned int dim>
template<unsigned int subdim>
arma::mat::fixed< dim+1, subdim+1 > RefElement< dim >::bary_coords ( unsigned int  sid)
inlinestatic
Parameters
sid- index of a sub-simplex in a simplex return an array of barycentric coordinates on <dim> simplex from <subdim> simplex for example: simplex<3> - ABCD and its subsubsimplex<1> AD (line index: 3) AD has barycoords for A (1,0), for D (0,1), but A in ABCD is (1,0,0,0) and D is (0,0,0,1) this method creates array ((1,0,0,0),(0,0,0,1))

Definition at line 450 of file ref_element.hh.

template<unsigned int dim>
auto RefElement< dim >::bary_to_local ( const BaryPoint bp)
static

Converts from barycentric to local coordinates.

Parameters
bppoint in barycentric coordinates
Returns
point in local coordinates

Definition at line 300 of file ref_element.cc.

Here is the caller graph for this function:

template<unsigned int dim>
auto RefElement< dim >::barycentric_on_face ( const BaryPoint barycentric,
unsigned int  i_face 
)
static

For given barycentric coordinates on the ref element returns barycentric coordinates on the ref. element of given face. Assumes that the input point is on the face. Barycentric order: (complanatory, local_coords )

Definition at line 378 of file ref_element.cc.

template<unsigned int dim>
auto RefElement< dim >::centers_of_subelements ( unsigned int  sub_dim)
static

Definition at line 465 of file ref_element.cc.

template<unsigned int dim>
auto RefElement< dim >::clip ( const BaryPoint barycentric)
static

If the given barycentric coordinate is in the ref. element, return unchanged. If the given barycentric coordinate is out of the ref. element, project it on the surface of the ref. element.

Definition at line 431 of file ref_element.cc.

Here is the caller graph for this function:

template<>
auto RefElement< 0 >::clip ( const BaryPoint barycentric) -> BaryPoint

Definition at line 411 of file ref_element.cc.

template<unsigned int dim>
template<unsigned int subdim>
static unsigned int RefElement< dim >::count ( )
static

Function returns number of subdim-simplices inside dim-simplex. The aim is covering all the n_**** members with a single function. TODO: think of generalization for n_****_per_**** members, like function interact: template<unsigned int subdimA, unsigned int subdimB> static unsigned int count();

template<>
unsigned int RefElement< 3 >::count< 0 > ( )
inline

Definition at line 504 of file ref_element.hh.

template<>
unsigned int RefElement< 3 >::count< 1 > ( )
inline

Definition at line 506 of file ref_element.hh.

template<>
unsigned int RefElement< 3 >::count< 2 > ( )
inline

Definition at line 508 of file ref_element.hh.

template<>
unsigned int RefElement< 3 >::count< 3 > ( )
inline

Definition at line 510 of file ref_element.hh.

template<>
unsigned int RefElement< 2 >::count< 0 > ( )
inline

Definition at line 512 of file ref_element.hh.

template<>
unsigned int RefElement< 2 >::count< 1 > ( )
inline

Definition at line 514 of file ref_element.hh.

template<>
unsigned int RefElement< 2 >::count< 2 > ( )
inline

Definition at line 516 of file ref_element.hh.

template<>
unsigned int RefElement< 2 >::count< 3 > ( )
inline

Definition at line 518 of file ref_element.hh.

template<>
unsigned int RefElement< 1 >::count< 0 > ( )
inline

Definition at line 520 of file ref_element.hh.

template<>
unsigned int RefElement< 1 >::count< 1 > ( )
inline

Definition at line 522 of file ref_element.hh.

template<>
unsigned int RefElement< 1 >::count< 2 > ( )
inline

Definition at line 524 of file ref_element.hh.

template<>
unsigned int RefElement< 1 >::count< 3 > ( )
inline

Definition at line 526 of file ref_element.hh.

template<unsigned int dim>
template<template< unsigned int OutDim, unsigned int InDim > class TInteraction, unsigned int OutDim, unsigned int InDim>
static const IdxVector< (InDim>OutDim ? InDim+1 : dim-InDim) > RefElement< dim >::interact ( TInteraction< OutDim, InDim >  interaction)
static

Usage: RefElement<3>::interact(Interaction<2,0>(1)) (means: In tetrahedron <3>, give indices of sides <2>, connected by node <0> with index 1) RefElement<3>::interact(Interaction<2,0>(1))[1] (as above, but give only the side with index 1)

Template usage: RefElement<dim>::interact(Interaction<OutDim, InDim>(i))[j] (means: on dim-dimensional reference element, go on InDim-dimensional subelement with index i, which connects OutDim-dimnesional subelements and select the one with index j)

This method serves as an interface to topology information of the reference element. It returns indices of OutDim-dimensional object of InDim-dimnesional object of given index in dim-dimnesional reference element.

Template Parameters
interaction- auxilliary object carying the index and the template arguments OutDim and InDim
OutDim- output dimension (give me node-0, line-1, side-2), <= dim
InDim- input dimension (for node-0, line-1, side-2), <= dim
Returns
vector of indices of OutDim-dimensional subelements represented by IdxVector object.

possible calls: dim OutDim InDim return 1,2,3 0 1 InDim+1 - give me indices of nodes of line of given index 3 0 2 InDim+1 - give me indices of nodes of a side (triangle) of given index 3 1 2 InDim+1 - give me indices of lines of side (triangle) of given index

1,2,3 1 0 dim-InDim - give me indices of lines with common node of given index 3 2 0 dim-InDim - give me indices of sides (triangles) with common node of given index 3 2 1 dim-InDim - give me indices of sides (triangles) with common line of given index

Here is the caller graph for this function:

template<unsigned int dim>
template<template< unsigned int OutDim, unsigned int InDim > class TInteraction, unsigned int OutDim, unsigned int InDim>
const IdxVector< (InDim>OutDim ? InDim+1 : dim-InDim) > RefElement< dim >::interact ( TInteraction< OutDim, InDim >  interaction)
inline

Definition at line 607 of file ref_element.hh.

template<unsigned int dim>
template<unsigned int OutDim, unsigned int InDim>
static const IdxVector< (InDim>OutDim ? InDim+1 : dim-InDim) > RefElement< dim >::interact_ ( unsigned int  index)
staticprivate

Internal part of the interact function.

template<>
const IdxVector< 1 > RefElement< 1 >::interact_< 0, 0 > ( unsigned int  i)
inlineprivate

This function is for "side_nodes" - for given side, give me nodes (0->0, 1->1).

Definition at line 542 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 1 >::interact_< 0, 1 > ( unsigned int  i)
inlineprivate

For line i {0}, give me indices of its nodes.

Definition at line 547 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 2 >::interact_< 0, 1 > ( unsigned int  i)
inlineprivate

For line i {0,1,2}, give me indices of its nodes.

Definition at line 552 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 3 >::interact_< 0, 1 > ( unsigned int  i)
inlineprivate

For line i {0,1,2,3,4,5}, give me indices of its nodes.

Definition at line 557 of file ref_element.hh.

template<>
const IdxVector< 1 > RefElement< 1 >::interact_< 1, 0 > ( unsigned int  i)
inlineprivate

For node i {0,1}, give me indices of lines.

Definition at line 562 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 2 >::interact_< 1, 0 > ( unsigned int  i)
inlineprivate

For node i {0,1,2}, give me indices of lines.

Definition at line 567 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::interact_< 1, 0 > ( unsigned int  i)
inlineprivate

For node i {0,1,2,3}, give me indices of lines.

Definition at line 572 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::interact_< 0, 2 > ( unsigned int  i)
inlineprivate

For side i {0,1,2}, give me indices of its nodes.

Definition at line 577 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::interact_< 2, 0 > ( unsigned int  i)
inlineprivate

For node i {0,1,2,3}, give me indices of sides.

Definition at line 582 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 3 >::interact_< 2, 1 > ( unsigned int  i)
inlineprivate

For line i {0,1,2,3}, give me indices of sides.

Definition at line 587 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::interact_< 1, 2 > ( unsigned int  i)
inlineprivate

For side i {0,1,2}, give me indices of its lines.

Definition at line 592 of file ref_element.hh.

template<unsigned int dim>
template<unsigned int OutDim, unsigned int InDim>
const IdxVector< (InDim>OutDim ? InDim+1 : dim-InDim) > RefElement< dim >::interact_ ( unsigned int  i)
inline

Definition at line 597 of file ref_element.hh.

template<unsigned int dim>
template<unsigned int subdim>
static BaryPoint RefElement< dim >::interpolate ( arma::vec::fixed< subdim+1 >  coord,
int  sub_simplex_idx 
)
static

Interpolate barycentric coords to a higher dimension of a simplex.

Parameters
coord- barycentric coords of a point on a sub-simplex
sub_simplex_idx- id of sub-simplex on a simplex
template<unsigned int dim>
template<unsigned int subdim>
auto RefElement< dim >::interpolate ( arma::vec::fixed< subdim+1 >  coord,
int  sub_simplex_idx 
) -> BaryPoint

Definition at line 481 of file ref_element.hh.

template<unsigned int dim>
arma::vec::fixed< dim+1 > RefElement< dim >::line_barycentric_interpolation ( BaryPoint  first_coords,
BaryPoint  second_coords,
double  first_theta,
double  second_theta,
double  theta 
)
static

Basic line interpolation.

Basic line interpolation

Definition at line 552 of file ref_element.cc.

Here is the caller graph for this function:

template<unsigned int dim>
static unsigned int RefElement< dim >::line_between_faces ( unsigned int  f1,
unsigned int  f2 
)
static

Return index of 1D line, shared by two faces f1 and f2 of the reference tetrahedron. Implemented only for dim == 3.

Here is the caller graph for this function:

template<>
unsigned int RefElement< 3 >::line_between_faces ( unsigned int  f1,
unsigned int  f2 
)

Definition at line 522 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 1 >::line_nodes_ ( )
private

Definition at line 50 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 2 >::line_nodes_ ( )
private

Definition at line 54 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 3 >::line_nodes_ ( )
private

Definition at line 60 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 1 >::line_nodes_ ( )
private

Definition at line 418 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 2 >::line_nodes_ ( )
private

Definition at line 419 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 3 >::line_nodes_ ( )
private

Definition at line 420 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 3 >::line_sides_ ( )
private

Definition at line 108 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 3 >::line_sides_ ( )
private

Definition at line 429 of file ref_element.hh.

template<unsigned int dim>
auto RefElement< dim >::local_to_bary ( const LocalPoint lp)
static

Converts from local to barycentric coordinates.

Parameters
lppoint in local coordinates (x,y)
Returns
point in barycentric coordinates (1-x-y, x, y)

Definition at line 286 of file ref_element.cc.

Here is the caller graph for this function:

template<unsigned int dim>
auto RefElement< dim >::make_bary_unit_vec ( )
static

Used in the clip method.

Definition at line 417 of file ref_element.cc.

template<unsigned int dim>
arma::vec::fixed< dim > RefElement< dim >::node_coords ( unsigned int  nid)
inlinestatic

Return coordinates of given node.

See also
the class documentation RefElement
Parameters
nidNode number. NOTE: Implementation is dependent on current node and side numbering.

Definition at line 465 of file ref_element.hh.

template<>
const IdxVector< 1 > RefElement< 1 >::node_lines_ ( )
private

Definition at line 70 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 2 >::node_lines_ ( )
private

Definition at line 75 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 3 >::node_lines_ ( )
private

Definition at line 82 of file ref_element.cc.

template<>
const IdxVector< 1 > RefElement< 1 >::node_lines_ ( )
private

Definition at line 422 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 2 >::node_lines_ ( )
private

Definition at line 423 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::node_lines_ ( )
private

Definition at line 424 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 3 >::node_sides_ ( )
private

Definition at line 100 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 3 >::node_sides_ ( )
private

Definition at line 427 of file ref_element.hh.

template<>
const std::vector< std::vector< std::vector< unsigned int > > > RefElement< 1 >::nodes_of_subelements ( )

Indexes of sides for each line - with right orientation Indexes of sides for each line - for Simplex<2>, with right orientation

Definition at line 266 of file ref_element.cc.

template<>
const std::vector< std::vector< std::vector< unsigned int > > > RefElement< 2 >::nodes_of_subelements ( )

Definition at line 271 of file ref_element.cc.

template<>
const std::vector< std::vector< std::vector< unsigned int > > > RefElement< 3 >::nodes_of_subelements ( )

Definition at line 277 of file ref_element.cc.

template<unsigned int dim>
unsigned int RefElement< dim >::normal_orientation ( unsigned int  sid)
static

Returns orientation of the normal of side sid. 0 -> OUT, 1 -> IN. NOTE: Implementation is dependent on current node and side numbering.

Definition at line 315 of file ref_element.cc.

Here is the caller graph for this function:

template<unsigned int dim>
static LocalPoint RefElement< dim >::normal_vector ( unsigned int  sid)
static

Compute normal vector to a given side.

Parameters
sidSide number.

Here is the caller graph for this function:

template<>
vec::fixed< 1 > RefElement< 1 >::normal_vector ( unsigned int  sid)

Definition at line 324 of file ref_element.cc.

template<>
vec::fixed< 2 > RefElement< 2 >::normal_vector ( unsigned int  sid)

Definition at line 332 of file ref_element.cc.

template<>
vec::fixed< 3 > RefElement< 3 >::normal_vector ( unsigned int  sid)

Definition at line 354 of file ref_element.cc.

template<unsigned int dim>
unsigned int RefElement< dim >::oposite_node ( unsigned int  sid)
inlinestatic

Returns index of the node that is oposite to side of given index sid. Note: It is dependent on current node and side numbering.

Parameters
sidSide number. NOTE: Implementation is dependent on current node and side numbering.

Definition at line 308 of file ref_element.cc.

Here is the caller graph for this function:

template<unsigned int dim>
unsigned int RefElement< dim >::permutation_index ( unsigned int  p[n_nodes_per_side])
static

For a given permutation p of nodes finds its index within side_permutations.

Parameters
pPermutation of nodes.

Definition at line 533 of file ref_element.cc.

Here is the caller graph for this function:

template<>
const IdxVector< 3 > RefElement< 3 >::side_lines_ ( )
private

Definition at line 118 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 3 >::side_lines_ ( )
private

Definition at line 431 of file ref_element.hh.

template<unsigned int dim>
static double RefElement< dim >::side_measure ( unsigned int  sid)
static
template<>
double RefElement< 1 >::side_measure ( unsigned int  sid)

Definition at line 493 of file ref_element.cc.

template<>
double RefElement< 2 >::side_measure ( unsigned int  sid)

Definition at line 502 of file ref_element.cc.

template<>
double RefElement< 3 >::side_measure ( unsigned int  sid)

Definition at line 511 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 3 >::side_nodes_ ( )
private

Definition at line 92 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 3 >::side_nodes_ ( )
private

Definition at line 426 of file ref_element.hh.

template<>
const unsigned int RefElement< 1 >::side_permutations ( )

Definition at line 127 of file ref_element.cc.

template<>
const unsigned int RefElement< 2 >::side_permutations ( )

Definition at line 129 of file ref_element.cc.

template<>
const unsigned int RefElement< 3 >::side_permutations ( )

Definition at line 131 of file ref_element.cc.

template<>
const unsigned int RefElement< 1 >::side_permutations ( )

Definition at line 433 of file ref_element.hh.

template<>
const unsigned int RefElement< 2 >::side_permutations ( )

Definition at line 434 of file ref_element.hh.

template<>
const unsigned int RefElement< 3 >::side_permutations ( )

Definition at line 435 of file ref_element.hh.

template<unsigned int dim>
template<unsigned int subdim>
unsigned int RefElement< dim >::topology_idx ( unsigned int  zeros_positions)
static

According to positions of zeros in barycentric coordinates, it gives the index of subdim-simplex in the reference element. Number of zeros must be equal to (3-subdim). e.g.: if 1 zeros, return index of side (subdim 2) if 2 zeros, return index of edge (subdim 1) if 3 zeros, return index of vertex (subdim 0)

Definition at line 531 of file ref_element.hh.

template<>
const IdxVector< 2 > RefElement< 1 >::topology_zeros_ ( )
private

Definition at line 141 of file ref_element.cc.

template<>
const IdxVector< 3 > RefElement< 2 >::topology_zeros_ ( )
private

Definition at line 149 of file ref_element.cc.

template<>
const IdxVector< 6 > RefElement< 3 >::topology_zeros_ ( )
private

Definition at line 161 of file ref_element.cc.

template<>
const IdxVector< 2 > RefElement< 1 >::topology_zeros_ ( )
private

Definition at line 437 of file ref_element.hh.

template<>
const IdxVector< 3 > RefElement< 2 >::topology_zeros_ ( )
private

Definition at line 438 of file ref_element.hh.

template<>
const IdxVector< 6 > RefElement< 3 >::topology_zeros_ ( )
private

Definition at line 439 of file ref_element.hh.

template<unsigned int dim>
std::pair< unsigned int, unsigned int > RefElement< dim >::zeros_positions ( const BaryPoint barycentric,
double  tolerance = std::numeric_limits<double>::epsilon()*2 
)
static

Return (1) number of zeros and (2) positions of zeros in barycentric coordinates. tolerance serves for testing zero values of barycentric coordinates.

Definition at line 393 of file ref_element.cc.

Here is the caller graph for this function:

Member Data Documentation

template<unsigned int dim>
const IdxVector<n_nodes_per_line> RefElement< dim >::line_nodes_[n_lines]
staticprivate

For given line, returns its nodes indices.

Definition at line 395 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<n_sides_per_line> RefElement< dim >::line_sides_[n_lines]
staticprivate

For given line, returns sides indices. For dim == 3.

Definition at line 399 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_lines = (unsigned int)((dim * (dim + 1)) / 2)
static

Number of lines, i.e. object of dimension dim-2 on the boundary of the reference element.

Definition at line 231 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_lines_per_node = dim
static

Number of lines with one common node.

Definition at line 222 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_lines_per_side = (unsigned int)((dim * (dim - 1)) / 2)
static

Number of lines on boundary of one side.

Definition at line 228 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_nodes = dim + 1
static

Number of nodes.

Definition at line 220 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_nodes_per_line = 2
static

Number of nodes in one line.

Definition at line 223 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_nodes_per_side = dim
static

Number of nodes on one side.

Definition at line 221 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_side_permutations = (dim+1)*(2*dim*dim-5*dim+6)/6
static

Number of permutations of nodes on sides.

dim value

1 1 2 2 3 6

Definition at line 265 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_sides = dim + 1
static

Number of sides.

Definition at line 219 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_sides_per_line = 2
static

Number of sides with one common line. dim == 3.

Definition at line 224 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::n_sides_per_node = dim
static

Number of sides with one common line.

Definition at line 225 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<n_lines_per_node> RefElement< dim >::node_lines_[n_nodes]
staticprivate

For given node, returns lines indices.

Definition at line 396 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<n_sides_per_node> RefElement< dim >::node_sides_[n_nodes]
staticprivate

For given node, returns sides indices. For dim == 3.

Definition at line 398 of file ref_element.hh.

template<unsigned int dim>
const std::vector< std::vector< std::vector<unsigned int> > > RefElement< dim >::nodes_of_subelements
static

Node numbers for each side. Indices of 1D lines of the 2D sides of an tetrahedron. Nonempty only for dim==3. Nodes of 1D lines of the tetrahedron. Indices of sides for each line. Nonempty only for dim==3 and dim==2.

Definition at line 255 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<n_lines_per_side> RefElement< dim >::side_lines_[n_sides]
staticprivate

For given side, returns lines indices. For dim == 3.

Definition at line 400 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<n_nodes_per_side> RefElement< dim >::side_nodes_[n_sides]
staticprivate

For given side, returns nodes indices. For dim == 3.

Definition at line 397 of file ref_element.hh.

template<unsigned int dim>
const unsigned int RefElement< dim >::side_permutations[n_side_permutations][n_nodes_per_side]
static

Permutations of nodes on sides.

Definition at line 270 of file ref_element.hh.

template<unsigned int dim>
const IdxVector<(n_lines > n_nodes) ? n_lines : n_nodes> RefElement< dim >::topology_zeros_[dim+1]
staticprivate

Consider an n-face (node, edge, face, bulk) with dimension subdim and index within subdimension idx. Barycentric coordinates of all points on the n-face have unique pattern of zero coordinates.

topology_zeros_[subdim][idx] is a bitfield with '1' where the pattern have zeros.

Definition at line 410 of file ref_element.hh.


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