Opm::TpsaNewtonMethod< TypeTag > Class Template Reference

Newton method solving for generic TPSA model. More...

#include <tpsanewtonmethod.hpp>

Public Member Functions
	TpsaNewtonMethod (Simulator &simulator)
	Constructor. More...
bool	apply ()
	Run the Newton method. More...
bool	converged () const
	Returns true if the error of the solution is below the tolerance. More...
Scalar	initialError () const
	Returns the error of the first linearization of the last apply(). More...
bool	initiallyConverged () const
	Returns true if the state handed to the last apply() already solved the system. More...
Problem &	problem ()
	Returns a reference to the object describing the current physical problem. More...
const Problem &	problem () const
	Returns a reference to the object describing the current physical problem. More...
Model &	model ()
	Returns a reference to the geomechanics model. More...
const Model &	model () const
	Returns a reference to the geomechanics model. More...
LinearSolverBackend &	linearSolver ()
	Returns the linear solver backend object for external use. More...
const LinearSolverBackend &	linearSolver () const
	Returns the linear solver backend object for external use. More...
int	numIterations () const
	Returns the number of iterations done since the Newton method was invoked. More...
int	numLinearizations () const
	Returns the number of linearizations that has done since the Newton method was invoked. More...
int	numTotLinearIterations () const
	Returns the number of linear solver iterations done since the Newton method was invoked. More...
Scalar	tolerance () const
	Return the current tolerance at which the Newton method considers itself to be converged. More...
Scalar	minIterations () const
	Returns minimum number of Newton iterations used. More...
const Timer &	prePostProcessTimer () const
	Return post-process timer. More...
const Timer &	linearizeTimer () const
	Return linearization timer. More...
const Timer &	solveTimer () const
	Return linear solver timer. More...
const Timer &	updateTimer () const
	Return solution update timer. More...

Static Public Member Functions
static void	registerParameters ()
	Register all run-time parameters for the Newton method. More...

Protected Member Functions
int	verbosity_ () const
	Verbosity level of Newton print messages. More...
bool	singleIteration_ () const
	Whether the Newton method is limited to a single iteration. More...
void	begin_ ()
	Called before the Newton method is applied to an non-linear system of equations. More...
void	beginIteration_ ()
	Calculations at the beginning of a Newton iteration. More...
void	linearizeDomain_ ()
	Linearize the global non-linear system of equations associated with the spatial domain. More...
void	preSolve_ (const SolutionVector &, const GlobalEqVector &currentResidual)
	Compute error before a Newton iteration. More...
void	update_ (SolutionVector &nextSolution, const SolutionVector &currentSolution, const GlobalEqVector &solutionUpdate, const GlobalEqVector &currentResidual)
	Update the current solution with a delta vector. More...
void	updatePrimaryVariables_ (unsigned, PrimaryVariables &nextValue, const PrimaryVariables &currentValue, const EqVector &update, const EqVector &)
	Update a single primary variables object. More...
void	endIteration_ ()
	Indicates that one Newton iteration was finished. More...
bool	proceed_ () const
	Returns true iff another Newton iteration should be done. More...
void	failed_ ()
	Called if the Newton method broke down. More...

Protected Attributes
Simulator &	simulator_
LinearSolverBackend	linearSolver_
Timer	prePostProcessTimer_
Timer	linearizeTimer_
Timer	solveTimer_
Timer	updateTimer_
Scalar	error_
Scalar	lastError_
Scalar	initialError_
TpsaNewtonMethodParams< Scalar >	params_
int	numIterations_
int	numLinearizations_
int	numTotLinearIterations_

Detailed Description

template<class TypeTag>
class Opm::TpsaNewtonMethod< TypeTag >

Newton method solving for generic TPSA model.

Generates the Jacobian matrix, J(u^n) and residual vector, R(u^n), with a solution vector, u^n, at iteration n. Subsequently the linear system J(u^n)\delta u^n = -R(u^n) is solved to get u^{n+1} = u^n + \Delta u^n.

By default only a single iteration is done to solve the linear elasticity system. Note that, in this case, only linear solver convergence is checked, and Newton error is only used to check if the linearization step produced a low enough residual to jump out before a linear solve.

