periodic_boundary_base.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2018 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 // Local Includes
19 #include "libmesh/libmesh_config.h"
20 
21 #ifdef LIBMESH_ENABLE_PERIODIC
22 
23 #include "libmesh/boundary_info.h" // BoundaryInfo::invalid_id
25 
26 namespace libMesh
27 {
28 
30  myboundary(BoundaryInfo::invalid_id),
31  pairedboundary(BoundaryInfo::invalid_id)
32 {
33 }
34 
35 
36 
38  myboundary(o.myboundary),
39  pairedboundary(o.pairedboundary),
40  variables(o.variables)
41 {
42  // Make a deep copy of _transformation_matrix, if it's not null
44  {
45  this->_transformation_matrix = libmesh_make_unique<DenseMatrix<Real>>();
47  }
48 }
49 
50 
51 
52 void PeriodicBoundaryBase::set_variable(unsigned int var)
53 {
54  variables.insert(var);
55 }
56 
57 
58 
60 {
61  variables.insert(pb.variables.begin(), pb.variables.end());
62 }
63 
64 
65 
66 bool PeriodicBoundaryBase::is_my_variable(unsigned int var_num) const
67 {
68  bool a = variables.empty() || (!variables.empty() && variables.find(var_num) != variables.end());
69  return a;
70 }
71 
72 
73 
75 {
76  return bool(_transformation_matrix);
77 }
78 
79 
80 
82 {
84  {
85  libmesh_error_msg("Transformation matrix is not defined");
86  }
87 
88  return *_transformation_matrix;
89 }
90 
91 
92 
94 {
95  // Make a deep copy of matrix
96  this->_transformation_matrix = libmesh_make_unique<DenseMatrix<Real>>();
97  *(this->_transformation_matrix) = matrix;
98 
99  // if _transformation_matrix is defined then it must be the same sie as variables.
100  libmesh_assert_equal_to(_transformation_matrix->m(), variables.size());
101  libmesh_assert_equal_to(_transformation_matrix->n(), variables.size());
102 }
103 
104 
105 
106 const std::set<unsigned int> & PeriodicBoundaryBase::get_variables() const
107 {
108  return variables;
109 }
110 
111 } // namespace libMesh
112 
113 #endif // LIBMESH_ENABLE_PERIODIC
const std::set< unsigned int > & get_variables() const
void merge(const PeriodicBoundaryBase &pb)
void set_variable(unsigned int var)
std::set< unsigned int > variables
void set_transformation_matrix(const DenseMatrix< Real > &matrix)
Used by the Mesh to keep track of boundary nodes and elements.
Definition: boundary_info.h:57
const DenseMatrix< Real > & get_transformation_matrix() const
bool is_my_variable(unsigned int var_num) const
Base class for all PeriodicBoundary implementations.
A matrix object used for finite element assembly and numerics.
Definition: dense_matrix.h:54
std::unique_ptr< DenseMatrix< Real > > _transformation_matrix