#include <deal.II/numerics/nonlinear.h>

Classes
class	AdditionalData

Public Member Functions
	NonlinearSolverSelector ()

	NonlinearSolverSelector (const AdditionalData &additional_data)

	NonlinearSolverSelector (const AdditionalData &additional_data, const MPI_Comm &mpi_communicator)

void	select (const typename AdditionalData::SolverType &type)

void	set_data (const AdditionalData &additional_data)

void	set_data (const typename TrilinosWrappers::NOXSolver< VectorType >::AdditionalData &additional_data, const Teuchos::RCP< Teuchos::ParameterList > &parameters=Teuchos::rcp(new Teuchos::ParameterList))

void	set_data (const typename SUNDIALS::KINSOL< VectorType >::AdditionalData &additional_data)

void	set_data (const typename PETScWrappers::NonlinearSolverData &additional_data)

void	solve (VectorType &initial_guess_and_solution)

Public Attributes
std::function< void(VectorType &)>	reinit_vector

std::function< void(const VectorType &src, VectorType &dst)>	residual

std::function< void(const VectorType &current_u)>	setup_jacobian

std::function< void(const VectorType &rhs, VectorType &dst, const double tolerance)>	solve_with_jacobian

Private Member Functions
void	solve_with_petsc (VectorType &initial_guess_and_solution)

void	solve_with_petsc (PETScWrappers::MPI::Vector &initial_guess_and_solution)

Private Attributes
AdditionalData	additional_data

MPI_Comm	mpi_communicator

SUNDIALS::KINSOL< VectorType >::AdditionalData	additional_data_kinsol

TrilinosWrappers::NOXSolver< VectorType >::AdditionalData	additional_data_nox

Teuchos::RCP< Teuchos::ParameterList >	parameters_nox

PETScWrappers::NonlinearSolverData	additional_data_petsc_snes

Detailed Description

template<typename VectorType = Vector<double>>
class NonlinearSolverSelector< VectorType >

This class offers a unified interface to several nonlinear solver implementations like KINSOL, NOX, and SNES.

KINSOL nonlinear solvers are part of the SUNDIALS package, while NOX is part of Trilinos and SNES is part of PETSc, respectively. If no solver is manually specified, this class will automaticlaly choose one of the available solvers depending on the enabled dependencies.

By calling the solve function of this NonlinearSolverSelector, it selects the solve function of that Solver that was specified in the constructor of this class, similar to the SolverSelector class.

Usage

An example of code one would run with this class is the following:

// Generate a @p NonlinearSolverSelector that uses @p KINSOL
NonlinearSolverSelector<Vector<double>>::AdditionalData additional_data;
additional_data.solver_type =
 NonlinearSolverSelector<Vector<double>>::AdditionalData::SolverType::kinsol;
 
NonlinearSolverSelector<Vector<double>> nonlinear_solver(additional_data);
 
// Functions that are required for solving a nonlinear problem,
// which are utilized in both @p KINSOL and @p NOX.
nonlinear_solver.reinit_vector = [&](Vector<double> &x) {...};
 
nonlinear_solver.residual =
          [&](const Vector<double> &evaluation_point,
              Vector<double> &      residual) {...};
 
nonlinear_solver.setup_jacobian =
          [&](const Vector<double> &current_u) {...};
 
nonlinear_solver.solve_with_jacobian =
          [&](const Vector<double> &rhs,
              Vector<double> &      dst,
              const double tolerance) {...};
 
// Calling the @p solve function with an initial guess.
nonlinear_solver.solve(current_solution);

Definition at line 75 of file nonlinear.h.

Constructor & Destructor Documentation

◆ NonlinearSolverSelector() [1/3]

template<typename VectorType >

NonlinearSolverSelector< VectorType >::NonlinearSolverSelector ( )

Constructor, filling in default values

Definition at line 488 of file nonlinear.h.

◆ NonlinearSolverSelector() [2/3]

template<typename VectorType >

NonlinearSolverSelector< VectorType >::NonlinearSolverSelector ( const AdditionalData & additional_data )

Constructor, selecting the solver and other parameters specified in additional_data.

Deprecated:: Use the other constructor with MPI_Comm instead.