Constructor & Destructor Documentation

TpsaNewtonMethod()

template<class TypeTag >

Opm::TpsaNewtonMethod< TypeTag >::TpsaNewtonMethod ( Simulator &  simulator )

inlineexplicit

Constructor.

Parameters

simulator

Simulator object

References Opm::TpsaNewtonMethod< TypeTag >::params_.

Member Function Documentation

apply()

template<class TypeTag >

bool Opm::TpsaNewtonMethod< TypeTag >::apply ( )

inline

Run the Newton method.

Returns

Bool indicating if Newton converged

References Opm::TpsaNewtonMethod< TypeTag >::begin_(), Opm::TpsaNewtonMethod< TypeTag >::beginIteration_(), Opm::TpsaNewtonMethod< TypeTag >::converged(), Opm::TpsaNewtonMethod< TypeTag >::endIteration_(), Opm::TpsaNewtonMethod< TypeTag >::error_, Opm::TpsaNewtonMethod< TypeTag >::failed_(), Opm::Timer::halt(), Opm::TpsaNewtonMethod< TypeTag >::linearizeDomain_(), Opm::TpsaNewtonMethod< TypeTag >::linearizeTimer_, Opm::TpsaNewtonMethod< TypeTag >::linearSolver_, Opm::TpsaNewtonMethod< TypeTag >::model(), Opm::TpsaNewtonMethod< TypeTag >::numIterations(), Opm::TpsaNewtonMethod< TypeTag >::numLinearizations(), Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations(), Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations_, Opm::TpsaNewtonMethod< TypeTag >::prePostProcessTimer_, Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), Opm::TpsaNewtonMethod< TypeTag >::proceed_(), Opm::Timer::realTimeElapsed(), Opm::TpsaNewtonMethod< TypeTag >::singleIteration_(), Opm::TpsaNewtonMethod< TypeTag >::solveTimer_, Opm::Timer::start(), Opm::Timer::stop(), Opm::TpsaNewtonMethod< TypeTag >::update_(), Opm::TpsaNewtonMethod< TypeTag >::updateTimer_, and Opm::TpsaNewtonMethod< TypeTag >::verbosity_().

begin_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::begin_ ( )

inlineprotected

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

References Opm::TpsaNewtonMethod< TypeTag >::error_, Opm::TpsaNewtonMethod< TypeTag >::initialError_, Opm::TpsaNewtonMethod< TypeTag >::numIterations_, Opm::TpsaNewtonMethod< TypeTag >::numLinearizations_, and Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

beginIteration_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::beginIteration_ ( )

inlineprotected

Calculations at the beginning of a Newton iteration.

References Opm::TpsaNewtonMethod< TypeTag >::error_, and Opm::TpsaNewtonMethod< TypeTag >::lastError_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

converged()

template<class TypeTag >

bool Opm::TpsaNewtonMethod< TypeTag >::converged ( ) const

inline

Returns true if the error of the solution is below the tolerance.

Returns

Bool indicating if convergence has been achived

Note

Decides whether another iteration is worth doing, but not whether apply() succeeded when only a single iteration is done, see singleIteration_()

References Opm::TpsaNewtonMethod< TypeTag >::error_, and Opm::TpsaNewtonMethod< TypeTag >::tolerance().

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::proceed_().

endIteration_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::endIteration_ ( )

inlineprotected

Indicates that one Newton iteration was finished.

References Opm::TpsaNewtonMethod< TypeTag >::error_, Opm::TpsaNewtonMethod< TypeTag >::numIterations_, Opm::to_string(), and Opm::TpsaNewtonMethod< TypeTag >::verbosity_().

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

failed_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::failed_ ( )

inlineprotected

Called if the Newton method broke down.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

initialError()

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::initialError ( ) const

inline

Returns the error of the first linearization of the last apply().

Returns

Weighted maximum norm of the residual the Newton method started from

References Opm::TpsaNewtonMethod< TypeTag >::initialError_.

initiallyConverged()

template<class TypeTag >

bool Opm::TpsaNewtonMethod< TypeTag >::initiallyConverged ( ) const

inline

Returns true if the state handed to the last apply() already solved the system.

