diff_solver.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2018 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 #include "libmesh/libmesh_common.h"
20 #include "libmesh/diff_solver.h"
21 #include "libmesh/newton_solver.h"
23 #include "libmesh/auto_ptr.h" // libmesh_make_unique
24 
25 namespace libMesh
26 {
27 
29  ParallelObject(s),
30  max_linear_iterations(1000),
31  max_nonlinear_iterations(100),
32  quiet(true),
33  verbose(false),
34  continue_after_max_iterations(true),
35  continue_after_backtrack_failure(false),
36  absolute_residual_tolerance(0.),
37  relative_residual_tolerance(0.),
38  absolute_step_tolerance(0.),
39  relative_step_tolerance(0.),
40  initial_linear_tolerance(1e-12),
41  minimum_linear_tolerance(TOLERANCE*TOLERANCE),
42  max_solution_norm(0.),
43  max_residual_norm(0.),
44  _outer_iterations(0),
45  _inner_iterations(0),
46  _system (s),
47  _solve_result(INVALID_SOLVE_RESULT)
48 {
49 }
50 
51 
52 
53 std::unique_ptr<DiffSolver> DiffSolver::build (sys_type & s)
54 {
55  return libmesh_make_unique<NewtonSolver>(s);
56 }
57 
58 
59 
61 {
62  // Reset the max_step_size and max_residual_norm for a new mesh
63  max_solution_norm = 0.;
64  max_residual_norm = 0.;
65 }
66 
67 
68 
70 {
71  // Reset the max_step_size and max_residual_norm for a new problem
72  max_solution_norm = 0.;
73  max_residual_norm = 0.;
74 }
75 
76 } // namespace libMesh
virtual void reinit()
Definition: diff_solver.C:60
static const Real TOLERANCE
static std::unique_ptr< DiffSolver > build(sys_type &s)
Definition: diff_solver.C:53
virtual void init()
Definition: diff_solver.C:69
An object whose state is distributed along a set of processors.
DiffSolver(sys_type &s)
Definition: diff_solver.C:28
Manages consistently variables, degrees of freedom, coefficient vectors, and matrices for implicit sy...