libMesh::DiffContext Class Reference

#include <diff_context.h>

Inheritance diagram for libMesh::DiffContext:

Public Types

typedef std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > >::iterator localized_vectors_iterator
 

Public Member Functions

 DiffContext (const System &)
 
virtual ~DiffContext ()
 
virtual void elem_reinit (Real)
 
virtual void elem_side_reinit (Real)
 
virtual void elem_edge_reinit (Real)
 
virtual void nonlocal_reinit (Real)
 
unsigned int n_vars () const
 
const Systemget_system () const
 
const DenseVector< Number > & get_elem_solution () const
 
DenseVector< Number > & get_elem_solution ()
 
const DenseSubVector< Number > & get_elem_solution (unsigned int var) const
 
DenseSubVector< Number > & get_elem_solution (unsigned int var)
 
const DenseVector< Number > & get_elem_solution_rate () const
 
DenseVector< Number > & get_elem_solution_rate ()
 
const DenseSubVector< Number > & get_elem_solution_rate (unsigned int var) const
 
DenseSubVector< Number > & get_elem_solution_rate (unsigned int var)
 
const DenseVector< Number > & get_elem_solution_accel () const
 
DenseVector< Number > & get_elem_solution_accel ()
 
const DenseSubVector< Number > & get_elem_solution_accel (unsigned int var) const
 
DenseSubVector< Number > & get_elem_solution_accel (unsigned int var)
 
const DenseVector< Number > & get_elem_fixed_solution () const
 
DenseVector< Number > & get_elem_fixed_solution ()
 
const DenseSubVector< Number > & get_elem_fixed_solution (unsigned int var) const
 
DenseSubVector< Number > & get_elem_fixed_solution (unsigned int var)
 
const DenseVector< Number > & get_elem_residual () const
 
DenseVector< Number > & get_elem_residual ()
 
const DenseSubVector< Number > & get_elem_residual (unsigned int var) const
 
DenseSubVector< Number > & get_elem_residual (unsigned int var)
 
const DenseMatrix< Number > & get_elem_jacobian () const
 
DenseMatrix< Number > & get_elem_jacobian ()
 
const DenseSubMatrix< Number > & get_elem_jacobian (unsigned int var1, unsigned int var2) const
 
DenseSubMatrix< Number > & get_elem_jacobian (unsigned int var1, unsigned int var2)
 
const std::vector< Number > & get_qois () const
 
std::vector< Number > & get_qois ()
 
const std::vector< DenseVector< Number > > & get_qoi_derivatives () const
 
std::vector< DenseVector< Number > > & get_qoi_derivatives ()
 
const DenseSubVector< Number > & get_qoi_derivatives (unsigned int qoi, unsigned int var) const
 
DenseSubVector< Number > & get_qoi_derivatives (unsigned int qoi, unsigned int var)
 
const std::vector< dof_id_type > & get_dof_indices () const
 
std::vector< dof_id_type > & get_dof_indices ()
 
const std::vector< dof_id_type > & get_dof_indices (unsigned int var) const
 
std::vector< dof_id_type > & get_dof_indices (unsigned int var)
 
Real get_system_time () const
 
Real get_time () const
 
void set_time (Real time_in)
 
Real get_elem_solution_derivative () const
 
Real get_elem_solution_rate_derivative () const
 
Real get_elem_solution_accel_derivative () const
 
Real get_fixed_solution_derivative () const
 
bool is_adjoint () const
 
bool & is_adjoint ()
 
void set_deltat_pointer (Real *dt)
 
Real get_deltat_value ()
 
void add_localized_vector (NumericVector< Number > &localized_vector, const System &sys)
 
DenseVector< Number > & get_localized_vector (const NumericVector< Number > &localized_vector)
 
const DenseVector< Number > & get_localized_vector (const NumericVector< Number > &localized_vector) const
 
DenseSubVector< Number > & get_localized_subvector (const NumericVector< Number > &localized_vector, unsigned int var)
 
const DenseSubVector< Number > & get_localized_subvector (const NumericVector< Number > &localized_vector, unsigned int var) const
 

Public Attributes

Real time
 
const Real system_time
 
Real elem_solution_derivative
 
Real elem_solution_rate_derivative
 
Real elem_solution_accel_derivative
 
Real fixed_solution_derivative
 

Protected Attributes

std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
 
DenseVector< Number_elem_solution
 
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolutions
 
DenseVector< Number_elem_solution_rate
 
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_rates
 
DenseVector< Number_elem_solution_accel
 
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_accels
 
DenseVector< Number_elem_fixed_solution
 
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_fixed_subsolutions
 
DenseVector< Number_elem_residual
 
DenseMatrix< Number_elem_jacobian
 
std::vector< Number_elem_qoi
 
std::vector< DenseVector< Number > > _elem_qoi_derivative
 
std::vector< std::vector< std::unique_ptr< DenseSubVector< Number > > > > _elem_qoi_subderivatives
 
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subresiduals
 
std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > > _elem_subjacobians
 
std::vector< dof_id_type_dof_indices
 
std::vector< std::vector< dof_id_type > > _dof_indices_var
 