Returns

Bool indicating if the first linearization was below the tolerance

References Opm::TpsaNewtonMethod< TypeTag >::initialError_, and Opm::TpsaNewtonMethod< TypeTag >::tolerance().

linearizeDomain_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::linearizeDomain_ ( )

inlineprotected

Linearize the global non-linear system of equations associated with the spatial domain.

References Opm::TpsaNewtonMethod< TypeTag >::model(), and Opm::TpsaNewtonMethod< TypeTag >::numLinearizations_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

linearizeTimer()

template<class TypeTag >

const Timer & Opm::TpsaNewtonMethod< TypeTag >::linearizeTimer ( ) const

inline

Return linearization timer.

Returns

Reference to linearization timer object

References Opm::TpsaNewtonMethod< TypeTag >::linearizeTimer_.

linearSolver() [1/2]

template<class TypeTag >

LinearSolverBackend & Opm::TpsaNewtonMethod< TypeTag >::linearSolver ( )

inline

Returns the linear solver backend object for external use.

Returns

Reference to linear solver object

References Opm::TpsaNewtonMethod< TypeTag >::linearSolver_.

linearSolver() [2/2]

template<class TypeTag >

const LinearSolverBackend & Opm::TpsaNewtonMethod< TypeTag >::linearSolver ( ) const

inline

Returns the linear solver backend object for external use.

Returns

Reference to linear solver object

References Opm::TpsaNewtonMethod< TypeTag >::linearSolver_.

minIterations()

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::minIterations ( ) const

inline

Returns minimum number of Newton iterations used.

Returns

Minimum number of Newton iterations

References Opm::TpsaNewtonMethod< TypeTag >::params_.

model() [1/2]

template<class TypeTag >

Model & Opm::TpsaNewtonMethod< TypeTag >::model ( )

inline

Returns a reference to the geomechanics model.

Returns

Reference to geomechanics model object

References Opm::TpsaNewtonMethod< TypeTag >::simulator_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), Opm::TpsaNewtonMethod< TypeTag >::linearizeDomain_(), Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), and Opm::TpsaNewtonMethod< TypeTag >::update_().

model() [2/2]

template<class TypeTag >

const Model & Opm::TpsaNewtonMethod< TypeTag >::model ( ) const

inline

Returns a reference to the geomechanics model.

Returns

Reference to geomechanics model object

References Opm::TpsaNewtonMethod< TypeTag >::simulator_.

numIterations()

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numIterations ( ) const

inline

Returns the number of iterations done since the Newton method was invoked.

Returns

Number of Newton iteratinos

References Opm::TpsaNewtonMethod< TypeTag >::numIterations_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::proceed_().

numLinearizations()

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numLinearizations ( ) const

inline

Returns the number of linearizations that has done since the Newton method was invoked.

Returns

Number of linearizations

References Opm::TpsaNewtonMethod< TypeTag >::numLinearizations_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

numTotLinearIterations()

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations ( ) const

inline

Returns the number of linear solver iterations done since the Newton method was invoked.

Returns

Number of linear iterations, summed over the Newton iterations

References Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

prePostProcessTimer()

template<class TypeTag >

const Timer & Opm::TpsaNewtonMethod< TypeTag >::prePostProcessTimer ( ) const

inline

Return post-process timer.

Returns

Reference to post-process timer object

References Opm::TpsaNewtonMethod< TypeTag >::prePostProcessTimer_.

preSolve_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::preSolve_ ( const SolutionVector &  , const GlobalEqVector &  currentResidual  )

inlineprotected

Compute error before a Newton iteration.

Parameters

currentSolution	Current solution vector
currentResidual	Current residual vector

References Opm::TpsaNewtonMethod< TypeTag >::error_, Opm::TpsaNewtonMethod< TypeTag >::initialError_, Opm::TpsaNewtonMethod< TypeTag >::lastError_, Opm::TpsaNewtonMethod< TypeTag >::model(), Opm::TpsaNewtonMethod< TypeTag >::numLinearizations_, Opm::TpsaNewtonMethod< TypeTag >::params_, Opm::TpsaNewtonMethod< TypeTag >::simulator_, and Opm::to_string().

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

problem() [1/2]

template<class TypeTag >

Problem & Opm::TpsaNewtonMethod< TypeTag >::problem ( )

inline

Returns a reference to the object describing the current physical problem.

Returns

Reference to problem object

References Opm::TpsaNewtonMethod< TypeTag >::simulator_.

problem() [2/2]

template<class TypeTag >

const Problem & Opm::TpsaNewtonMethod< TypeTag >::problem ( ) const

inline

Returns a reference to the object describing the current physical problem.

Returns

Reference to problem object

References Opm::TpsaNewtonMethod< TypeTag >::simulator_.

proceed_()

template<class TypeTag >

bool Opm::TpsaNewtonMethod< TypeTag >::proceed_ ( ) const

inlineprotected

Returns true iff another Newton iteration should be done.

Returns

Bool indicating if Newton iterations should continue

References Opm::TpsaNewtonMethod< TypeTag >::converged(), Opm::TpsaNewtonMethod< TypeTag >::error_, Opm::TpsaNewtonMethod< TypeTag >::lastError_, Opm::TpsaNewtonMethod< TypeTag >::numIterations(), Opm::TpsaNewtonMethod< TypeTag >::params_, and Opm::TpsaNewtonMethod< TypeTag >::singleIteration_().

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

registerParameters()

template<class TypeTag >

static void Opm::TpsaNewtonMethod< TypeTag >::registerParameters ( )

inlinestatic

Register all run-time parameters for the Newton method.

References Opm::TpsaNewtonMethodParams< Scalar >::registerParameters().

singleIteration_()

template<class TypeTag >

bool Opm::TpsaNewtonMethod< TypeTag >::singleIteration_ ( ) const

inlineprotected

Whether the Newton method is limited to a single iteration.

Returns

Bool indicating if at most one linearization and linear solve is done

If the elasticity system is linear, a single solve is exact and more than one Newton error is unnecessary.

References Opm::TpsaNewtonMethod< TypeTag >::params_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::proceed_().

solveTimer()

template<class TypeTag >

const Timer & Opm::TpsaNewtonMethod< TypeTag >::solveTimer ( ) const

inline

Return linear solver timer.

Returns

Reference to linear solver timer object

References Opm::TpsaNewtonMethod< TypeTag >::solveTimer_.

tolerance()

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::tolerance ( ) const

inline

Return the current tolerance at which the Newton method considers itself to be converged.

Returns

Tolerance for Newton error

References Opm::TpsaNewtonMethod< TypeTag >::params_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::converged(), and Opm::TpsaNewtonMethod< TypeTag >::initiallyConverged().

update_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::update_ ( SolutionVector &  nextSolution, const SolutionVector &  currentSolution, const GlobalEqVector &  solutionUpdate, const GlobalEqVector &  currentResidual  )

inlineprotected

Update the current solution with a delta vector.

Parameters

nextSolution	The solution vector after the current iteration
currentSolution	The solution vector after the last iteration
solutionUpdate	The delta vector as calculated by solving the linear system of equations
currentResidual	The residual vector of the current Newton-Raphson iteraton

Different update strategies, such as chopped updates can be implemented by overriding this method. The default behavior is use the standard Newton-Raphson update strategy, i.e. u^{n+1} = u^n - \Delta u^n

References Opm::TpsaNewtonMethod< TypeTag >::model(), and Opm::TpsaNewtonMethod< TypeTag >::updatePrimaryVariables_().

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply().

updatePrimaryVariables_()

template<class TypeTag >

void Opm::TpsaNewtonMethod< TypeTag >::updatePrimaryVariables_ ( unsigned  , PrimaryVariables &  nextValue, const PrimaryVariables &  currentValue, const EqVector &  update, const EqVector &    )

inlineprotected

Update a single primary variables object.

Parameters

nextValue	The solution vector after the current iteration
currentValue	The solution vector after the last iteration
update	The delta vector as calculated by solving the linear system of equations
currentResidual	The residual vector of the current Newton-Raphson iteraton

Referenced by Opm::TpsaNewtonMethod< TypeTag >::update_().

updateTimer()

template<class TypeTag >

