libMesh
Public Member Functions | Static Public Member Functions | Public Attributes | Protected Types | Protected Member Functions | Static Protected Attributes | List of all members
AssemblyA1 Struct Reference

#include <assembly.h>

Inheritance diagram for AssemblyA1:
[legend]

Public Member Functions

 AssemblyA1 (ElasticityRBConstruction &rb_sys_in)
 
virtual void interior_assembly (FEMContext &c)
 Perform the element interior assembly. More...
 
virtual void boundary_assembly (FEMContext &c)
 Perform the element boundary assembly. More...
 
virtual void get_nodal_values (std::vector< dof_id_type > &, DenseMatrix< Number > &, DenseVector< Number > &, const System &, const Node &)
 Get values to add to the matrix or rhs vector based on node. More...
 
virtual void get_nodal_rhs_values (std::map< numeric_index_type, Number > &values, const System &, const Node &)
 Get values to add to the RHS vector based on node. More...
 

Static Public Member Functions

static std::string get_info ()
 Gets a string containing the reference information. More...
 
static void print_info (std::ostream &out=libMesh::out)
 Prints the reference information, by default to libMesh::out. More...
 
static unsigned int n_objects ()
 Prints the number of outstanding (created, but not yet destroyed) objects. More...
 
static void enable_print_counter_info ()
 Methods to enable/disable the reference counter output from print_info() More...
 
static void disable_print_counter_info ()
 
static std::string get_info ()
 Gets a string containing the reference information. More...
 
static void print_info (std::ostream &out=libMesh::out)
 Prints the reference information, by default to libMesh::out. More...
 
static unsigned int n_objects ()
 Prints the number of outstanding (created, but not yet destroyed) objects. More...
 
static void enable_print_counter_info ()
 Methods to enable/disable the reference counter output from print_info() More...
 
static void disable_print_counter_info ()
 

Public Attributes

RBConstructionrb_con
 
bool is_nodal_rhs_values_overriden
 Temporary flag to help us figure out if we should call the deprecated get_nodal_rhs_values method or not. More...
 
ElasticityRBConstructionrb_sys
 The ElasticityRBConstruction object that will use this assembly. More...
 
bool is_nodal_rhs_values_overriden
 Temporary flag to help us figure out if we should call the deprecated get_nodal_rhs_values method or not. More...
 

Protected Types

typedef std::map< std::string, std::pair< unsigned int, unsigned int > > Counts
 Data structure to log the information. More...
 
typedef std::map< std::string, std::pair< unsigned int, unsigned int > > Counts
 Data structure to log the information. More...
 

Protected Member Functions

void increment_constructor_count (const std::string &name)
 Increments the construction counter. More...
 
void increment_destructor_count (const std::string &name)
 Increments the destruction counter. More...
 
void increment_constructor_count (const std::string &name)
 Increments the construction counter. More...
 
void increment_destructor_count (const std::string &name)
 Increments the destruction counter. More...
 

Static Protected Attributes

static Counts _counts
 Actually holds the data. More...
 
static Threads::atomic< unsigned int_n_objects
 The number of objects. More...
 
static Threads::spin_mutex _mutex
 Mutual exclusion object to enable thread-safe reference counting. More...
 
static bool _enable_print_counter = true
 Flag to control whether reference count information is printed when print_info is called. More...
 
static Counts _counts
 Actually holds the data. More...
 
static Threads::atomic< unsigned int_n_objects
 The number of objects. More...
 
static Threads::spin_mutex _mutex
 Mutual exclusion object to enable thread-safe reference counting. More...
 
static bool _enable_print_counter = true
 Flag to control whether reference count information is printed when print_info is called. More...
 

Detailed Description

Definition at line 86 of file assembly.h.

Member Typedef Documentation

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts
protectedinherited

Data structure to log the information.

The log is identified by the class name.

Definition at line 119 of file reference_counter.h.

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts
protectedinherited

Data structure to log the information.

The log is identified by the class name.

Definition at line 119 of file reference_counter.h.

Constructor & Destructor Documentation

