 |
version 3.10-dev
|
|
version 3.10-dev
|
A Richards model specific Newton solver. More...
#include <dumux/porousmediumflow/richards/newtonsolver.hh>
This solver 'knows' what a 'physically meaningful' solution is and can thus do update smarter than the plain Newton solver.
Public Types | |
| using | Communication = Dune::Communication< Dune::MPIHelper::MPICommunicator > |
| using | Assembler = Assembler |
| export the assembler and linear solver types More... | |
| using | LinearSolver = LinearSolver |
| using | Variables = Detail::PDESolver::AssemblerVariables< Assembler > |
| export the type of variables that represent a numerical solution More... | |
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 (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... | |
| bool | apply (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 (ResidualVector &deltaU) |
| Solve the linear system of equations \(\mathbf{A}x - b = 0\). More... | |
| void | newtonUpdate (Variables &vars, const SolutionVector &uLastIter, const ResidualVector &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 bool | newtonConverged () const |
| Returns true if the error of the solution is below the tolerance. More... | |
| virtual void | newtonFail (Variables &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) |
| Specify the verbosity level. More... | |
| int | verbosity () const |
| Return the verbosity level. More... | |
| void | setUseLineSearch (bool val=true) |
| Specify whether line search is enabled or not. More... | |
| bool | useLineSearch () const |
| Return whether line search is enabled or not. 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 | LinearAlgebraNativeBackend = VariablesBackend< ResidualVector > |
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... | |
| bool | checkSizesOfSubMatrices (const M &) const |
| Default implementation for any matrix type. 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... | |
|
inherited |
|
inherited |
|
protectedinherited |
|
inherited |
|
inherited |
|
inlineoverridevirtualinherited |
| vars | The variables object representing the current state of the numerical solution (primary and possibly secondary variables). |
Variables will represent the solution. If convergence fails, they are in some intermediate, undefined state. Implements Dumux::PDESolver< Assembler, LinearSolver >.
|
inlinevirtualinherited |
| vars | The current iteration's variables |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inlineinherited |
|
inlineprotectedinherited |
|
inlineinherited |
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inlineinherited |
|
inlinevirtualinherited |
| initVars | The variables representing the initial solution |
|
inlinevirtualinherited |
|
inlinevirtualinherited |
|
inlinevirtualinherited |
|
inlinevirtualinherited |
| vars | The variables after the current Newton iteration |
| uLastIter | The solution at the beginning of the current Newton iteration |
|
inlinevirtualinherited |
|
inlinevirtualinherited |
| varsCurrentIter | The variables of the current Newton iteration |
| converged | if the Newton method's convergence criterion was met in this step |
|
inlinevirtualinherited |
|
inlineinherited |
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 \]
| 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. |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
| tolerance | The maximum absolute residual at which the scheme is considered finished |
|
inlineinherited |
| tolerance | The maximum relative shift between two Newton iterations at which the scheme is considered finished |
|
inlineinherited |
| maxSteps | Number of iterations after we give up |
|
inlineinherited |
| minSteps | Minimum number of iterations |
|
inlineinherited |
| tolerance | The maximum reduction of the residual norm at which the scheme is considered finished |
|
inlineinherited |
|
inlineinherited |
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.
| targetSteps | Number of iterations which are considered "optimal" |
|
inlineinherited |
|
inlineinherited |
|
inlineprotectedvirtualinherited |
|
inlineoverridevirtualinherited |
| vars | The variables object representing the current state of the numerical solution (primary and possibly secondary variables). |
Implements Dumux::PDESolver< Assembler, LinearSolver >.
|
inlineoverridevirtualinherited |
| 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 >.
|
inlineinherited |
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.
| deltaU | The difference between the current and the next solution |
|
inlineinherited |
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.
|
inlineinherited |
|
inlineinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
3.10-dev
Generated by 
1.9.3