libMesh
transient_rb_assembly_expansion.C
Go to the documentation of this file.
1 // rbOOmit: An implementation of the Certified Reduced Basis method.
2 // Copyright (C) 2009, 2010 David J. Knezevic
3 
4 // This file is part of rbOOmit.
5 
6 // rbOOmit is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 
11 // rbOOmit is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
15 
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
19 
20 // rbOOmit includes
21 #include "libmesh/transient_rb_assembly_expansion.h"
22 #include "libmesh/elem_assembly.h"
23 
24 namespace libMesh
25 {
26 
27 // ------------------------------------------------------------
28 // TransientRBAssemblyExpansion implementation
29 
31 {
32 }
33 
35  FEMContext & context)
36 {
37  if (q >= get_n_M_terms())
38  libmesh_error_msg("Error: We must have q < get_n_M_terms in perform_M_interior_assembly.");
39 
41 
42  return _M_assembly_vector[q]->interior_assembly( context );
43 }
44 
46  FEMContext & context)
47 {
48  if (q >= get_n_M_terms())
49  libmesh_error_msg("Error: We must have q < get_n_M_terms in perform_M_boundary_assembly.");
50 
52 
53  return _M_assembly_vector[q]->boundary_assembly( context );
54 }
55 
57 {
58  return cast_int<unsigned int>(_M_assembly_vector.size());
59 }
60 
62 {
63  _M_assembly_vector.push_back(M_q_assembly);
64 }
65 
67 {
68  if (q >= get_n_M_terms())
69  libmesh_error_msg("Error: We must have q < get_n_M_terms in get_M_assembly.");
70 
71  return *_M_assembly_vector[q];
72 }
73 
74 }
void attach_M_assembly(ElemAssembly *A_q_assembly)
Attach ElemAssembly object for the time-derivative (both interior and boundary assembly).
The libMesh namespace provides an interface to certain functionality in the library.
libmesh_assert(j)
void perform_M_boundary_assembly(unsigned int q, FEMContext &context)
Perform the specified M boundary assembly.
ElemAssembly & get_M_assembly(unsigned int q)
Return a reference to the specified M_assembly object.
std::vector< ElemAssembly * > _M_assembly_vector
Vectors storing the function pointers to the assembly routines for the time-derivative operators...
This class provides all data required for a physics package (e.g.
Definition: fem_context.h:61
ElemAssembly provides a per-element (interior and boundary) assembly functionality.
Definition: elem_assembly.h:38
unsigned int get_n_M_terms() const
Get Q_m, the number of terms in the affine expansion for the bilinear form.
void perform_M_interior_assembly(unsigned int q, FEMContext &context)
Perform the specified M interior assembly.