Definition at line 495 of file nonlinear.h.

◆ NonlinearSolverSelector() [3/3]

template<typename VectorType >

NonlinearSolverSelector< VectorType >::NonlinearSolverSelector	(	const AdditionalData &	additional_data,
		const MPI_Comm &	mpi_communicator
	)

Constructor.

Parameters

additional_data	NonlinearSolverSelector configuration data
mpi_communicator	MPI communicator used by the nonlinear solver.

Definition at line 506 of file nonlinear.h.

Member Function Documentation

◆ select()

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::select ( const typename AdditionalData::SolverType & type )

Select a new nonlinear solver. All solver names used in this class are all lower case.

Definition at line 519 of file nonlinear.h.

◆ set_data() [1/4]

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::set_data ( const AdditionalData & additional_data )

Set the generic additional data for all nonlinear solvers.

Definition at line 407 of file nonlinear.h.

◆ set_data() [2/4]

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::set_data	(	const typename TrilinosWrappers::NOXSolver< VectorType >::AdditionalData &	additional_data,
		const Teuchos::RCP< Teuchos::ParameterList > &	parameters = `Teuchos::rcp(new Teuchos::ParameterList)`
	)

Set the additional data for NOX. See TrilinosWrappers::NOXSolver for more information.

Definition at line 550 of file nonlinear.h.

◆ set_data() [3/4]

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::set_data ( const typename SUNDIALS::KINSOL< VectorType >::AdditionalData & additional_data )

Set the additional data for KINSOL. See SUNDIALS::KINSOL for more information.

Definition at line 565 of file nonlinear.h.

◆ set_data() [4/4]

template<typename VectorType = Vector<double>>

void NonlinearSolverSelector< VectorType >::set_data ( const typename PETScWrappers::NonlinearSolverData & additional_data )

Set the additional data for SNES. See PETScWrappers::NonlinearSolverData for more information.

◆ solve()

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::solve ( VectorType & initial_guess_and_solution )

Solve the nonlinear system.

The content of initial_guess_and_solution is used as an initial guess and the final solution is stored in the same vector.

Definition at line 615 of file nonlinear.h.

◆ solve_with_petsc() [1/2]

template<typename VectorType >

void NonlinearSolverSelector< VectorType >::solve_with_petsc ( VectorType & initial_guess_and_solution )

private

Solve with PETSc SNES. Internal functions specialized for PETSc Vectors. This is necessary to only instantiate the SNES solvers with PETSc Vectors, as this is the only vector type currently supported.

Definition at line 576 of file nonlinear.h.

◆ solve_with_petsc() [2/2]

void NonlinearSolverSelector< PETScWrappers::MPI::Vector >::solve_with_petsc ( PETScWrappers::MPI::Vector & initial_guess_and_solution )

private

Definition at line 588 of file nonlinear.h.

Member Data Documentation

◆ reinit_vector

template<typename VectorType = Vector<double>>

std::function<void(VectorType &)> NonlinearSolverSelector< VectorType >::reinit_vector

A function object that users need to supply and that is intended to reinitize the given vector to its correct size, and block structure (if block vectors are used), along with any other properties necessary.

Definition at line 287 of file nonlinear.h.

◆ residual

template<typename VectorType = Vector<double>>

std::function<void(const VectorType &src, VectorType &dst)> NonlinearSolverSelector< VectorType >::residual

A function object that users should supply and that is intended to compute the residual dst = F(src).

Note: This variable represents a user provided callback. See there for a description of how to deal with errors and other requirements and conventions. Some of the underlying packages used by this class can deal with recoverable exceptions, whereas others cannot. As a consequence, if a callback throws an exception of type RecoverableUserCallbackError, then this exception may or may not be treated like any other exception.

Definition at line 302 of file nonlinear.h.

◆ setup_jacobian

template<typename VectorType = Vector<double>>

std::function<void(const VectorType &current_u)> NonlinearSolverSelector< VectorType >::setup_jacobian

A function object that users may supply and that is intended to prepare the linear solver for subsequent calls to solve_with_jacobian).

The job of setup_jacobian() is to prepare the linear solver for subsequent calls to solve_with_jacobian(), in the solution of linear systems \(Ax = b\). The exact nature of this system depends on the SolutionStrategy that has been selected.

