Inheritance diagram for aspect::StokesMatrixFreeHandlerImplementation< dim, velocity_degree >:

Public Member Functions
	StokesMatrixFreeHandlerImplementation (Simulator< dim > &, ParameterHandler &prm)

	~StokesMatrixFreeHandlerImplementation () override=default

std::pair< double, double >	solve () override

void	setup_dofs () override

void	assemble () override

void	build_preconditioner () override

const DoFHandler< dim > &	get_dof_handler_v () const override

const DoFHandler< dim > &	get_dof_handler_p () const override

const DoFHandler< dim > &	get_dof_handler_projection () const override

const AffineConstraints< double > &	get_constraints_v () const override

const AffineConstraints< double > &	get_constraints_p () const override

const MGTransferMF< dim, GMGNumberType > &	get_mg_transfer_A () const override

const MGTransferMF< dim, GMGNumberType > &	get_mg_transfer_S () const override

std::size_t	get_cell_data_memory_consumption () const override

Public Member Functions inherited from aspect::StokesMatrixFreeHandler< dim >
virtual	~StokesMatrixFreeHandler ()=default

Static Public Member Functions
static void	declare_parameters (ParameterHandler &prm)

Static Public Member Functions inherited from aspect::StokesMatrixFreeHandler< dim >
static void	declare_parameters (ParameterHandler &prm)

Private Types
using	StokesMatrixType = MatrixFreeStokesOperators::StokesOperator< dim, velocity_degree, double >

using	SchurComplementMatrixType = MatrixFreeStokesOperators::MassMatrixOperator< dim, velocity_degree-1, double >

using	ABlockMatrixType = MatrixFreeStokesOperators::ABlockOperator< dim, velocity_degree, double >

using	GMGSchurComplementMatrixType = MatrixFreeStokesOperators::MassMatrixOperator< dim, velocity_degree-1, GMGNumberType >

using	GMGABlockMatrixType = MatrixFreeStokesOperators::ABlockOperator< dim, velocity_degree, GMGNumberType >

Private Member Functions
void	parse_parameters (ParameterHandler &prm)

void	evaluate_material_model ()

void	correct_stokes_rhs ()

Private Attributes
Simulator< dim > &	sim

bool	print_details

bool	do_timings

double	minimum_viscosity

double	maximum_viscosity

DoFHandler< dim >	dof_handler_v

DoFHandler< dim >	dof_handler_p

DoFHandler< dim >	dof_handler_projection

FESystem< dim >	fe_v

FESystem< dim >	fe_p

FESystem< dim >	fe_projection

MatrixFreeStokesOperators::OperatorCellData< dim, GMGNumberType >	active_cell_data

MGLevelObject< MatrixFreeStokesOperators::OperatorCellData< dim, GMGNumberType > >	level_cell_data

StokesMatrixType	stokes_matrix

ABlockMatrixType	A_block_matrix

SchurComplementMatrixType	Schur_complement_block_matrix

AffineConstraints< double >	constraints_v

AffineConstraints< double >	constraints_p

MGLevelObject< GMGABlockMatrixType >	mg_matrices_A_block

MGLevelObject< GMGSchurComplementMatrixType >	mg_matrices_Schur_complement

MGConstrainedDoFs	mg_constrained_dofs_A_block

MGConstrainedDoFs	mg_constrained_dofs_Schur_complement

MGConstrainedDoFs	mg_constrained_dofs_projection

MGTransferMF< dim, GMGNumberType >	mg_transfer_A_block

MGTransferMF< dim, GMGNumberType >	mg_transfer_Schur_complement

std::vector< std::shared_ptr< MatrixFree< dim, double > > >	matrix_free_objects

Detailed Description

template<int dim, int velocity_degree>
class aspect::StokesMatrixFreeHandlerImplementation< dim, velocity_degree >

Main class of the Matrix-free method. Here are all the functions for setup, assembly and solving the Stokes system.

We need to derive from StokesMatrixFreeHandler to be able to introduce a second template argument for the degree of the Stokes finite element. This way, the main simulator does not need to know about the degree by using a pointer to the base class and we can pick the desired velocity degree at runtime.

Definition at line 90 of file simulator.h.