AssemblyA1::AssemblyA1 ( ElasticityRBConstruction rb_sys_in)

Definition at line 88 of file assembly.h.

References libMesh::ElemAssembly::interior_assembly().

88  :
89  ElasticityAssembly(rb_sys_in)
90  {}
ElasticityAssembly(ElasticityRBConstruction &rb_sys_in)
Definition: assembly.h:49

Member Function Documentation

virtual void AssemblyA1::boundary_assembly ( FEMContext )
virtual

Perform the element boundary assembly.

Reimplemented from libMesh::ElemAssembly.

Definition at line 140 of file assembly.h.

References libMesh::BoundaryInfo::boundary_ids(), libMesh::MeshBase::get_boundary_info(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem(), libMesh::DiffContext::get_elem_jacobian(), libMesh::System::get_mesh(), libMesh::FEMContext::get_side_fe(), libMesh::FEMContext::get_side_qrule(), libmesh_nullptr, libMesh::QBase::n_points(), ElemAssemblyWithConstruction::rb_con, libMesh::Real, and libMesh::FEMContext::side.

141  {
142  std::vector<boundary_id_type> bc_ids;
143  rb_con->get_mesh().get_boundary_info().boundary_ids (&c.get_elem(), c.side, bc_ids);
144  for (std::vector<boundary_id_type>::const_iterator b =
145  bc_ids.begin(); b != bc_ids.end(); ++b)
146  if (*b == 1 || *b == 3) // y == -0.2, y == 0.2
147  {
148  const unsigned int u_var = 0;
149 
150  FEBase * side_fe = libmesh_nullptr;
151  c.get_side_fe(u_var, side_fe);
152 
153  const std::vector<Real> & JxW_side = side_fe->get_JxW();
154 
155  const std::vector<std::vector<Real>> & phi_side = side_fe->get_phi();
156 
157  const std::vector<Point> & xyz = side_fe->get_xyz();
158 
159  // The number of local degrees of freedom in each variable
160  const unsigned int n_u_dofs = c.get_dof_indices(u_var).size();
161 
162  // Now we will build the affine operator
163  unsigned int n_sidepoints = c.get_side_qrule().n_points();
164 
165  for (unsigned int qp=0; qp != n_sidepoints; qp++)
166  {
167  Real x_hat = xyz[qp](0);
168 
169  for (unsigned int i=0; i != n_u_dofs; i++)
170  for (unsigned int j=0; j != n_u_dofs; j++)
171  c.get_elem_jacobian()(i,j) += JxW_side[qp] * x_hat * phi_side[j][qp]*phi_side[i][qp];
172  }
173 
174  break;
175  }
176  }
const BoundaryInfo & get_boundary_info() const
The information about boundary ids on the mesh.
Definition: mesh_base.h:117
const class libmesh_nullptr_t libmesh_nullptr
std::vector< boundary_id_type > boundary_ids(const Node *node) const
const MeshBase & get_mesh() const
Definition: system.h:2014
FEGenericBase< Real > FEBase
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
RBConstruction * rb_con
Definition: assembly.h:48
void libMesh::ReferenceCounter::disable_print_counter_info ( )
staticinherited

Definition at line 107 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit(), and libMesh::ReferenceCounter::n_objects().

108 {
109  _enable_print_counter = false;
110  return;
111 }
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...
void libMesh::ReferenceCounter::disable_print_counter_info ( )
staticinherited

Definition at line 107 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit(), and libMesh::ReferenceCounter::n_objects().

108 {
109  _enable_print_counter = false;
110  return;
111 }
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...
void libMesh::ReferenceCounter::enable_print_counter_info ( )
staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 101 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::ReferenceCounter::n_objects().

102 {
103  _enable_print_counter = true;
104  return;
105 }
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...
void libMesh::ReferenceCounter::enable_print_counter_info ( )
staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 101 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::ReferenceCounter::n_objects().

102 {
103  _enable_print_counter = true;
104  return;
105 }
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...
std::string libMesh::ReferenceCounter::get_info ( )
staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

48 {
49 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
50 
51  std::ostringstream oss;
52 
53  oss << '\n'
54  << " ---------------------------------------------------------------------------- \n"
55  << "| Reference count information |\n"
56  << " ---------------------------------------------------------------------------- \n";
57 
58  for (Counts::iterator it = _counts.begin();
59  it != _counts.end(); ++it)
60  {
61  const std::string name(it->first);
62  const unsigned int creations = it->second.first;
63  const unsigned int destructions = it->second.second;
64 
65  oss << "| " << name << " reference count information:\n"
66  << "| Creations: " << creations << '\n'
67  << "| Destructions: " << destructions << '\n';
68  }
69 
70  oss << " ---------------------------------------------------------------------------- \n";
71 
72  return oss.str();
73 
74 #else
75 
76  return "";
77 
78 #endif
79 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
static Counts _counts
Actually holds the data.
std::string libMesh::ReferenceCounter::get_info ( )
staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

48 {
49 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
50 
51  std::ostringstream oss;
52 
53  oss << '\n'
54  << " ---------------------------------------------------------------------------- \n"
55  << "| Reference count information |\n"
56  << " ---------------------------------------------------------------------------- \n";
57 
58  for (Counts::iterator it = _counts.begin();
59  it != _counts.end(); ++it)
60  {
61  const std::string name(it->first);
62  const unsigned int creations = it->second.first;
63  const unsigned int destructions = it->second.second;
64 
65  oss << "| " << name << " reference count information:\n"
66  << "| Creations: " << creations << '\n'
67  << "| Destructions: " << destructions << '\n';
68  }
69 
70  oss << " ---------------------------------------------------------------------------- \n";
71 
72  return oss.str();
73 
74 #else
75 
76  return "";
77 
78 #endif
79 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
static Counts _counts
Actually holds the data.
virtual void libMesh::ElemAssembly::get_nodal_rhs_values ( std::map< numeric_index_type, Number > &  values,
const System ,
const Node  
)
virtualinherited

Get values to add to the RHS vector based on node.

This allows one to impose point loads, for example. Get values to add to the matrix or rhs vector based on node. This allows one to impose point loads or springs, for example.

Definition at line 88 of file elem_assembly.h.

References libMesh::ElemAssembly::is_nodal_rhs_values_overriden.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), AssemblyPointLoadX::AssemblyPointLoadX(), AssemblyPointLoadY::AssemblyPointLoadY(), and AssemblyPointLoadZ::AssemblyPointLoadZ().

91  {
92  // use get_nodal_values instead, since that enables matrix and vector assembly
93  libmesh_deprecated();
94 
95  // By default, just clear the values map
96  values.clear();
97 
98  // Set flag so that we know this is the default implementation
100  }
bool is_nodal_rhs_values_overriden
Temporary flag to help us figure out if we should call the deprecated get_nodal_rhs_values method or ...
virtual void libMesh::ElemAssembly::get_nodal_values ( std::vector< dof_id_type > &  ,
DenseMatrix< Number > &  ,
DenseVector< Number > &  ,
const System ,
const Node  
)
virtualinherited

Get values to add to the matrix or rhs vector based on node.

This allows one to impose point loads or springs, for example.

Definition at line 71 of file elem_assembly.h.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector().

76  {
77  // Do nothing by default
78  }
void libMesh::ReferenceCounter::increment_constructor_count ( const std::string &  name)
protectedinherited

Increments the construction counter.

Should be called in the constructor of any derived class that will be reference counted.

Definition at line 185 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

186 {
187  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
188  std::pair<unsigned int, unsigned int> & p = _counts[name];
189 
190  p.first++;
191 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
spin_mutex spin_mtx
A convenient spin mutex object which can be used for obtaining locks.
Definition: threads.C:29
static Counts _counts
Actually holds the data.
void libMesh::ReferenceCounter::increment_constructor_count ( const std::string &  name)
protectedinherited

Increments the construction counter.

Should be called in the constructor of any derived class that will be reference counted.

Definition at line 185 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

186 {
187  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
188  std::pair<unsigned int, unsigned int> & p = _counts[name];
189 
190  p.first++;
191 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
spin_mutex spin_mtx
A convenient spin mutex object which can be used for obtaining locks.
Definition: threads.C:29
static Counts _counts
Actually holds the data.
void libMesh::ReferenceCounter::increment_destructor_count ( const std::string &  name)
protectedinherited

Increments the destruction counter.

Should be called in the destructor of any derived class that will be reference counted.

Definition at line 198 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

199 {
200  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
201  std::pair<unsigned int, unsigned int> & p = _counts[name];
202 
203  p.second++;
204 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
spin_mutex spin_mtx
A convenient spin mutex object which can be used for obtaining locks.
Definition: threads.C:29
static Counts _counts
Actually holds the data.
void libMesh::ReferenceCounter::increment_destructor_count ( const std::string &  name)
protectedinherited

Increments the destruction counter.

Should be called in the destructor of any derived class that will be reference counted.

Definition at line 198 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

199 {
200  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
201  std::pair<unsigned int, unsigned int> & p = _counts[name];
202 
203  p.second++;
204 }
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:39
spin_mutex spin_mtx
A convenient spin mutex object which can be used for obtaining locks.
Definition: threads.C:29
static Counts _counts
Actually holds the data.
void AssemblyA1::interior_assembly ( FEMContext )
virtual

Perform the element interior assembly.

Reimplemented from libMesh::ElemAssembly.

Definition at line 114 of file assembly.C.

References elasticity_tensor(), libMesh::DiffContext::get_dof_indices(), libMesh::DiffContext::get_elem_jacobian(), libMesh::FEMContext::get_element_fe(), libMesh::FEMContext::get_element_qrule(), libmesh_nullptr, libMesh::QBase::n_points(), libMesh::System::n_vars(), ElasticityAssembly::rb_sys, libMesh::Real, ElasticityRBConstruction::u_var, ElasticityRBConstruction::v_var, and ElasticityRBConstruction::w_var.

115 {
116  const unsigned int n_components = rb_sys.n_vars();
117 
118  // make sure we have three components
119  libmesh_assert_equal_to (n_components, 3);
120 
121  const unsigned int u_var = rb_sys.u_var;
122  const unsigned int v_var = rb_sys.v_var;
123  const unsigned int w_var = rb_sys.w_var;
124 
125  FEBase * elem_fe = libmesh_nullptr;
126  c.get_element_fe(u_var, elem_fe);
127 
128  const std::vector<Real> & JxW = elem_fe->get_JxW();
129 
130  // The velocity shape function gradients at interior
131  // quadrature points.
132  const std::vector<std::vector<RealGradient>> & dphi = elem_fe->get_dphi();
133 
134  // Now we will build the affine operator
135  unsigned int n_qpoints = c.get_element_qrule().n_points();
136 
137  std::vector<unsigned int> n_var_dofs(n_components);
138  n_var_dofs[u_var] = c.get_dof_indices(u_var).size();
139  n_var_dofs[v_var] = c.get_dof_indices(v_var).size();
140  n_var_dofs[w_var] = c.get_dof_indices(w_var).size();
141 
142  for (unsigned int C_i = 0; C_i < n_components; C_i++)
143  for (unsigned int C_j = 1; C_j < n_components; C_j++)
144  for (unsigned int C_k = 0; C_k < n_components; C_k++)
145  for (unsigned int C_l = 1; C_l < n_components; C_l++)
146  {
147  Real C_ijkl = elasticity_tensor(C_i,C_j,C_k,C_l);
148  for (unsigned int qp=0; qp<n_qpoints; qp++)
149  for (unsigned int i=0; i<n_var_dofs[C_i]; i++)
150  for (unsigned int j=0; j<n_var_dofs[C_k]; j++)
151  (c.get_elem_jacobian(C_i,C_k))(i,j) +=
152  JxW[qp]*(C_ijkl * dphi[i][qp](C_j)*dphi[j][qp](C_l));
153  }
154 }
unsigned int u_var
Variable numbers.
Definition: rb_classes.h:121
ElasticityRBConstruction & rb_sys
The ElasticityRBConstruction object that will use this assembly.
Definition: assembly.h:56
Real elasticity_tensor(unsigned int i, unsigned int j, unsigned int k, unsigned int l)
Definition: assembly.C:38
const class libmesh_nullptr_t libmesh_nullptr
FEGenericBase< Real > FEBase
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
unsigned int n_vars() const
Definition: system.h:2086
static unsigned int libMesh::ReferenceCounter::n_objects ( )
staticinherited
static unsigned int libMesh::ReferenceCounter::n_objects ( )
staticinherited
void libMesh::ReferenceCounter::print_info ( std::ostream &  out = libMesh::out)
staticinherited

Prints the reference information, by default to libMesh::out.

Definition at line 88 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter, and libMesh::ReferenceCounter::get_info().

Referenced by libMesh::LibMeshInit::LibMeshInit().

89 {
91  out_stream << ReferenceCounter::get_info();
92 }
static std::string get_info()
Gets a string containing the reference information.
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...
void libMesh::ReferenceCounter::print_info ( std::ostream &  out = libMesh::out)
staticinherited

Prints the reference information, by default to libMesh::out.

Definition at line 88 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter, and libMesh::ReferenceCounter::get_info().

Referenced by libMesh::LibMeshInit::LibMeshInit().

89 {
91  out_stream << ReferenceCounter::get_info();
92 }
static std::string get_info()
Gets a string containing the reference information.
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called...

Member Data Documentation

ReferenceCounter::Counts libMesh::ReferenceCounter::_counts
staticprotectedinherited
ReferenceCounter::Counts libMesh::ReferenceCounter::_counts
staticprotectedinherited
bool libMesh::ReferenceCounter::_enable_print_counter = true
staticprotectedinherited

Flag to control whether reference count information is printed when print_info is called.

Definition at line 143 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), and libMesh::ReferenceCounter::print_info().

bool libMesh::ReferenceCounter::_enable_print_counter = true
staticprotectedinherited

Flag to control whether reference count information is printed when print_info is called.

Definition at line 143 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), and libMesh::ReferenceCounter::print_info().

Threads::spin_mutex libMesh::ReferenceCounter::_mutex
staticprotectedinherited

Mutual exclusion object to enable thread-safe reference counting.

Definition at line 137 of file reference_counter.h.

Threads::spin_mutex libMesh::ReferenceCounter::_mutex
staticprotectedinherited

Mutual exclusion object to enable thread-safe reference counting.

Definition at line 137 of file reference_counter.h.

Threads::atomic< unsigned int > libMesh::ReferenceCounter::_n_objects
staticprotectedinherited

The number of objects.

Print the reference count information when the number returns to 0.

Definition at line 132 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::n_objects(), libMesh::ReferenceCounter::ReferenceCounter(), and libMesh::ReferenceCounter::~ReferenceCounter().

Threads::atomic< unsigned int > libMesh::ReferenceCounter::_n_objects
staticprotectedinherited

The number of objects.

Print the reference count information when the number returns to 0.

Definition at line 132 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::n_objects(), libMesh::ReferenceCounter::ReferenceCounter(), and libMesh::ReferenceCounter::~ReferenceCounter().

bool libMesh::ElemAssembly::is_nodal_rhs_values_overriden
inherited

Temporary flag to help us figure out if we should call the deprecated get_nodal_rhs_values method or not.

Definition at line 107 of file elem_assembly.h.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and libMesh::ElemAssembly::get_nodal_rhs_values().

bool libMesh::ElemAssembly::is_nodal_rhs_values_overriden
inherited

Temporary flag to help us figure out if we should call the deprecated get_nodal_rhs_values method or not.

Definition at line 107 of file elem_assembly.h.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and libMesh::ElemAssembly::get_nodal_rhs_values().

RBConstruction* ElemAssemblyWithConstruction::rb_con
inherited
ElasticityRBConstruction& ElasticityAssembly::rb_sys
inherited

The documentation for this struct was generated from the following files: