A Richards model specific Newton solver. More...

#include <dumux/porousmediumflow/richards/newtonsolver.hh>

Inheritance diagram for Dumux::RichardsNewtonSolver< Assembler, LinearSolver >:

Description

template<class Assembler, class LinearSolver>
class Dumux::RichardsNewtonSolver< Assembler, LinearSolver >

A Richards model specific Newton solver.

This solver 'knows' what a 'physically meaningful' solution is and can thus do update smarter than the plain Newton solver.

Public Types
using	Communication = Comm

Public Member Functions
const Communication &	comm () const
	the communicator for parallel runs More...

void	setMaxRelativeShift (Scalar tolerance)
	Set the maximum acceptable difference of any primary variable between two iterations for declaring convergence. More...

void	setMaxAbsoluteResidual (Scalar tolerance)
	Set the maximum acceptable absolute residual for declaring convergence. More...

void	setResidualReduction (Scalar tolerance)
	Set the maximum acceptable residual norm reduction. More...

void	setTargetSteps (int targetSteps)
	Set the number of iterations at which the Newton method should aim at. More...

void	setMinSteps (int minSteps)
	Set the number of minimum iterations for the Newton method. More...

void	setMaxSteps (int maxSteps)
	Set the number of iterations after which the Newton method gives up. More...

void	solve (SolutionVector &uCurrentIter, TimeLoop &timeLoop) override
	Run the Newton method to solve a non-linear system. Does time step control when the Newton fails to converge. More...

void	solve (SolutionVector &uCurrentIter) override
	Run the Newton method to solve a non-linear system. The solver is responsible for all the strategic decisions. More...

virtual void	newtonBegin (SolutionVector &u)
	Called before the Newton method is applied to an non-linear system of equations. More...

virtual bool	newtonProceed (const SolutionVector &uCurrentIter, bool converged)
	Returns true if another iteration should be done. More...

virtual void	newtonBeginStep (const SolutionVector &u)
	Indicates the beginning of a Newton iteration. More...

virtual void	assembleLinearSystem (const SolutionVector &uCurrentIter)
	Assemble the linear system of equations \(\mathbf{A}x - b = 0\). More...

void	solveLinearSystem (SolutionVector &deltaU)
	Solve the linear system of equations \(\mathbf{A}x - b = 0\). More...

void	newtonUpdate (SolutionVector &uCurrentIter, const SolutionVector &uLastIter, const SolutionVector &deltaU)
	Update the current solution with a delta vector. More...

virtual void	newtonEndStep (SolutionVector &uCurrentIter, const SolutionVector &uLastIter)
	Indicates that one Newton iteration was finished. More...

virtual void	newtonEnd ()
	Called if the Newton method ended (not known yet if we failed or succeeded) More...

virtual bool	newtonConverged () const
	Returns true if the error of the solution is below the tolerance. More...

virtual void	newtonFail (SolutionVector &u)
	Called if the Newton method broke down. This method is called after newtonEnd() More...

virtual void	newtonSucceed ()
	Called if the Newton method ended successfully This method is called after newtonEnd() More...

void	report (std::ostream &sout=std::cout) const
	output statistics / report More...

void	resetReport ()
	reset the statistics More...

void	reportParams (std::ostream &sout=std::cout) const
	Report the options and parameters this Newton is configured with. More...

Scalar	suggestTimeStepSize (Scalar oldTimeStep) const
	Suggest a new time-step size based on the old time-step size. More...

void	setVerbosity (int val)
	Specifies the verbosity level. More...

int	verbosity () const
	Return the verbosity level. More...

const std::string &	paramGroup () const
	Returns the parameter group. More...

void	attachConvergenceWriter (std::shared_ptr< ConvergenceWriter > convWriter)
	Attach a convergence writer to write out intermediate results after each iteration. More...

void	detachConvergenceWriter ()
	Detach the convergence writer to stop the output. More...

Protected Member Functions
virtual void	solutionChanged_ (const SolutionVector &uCurrentIter)
	Update solution-depended quantities like grid variables after the solution has changed. More...

void	computeResidualReduction_ (const SolutionVector &uCurrentIter)

bool	enableResidualCriterion () const

void	initPriVarSwitch_ (SolutionVector &, std::false_type)
	Initialize the privar switch, noop if there is no priVarSwitch. More...

void	initPriVarSwitch_ (SolutionVector &sol, std::true_type)
	Initialize the privar switch. More...

void	invokePriVarSwitch_ (SolutionVector &, std::false_type)
	Switch primary variables if necessary, noop if there is no priVarSwitch. More...

