ASPECT
|
Public Member Functions | |
AdvectionSystem (const FiniteElement< dim > &finite_element, const bool field_is_discontinuous) | |
Public Member Functions inherited from aspect::internal::Assembly::CopyData::CopyDataBase< dim > | |
virtual | ~CopyDataBase ()=default |
Public Attributes | |
FullMatrix< double > | local_matrix |
std::vector< FullMatrix< double > > | local_matrices_int_ext |
std::vector< FullMatrix< double > > | local_matrices_ext_int |
std::vector< FullMatrix< double > > | local_matrices_ext_ext |
Vector< double > | local_rhs |
std::vector< bool > | assembled_matrices |
std::vector< types::global_dof_index > | local_dof_indices |
std::vector< std::vector< types::global_dof_index > > | neighbor_dof_indices |
Additionally to the Stokes system the Advection system copy data object also needs to keep track of contributions across faces (mostly for discontinuous elements that contain DG terms).
Definition at line 116 of file simulator.h.
aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::AdvectionSystem | ( | const FiniteElement< dim > & | finite_element, |
const bool | field_is_discontinuous | ||
) |
Constructor.
finite_element | The element that describes the field for which we are trying to assemble a linear system. Not the global finite element. |
field_is_discontinuous | If true, the field is a DG element. |
FullMatrix<double> aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_matrix |
Local contributions to the global matrix that correspond only to the variables listed in local_dof_indices
Definition at line 430 of file interface.h.
std::vector<FullMatrix<double> > aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_matrices_int_ext |
Local contributions to the global matrix from the face terms in the discontinuous Galerkin method. These arrays are of a length sufficient to hold one matrix for each possible face or subface of the cell. The discontinuous Galerkin bilinear form contains terms arising from internal (to the cell) values and external (to the cell) values. _int_ext
and _ext_int
hold the terms arising from the pairing between a cell and its neighbor, while _ext_ext
is the pairing of the neighbor's dofs with themselves. In the continuous Galerkin case, these are unused, and set to size zero.
Definition at line 443 of file interface.h.
std::vector<FullMatrix<double> > aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_matrices_ext_int |
Definition at line 444 of file interface.h.
std::vector<FullMatrix<double> > aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_matrices_ext_ext |
Definition at line 445 of file interface.h.
Vector<double> aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_rhs |
Local contributions to the right hand side that correspond only to the variables listed in local_dof_indices
Definition at line 451 of file interface.h.
std::vector<bool> aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::assembled_matrices |
Denotes which face matrices have actually been assembled in the DG field assembly. Entries for matrices not used (for example, those corresponding to non-existent subfaces; or faces being assembled by the neighboring cell) are set to false.
Definition at line 459 of file interface.h.
std::vector<types::global_dof_index> aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::local_dof_indices |
Indices of those degrees of freedom that actually correspond to the temperature or compositional field. since this structure is used to represent just contributions to the advection systems, there will be no contributions to other parts of the system and consequently, we do not need to list here indices that correspond to velocity or pressure degrees (or, in fact any other variable outside the block we are currently considering)
Definition at line 470 of file interface.h.
std::vector<std::vector<types::global_dof_index> > aspect::internal::Assembly::CopyData::AdvectionSystem< dim >::neighbor_dof_indices |
Indices of the degrees of freedom corresponding to the temperature or composition field on all possible neighboring cells. This is used in the discontinuous Galerkin method. The outer array has a length sufficient to hold one element for each possible face and sub-face of the current cell. The object is not used and has size zero if in the continuous Galerkin case.
Definition at line 480 of file interface.h.