ref_element.hh File Reference

Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc. More...

#include <vector>
#include <array>
#include <armadillo>
#include "system/armor.hh"
#include "system/asserts.hh"

Go to the source code of this file.

Classes
struct	Interaction< OutDim, InDim >
class	RefElement< dim >

Typedefs
template<std::size_t Size>
using	IdxVector = std::array< unsigned int, Size >

Detailed Description

Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc.

 * Copyright (C) 2015 Technical University of Liberec. All rights reserved.

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License version 3 as published by the Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

Author

Jan Stebel

Todo:

TODO: reconsider following whether it is actual...

• design interface in such a way, that we can change numbering

• design numbering and orientations on ref element that is consistent (orientation and numbering od 2d el. match sides of 3d), and possibly allows permute vertices of elements so that sides sharing an edge match numbering and orientation (we dos'nt need permute faces)

  Proposal(prefers combinatoric order) : 1D - orientation V0 -> V1

  2D - edges: E0: V0 -> V1, E1: V0 -> V2 E2: V1 -> V2 This maximize number of edge orientations achievable by edge permutations edge numbering edge orientation( in original numbering) 0 1 2 + + + 0 2 1 - + + 1 0 2 + + - 1 2 0 - - + 2 0 1 + - - 2 1 0 - - -

        vertices   edges       normal (out = +)
  

  3D - sides: S0: 0 1 2 E0 E1 E3 - S1: 0 1 3 E0 E2 E4 + S2: 0 2 3 E1 E2 E5 - S3: 1 2 3 E3 E4 E5 -

  edges: E0: V0 -> V1 x direction E1: V0 -> V2 y direction E2: V0 -> V3 z direction E3: V1 -> V2 E4: V1 -> V3 E5: V2 -> V3

• functions from DEAL.ii: bool is_inside_unit_cell( point ) line_to_cell_vertices(line, vertex) vertex index on line to index on whole element face_to_cell_vertices(face, vertex, Orientation), Orientation should be some class describing permutation of shared face to element's face/side face_to_cel_lines standard_to_real_face_vertex(vertex, Orientation), maps vertex to permuted face real_to_standard_face_vertex - inverse ... same for line; we should need also something like standard_to_real_line_vertex; Deal dosn;t support Orientation changes for 2D element faces Point unit_cell_vertex(vertex) - coordinates project_to_unit_cell distance_to_unit_cell d_linear_shape_function
• can not change numbering of element sides due to DarcyFlow, which use hardwired side numbering in construction of basis functions
• any change of side numbering requires also change in flow/old_bcd.cc

Definition in file ref_element.hh.

Typedef Documentation

IdxVector

template<std::size_t Size>

using IdxVector = std::array<unsigned int, Size>

Auxilliary class representing vector of indices (unsigned int).

Template Parameters

Size	is the fixed size of the vector.

Definition at line 147 of file ref_element.hh.