15#ifndef dealii_nonlinear_h
16#define dealii_nonlinear_h
79template <
typename VectorType = Vector<
double>>
249#ifdef DEAL_II_TRILINOS_WITH_NOX
254 const Teuchos::RCP<Teuchos::ParameterList> ¶meters =
255 Teuchos::rcp(
new Teuchos::ParameterList));
262#ifdef DEAL_II_WITH_SUNDIALS
272#ifdef DEAL_II_WITH_PETSC
284 solve(VectorType &initial_guess_and_solution);
307 std::function<void(
const VectorType &src, VectorType &dst)>
residual;
359 void(
const VectorType &rhs, VectorType &dst,
const double tolerance)>
377#ifdef DEAL_II_WITH_SUNDIALS
384#ifdef DEAL_II_TRILINOS_WITH_NOX
388 Teuchos::rcp(
new Teuchos::ParameterList);
394#ifdef DEAL_II_WITH_PETSC
410template <
typename VectorType>
415 (void)additional_data;
417#ifdef DEAL_II_WITH_SUNDIALS
423 else if (additional_data.
strategy ==
425 additional_data_kinsol.strategy =
427 else if (additional_data.
strategy ==
429 additional_data_kinsol.strategy =
434 additional_data_kinsol.maximum_non_linear_iterations =
436 additional_data_kinsol.function_tolerance =
438 additional_data_kinsol.step_tolerance = additional_data.
step_tolerance;
439 additional_data_kinsol.anderson_subspace_size =
444#ifdef DEAL_II_TRILINOS_WITH_NOX
445 parameters_nox->set(
"Nonlinear Solver",
"Line Search Based");
446 Teuchos::ParameterList &Line_Search = parameters_nox->sublist(
"Line Search");
447 Line_Search.set(
"Method",
"Full Step");
454#ifdef DEAL_II_WITH_PETSC
455 additional_data_petsc_snes.options_prefix =
"";
461 additional_data_petsc_snes.snes_type =
"anderson";
464 else if (additional_data.
strategy ==
467 additional_data_petsc_snes.snes_type =
"newtonls";
468 additional_data_petsc_snes.snes_linesearch_type =
"bt";
470 else if (additional_data.
strategy ==
473 additional_data_petsc_snes.snes_linesearch_type =
"newtonls";
474 additional_data_petsc_snes.snes_linesearch_type =
"basic";
476 else if (additional_data.
strategy ==
478 additional_data_petsc_snes.snes_type =
"nrichardson";
480 additional_data_petsc_snes.absolute_tolerance =
482 additional_data_petsc_snes.relative_tolerance =
484 additional_data_petsc_snes.step_tolerance = additional_data.
step_tolerance;
485 additional_data_petsc_snes.maximum_non_linear_iterations =
487 additional_data_petsc_snes.max_n_function_evaluations = -1;
494template <
typename VectorType>
496 : mpi_communicator(MPI_COMM_SELF)
501template <
typename VectorType>
504 : additional_data(additional_data)
505 , mpi_communicator(MPI_COMM_SELF)
512template <
typename VectorType>
516 : additional_data(additional_data)
517 , mpi_communicator(mpi_communicator)
524template <
typename VectorType>
529 additional_data.solver_type = type;
534template <
typename VectorType>
538 const unsigned int maximum_non_linear_iterations,
539 const double function_tolerance,
540 const double relative_tolerance,
541 const double step_tolerance,
542 const unsigned int anderson_subspace_size)
543 : solver_type(solver_type)
545 , maximum_non_linear_iterations(maximum_non_linear_iterations)
546 , function_tolerance(function_tolerance)
547 , relative_tolerance(relative_tolerance)
548 , step_tolerance(step_tolerance)
549 , anderson_subspace_size(anderson_subspace_size)
554#ifdef DEAL_II_TRILINOS_WITH_NOX
555template <
typename VectorType>
560 const Teuchos::RCP<Teuchos::ParameterList> ¶meters)
569#ifdef DEAL_II_WITH_SUNDIALS
570template <
typename VectorType>
581template <
typename VectorType>
587 ExcMessage(
"PETSc SNES requires you to use PETSc vectors."));
592#ifdef DEAL_II_WITH_PETSC
610 const double tolerance = 1e-6;
614 nonlinear_solver.
solve(initial_guess_and_solution);
620template <
typename VectorType>
623 VectorType &initial_guess_and_solution)
629#ifdef DEAL_II_WITH_PETSC
632#ifdef DEAL_II_TRILINOS_WITH_NOX
635#ifdef DEAL_II_WITH_SUNDIALS
646#ifdef DEAL_II_WITH_SUNDIALS
655 nonlinear_solver.
setup_jacobian = [&](
const VectorType ¤t_u,
656 const VectorType & ) {
662 nonlinear_solver.
solve(initial_guess_and_solution);
665 false,
ExcMessage(
"You do not have SUNDIALS configured with deal.II!"));
670#ifdef DEAL_II_TRILINOS_WITH_NOX
684 nonlinear_solver.
solve(initial_guess_and_solution);
687 false,
ExcMessage(
"You do not have Trilinos configured with deal.II"));
699 const std::string solvers =
""
700#ifdef DEAL_II_WITH_SUNDIALS
703#ifdef DEAL_II_TRILINOS_WITH_NOX
706#ifdef DEAL_II_WITH_PETSC
713 "Invalid nonlinear solver specified. "
714 "The solvers available in your installation are:\n" +
unsigned int anderson_subspace_size
SolutionStrategy strategy
AdditionalData(const SolverType &solver_type=automatic, const SolutionStrategy &strategy=linesearch, const unsigned int maximum_non_linear_iterations=200, const double function_tolerance=1e-8, const double relative_tolerance=1e-5, const double step_tolerance=0.0, const unsigned int anderson_subspace_size=0)
double function_tolerance
unsigned int maximum_non_linear_iterations
const double relative_tolerance
TrilinosWrappers::NOXSolver< VectorType >::AdditionalData additional_data_nox
void set_data(const AdditionalData &additional_data)
void solve_with_petsc(VectorType &initial_guess_and_solution)
void set_data(const typename PETScWrappers::NonlinearSolverData &additional_data)
SUNDIALS::KINSOL< VectorType >::AdditionalData additional_data_kinsol
std::function< void(const VectorType &src, VectorType &dst)> residual
std::function< void(VectorType &)> reinit_vector
AdditionalData additional_data
Teuchos::RCP< Teuchos::ParameterList > parameters_nox
PETScWrappers::NonlinearSolverData additional_data_petsc_snes
MPI_Comm mpi_communicator
void solve(VectorType &initial_guess_and_solution)
void select(const typename AdditionalData::SolverType &type)
NonlinearSolverSelector()
std::function< void(const VectorType ¤t_u)> setup_jacobian
std::function< void(const VectorType &rhs, VectorType &dst, const double tolerance)> solve_with_jacobian
std::function< void(const VectorType &src, VectorType &dst)> solve_with_jacobian
std::function< void(const VectorType &x, VectorType &res)> residual
unsigned int solve(VectorType &x)
std::function< void(const VectorType &x)> setup_jacobian
std::function< void(const VectorType ¤t_u, const VectorType ¤t_f)> setup_jacobian
unsigned int solve(VectorType &initial_guess_and_solution)
std::function< void(const VectorType &src, VectorType &dst)> residual
std::function< void(const VectorType &rhs, VectorType &dst, const double tolerance)> solve_with_jacobian
std::function< void(VectorType &)> reinit_vector
std::function< void(const VectorType &y, VectorType &x, const double tolerance)> solve_with_jacobian
std::function< void(const VectorType &u, VectorType &F)> residual
unsigned int solve(VectorType &solution)
std::function< void(const VectorType ¤t_u)> setup_jacobian
#define DEAL_II_DEPRECATED
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)