In the cases strategy = SolutionStrategy::newton or SolutionStrategy::linesearch, \(A\) is the Jacobian \(J = \partial F/\partial u\). If strategy = SolutionStrategy::picard, \(A\) is the approximate Jacobian matrix \(L\).

Parameters

current_u Current value of \(u\)

Note: This variable represents a user provided callback. See there for a description of how to deal with errors and other requirements and conventions. Some of the underlying packages used by this class can deal with recoverable exceptions, whereas others cannot. As a consequence, if a callback throws an exception of type RecoverableUserCallbackError, then this exception may or may not be treated like any other exception.

Definition at line 330 of file nonlinear.h.

◆ solve_with_jacobian

template<typename VectorType = Vector<double>>

std::function< void(const VectorType &rhs, VectorType &dst, const double tolerance)> NonlinearSolverSelector< VectorType >::solve_with_jacobian

A function object that users may supply and that is intended to solve a linear system with the Jacobian matrix.

Note: PETSc SNES does not provide us with a linear tolerance to solve linear system with. We will provide a value of 1e-6 in that case.

Parameters

[in]	rhs	The system right hand side to solve for.
[out]	dst	The solution of \(J^{-1} * \texttt{src}\).
[in]	tolerance	The tolerance with which to solve the linear system of equations.

Note: This variable represents a user provided callback. See there for a description of how to deal with errors and other requirements and conventions. Some of the underlying packages used by this class can deal with recoverable exceptions, whereas others cannot. As a consequence, if a callback throws an exception of type RecoverableUserCallbackError, then this exception may or may not be treated like any other exception.

Definition at line 355 of file nonlinear.h.

◆ additional_data

template<typename VectorType = Vector<double>>

AdditionalData NonlinearSolverSelector< VectorType >::additional_data

private

NonlinearSolverSelector configuration data.

Definition at line 361 of file nonlinear.h.

◆ mpi_communicator

template<typename VectorType = Vector<double>>

MPI_Comm NonlinearSolverSelector< VectorType >::mpi_communicator

private

The MPI communicator to be used by this solver, if any.

Definition at line 367 of file nonlinear.h.

◆ additional_data_kinsol

template<typename VectorType = Vector<double>>

SUNDIALS::KINSOL<VectorType>::AdditionalData NonlinearSolverSelector< VectorType >::additional_data_kinsol

private

KINSOL configuration data

Definition at line 373 of file nonlinear.h.

◆ additional_data_nox

template<typename VectorType = Vector<double>>

TrilinosWrappers::NOXSolver<VectorType>::AdditionalData NonlinearSolverSelector< VectorType >::additional_data_nox

private

NOX configuration data

Definition at line 381 of file nonlinear.h.

◆ parameters_nox

template<typename VectorType = Vector<double>>

Teuchos::RCP<Teuchos::ParameterList> NonlinearSolverSelector< VectorType >::parameters_nox

private

Initial value:

=

Teuchos::rcp(new Teuchos::ParameterList)

Definition at line 382 of file nonlinear.h.

◆ additional_data_petsc_snes

template<typename VectorType = Vector<double>>

PETScWrappers::NonlinearSolverData NonlinearSolverSelector< VectorType >::additional_data_petsc_snes

private

PETSc SNES configuration data

Definition at line 390 of file nonlinear.h.

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

include/deal.II/numerics/nonlinear.h

Classes

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Usage

Constructor & Destructor Documentation

◆ NonlinearSolverSelector() [1/3]

◆ NonlinearSolverSelector() [2/3]

◆ NonlinearSolverSelector() [3/3]

Member Function Documentation

◆ select()

◆ set_data() [1/4]

◆ set_data() [2/4]

◆ set_data() [3/4]

◆ set_data() [4/4]

◆ solve()

◆ solve_with_petsc() [1/2]

◆ solve_with_petsc() [2/2]

Member Data Documentation

◆ reinit_vector

◆ residual

◆ setup_jacobian

◆ solve_with_jacobian

◆ additional_data

◆ mpi_communicator

◆ additional_data_kinsol

◆ additional_data_nox

◆ parameters_nox

◆ additional_data_petsc_snes