Flow123d  JS_before_hm-1002-gafa1d04
elements.cc
Go to the documentation of this file.
1 /*!
2  *
3  * Copyright (C) 2015 Technical University of Liberec. All rights reserved.
4  *
5  * This program is free software; you can redistribute it and/or modify it under
6  * the terms of the GNU General Public License version 3 as published by the
7  * Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
8  *
9  * This program is distributed in the hope that it will be useful, but WITHOUT
10  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11  * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
12  *
13  *
14  * @file elements.cc
15  * @ingroup mesh
16  * @brief Various element oriented stuff, should be restricted to purely geometric functions
17  */
18 
19 #include <vector>
20 #include <string>
21 
22 #include "system/system.hh"
23 #include "mesh/accessors.hh"
24 #include "elements.h"
25 #include "mesh/mesh.h"
26 #include "mesh/ref_element.hh"
27 
28 // following deps. should be removed
29 //#include "materials.hh"
30 #include "mesh/accessors.hh"
31 #include "la/distribution.hh"
32 
33 
34 
36 : boundary_idx_(NULL),
37  neigh_vb(NULL),
38  pid_(0),
39  n_neighs_vb_(0),
40  dim_(0)
41 {
42 }
43 
44 
45 Element::Element(unsigned int dim, RegionIdx reg)
46 {
47  init(dim, reg);
48 }
49 
50 
51 
52 void Element::init(unsigned int dim, RegionIdx reg) {
53  pid_=0;
54  n_neighs_vb_=0;
55  neigh_vb=NULL;
56  dim_=dim;
57  region_idx_=reg;
58 
59  edge_idx_.resize( n_sides() );
60  permutation_idx_.resize( n_sides() );
61  boundary_idx_ = NULL;
62 
63  for (unsigned int si=0; si<this->n_sides(); si++) {
66  }
67 }
68 
69 
71  // Can not make deallocation here since then resize of
72  // vectors of elements deallocates what should be keeped.
73 }
74 
75 
76 /**
77  * Count element sides of the space dimension @p side_dim.
78  */
79 
80 /* If we use this method, it will be moved to mesh accessor class.
81 unsigned int Element::n_sides_by_dim(unsigned int side_dim)
82 {
83  if (side_dim == dim()) return 1;
84 
85  unsigned int n = 0;
86  for (unsigned int i=0; i<n_sides(); i++)
87  if (side(i)->dim() == side_dim) n++;
88  return n;
89 }*/
90 
91 
92 
93 //-----------------------------------------------------------------------------
94 // vim: set cindent:
unsigned int * boundary_idx_
Definition: elements.h:80
static const unsigned int undef_idx
Definition: mesh.h:104
Element()
Definition: elements.cc:35
void init(unsigned int dim, RegionIdx reg)
Definition: elements.cc:52
unsigned int dim() const
Definition: elements.h:121
~Element()
Definition: elements.cc:70
Neighbour ** neigh_vb
Definition: elements.h:84
unsigned int n_sides() const
Definition: elements.h:132
std::vector< unsigned int > edge_idx_
Edges on sides.
Definition: elements.h:90
std::vector< unsigned int > permutation_idx_
Definition: elements.h:104
unsigned int dim_
Definition: elements.h:108
Support classes for parallel programing.
RegionIdx region_idx_
Definition: elements.h:107
int pid_
Id # of mesh partition.
Definition: elements.h:89
Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc...
unsigned int n_neighs_vb_
of neighbours, V-B type (comp.)
Definition: elements.h:91