Private Attributes

Real_deltat
 
const System_system
 
bool _is_adjoint
 

Detailed Description

This class provides all data required for a physics package (e.g. a DifferentiableSystem subclass) to perform local element residual and jacobian integrations.

This class is part of the new DifferentiableSystem framework, which is still experimental. Users of this framework should beware of bugs and future API changes.

Author
Roy H. Stogner
Date
2009

Definition at line 55 of file diff_context.h.

Member Typedef Documentation

typedef std::map<const NumericVector<Number> *, std::pair<DenseVector<Number>, std::vector<std::unique_ptr<DenseSubVector<Number> > > > >::iterator libMesh::DiffContext::localized_vectors_iterator

Typedef for the localized_vectors iterator

Definition at line 528 of file diff_context.h.

Constructor & Destructor Documentation

libMesh::DiffContext::DiffContext ( const System sys)
explicit

Constructor. Optionally initializes required data structures.

Definition at line 30 of file diff_context.C.

References _elem_fixed_solution, _elem_fixed_subsolutions, _elem_jacobian, _elem_qoi, _elem_qoi_derivative, _elem_qoi_subderivatives, _elem_residual, _elem_solution, _elem_solution_accel, _elem_solution_rate, _elem_subjacobians, _elem_subresiduals, _elem_subsolution_accels, _elem_subsolution_rates, _elem_subsolutions, libMesh::DifferentiablePhysics::get_second_order_vars(), libMesh::DifferentiableSystem::get_time_solver(), libMesh::TimeSolver::is_steady(), libMesh::System::n_vars(), libMesh::System::qoi, libMesh::UnsteadySolver::time_order(), and libMesh::System::use_fixed_solution.

30  :
31  time(sys.time),
32  system_time(sys.time),
37  _dof_indices_var(sys.n_vars()),
39  _system(sys),
40  _is_adjoint(false)
41 {
42  // Finally initialize solution/residual/jacobian data structures
43  unsigned int nv = sys.n_vars();
44 
45  _elem_subsolutions.reserve(nv);
46  _elem_subresiduals.reserve(nv);
47  _elem_subjacobians.resize(nv);
48  _elem_subsolution_rates.reserve(nv);
49  _elem_subsolution_accels.reserve(nv);
50  if (sys.use_fixed_solution)
51  _elem_fixed_subsolutions.reserve(nv);
52 
53  // If the user resizes sys.qoi, it will invalidate us
54 
55  std::size_t n_qoi = sys.qoi.size();
56  _elem_qoi.resize(n_qoi);
57  _elem_qoi_derivative.resize(n_qoi);
58  _elem_qoi_subderivatives.resize(n_qoi);
59  for (std::size_t q=0; q != n_qoi; ++q)
60  _elem_qoi_subderivatives[q].reserve(nv);
61 
62  for (unsigned int i=0; i != nv; ++i)
63  {
64  _elem_subsolutions.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_solution));
65  _elem_subresiduals.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_residual));
66  for (std::size_t q=0; q != n_qoi; ++q)
67  _elem_qoi_subderivatives[q].emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_qoi_derivative[q]));
68  _elem_subjacobians[i].reserve(nv);
69 
70  // Only make space for these if we're using DiffSystem
71  // This is assuming *only* DiffSystem is using elem_solution_rate/accel
72  const DifferentiableSystem * diff_system = dynamic_cast<const DifferentiableSystem *>(&sys);
73  if (diff_system)
74  {
75  // Now, we only need these if the solver is unsteady
76  if (!diff_system->get_time_solver().is_steady())
77  {
78  _elem_subsolution_rates.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_solution_rate));
79 
80  // We only need accel space if the TimeSolver is second order
81  const UnsteadySolver & time_solver = cast_ref<const UnsteadySolver &>(diff_system->get_time_solver());
82 
83  if (time_solver.time_order() >= 2 || !diff_system->get_second_order_vars().empty())
84  _elem_subsolution_accels.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_solution_accel));
85  }
86  }
87 
88  if (sys.use_fixed_solution)
89  _elem_fixed_subsolutions.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_elem_fixed_solution));
90 
91  for (unsigned int j=0; j != nv; ++j)
92  _elem_subjacobians[i].emplace_back(libmesh_make_unique<DenseSubMatrix<Number>>(_elem_jacobian));
93  }
94 }
std::vector< std::vector< dof_id_type > > _dof_indices_var
Definition: diff_context.h:622
DenseVector< Number > _elem_solution
Definition: diff_context.h:565
std::vector< Number > _elem_qoi
Definition: diff_context.h:604
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_rates
Definition: diff_context.h:573
DenseVector< Number > _elem_fixed_solution
Definition: diff_context.h:587
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subresiduals
Definition: diff_context.h:615
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_accels
Definition: diff_context.h:580
std::vector< DenseVector< Number > > _elem_qoi_derivative
Definition: diff_context.h:609
DenseMatrix< Number > _elem_jacobian
Definition: diff_context.h:599
std::vector< std::vector< std::unique_ptr< DenseSubVector< Number > > > > _elem_qoi_subderivatives
Definition: diff_context.h:610
const class libmesh_nullptr_t libmesh_nullptr
Real elem_solution_rate_derivative
Definition: diff_context.h:491
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolutions
Definition: diff_context.h:566
const System & _system
Definition: diff_context.h:642
std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > > _elem_subjacobians
Definition: diff_context.h:616
const Real system_time
Definition: diff_context.h:478
DenseVector< Number > _elem_residual
Definition: diff_context.h:593
DenseVector< Number > _elem_solution_accel
Definition: diff_context.h:579
Real elem_solution_accel_derivative
Definition: diff_context.h:498
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_fixed_subsolutions
Definition: diff_context.h:588
DenseVector< Number > _elem_solution_rate
Definition: diff_context.h:572
libMesh::DiffContext::~DiffContext ( )
virtual