void	invokePriVarSwitch_ (SolutionVector &uCurrentIter, std::true_type)
	Switch primary variables if necessary. More...

const Assembler &	assembler () const
	Access the assembler. More...

Assembler &	assembler ()
	Access the assembler. More...

const LinearSolver &	linearSolver () const
	Access the linear solver. More...

LinearSolver &	linearSolver ()
	Access the linear solver. More...

template<class FirstRow , class ... Args>
bool	checkSizesOfSubMatrices (const Dune::MultiTypeBlockMatrix< FirstRow, Args... > &matrix) const
	Helper function to assure the MultiTypeBlockMatrix's sub-blocks have the correct sizes. More...

Protected Attributes
int	targetSteps_
	optimal number of iterations we want to achieve More...

int	minSteps_
	minimum number of iterations we do More...

int	maxSteps_
	maximum number of iterations we do before giving up More...

int	numSteps_
	actual number of steps done so far More...

Scalar	reduction_

Scalar	residualNorm_

Scalar	lastReduction_

Scalar	initialResidual_

Scalar	shift_

Scalar	lastShift_

std::ostringstream	endIterMsgStream_
	message stream to be displayed at the end of iterations More...

Member Typedef Documentation

◆ Communication

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

using Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::Communication = Comm

inherited

Member Function Documentation

◆ assembleLinearSystem()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::assembleLinearSystem ( const SolutionVector & uCurrentIter )

inlinevirtualinherited

Assemble the linear system of equations \(\mathbf{A}x - b = 0\).

Parameters

uCurrentIter The current iteration's solution vector

◆ assembler() [1/2]

template<class Assembler , class LinearSolver >

Assembler & Dumux::PDESolver< Assembler, LinearSolver >::assembler ( )

inlineprotectedinherited

Access the assembler.

◆ assembler() [2/2]

template<class Assembler , class LinearSolver >

const Assembler & Dumux::PDESolver< Assembler, LinearSolver >::assembler ( ) const

inlineprotectedinherited

Access the assembler.

◆ attachConvergenceWriter()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::attachConvergenceWriter ( std::shared_ptr< ConvergenceWriter > convWriter )

inlineinherited

Attach a convergence writer to write out intermediate results after each iteration.

◆ checkSizesOfSubMatrices()

template<class Assembler , class LinearSolver >

template<class FirstRow , class ... Args>

bool Dumux::PDESolver< Assembler, LinearSolver >::checkSizesOfSubMatrices ( const Dune::MultiTypeBlockMatrix< FirstRow, Args... > & matrix ) const

inlineprotectedinherited

Helper function to assure the MultiTypeBlockMatrix's sub-blocks have the correct sizes.

◆ comm()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

const Communication & Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::comm ( ) const

inlineinherited

the communicator for parallel runs

◆ computeResidualReduction_()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::computeResidualReduction_ ( const SolutionVector & uCurrentIter )

inlineprotectedinherited

◆ detachConvergenceWriter()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::detachConvergenceWriter ( )

inlineinherited

Detach the convergence writer to stop the output.

◆ enableResidualCriterion()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

bool Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::enableResidualCriterion ( ) const

inlineprotectedinherited