const Timer & Opm::TpsaNewtonMethod< TypeTag >::updateTimer ( ) const

inline

Return solution update timer.

Returns

Reference to solution update timer object

References Opm::TpsaNewtonMethod< TypeTag >::updateTimer_.

verbosity_()

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::verbosity_ ( ) const

inlineprotected

Verbosity level of Newton print messages.

Returns

Level indicating verbosity

Newton procedure verbosity: 0=none, 1=basic, 2=all

References Opm::TpsaNewtonMethod< TypeTag >::params_, and Opm::TpsaNewtonMethod< TypeTag >::simulator_.

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::endIteration_().

Member Data Documentation

error_

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::error_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), Opm::TpsaNewtonMethod< TypeTag >::begin_(), Opm::TpsaNewtonMethod< TypeTag >::beginIteration_(), Opm::TpsaNewtonMethod< TypeTag >::converged(), Opm::TpsaNewtonMethod< TypeTag >::endIteration_(), Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), and Opm::TpsaNewtonMethod< TypeTag >::proceed_().

initialError_

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::initialError_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::begin_(), Opm::TpsaNewtonMethod< TypeTag >::initialError(), Opm::TpsaNewtonMethod< TypeTag >::initiallyConverged(), and Opm::TpsaNewtonMethod< TypeTag >::preSolve_().

lastError_

template<class TypeTag >

Scalar Opm::TpsaNewtonMethod< TypeTag >::lastError_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::beginIteration_(), Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), and Opm::TpsaNewtonMethod< TypeTag >::proceed_().

linearizeTimer_

template<class TypeTag >

Timer Opm::TpsaNewtonMethod< TypeTag >::linearizeTimer_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::linearizeTimer().

linearSolver_

template<class TypeTag >

LinearSolverBackend Opm::TpsaNewtonMethod< TypeTag >::linearSolver_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::linearSolver().

numIterations_

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numIterations_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::begin_(), Opm::TpsaNewtonMethod< TypeTag >::endIteration_(), and Opm::TpsaNewtonMethod< TypeTag >::numIterations().

numLinearizations_

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numLinearizations_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::begin_(), Opm::TpsaNewtonMethod< TypeTag >::linearizeDomain_(), Opm::TpsaNewtonMethod< TypeTag >::numLinearizations(), and Opm::TpsaNewtonMethod< TypeTag >::preSolve_().

numTotLinearIterations_

template<class TypeTag >

int Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), Opm::TpsaNewtonMethod< TypeTag >::begin_(), and Opm::TpsaNewtonMethod< TypeTag >::numTotLinearIterations().

params_

template<class TypeTag >

TpsaNewtonMethodParams<Scalar> Opm::TpsaNewtonMethod< TypeTag >::params_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::minIterations(), Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), Opm::TpsaNewtonMethod< TypeTag >::proceed_(), Opm::TpsaNewtonMethod< TypeTag >::singleIteration_(), Opm::TpsaNewtonMethod< TypeTag >::tolerance(), Opm::TpsaNewtonMethod< TypeTag >::TpsaNewtonMethod(), and Opm::TpsaNewtonMethod< TypeTag >::verbosity_().

prePostProcessTimer_

template<class TypeTag >

Timer Opm::TpsaNewtonMethod< TypeTag >::prePostProcessTimer_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::prePostProcessTimer().

simulator_

template<class TypeTag >

Simulator& Opm::TpsaNewtonMethod< TypeTag >::simulator_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::model(), Opm::TpsaNewtonMethod< TypeTag >::preSolve_(), Opm::TpsaNewtonMethod< TypeTag >::problem(), and Opm::TpsaNewtonMethod< TypeTag >::verbosity_().

solveTimer_

template<class TypeTag >

Timer Opm::TpsaNewtonMethod< TypeTag >::solveTimer_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::solveTimer().

updateTimer_

template<class TypeTag >

Timer Opm::TpsaNewtonMethod< TypeTag >::updateTimer_

protected

Referenced by Opm::TpsaNewtonMethod< TypeTag >::apply(), and Opm::TpsaNewtonMethod< TypeTag >::updateTimer().

---

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

• tpsanewtonmethod.hpp