Destructor.

Definition at line 98 of file diff_context.C.

99 {
100 }

Member Function Documentation

void libMesh::DiffContext::add_localized_vector ( NumericVector< Number > &  localized_vector,
const System sys 
)

Adds a vector to the map of localized vectors. We can later evaluate interior_values, interior_gradients and side_values for these fields these vectors represent.

Definition at line 119 of file diff_context.C.

References _localized_vectors, and libMesh::System::n_vars().

120 {
121  // Make an empty pair keyed with a reference to this _localized_vector
122  _localized_vectors[&localized_vector] = std::make_pair(DenseVector<Number>(), std::vector<std::unique_ptr<DenseSubVector<Number>>>());
123 
124  unsigned int nv = sys.n_vars();
125 
126  _localized_vectors[&localized_vector].second.reserve(nv);
127 
128  // Fill the DenseSubVector with nv copies of DenseVector
129  for (unsigned int i=0; i != nv; ++i)
130  _localized_vectors[&localized_vector].second.emplace_back(libmesh_make_unique<DenseSubVector<Number>>(_localized_vectors[&localized_vector].first));
131 }
std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
Definition: diff_context.h:559
virtual void libMesh::DiffContext::elem_edge_reinit ( Real  )
inlinevirtual

Gives derived classes the opportunity to reinitialize data needed for an edge integration at a new point within a timestep

Reimplemented in libMesh::FEMContext.

Definition at line 88 of file diff_context.h.

88 {}
virtual void libMesh::DiffContext::elem_reinit ( Real  )
inlinevirtual

Gives derived classes the opportunity to reinitialize data (FE objects in FEMSystem, for example) needed for an interior integration at a new point within a timestep

Reimplemented in libMesh::FEMContext.

Definition at line 76 of file diff_context.h.

Referenced by libMesh::EulerSolver::element_residual(), libMesh::Euler2Solver::element_residual(), and libMesh::NewmarkSolver::element_residual().

76 {}
virtual void libMesh::DiffContext::elem_side_reinit ( Real  )
inlinevirtual

Gives derived classes the opportunity to reinitialize data needed for a side integration at a new point within a timestep

Reimplemented in libMesh::FEMContext.

Definition at line 82 of file diff_context.h.

Referenced by libMesh::EulerSolver::side_residual(), libMesh::Euler2Solver::side_residual(), and libMesh::NewmarkSolver::side_residual().

82 {}
Real libMesh::DiffContext::get_deltat_value ( )
Returns
The value currently pointed to by this class's _deltat member

Definition at line 111 of file diff_context.C.

References _deltat.

Referenced by libMesh::FEMContext::_update_time_from_system().

112 {
113  libmesh_assert(_deltat);
114 
115  return *_deltat;
116 }
const std::vector<dof_id_type>& libMesh::DiffContext::get_dof_indices ( ) const
inline

Accessor for element dof indices

Definition at line 367 of file diff_context.h.

References _dof_indices.

Referenced by libMesh::EulerSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::NewmarkSolver::_general_residual(), libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), libMesh::FEMContext::fixed_point_gradient(), libMesh::FEMContext::fixed_point_hessian(), libMesh::FEMContext::fixed_point_value(), libMesh::VectorSetAction< Val >::insert(), libMesh::MatrixFillAction< ValIn, ValOut >::insert(), libMesh::FEMContext::interior_curl(), libMesh::FEMContext::interior_div(), libMesh::FEMContext::interior_gradients(), libMesh::FEMContext::interior_hessians(), libMesh::FEMContext::interior_values(), libMesh::FEMSystem::mesh_position_get(), libMesh::DGFEMContext::neighbor_side_fe_reinit(), libMesh::FEMSystem::numerical_jacobian(), libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::operator()(), libMesh::FEMContext::point_curl(), libMesh::FEMContext::point_gradient(), libMesh::FEMContext::point_hessian(), libMesh::FEMContext::point_value(), libMesh::FEMContext::pre_fe_reinit(), libMesh::FEMContext::side_gradient(), libMesh::FEMContext::side_gradients(), libMesh::FEMContext::side_hessians(), libMesh::FEMContext::side_values(), libMesh::FEMContext::some_gradient(), libMesh::FEMContext::some_hessian(), and libMesh::FEMContext::some_value().

