Flow123d  master-f44eb46
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  inverted(false),
39  pid_(0),
40  n_neighs_vb_(0),
41  dim_(0)
42 {
43 }
44 
45 
46 Element::Element(unsigned int dim, RegionIdx reg)
47 {
48  init(dim, reg);
49 }
50 
51 
52 
53 void Element::init(unsigned int dim, RegionIdx reg) {
54  inverted = false;
55  pid_=0;
56  n_neighs_vb_=0;
57  neigh_vb=NULL;
58  dim_=dim;
59  region_idx_=reg;
60 
61  edge_idx_.resize( n_sides() );
62  boundary_idx_ = NULL;
63 
64  for (unsigned int si=0; si<this->n_sides(); si++) {
65  edge_idx_[ si ]=undef_idx;
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:
~Element()
Definition: elements.cc:70
unsigned int * boundary_idx_
Definition: elements.h:79
unsigned int dim_
Definition: elements.h:105
int pid_
Id # of mesh partition.
Definition: elements.h:95
Neighbour ** neigh_vb
Definition: elements.h:83
bool inverted
Inverted permutation of element nodes, negative Jacobian.
Definition: elements.h:88
unsigned int n_neighs_vb_
Definition: elements.h:97
void init(unsigned int dim, RegionIdx reg)
Definition: elements.cc:53
RegionIdx region_idx_
Definition: elements.h:104
Element()
Definition: elements.cc:35
unsigned int n_sides() const
Definition: elements.h:131
unsigned int dim() const
Definition: elements.h:120
std::vector< unsigned int > edge_idx_
Edges on sides.
Definition: elements.h:96
Support classes for parallel programing.
const unsigned int undef_idx
Definition: index_types.hh:32
Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc.