Peano

Searching...

No Matches

Introduce multiplicative CC Coupling. More...

Inheritance diagram for solvers.api.actionsets.DGCGCoupling.MultiplicativeDGCGCoupling:

[legend]

Collaboration diagram for solvers.api.actionsets.DGCGCoupling.MultiplicativeDGCGCoupling:

[legend]

Public Member Functions
	__init__ (self, dg_solver, cg_solver, prolongation_matrix, prolongation_matrix_scaling, restriction_matrix, restriction_matrix_scaling, injection_matrix, injection_matrix_scaling, use_fas, smoothing_steps_DG=4, smoothing_steps_CG=-1, vcycles=1)
	Construct the solver.

	get_attributes (self)
	Add a new static attribute to the class.

	get_static_initialisations (self, full_qualified_classname)
	Initialise the new counter.

	get_body_of_prepareTraversal (self)

Public Member Functions inherited from solvers.api.actionsets.DGCGCoupling.AbstractDGCGCoupling
	get_body_of_operation (self, operation_name)
	Return actual C++ code snippets to be inserted into C++ code.

	get_action_set_name (self)
	Configure name of generated C++ action set.

	user_should_modify_template (self)
	The action set that Peano will generate that corresponds to this class should not be modified by users and can safely be overwritten every time we run the Python toolkit.

	get_includes (self)
	We need the solver repository in this action set, as we directly access the solver object.

Public Member Functions inherited from peano4.solversteps.ActionSet.ActionSet
	get_constructor_body (self)
	Define a tailored constructor body.

	get_destructor_body (self)

	get_body_of_getGridControlEvents (self)

	get_body_of_unprepareTraversal (self)

Data Fields
	smoothing_steps_DG

	smoothing_steps_CG

	vcycles

	cg_solver

	d

Data Fields inherited from solvers.api.actionsets.DGCGCoupling.AbstractDGCGCoupling
	d

Data Fields inherited from peano4.solversteps.ActionSet.ActionSet
	descend_invocation_order

	parallel

Additional Inherited Members
Static Public Attributes inherited from solvers.api.actionsets.DGCGCoupling.AbstractDGCGCoupling
str	templateTouchCellFirstTime

Static Public Attributes inherited from peano4.solversteps.ActionSet.ActionSet
str	OPERATION_BEGIN_TRAVERSAL = "beginTraversal"

str	OPERATION_END_TRAVERSAL = "endTraversal"

str	OPERATION_CREATE_PERSISTENT_VERTEX = "createPersistentVertex"

str	OPERATION_DESTROY_PERSISTENT_VERTEX = "destroyPersistentVertex"

str	OPERATION_CREATE_HANGING_VERTEX = "createHangingVertex"

str	OPERATION_DESTROY_HANGING_VERTEX = "destroyHangingVertex"

str	OPERATION_CREATE_PERSISTENT_FACE = "createPersistentFace"

str	OPERATION_DESTROY_PERSISTENT_FACE = "destroyPersistentFace"

str	OPERATION_CREATE_HANGING_FACE = "createHangingFace"

str	OPERATION_DESTROY_HANGING_FACE = "destroyHangingFace"

str	OPERATION_CREATE_CELL = "createCell"

str	OPERATION_DESTROY_CELL = "destroyCell"

str	OPERATION_TOUCH_VERTEX_FIRST_TIME = "touchVertexFirstTime"

str	OPERATION_TOUCH_VERTEX_LAST_TIME = "touchVertexLastTime"

str	OPERATION_TOUCH_FACE_FIRST_TIME = "touchFaceFirstTime"

str	OPERATION_TOUCH_FACE_LAST_TIME = "touchFaceLastTime"

str	OPERATION_TOUCH_CELL_FIRST_TIME = "touchCellFirstTime"

str	OPERATION_TOUCH_CELL_LAST_TIME = "touchCellLastTime"

Protected Attributes inherited from solvers.api.actionsets.DGCGCoupling.AbstractDGCGCoupling
	_use_fas

Detailed Description

Introduce multiplicative CC Coupling.

The control logic here is simplistic and prescripted. We follow the following steps for the additive solver with one iteration:

Pre First smoothing step
- When we touchVertexFirstTime, we pass over rollOverAndPrepareCGVertex, which has the effect of taking whatever is in the CG solution and writing it into the Delta. This is then interpolated during touchCellFirstTime.
- Interpolate correction, i.e. sum all the hierarchies up. This happens in touchCellFirstTime().
- Suspend CG solver
- Suspend DG solver
- Let DG solver project data onto faces
DG smoothing step
- Suspend CG solver
- Run DG solver
Nth DG smoothing step
- Still run DG solver (finish smoothing) During touchCellLastTime, we:
- Compute hierarchical residual
- Restrict
- Suspend CG solver These steps will set newRestrictedRhs correctly
Any CG solver On the first CG smoothing step, we should, during touchVertexFirstTime, pass over rollOverAndPrepareCGVertex which takes newRestrictedRhs and writes it into rhs, and we reset value to 0.
- Suspend DG solver

Realisation

This action set will map onto a class of its own. We give the class a class attribute which we increment by one in each step. We then make the actual logic depend upon this counter, as we inject the predicates into the superclass: That is, the superclass has some boolean expressions which allow us to switch its features on and off. We make those guys depends upon the new counter.

Definition at line 636 of file DGCGCoupling.py.

Constructor & Destructor Documentation

◆ init()

solvers.api.actionsets.DGCGCoupling.MultiplicativeDGCGCoupling.__init__	(	self,
		dg_solver,
		cg_solver,
		prolongation_matrix,
		prolongation_matrix_scaling,
		restriction_matrix,
		restriction_matrix_scaling,
		injection_matrix,
		injection_matrix_scaling,
		use_fas,
		smoothing_steps_DG = 4,
		smoothing_steps_CG = -1,
		vcycles = 1 )

Construct the solver.

Parameters

smoothing_steps_DG	positive integer
smoothing_steps_CG	positive integer or -1 The -1 indicates that we wait for full convergence on the coarse mesh.
vcycles	positive integer or -1. Terminate after the given number of full V-cycles. The -1 indicates that we wait for full convergence on the coarse mesh. If a positive integer is set, the correction solver tolerance needs to be set to a small value to ensure the given number of cycles is completed before the tolerance is reached.

Reimplemented from solvers.api.actionsets.DGCGCoupling.AbstractDGCGCoupling.

Definition at line 681 of file DGCGCoupling.py.