368  { return _dof_indices; }
std::vector< dof_id_type > _dof_indices
Definition: diff_context.h:621
std::vector<dof_id_type>& libMesh::DiffContext::get_dof_indices ( )
inline

Non-const accessor for element dof indices

Definition at line 373 of file diff_context.h.

References _dof_indices.

374  { return _dof_indices; }
std::vector< dof_id_type > _dof_indices
Definition: diff_context.h:621
const std::vector<dof_id_type>& libMesh::DiffContext::get_dof_indices ( unsigned int  var) const
inline

Accessor for element dof indices of a particular variable corresponding to the index argument.

Definition at line 380 of file diff_context.h.

References _dof_indices_var.

381  {
382  libmesh_assert_greater(_dof_indices_var.size(), var);
383  return _dof_indices_var[var];
384  }
std::vector< std::vector< dof_id_type > > _dof_indices_var
Definition: diff_context.h:622
std::vector<dof_id_type>& libMesh::DiffContext::get_dof_indices ( unsigned int  var)
inline

Accessor for element dof indices of a particular variable corresponding to the index argument.

Definition at line 390 of file diff_context.h.

References _dof_indices_var.

391  {
392  libmesh_assert_greater(_dof_indices_var.size(), var);
393  return _dof_indices_var[var];
394  }
std::vector< std::vector< dof_id_type > > _dof_indices_var
Definition: diff_context.h:622
DenseVector<Number>& libMesh::DiffContext::get_elem_fixed_solution ( )
inline

Non-const accessor for element fixed solution.

Definition at line 221 of file diff_context.h.

References _elem_fixed_solution.

222  { return _elem_fixed_solution; }
DenseVector< Number > _elem_fixed_solution
Definition: diff_context.h:587
const DenseSubVector<Number>& libMesh::DiffContext::get_elem_fixed_solution ( unsigned int  var) const
inline

Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

Definition at line 228 of file diff_context.h.

References _elem_fixed_subsolutions.