◆ initPriVarSwitch_() [1/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::initPriVarSwitch_	(	SolutionVector &	,
		std::false_type
	)

inlineprotectedinherited

Initialize the privar switch, noop if there is no priVarSwitch.

◆ initPriVarSwitch_() [2/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::initPriVarSwitch_	(	SolutionVector &	sol,
		std::true_type
	)

inlineprotectedinherited

Initialize the privar switch.

◆ invokePriVarSwitch_() [1/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::invokePriVarSwitch_	(	SolutionVector &	,
		std::false_type
	)

inlineprotectedinherited

Switch primary variables if necessary, noop if there is no priVarSwitch.

◆ invokePriVarSwitch_() [2/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::invokePriVarSwitch_	(	SolutionVector &	uCurrentIter,
		std::true_type
	)

inlineprotectedinherited

Switch primary variables if necessary.

◆ linearSolver() [1/2]

template<class Assembler , class LinearSolver >

LinearSolver & Dumux::PDESolver< Assembler, LinearSolver >::linearSolver ( )

inlineprotectedinherited

Access the linear solver.

◆ linearSolver() [2/2]

template<class Assembler , class LinearSolver >

const LinearSolver & Dumux::PDESolver< Assembler, LinearSolver >::linearSolver ( ) const

inlineprotectedinherited

Access the linear solver.

◆ newtonBegin()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonBegin ( SolutionVector & u )

inlinevirtualinherited

Called before the Newton method is applied to an non-linear system of equations.

Parameters

u	The initial solution

Reimplemented in Dumux::MultiDomainNewtonSolver< Assembler, LinearSolver, CouplingManager, Reassembler, Comm >.

◆ newtonBeginStep()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonBeginStep ( const SolutionVector & u )

inlinevirtualinherited

Indicates the beginning of a Newton iteration.

Reimplemented in Dumux::MultiDomainNewtonSolver< Assembler, LinearSolver, CouplingManager, Reassembler, Comm >.

◆ newtonConverged()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual bool Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonConverged ( ) const

inlinevirtualinherited

Returns true if the error of the solution is below the tolerance.

Reimplemented in Dumux::MultiDomainNewtonSolver< Assembler, LinearSolver, CouplingManager, Reassembler, Comm >.

◆ newtonEnd()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonEnd ( )

inlinevirtualinherited

Called if the Newton method ended (not known yet if we failed or succeeded)

◆ newtonEndStep()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonEndStep	(	SolutionVector &	uCurrentIter,
		const SolutionVector &	uLastIter
	)

inlinevirtualinherited

Indicates that one Newton iteration was finished.

Parameters

uCurrentIter	The solution after the current Newton iteration
uLastIter	The solution at the beginning of the current Newton iteration

Reimplemented in Dumux::NonEquilibriumNewtonSolver< Assembler, LinearSolver >, and Dumux::MultiDomainNewtonSolver< Assembler, LinearSolver, CouplingManager, Reassembler, Comm >.

◆ newtonFail()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonFail ( SolutionVector & u )

inlinevirtualinherited

Called if the Newton method broke down. This method is called after newtonEnd()

◆ newtonProceed()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual bool Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonProceed	(	const SolutionVector &	uCurrentIter,
		bool	converged
	)

inlinevirtualinherited

Returns true if another iteration should be done.

Parameters

uCurrentIter	The solution of the current Newton iteration
converged	if the Newton method's convergence criterion was met in this step

◆ newtonSucceed()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonSucceed ( )

inlinevirtualinherited

Called if the Newton method ended successfully This method is called after newtonEnd()

◆ newtonUpdate()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonUpdate	(	SolutionVector &	uCurrentIter,
		const SolutionVector &	uLastIter,
		const SolutionVector &	deltaU
	)

inlineinherited

Update the current solution with a delta vector.

The error estimates required for the newtonConverged() and newtonProceed() methods should be updated inside this method.

Different update strategies, such as line search and chopped updates can be implemented. The default behavior is just to subtract deltaU from uLastIter, i.e.

\[ u^{k+1} = u^k - \Delta u^k \]

Parameters

uCurrentIter	The solution vector after the current iteration
uLastIter	The solution vector after the last iteration
deltaU	The delta as calculated from solving the linear system of equations. This parameter also stores the updated solution.

◆ paramGroup()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

const std::string & Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::paramGroup ( ) const

inlineinherited

Returns the parameter group.

◆ report()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::report ( std::ostream & sout = std::cout ) const

inlineinherited

output statistics / report

◆ reportParams()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::reportParams ( std::ostream & sout = std::cout ) const

inlineinherited

Report the options and parameters this Newton is configured with.

◆ resetReport()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::resetReport ( )

inlineinherited

reset the statistics

◆ setMaxAbsoluteResidual()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setMaxAbsoluteResidual ( Scalar tolerance )

inlineinherited

Set the maximum acceptable absolute residual for declaring convergence.

Parameters

tolerance The maximum absolute residual at which the scheme is considered finished

◆ setMaxRelativeShift()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setMaxRelativeShift ( Scalar tolerance )

inlineinherited

Set the maximum acceptable difference of any primary variable between two iterations for declaring convergence.

Parameters

tolerance The maximum relative shift between two Newton iterations at which the scheme is considered finished

◆ setMaxSteps()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setMaxSteps ( int maxSteps )

inlineinherited

Set the number of iterations after which the Newton method gives up.

Parameters

maxSteps Number of iterations after we give up

◆ setMinSteps()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setMinSteps ( int minSteps )

inlineinherited

Set the number of minimum iterations for the Newton method.

Parameters

minSteps Minimum number of iterations

◆ setResidualReduction()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setResidualReduction ( Scalar tolerance )

inlineinherited

Set the maximum acceptable residual norm reduction.

Parameters

tolerance The maximum reduction of the residual norm at which the scheme is considered finished

◆ setTargetSteps()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setTargetSteps ( int targetSteps )

inlineinherited

Set the number of iterations at which the Newton method should aim at.

This is used to control the time-step size. The heuristic used is to scale the last time-step size by the deviation of the number of iterations used from the target steps.

Parameters

targetSteps Number of iterations which are considered "optimal"

◆ setVerbosity()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setVerbosity ( int val )

inlineinherited

Specifies the verbosity level.

◆ solutionChanged_()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::solutionChanged_ ( const SolutionVector & uCurrentIter )

inlineprotectedvirtualinherited

Update solution-depended quantities like grid variables after the solution has changed.

◆ solve() [1/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::solve ( SolutionVector & uCurrentIter )

inlineoverridevirtualinherited

Run the Newton method to solve a non-linear system. The solver is responsible for all the strategic decisions.

Implements Dumux::PDESolver< Assembler, LinearSolver >.

◆ solve() [2/2]

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::solve	(	SolutionVector &	uCurrentIter,
		TimeLoop &	timeLoop
	)

inlineoverridevirtualinherited

Run the Newton method to solve a non-linear system. Does time step control when the Newton fails to converge.

Reimplemented from Dumux::PDESolver< Assembler, LinearSolver >.

◆ solveLinearSystem()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::solveLinearSystem ( SolutionVector & deltaU )

inlineinherited

Solve the linear system of equations \(\mathbf{A}x - b = 0\).

Throws Dumux::NumericalProblem if the linear solver didn't converge.

If the linear solver doesn't accept multitype matrices we copy the matrix into a 1x1 block BCRS matrix for solving.

Parameters

deltaU The difference between the current and the next solution

◆ suggestTimeStepSize()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::suggestTimeStepSize ( Scalar oldTimeStep ) const

inlineinherited

Suggest a new time-step size based on the old time-step size.

The default behavior is to suggest the old time-step size scaled by the ratio between the target iterations and the iterations required to actually solve the last time-step.

◆ verbosity()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

int Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::verbosity ( ) const

inlineinherited

Return the verbosity level.

Member Data Documentation

◆ endIterMsgStream_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

std::ostringstream Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::endIterMsgStream_

protectedinherited

message stream to be displayed at the end of iterations

◆ initialResidual_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::initialResidual_

protectedinherited

◆ lastReduction_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::lastReduction_

protectedinherited

◆ lastShift_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::lastShift_

protectedinherited

◆ maxSteps_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

int Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::maxSteps_

protectedinherited

maximum number of iterations we do before giving up

◆ minSteps_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

int Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::minSteps_

protectedinherited

minimum number of iterations we do

◆ numSteps_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

int Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::numSteps_

protectedinherited

actual number of steps done so far

◆ reduction_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::reduction_

protectedinherited

◆ residualNorm_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::residualNorm_

protectedinherited

◆ shift_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

Scalar Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::shift_

protectedinherited

◆ targetSteps_

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::CollectiveCommunication<Dune::MPIHelper::MPICommunicator>>

int Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::targetSteps_

protectedinherited

optimal number of iterations we want to achieve

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

porousmediumflow/richards/newtonsolver.hh

Description

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Member Typedef Documentation

◆ Communication

Member Function Documentation

◆ assembleLinearSystem()

◆ assembler() [1/2]

◆ assembler() [2/2]

◆ attachConvergenceWriter()

◆ checkSizesOfSubMatrices()

◆ comm()

◆ computeResidualReduction_()

◆ detachConvergenceWriter()

◆ enableResidualCriterion()

◆ initPriVarSwitch_() [1/2]

◆ initPriVarSwitch_() [2/2]

◆ invokePriVarSwitch_() [1/2]

◆ invokePriVarSwitch_() [2/2]

◆ linearSolver() [1/2]

◆ linearSolver() [2/2]

◆ newtonBegin()

◆ newtonBeginStep()

◆ newtonConverged()

◆ newtonEnd()

◆ newtonEndStep()

◆ newtonFail()

◆ newtonProceed()

◆ newtonSucceed()

◆ newtonUpdate()

◆ paramGroup()

◆ report()

◆ reportParams()

◆ resetReport()

◆ setMaxAbsoluteResidual()

◆ setMaxRelativeShift()

◆ setMaxSteps()

◆ setMinSteps()

◆ setResidualReduction()

◆ setTargetSteps()

◆ setVerbosity()

◆ solutionChanged_()

◆ solve() [1/2]

◆ solve() [2/2]

◆ solveLinearSystem()

◆ suggestTimeStepSize()

◆ verbosity()

Member Data Documentation

◆ endIterMsgStream_

◆ initialResidual_

◆ lastReduction_

◆ lastShift_

◆ maxSteps_

◆ minSteps_

◆ numSteps_

◆ reduction_

◆ residualNorm_

◆ shift_

◆ targetSteps_