A two-phase PNM specific newton solver. More...

#include <dumux/porenetwork/2p/newtonsolver.hh>

Inheritance diagram for Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >:

Description

template<class Assembler, class LinearSolver, template< class, class > class NewtonConsistencyChecks = TwoPNewtonConsistencyChecks>
class Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >

A two-phase PNM specific newton solver.

This controller 'knows' what a 'physically meaningful' solution is which allows the newton method to abort quicker if the solution is way out of bounds.

Public Types
using	Communication = Comm

using	Assembler = Assembler
	export the assembler and linear solver types More...

using	LinearSolver = LinearSolver

using	Variables = Detail::AssemblerVariables< Assembler >
	export the type of variables that represent a numerical solution More...

Public Member Functions
void	newtonEndStep (SolutionVector &uCurrentIter, const SolutionVector &uLastIter) final
	Called after each Newton update. More...

bool	newtonConverged () const final
	Returns true if the current solution can be considered to be accurate enough. We enforce an additional Newton iteration if the invasion switch was triggered in the last iteration. More...

void	newtonFail (SolutionVector &u) final
	Called if the Newton method broke down. This method is called after newtonEnd() and resets the invasion state. More...

void	newtonSucceed () final
	Called if the Newton method ended successfully This method is called after newtonEnd() and advances the invasion state. More...

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 (Variables &vars, 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 (Variables &vars) override
	Run the Newton method to solve a non-linear system. The solver is responsible for all the strategic decisions. More...

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

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

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

virtual void	assembleLinearSystem (const Variables &vars)
	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 (Variables &vars, const SolutionVector &uLastIter, const SolutionVector &deltaU)
	Update the current solution with a delta vector. More...

virtual void	newtonEndStep (Variables &vars, 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 void	newtonFail (Variables &u)
	Called if the Newton method broke down. 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...

Scalar	retryTimeStepReductionFactor () const
	Return the factor for reducing the time step after a Newton iteration has failed. More...

void	setRetryTimeStepReductionFactor (const Scalar factor)
	Set the factor for reducing the time step after a Newton iteration has failed. More...

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

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

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

Protected Types
using	Backend = VariablesBackend< typename ParentType::Variables >

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

void	computeResidualReduction_ (const Variables &vars)

bool	enableResidualCriterion () const

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

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

◆ Assembler

using Dumux::PDESolver< Assembler , LinearSolver >::Assembler = Assembler

inherited

export the assembler and linear solver types

◆ Backend

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

using Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::Backend = VariablesBackend<typename ParentType::Variables>

protectedinherited

◆ Communication

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

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

inherited

◆ LinearSolver

using Dumux::PDESolver< Assembler , LinearSolver >::LinearSolver = LinearSolver

inherited

◆ Variables

using Dumux::PDESolver< Assembler , LinearSolver >::Variables = Detail::AssemblerVariables<Assembler>

inherited

export the type of variables that represent a numerical solution

Member Function Documentation

◆ assembleLinearSystem()

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

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

inlinevirtualinherited

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

Parameters

vars	The current iteration's variables

◆ assembler() [1/2]

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

inlineinherited

Access the assembler.

◆ assembler() [2/2]

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

inlineinherited

Access the assembler.

◆ attachConvergenceWriter()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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()

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::Communication<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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlineprotectedinherited

◆ detachConvergenceWriter()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlineprotectedinherited

◆ linearSolver() [1/2]

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

inlineprotectedinherited

Access the linear solver.

◆ linearSolver() [2/2]

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

inlineinherited

Access the linear solver.

◆ newtonBegin()

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

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

inlinevirtualinherited

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

Parameters

initVars The variables representing 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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlinevirtualinherited

Indicates the beginning of a Newton iteration.

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

◆ newtonConverged()

template<class Assembler , class LinearSolver , template< class, class > class NewtonConsistencyChecks = TwoPNewtonConsistencyChecks>

bool Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >::newtonConverged ( ) const

inlinefinalvirtual

Returns true if the current solution can be considered to be accurate enough. We enforce an additional Newton iteration if the invasion switch was triggered in the last iteration.

Reimplemented from Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >.

◆ newtonEnd()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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() [1/2]

template<class Assembler , class LinearSolver , template< class, class > class NewtonConsistencyChecks = TwoPNewtonConsistencyChecks>

void Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >::newtonEndStep	(	SolutionVector &	uCurrentIter,
		const SolutionVector &	uLastIter
	)

inlinefinal

Called after each Newton update.

Parameters

uCurrentIter	The current global solution vector
uLastIter	The previous global solution vector

◆ newtonEndStep() [2/2]

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

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

inlinevirtualinherited

Indicates that one Newton iteration was finished.

Parameters

vars	The variables 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() [1/2]

template<class Assembler , class LinearSolver , template< class, class > class NewtonConsistencyChecks = TwoPNewtonConsistencyChecks>

void Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >::newtonFail ( SolutionVector & u )

inlinefinal

Called if the Newton method broke down. This method is called after newtonEnd() and resets the invasion state.

◆ newtonFail() [2/2]

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

virtual void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonFail ( Variables & 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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlinevirtualinherited

Returns true if another iteration should be done.

Parameters

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

◆ newtonSucceed()

template<class Assembler , class LinearSolver , template< class, class > class NewtonConsistencyChecks = TwoPNewtonConsistencyChecks>

void Dumux::PoreNetwork::TwoPNewtonSolver< Assembler, LinearSolver, NewtonConsistencyChecks >::newtonSucceed ( )

inlinefinalvirtual

Called if the Newton method ended successfully This method is called after newtonEnd() and advances the invasion state.

Reimplemented from Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >.

◆ newtonUpdate()

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

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::newtonUpdate	(	Variables &	vars,
		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

vars	The variables 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::Communication<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::Communication<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::Communication<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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlineinherited

reset the statistics

◆ retryTimeStepReductionFactor()

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

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

inlineinherited

Return the factor for reducing the time step after a Newton iteration has failed.

◆ setMaxAbsoluteResidual()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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

◆ setRetryTimeStepReductionFactor()

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

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::setRetryTimeStepReductionFactor ( const Scalar factor )

inlineinherited

Set the factor for reducing the time step after a Newton iteration has failed.

◆ setTargetSteps()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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::Communication<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::Communication<Dune::MPIHelper::MPICommunicator>>

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

inlineprotectedvirtualinherited

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

Todo:: TODO: In case we stop support for old-style grid variables / assemblers at one point, this would become obsolete as only the update call to the backend would remain.

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

◆ solve() [1/2]

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

void Dumux::NewtonSolver< Assembler, LinearSolver, Reassembler, Comm >::solve ( Variables & vars )

inlineoverridevirtualinherited

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

Parameters

vars	The variables object representing the current state of the numerical solution (primary and possibly secondary variables).

Implements Dumux::PDESolver< Assembler, LinearSolver >.

◆ solve() [2/2]

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

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

inlineoverridevirtualinherited

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

Parameters

vars	The variables object representing the current state of the numerical solution (primary and possibly secondary variables).
timeLoop	The time loop.

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

◆ solveLinearSystem()

template<class Assembler , class LinearSolver , class Reassembler = PartialReassembler<Assembler>, class Comm = Dune::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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::Communication<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:

porenetwork/2p/newtonsolver.hh

Description

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Member Typedef Documentation

◆ Assembler

◆ Backend

◆ Communication

◆ LinearSolver

◆ Variables

Member Function Documentation

◆ assembleLinearSystem()

◆ assembler() [1/2]

◆ assembler() [2/2]

◆ attachConvergenceWriter()

◆ checkSizesOfSubMatrices()

◆ comm()

◆ computeResidualReduction_()

◆ detachConvergenceWriter()

◆ enableResidualCriterion()

◆ linearSolver() [1/2]

◆ linearSolver() [2/2]

◆ newtonBegin()

◆ newtonBeginStep()

◆ newtonConverged()

◆ newtonEnd()

◆ newtonEndStep() [1/2]

◆ newtonEndStep() [2/2]

◆ newtonFail() [1/2]

◆ newtonFail() [2/2]

◆ newtonProceed()

◆ newtonSucceed()

◆ newtonUpdate()

◆ paramGroup()

◆ report()

◆ reportParams()

◆ resetReport()

◆ retryTimeStepReductionFactor()

◆ setMaxAbsoluteResidual()

◆ setMaxRelativeShift()

◆ setMaxSteps()

◆ setMinSteps()

◆ setResidualReduction()

◆ setRetryTimeStepReductionFactor()

◆ 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_