ug4
assemble_interface.h
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2010-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__ASSEMBLE__
34 #define __H__UG__LIB_DISC__ASSEMBLE__
35 
39 
40 namespace ug{
41 
42 //predeclaration
43 template <typename TAlgebra>
44 class IConstraint;
45 
46 
58 
61 
108 template <typename TAlgebra>
110 {
111  public:
113  typedef TAlgebra algebra_type;
114 
116  typedef typename TAlgebra::matrix_type matrix_type;
117 
119  typedef typename TAlgebra::vector_type vector_type;
120 
121  public:
123 
130  virtual void assemble_jacobian(matrix_type& J, const vector_type& u, const GridLevel& gl) = 0;
132  {assemble_jacobian(J,u,GridLevel());}
133 
135 
142  virtual void assemble_defect(vector_type& d, const vector_type& u, const GridLevel& gl) = 0;
144  {assemble_defect(d,u, GridLevel());}
145 
147 
154  virtual void assemble_linear(matrix_type& A, vector_type& b, const GridLevel& gl) = 0;
156  {assemble_linear(A,b, GridLevel());}
157 
159  virtual void assemble_rhs(vector_type& rhs, const vector_type& u, const GridLevel& gl) = 0;
160  virtual void assemble_rhs(vector_type& rhs, const vector_type& u)
161  {assemble_rhs(rhs, u, GridLevel());}
162 
164 
170  virtual void assemble_rhs(vector_type& b, const GridLevel& gl) = 0;
172  {assemble_rhs(b, GridLevel());}
173 
175 
182  virtual void adjust_solution(vector_type& u, const GridLevel& gl) = 0;
184  {adjust_solution(u, GridLevel());}
185 
187  virtual void assemble_mass_matrix(matrix_type& M, const vector_type& u, const GridLevel& gl)
188  {UG_THROW("IAssemble: assemble_mass_matrix not implemented.");}
191 
193  virtual void assemble_stiffness_matrix(matrix_type& A, const vector_type& u, const GridLevel& gl)
194  {UG_THROW("IAssemble: assemble_stiffness_matrix not implemented.");}
197 
202 
204  virtual size_t num_constraints() const = 0;
205 
208 
210  virtual ~IAssemble(){};
211 
212 };
213 
215 
216 }; // name space ug
217 
218 #endif /* __H__UG__LIB_DISC__ASSEMBLE__ */
Definition: smart_pointer.h:296
Definition: smart_pointer.h:108
Definition: grid_level.h:42
Interface providing Jacobian and Defect of a discretization.
Definition: assemble_interface.h:110
TAlgebra algebra_type
Algebra type.
Definition: assemble_interface.h:113
virtual void adjust_solution(vector_type &u, const GridLevel &gl)=0
sets dirichlet values in solution vector
virtual SmartPtr< IConstraint< TAlgebra > > constraint(size_t i)=0
returns the i'th constraint
void assemble_stiffness_matrix(matrix_type &A, const vector_type &u)
Definition: assemble_interface.h:195
virtual ~IAssemble()
Virtual Destructor.
Definition: assemble_interface.h:210
TAlgebra::vector_type vector_type
Type of algebra vector.
Definition: assemble_interface.h:119
virtual void assemble_stiffness_matrix(matrix_type &A, const vector_type &u, const GridLevel &gl)
assembles stiffness matrix
Definition: assemble_interface.h:193
virtual void assemble_defect(vector_type &d, const vector_type &u, const GridLevel &gl)=0
assembles Defect
void assemble_linear(matrix_type &A, vector_type &b)
Definition: assemble_interface.h:155
virtual void assemble_linear(matrix_type &A, vector_type &b, const GridLevel &gl)=0
Assembles Matrix and Right-Hand-Side for a linear problem.
void assemble_defect(vector_type &d, const vector_type &u)
Definition: assemble_interface.h:143
void assemble_rhs(vector_type &b)
Definition: assemble_interface.h:171
void assemble_mass_matrix(matrix_type &M, const vector_type &u)
Definition: assemble_interface.h:189
void adjust_solution(vector_type &u)
Definition: assemble_interface.h:183
virtual SmartPtr< AssemblingTuner< TAlgebra > > ass_tuner()=0
virtual ConstSmartPtr< AssemblingTuner< TAlgebra > > ass_tuner() const =0
virtual void assemble_mass_matrix(matrix_type &M, const vector_type &u, const GridLevel &gl)
assembles mass matrix
Definition: assemble_interface.h:187
virtual void assemble_rhs(vector_type &rhs, const vector_type &u)
Definition: assemble_interface.h:160
virtual void assemble_rhs(vector_type &rhs, const vector_type &u, const GridLevel &gl)=0
assembles rhs
virtual void assemble_rhs(vector_type &b, const GridLevel &gl)=0
Assembles Right-Hand-Side for a linear problem.
TAlgebra::matrix_type matrix_type
Type of algebra matrix.
Definition: assemble_interface.h:116
void assemble_jacobian(matrix_type &J, const vector_type &u)
Definition: assemble_interface.h:131
virtual size_t num_constraints() const =0
returns the number of constraints
virtual void assemble_jacobian(matrix_type &J, const vector_type &u, const GridLevel &gl)=0
assembles Jacobian (or Approximation of Jacobian)
#define UG_THROW(msg)
Definition: error.h:57
the ug namespace