Loading [MathJax]/extensions/tex2jax.js
ug4
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
neumann_boundary_fe.h
Go to the documentation of this file.
1/*
2 * Copyright (c) 2013-2015: G-CSC, Goethe University Frankfurt
3 * Author: Andreas Vogel
4 *
5 * This file is part of UG4.
6 *
7 * UG4 is free software: you can redistribute it and/or modify it under the
8 * terms of the GNU Lesser General Public License version 3 (as published by the
9 * Free Software Foundation) with the following additional attribution
10 * requirements (according to LGPL/GPL v3 §7):
11 *
12 * (1) The following notice must be displayed in the Appropriate Legal Notices
13 * of covered and combined works: "Based on UG4 (www.ug4.org/license)".
14 *
15 * (2) The following notice must be displayed at a prominent place in the
16 * terminal output of covered works: "Based on UG4 (www.ug4.org/license)".
17 *
18 * (3) The following bibliography is recommended for citation and must be
19 * preserved in all covered files:
20 * "Reiter, S., Vogel, A., Heppner, I., Rupp, M., and Wittum, G. A massively
21 * parallel geometric multigrid solver on hierarchically distributed grids.
22 * Computing and visualization in science 16, 4 (2013), 151-164"
23 * "Vogel, A., Reiter, S., Rupp, M., Nägel, A., and Wittum, G. UG4 -- a novel
24 * flexible software system for simulating pde based models on high performance
25 * computers. Computing and visualization in science 16, 4 (2013), 165-179"
26 *
27 * This program is distributed in the hope that it will be useful,
28 * but WITHOUT ANY WARRANTY; without even the implied warranty of
29 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
30 * GNU Lesser General Public License for more details.
31 */
32
33#ifndef __H__UG__LIB_DISC__SPATIAL_DISC__ELEM_DISC__NEUMANN_BOUNDARY___NEUMANN_BOUNDARY_FE__
34#define __H__UG__LIB_DISC__SPATIAL_DISC__ELEM_DISC__NEUMANN_BOUNDARY___NEUMANN_BOUNDARY_FE__
35
36// other ug4 modules
37#include "common/common.h"
38
39// library intern headers
40#include "../neumann_boundary_base.h"
41
42namespace ug{
43
44template<typename TDomain>
46 : public NeumannBoundaryBase<TDomain>
47{
48 private:
51
54
55 public:
57 static const int dim = base_type::dim;
58
59 public:
61 NeumannBoundaryFE(const char* function);
62
65 void add(SmartPtr<CplUserData<number, dim> > data, const char* BndSubsets, const char* InnerSubsets);
66 void add(SmartPtr<CplUserData<number, dim, bool> > user, const char* BndSubsets, const char* InnerSubsets);
67 void add(SmartPtr<CplUserData<MathVector<dim>, dim> > user, const char* BndSubsets, const char* InnerSubsets);
69
70 protected:
71 using typename base_type::Data;
72
75 {
77 std::string BndSubsets, std::string InnerSubsets,
78 NeumannBoundaryFE* this_)
79 : base_type::Data(BndSubsets, InnerSubsets), This(this_)
80 {
81 import.set_data(data);
82 }
83
84 template<typename TElem, typename TFEGeom>
85 void extract_bip(const TFEGeom& geo);
86
87 template <typename TElem, typename TFEGeom>
88 void lin_def(const LocalVector& u,
89 std::vector<std::vector<number> > vvvLinDef[],
90 const size_t nip);
91
93 std::vector<MathVector<dim> > vLocIP;
94 std::vector<MathVector<dim> > vGloIP;
96 };
97 friend struct NumberData;
98
101 {
103 std::string BndSubsets, std::string InnerSubsets)
104 : base_type::Data(BndSubsets, InnerSubsets), functor(functor_) {}
105
107 };
108
111 {
113 std::string BndSubsets, std::string InnerSubsets)
114 : base_type::Data(BndSubsets, InnerSubsets), functor(functor_) {}
115
117 };
118
119 std::vector<NumberData> m_vNumberData;
120 std::vector<BNDNumberData> m_vBNDNumberData;
121 std::vector<VectorData> m_vVectorData;
122
124
125 public:
127 virtual void prepare_setting(const std::vector<LFEID>& vLfeID, bool bNonRegularGrid);
128
129 protected:
133
136
138 int m_si;
139
140 protected:
143 template<typename TElem, typename TFEGeom>
144 void prep_elem_loop(const ReferenceObjectID roid, const int si);
145 template<typename TElem, typename TFEGeom>
146 void prep_elem(const LocalVector& u, GridObject* elem, const ReferenceObjectID roid, const MathVector<dim> vCornerCoords[]);
147 template<typename TElem, typename TFEGeom>
148 void finish_elem_loop();
149 template<typename TElem, typename TFEGeom>
150 void add_rhs_elem(LocalVector& d, GridObject* elem, const MathVector<dim> vCornerCoords[]);
152
153 static const int _C_ = 0;
154
155 protected:
156 void register_all_funcs(int order);
157 template<typename TElem, typename TFEGeom> void register_func();
158};
159
160} // end namespac ug
161
162#endif /*__H__UG__LIB_DISC__SPATIAL_DISC__ELEM_DISC__NEUMANN_BOUNDARY___NEUMANN_BOUNDARY_FE1__*/
Definition smart_pointer.h:108
Type based UserData.
Definition user_data.h:501
Data import.
Definition data_import.h:180
The base class for all geometric objects, such as vertices, edges, faces, volumes,...
Definition grid_base_objects.h:157
Identifier for Local Finite Elements.
Definition local_finite_element_id.h:98
Definition local_algebra.h:198
a mathematical Vector with N entries.
Definition math_vector.h:97
Definition neumann_boundary_base.h:47
static const int dim
World dimension.
Definition neumann_boundary_base.h:57
Definition neumann_boundary_fe.h:47
virtual void prepare_setting(const std::vector< LFEID > &vLfeID, bool bNonRegularGrid)
type of trial space for each function used
Definition neumann_boundary_fe.cpp:53
void prep_elem(const LocalVector &u, GridObject *elem, const ReferenceObjectID roid, const MathVector< dim > vCornerCoords[])
Definition neumann_boundary_fe.cpp:174
void register_all_funcs(int order)
std::vector< NumberData > m_vNumberData
Definition neumann_boundary_fe.h:119
int m_si
current inner subset
Definition neumann_boundary_fe.h:138
int m_quadOrder
Definition neumann_boundary_fe.h:132
static const int _C_
Definition neumann_boundary_fe.h:153
int m_order
current order of disc scheme
Definition neumann_boundary_fe.h:131
LFEID m_lfeID
current shape function set
Definition neumann_boundary_fe.h:135
void add_rhs_elem(LocalVector &d, GridObject *elem, const MathVector< dim > vCornerCoords[])
Definition neumann_boundary_fe.cpp:194
void prep_elem_loop(const ReferenceObjectID roid, const int si)
Definition neumann_boundary_fe.cpp:114
NeumannBoundaryFE< TDomain > this_type
Base class type.
Definition neumann_boundary_fe.h:53
void add(SmartPtr< CplUserData< number, dim > > data, const char *BndSubsets, const char *InnerSubsets)
Definition neumann_boundary_fe.cpp:73
std::vector< VectorData > m_vVectorData
Definition neumann_boundary_fe.h:121
NeumannBoundaryBase< TDomain > base_type
Base class type.
Definition neumann_boundary_fe.h:50
static const int dim
World dimension.
Definition neumann_boundary_fe.h:57
void update_subset_groups()
Definition neumann_boundary_fe.cpp:96
void finish_elem_loop()
Definition neumann_boundary_fe.cpp:260
std::vector< BNDNumberData > m_vBNDNumberData
Definition neumann_boundary_fe.h:120
void register_func()
Definition neumann_boundary_fe.cpp:381
the ug namespace
ReferenceObjectID
these ids are used to identify the shape of a geometric object.
Definition grid_base_objects.h:74
base class for user data
Definition neumann_boundary_base.h:87
Conditional scalar user data.
Definition neumann_boundary_fe.h:101
SmartPtr< CplUserData< number, dim, bool > > functor
Definition neumann_boundary_fe.h:106
BNDNumberData(SmartPtr< CplUserData< number, dim, bool > > functor_, std::string BndSubsets, std::string InnerSubsets)
Definition neumann_boundary_fe.h:102
Unconditional scalar user data.
Definition neumann_boundary_fe.h:75
std::vector< MathVector< dim > > vLocIP
Definition neumann_boundary_fe.h:93
void extract_bip(const TFEGeom &geo)
Definition neumann_boundary_fe.cpp:324
std::vector< MathVector< dim > > vGloIP
Definition neumann_boundary_fe.h:94
NeumannBoundaryFE * This
Definition neumann_boundary_fe.h:95
NumberData(SmartPtr< CplUserData< number, dim > > data, std::string BndSubsets, std::string InnerSubsets, NeumannBoundaryFE *this_)
Definition neumann_boundary_fe.h:76
void lin_def(const LocalVector &u, std::vector< std::vector< number > > vvvLinDef[], const size_t nip)
Definition neumann_boundary_fe.cpp:300
Unconditional vector user data.
Definition neumann_boundary_fe.h:111
SmartPtr< CplUserData< MathVector< dim >, dim > > functor
Definition neumann_boundary_fe.h:116
VectorData(SmartPtr< CplUserData< MathVector< dim >, dim > > functor_, std::string BndSubsets, std::string InnerSubsets)
Definition neumann_boundary_fe.h:112