229  {
230  libmesh_assert_greater(_elem_fixed_subsolutions.size(), var);
231  libmesh_assert(_elem_fixed_subsolutions[var]);
232  return *(_elem_fixed_subsolutions[var]);
233  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_fixed_subsolutions
Definition: diff_context.h:588
DenseSubVector<Number>& libMesh::DiffContext::get_elem_fixed_solution ( unsigned int  var)
inline

Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

Definition at line 239 of file diff_context.h.

References _elem_fixed_subsolutions.

240  {
241  libmesh_assert_greater(_elem_fixed_subsolutions.size(), var);
242  libmesh_assert(_elem_fixed_subsolutions[var]);
243  return *(_elem_fixed_subsolutions[var]);
244  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_fixed_subsolutions
Definition: diff_context.h:588
DenseMatrix<Number>& libMesh::DiffContext::get_elem_jacobian ( )
inline

Non-const accessor for element Jacobian.

Definition at line 289 of file diff_context.h.

References _elem_jacobian.

290  { return _elem_jacobian; }
DenseMatrix< Number > _elem_jacobian
Definition: diff_context.h:599
const DenseSubMatrix<Number>& libMesh::DiffContext::get_elem_jacobian ( unsigned int  var1,
unsigned int  var2 
) const
inline

Const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 296 of file diff_context.h.

References _elem_subjacobians.

297  {
298  libmesh_assert_greater(_elem_subjacobians.size(), var1);
299  libmesh_assert_greater(_elem_subjacobians[var1].size(), var2);
300  libmesh_assert(_elem_subjacobians[var1][var2]);
301  return *(_elem_subjacobians[var1][var2]);
302  }
std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > > _elem_subjacobians
Definition: diff_context.h:616
DenseSubMatrix<Number>& libMesh::DiffContext::get_elem_jacobian ( unsigned int  var1,
unsigned int  var2 
)
inline

Non-const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 308 of file diff_context.h.

References _elem_subjacobians.

309  {
310  libmesh_assert_greater(_elem_subjacobians.size(), var1);
311  libmesh_assert_greater(_elem_subjacobians[var1].size(), var2);
312  libmesh_assert(_elem_subjacobians[var1][var2]);
313  return *(_elem_subjacobians[var1][var2]);
314  }
std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > > _elem_subjacobians
Definition: diff_context.h:616
const DenseVector<Number>& libMesh::DiffContext::get_elem_residual ( ) const
inline
DenseVector<Number>& libMesh::DiffContext::get_elem_residual ( )
inline

Non-const accessor for element residual.

Definition at line 255 of file diff_context.h.

References _elem_residual.

256  { return _elem_residual; }
DenseVector< Number > _elem_residual
Definition: diff_context.h:593
const DenseSubVector<Number>& libMesh::DiffContext::get_elem_residual ( unsigned int  var) const
inline

Const accessor for element residual of a particular variable corresponding to the variable index argument.

Definition at line 262 of file diff_context.h.

References _elem_subresiduals.

263  {
264  libmesh_assert_greater(_elem_subresiduals.size(), var);
265  libmesh_assert(_elem_subresiduals[var]);
266  return *(_elem_subresiduals[var]);
267  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subresiduals
Definition: diff_context.h:615
DenseSubVector<Number>& libMesh::DiffContext::get_elem_residual ( unsigned int  var)
inline

Non-const accessor for element residual of a particular variable corresponding to the variable index argument.

Definition at line 273 of file diff_context.h.

References _elem_subresiduals.

274  {
275  libmesh_assert_greater(_elem_subresiduals.size(), var);
276  libmesh_assert(_elem_subresiduals[var]);
277  return *(_elem_subresiduals[var]);
278  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subresiduals
Definition: diff_context.h:615
DenseVector<Number>& libMesh::DiffContext::get_elem_solution ( )
inline

Non-const accessor for element solution.

Definition at line 117 of file diff_context.h.

References _elem_solution.

118  { return _elem_solution; }
DenseVector< Number > _elem_solution
Definition: diff_context.h:565
const DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution ( unsigned int  var) const
inline

Accessor for element solution of a particular variable corresponding to the variable index argument.

Definition at line 124 of file diff_context.h.

References _elem_subsolutions.

125  {
126  libmesh_assert_greater(_elem_subsolutions.size(), var);
127  libmesh_assert(_elem_subsolutions[var]);
128  return *(_elem_subsolutions[var]);
129  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolutions
Definition: diff_context.h:566
DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution ( unsigned int  var)
inline

Accessor for element solution of a particular variable corresponding to the variable index argument.

Definition at line 135 of file diff_context.h.

References _elem_subsolutions.

136  {
137  libmesh_assert_greater(_elem_subsolutions.size(), var);
138  libmesh_assert(_elem_subsolutions[var]);
139  return *(_elem_subsolutions[var]);
140  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolutions
Definition: diff_context.h:566
const DenseVector<Number>& libMesh::DiffContext::get_elem_solution_accel ( ) const
inline
DenseVector<Number>& libMesh::DiffContext::get_elem_solution_accel ( )
inline

Non-const accessor for element solution accel of change w.r.t. time.

Definition at line 187 of file diff_context.h.

References _elem_solution_accel.

188  { return _elem_solution_accel; }
DenseVector< Number > _elem_solution_accel
Definition: diff_context.h:579
const DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution_accel ( unsigned int  var) const
inline

Accessor for element solution accel for a particular variable corresponding to the variable index argument.

Definition at line 194 of file diff_context.h.

References _elem_subsolution_accels.

195  {
196  libmesh_assert_greater(_elem_subsolution_accels.size(), var);
197  libmesh_assert(_elem_subsolution_accels[var]);
198  return *(_elem_subsolution_accels[var]);
199  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_accels
Definition: diff_context.h:580
DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution_accel ( unsigned int  var)
inline

Accessor for element solution accel for a particular variable corresponding to the variable index argument.

Definition at line 205 of file diff_context.h.

References _elem_subsolution_accels.

206  {
207  libmesh_assert_greater(_elem_subsolution_accels.size(), var);
208  libmesh_assert(_elem_subsolution_accels[var]);
209  return *(_elem_subsolution_accels[var]);
210  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_accels
Definition: diff_context.h:580
Real libMesh::DiffContext::get_elem_solution_accel_derivative ( ) const
inline

The derivative of the current elem_solution_accel w.r.t. the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_accel_derivative() is 0.

Definition at line 438 of file diff_context.h.

References elem_solution_accel_derivative.

Real elem_solution_accel_derivative
Definition: diff_context.h:498
Real libMesh::DiffContext::get_elem_solution_derivative ( ) const
inline

The derivative of the current elem_solution w.r.t. the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_derivative() is 0.

Definition at line 420 of file diff_context.h.

References elem_solution_derivative.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns().

421  { return elem_solution_derivative; }
const DenseVector<Number>& libMesh::DiffContext::get_elem_solution_rate ( ) const
inline
DenseVector<Number>& libMesh::DiffContext::get_elem_solution_rate ( )
inline

Non-const accessor for element solution rate of change w.r.t. time.

Definition at line 152 of file diff_context.h.

References _elem_solution_rate.

153  { return _elem_solution_rate; }
DenseVector< Number > _elem_solution_rate
Definition: diff_context.h:572
const DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution_rate ( unsigned int  var) const
inline

Accessor for element solution rate for a particular variable corresponding to the variable index argument.

Definition at line 159 of file diff_context.h.

References _elem_subsolution_rates.

160  {
161  libmesh_assert_greater(_elem_subsolution_rates.size(), var);
162  libmesh_assert(_elem_subsolution_rates[var]);
163  return *(_elem_subsolution_rates[var]);
164  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_rates
Definition: diff_context.h:573
DenseSubVector<Number>& libMesh::DiffContext::get_elem_solution_rate ( unsigned int  var)
inline

Accessor for element solution rate for a particular variable corresponding to the variable index argument.

Definition at line 170 of file diff_context.h.

References _elem_subsolution_rates.

171  {
172  libmesh_assert_greater(_elem_subsolution_rates.size(), var);
173  libmesh_assert(_elem_subsolution_rates[var]);
174  return *(_elem_subsolution_rates[var]);
175  }
std::vector< std::unique_ptr< DenseSubVector< Number > > > _elem_subsolution_rates
Definition: diff_context.h:573
Real libMesh::DiffContext::get_elem_solution_rate_derivative ( ) const
inline

The derivative of the current elem_solution_rate w.r.t. the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_rate_derivative() is 0.

Definition at line 429 of file diff_context.h.

References elem_solution_rate_derivative.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), and libMesh::FirstOrderUnsteadySolver::prepare_accel().

Real elem_solution_rate_derivative
Definition: diff_context.h:491
Real libMesh::DiffContext::get_fixed_solution_derivative ( ) const
inline

The derivative of the current fixed_elem_solution w.r.t. the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_fixed_elem_solution_derivative() is 0.

Definition at line 447 of file diff_context.h.

References fixed_solution_derivative.

448  { return fixed_solution_derivative; }
DenseSubVector< Number > & libMesh::DiffContext::get_localized_subvector ( const NumericVector< Number > &  localized_vector,
unsigned int  var 
)

Return a reference to DenseSubVector localization of localized_vector at variable var contained in the _localized_vectors map

Definition at line 148 of file diff_context.C.

References _localized_vectors.

Referenced by libMesh::FEMContext::interior_gradients(), libMesh::FEMContext::interior_hessians(), libMesh::FEMContext::interior_values(), libMesh::FEMContext::side_gradients(), libMesh::FEMContext::side_hessians(), and libMesh::FEMContext::side_values().

149 {
150  return *_localized_vectors[&localized_vector].second[var];
151 }
std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
Definition: diff_context.h:559
const DenseSubVector< Number > & libMesh::DiffContext::get_localized_subvector ( const NumericVector< Number > &  localized_vector,
unsigned int  var 
) const

const accessible version of get_localized_subvector function

Definition at line 154 of file diff_context.C.

References _localized_vectors.

155 {
156  auto localized_vectors_it = _localized_vectors.find(&localized_vector);
157  libmesh_assert(localized_vectors_it != _localized_vectors.end());
158  return *localized_vectors_it->second.second[var];
159 }
std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
Definition: diff_context.h:559
DenseVector< Number > & libMesh::DiffContext::get_localized_vector ( const NumericVector< Number > &  localized_vector)

Return a reference to DenseVector localization of localized_vector contained in the _localized_vectors map

Definition at line 134 of file diff_context.C.

References _localized_vectors.

135 {
136  return _localized_vectors[&localized_vector].first;
137 }
std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
Definition: diff_context.h:559
const DenseVector< Number > & libMesh::DiffContext::get_localized_vector ( const NumericVector< Number > &  localized_vector) const

const accessible version of get_localized_vector function

Definition at line 140 of file diff_context.C.

References _localized_vectors.

141 {
142  auto localized_vectors_it = _localized_vectors.find(&localized_vector);
143  libmesh_assert(localized_vectors_it != _localized_vectors.end());
144  return localized_vectors_it->second.first;
145 }
std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< std::unique_ptr< DenseSubVector< Number > > > > > _localized_vectors
Definition: diff_context.h:559
const std::vector<DenseVector<Number> >& libMesh::DiffContext::get_qoi_derivatives ( ) const
inline

Const accessor for QoI derivatives.

Definition at line 331 of file diff_context.h.

References _elem_qoi_derivative.

Referenced by libMesh::FEMContext::pre_fe_reinit().

332  { return _elem_qoi_derivative; }
std::vector< DenseVector< Number > > _elem_qoi_derivative
Definition: diff_context.h:609
std::vector<DenseVector<Number> >& libMesh::DiffContext::get_qoi_derivatives ( )
inline

Non-const accessor for QoI derivatives.

Definition at line 337 of file diff_context.h.

References _elem_qoi_derivative.

338  { return _elem_qoi_derivative; }
std::vector< DenseVector< Number > > _elem_qoi_derivative
Definition: diff_context.h:609
const DenseSubVector<Number>& libMesh::DiffContext::get_qoi_derivatives ( unsigned int  qoi,
unsigned int  var 
) const
inline

Const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

Definition at line 344 of file diff_context.h.

References _elem_qoi_subderivatives.

345  {
346  libmesh_assert_greater(_elem_qoi_subderivatives.size(), qoi);
347  libmesh_assert_greater(_elem_qoi_subderivatives[qoi].size(), var);
348  libmesh_assert(_elem_qoi_subderivatives[qoi][var]);
349  return *(_elem_qoi_subderivatives[qoi][var]);
350  }
std::vector< std::vector< std::unique_ptr< DenseSubVector< Number > > > > _elem_qoi_subderivatives
Definition: diff_context.h:610
DenseSubVector<Number>& libMesh::DiffContext::get_qoi_derivatives ( unsigned int  qoi,
unsigned int  var 
)
inline

Non-const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

Definition at line 356 of file diff_context.h.

References _elem_qoi_subderivatives.

357  {
358  libmesh_assert_greater(_elem_qoi_subderivatives.size(), qoi);
359  libmesh_assert_greater(_elem_qoi_subderivatives[qoi].size(), var);
360  libmesh_assert(_elem_qoi_subderivatives[qoi][var]);
361  return *(_elem_qoi_subderivatives[qoi][var]);
362  }
std::vector< std::vector< std::unique_ptr< DenseSubVector< Number > > > > _elem_qoi_subderivatives
Definition: diff_context.h:610
const std::vector<Number>& libMesh::DiffContext::get_qois ( ) const
inline

Const accessor for QoI vector.

Definition at line 319 of file diff_context.h.

References _elem_qoi.

320  { return _elem_qoi; }
std::vector< Number > _elem_qoi
Definition: diff_context.h:604
std::vector<Number>& libMesh::DiffContext::get_qois ( )
inline

Non-const accessor for QoI vector.

Definition at line 325 of file diff_context.h.

References _elem_qoi.

326  { return _elem_qoi; }
std::vector< Number > _elem_qoi
Definition: diff_context.h:604
const System& libMesh::DiffContext::get_system ( ) const
inline

Accessor for associated system.

Definition at line 105 of file diff_context.h.

References _system.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), and libMesh::DGFEMContext::neighbor_side_fe_reinit().

106  { return _system; }
const System & _system
Definition: diff_context.h:642
Real libMesh::DiffContext::get_system_time ( ) const
inline

Accessor for the time variable stored in the system class.

Definition at line 399 of file diff_context.h.

References system_time.

Referenced by libMesh::FEMContext::_update_time_from_system().

400  { return system_time; }
const Real system_time
Definition: diff_context.h:478
Real libMesh::DiffContext::get_time ( ) const
inline

Accessor for the time for which the current nonlinear_solution is defined.

Definition at line 405 of file diff_context.h.

References time.

406  { return time; }
bool libMesh::DiffContext::is_adjoint ( ) const
inline

Accessor for querying whether we need to do a primal or adjoint solve

Definition at line 454 of file diff_context.h.

References _is_adjoint.

Referenced by libMesh::FEMSystem::build_context().

455  { return _is_adjoint; }
bool& libMesh::DiffContext::is_adjoint ( )
inline

Accessor for setting whether we need to do a primal or adjoint solve

Definition at line 461 of file diff_context.h.

References _is_adjoint.

462  { return _is_adjoint; }
unsigned int libMesh::DiffContext::n_vars ( ) const
inline
virtual void libMesh::DiffContext::nonlocal_reinit ( Real  )
inlinevirtual

Gives derived classes the opportunity to reinitialize data needed for nonlocal calculations at a new point within a timestep

Reimplemented in libMesh::FEMContext.

Definition at line 94 of file diff_context.h.

Referenced by libMesh::EulerSolver::nonlocal_residual(), libMesh::Euler2Solver::nonlocal_residual(), and libMesh::NewmarkSolver::nonlocal_residual().

94 {}
void libMesh::DiffContext::set_deltat_pointer ( Real dt)

Points the _deltat member of this class at a timestep value stored in the creating System, for example DiffSystem::deltat

Definition at line 103 of file diff_context.C.

References _deltat.

Referenced by libMesh::FEMSystem::build_context(), libMesh::DifferentiableSystem::build_context(), and libMesh::FEMSystem::init_context().

104 {
105  // We may actually want to be able to set this pointer to NULL, so
106  // don't report an error for that.
107  _deltat = dt;
108 }
void libMesh::DiffContext::set_time ( Real  time_in)
inline

Set the time for which the current nonlinear_solution is defined.

Definition at line 411 of file diff_context.h.

References time.

Referenced by libMesh::FEMContext::_update_time_from_system().

412  { time = time_in; }

Member Data Documentation

Real* libMesh::DiffContext::_deltat
private

Default NULL, can optionally be used to point to a timestep value in the System-derived class responsible for creating this DiffContext.

In DiffSystem's build_context() function, is assigned to point to the deltat member of that class.

Accessible via public get_deltat()/set_deltat() methods of this class.

Always test for NULL before using!

Definition at line 637 of file diff_context.h.

Referenced by get_deltat_value(), and set_deltat_pointer().

std::vector<dof_id_type> libMesh::DiffContext::_dof_indices
protected

Global Degree of freedom index lists

Definition at line 621 of file diff_context.h.

Referenced by get_dof_indices().

std::vector<std::vector<dof_id_type> > libMesh::DiffContext::_dof_indices_var
protected

Definition at line 622 of file diff_context.h.

Referenced by get_dof_indices(), and n_vars().

DenseVector<Number> libMesh::DiffContext::_elem_fixed_solution
protected

Element by element components of nonlinear_solution at a fixed point in a timestep, for optional use by e.g. stabilized methods

Definition at line 587 of file diff_context.h.

Referenced by DiffContext(), and get_elem_fixed_solution().

std::vector<std::unique_ptr<DenseSubVector<Number> > > libMesh::DiffContext::_elem_fixed_subsolutions
protected
DenseMatrix<Number> libMesh::DiffContext::_elem_jacobian
protected

Element jacobian: derivatives of elem_residual with respect to elem_solution

Definition at line 599 of file diff_context.h.

Referenced by DiffContext(), and get_elem_jacobian().

std::vector<Number> libMesh::DiffContext::_elem_qoi
protected

Element quantity of interest contributions

Definition at line 604 of file diff_context.h.

Referenced by DiffContext(), and get_qois().

std::vector<DenseVector<Number> > libMesh::DiffContext::_elem_qoi_derivative
protected

Element quantity of interest derivative contributions

Definition at line 609 of file diff_context.h.

Referenced by DiffContext(), and get_qoi_derivatives().

std::vector<std::vector<std::unique_ptr<DenseSubVector<Number> > > > libMesh::DiffContext::_elem_qoi_subderivatives
protected
DenseVector<Number> libMesh::DiffContext::_elem_residual
protected

Element residual vector

Definition at line 593 of file diff_context.h.

Referenced by DiffContext(), and get_elem_residual().

DenseVector<Number> libMesh::DiffContext::_elem_solution
protected

Element by element components of nonlinear_solution as adjusted by a time_solver

Definition at line 565 of file diff_context.h.

Referenced by DiffContext(), and get_elem_solution().

DenseVector<Number> libMesh::DiffContext::_elem_solution_accel
protected

Element by element components of du/dt as adjusted by a time_solver

Definition at line 579 of file diff_context.h.

Referenced by DiffContext(), and get_elem_solution_accel().

DenseVector<Number> libMesh::DiffContext::_elem_solution_rate
protected

Element by element components of du/dt as adjusted by a time_solver

Definition at line 572 of file diff_context.h.

Referenced by DiffContext(), and get_elem_solution_rate().

std::vector<std::vector<std::unique_ptr<DenseSubMatrix<Number> > > > libMesh::DiffContext::_elem_subjacobians
protected

Definition at line 616 of file diff_context.h.

Referenced by DiffContext(), and get_elem_jacobian().

std::vector<std::unique_ptr<DenseSubVector<Number> > > libMesh::DiffContext::_elem_subresiduals
protected

Element residual subvectors and Jacobian submatrices

Definition at line 615 of file diff_context.h.

Referenced by DiffContext(), and get_elem_residual().

std::vector<std::unique_ptr<DenseSubVector<Number> > > libMesh::DiffContext::_elem_subsolution_accels
protected

Definition at line 580 of file diff_context.h.

Referenced by DiffContext(), and get_elem_solution_accel().

std::vector<std::unique_ptr<DenseSubVector<Number> > > libMesh::DiffContext::_elem_subsolution_rates
protected

Definition at line 573 of file diff_context.h.

Referenced by DiffContext(), and get_elem_solution_rate().

bool libMesh::DiffContext::_is_adjoint
private

Is this context to be used for a primal or adjoint solve?

Definition at line 647 of file diff_context.h.

Referenced by is_adjoint().

std::map<const NumericVector<Number> *, std::pair<DenseVector<Number>, std::vector<std::unique_ptr<DenseSubVector<Number> > > > > libMesh::DiffContext::_localized_vectors
protected

Contains pointers to vectors the user has asked to be localized, keyed with pairs of element localized versions of that vector and per variable views

Definition at line 559 of file diff_context.h.

Referenced by add_localized_vector(), get_localized_subvector(), get_localized_vector(), and libMesh::FEMContext::pre_fe_reinit().

const System& libMesh::DiffContext::_system
private

A reference to the system this context is constructed with

Definition at line 642 of file diff_context.h.

Referenced by get_system().

Real libMesh::DiffContext::elem_solution_accel_derivative

The derivative of elem_solution_accel with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Definition at line 498 of file diff_context.h.

Referenced by libMesh::NewmarkSolver::_general_residual(), get_elem_solution_accel_derivative(), and libMesh::FirstOrderUnsteadySolver::prepare_accel().

Real libMesh::DiffContext::elem_solution_derivative
Real libMesh::DiffContext::elem_solution_rate_derivative

The derivative of elem_solution_rate with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Definition at line 491 of file diff_context.h.

Referenced by libMesh::EulerSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::NewmarkSolver::_general_residual(), libMesh::EigenTimeSolver::element_residual(), get_elem_solution_rate_derivative(), and libMesh::EigenTimeSolver::side_residual().

Real libMesh::DiffContext::fixed_solution_derivative

The derivative of elem_fixed_solution with respect to the nonlinear solution, for use by systems constructing jacobians with elem_fixed_solution based methods

Definition at line 505 of file diff_context.h.

Referenced by libMesh::EulerSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::SteadySolver::_general_residual(), and get_fixed_solution_derivative().

const Real libMesh::DiffContext::system_time

This is the time stored in the System class at the time this context was created, i.e. the time at the beginning of the current timestep. This value gets set in the constructor and unlike DiffContext::time, is not tweaked mid-timestep by transient solvers: it remains equal to the value it was assigned at construction.

Definition at line 478 of file diff_context.h.

Referenced by get_system_time().

Real libMesh::DiffContext::time

For time-dependent problems, this is the time t for which the current nonlinear_solution is defined. FIXME - this needs to be tweaked mid-timestep by all transient solvers!

Definition at line 469 of file diff_context.h.

Referenced by get_time(), libMesh::WeightedPatchRecoveryErrorEstimator::EstimateError::operator()(), and set_time().


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