A convergence test based on the weighted root-mean-square norm of the update to the null vector component for turning point location. More...

#include <LOCA_Bifurcation_TPBord_StatusTest_NullVectorNormWRMS.H>

Inheritance diagram for LOCA::Bifurcation::TPBord::StatusTest::NullVectorNormWRMS:

[legend]

Collaboration diagram for LOCA::Bifurcation::TPBord::StatusTest::NullVectorNormWRMS:

[legend]

Public Member Functions
	NullVectorNormWRMS (double rtol, double atol, double tol)
	Constructor. More...

virtual	~NullVectorNormWRMS ()
	Destructor.

virtual NOX::StatusTest::StatusType	checkStatus (const NOX::Solver::Generic &problem)
	Evaluates convergence criteria specified above.

virtual NOX::StatusTest::StatusType	getStatus () const
	Returns status as defined above.

virtual std::ostream &	print (ostream &stream, int indent=0) const
	Prints current status.


double	getNullVectorNormWRMS () const
	Returns the value of weighted parameter update norm.

double	getRTOL () const
	Returns the realative tolerance set in the constructor.

double	getATOL () const
	Returns the absolute tolerance set in the constructor.

double	getTOL () const
	Returns the tolerance set in the constructor.

Public Member Functions inherited from NOX::StatusTest::Generic
	Generic ()
	Constructor.

virtual	~Generic ()
	Destructor.

virtual NOX::StatusTest::StatusType	checkStatus (const NOX::Solver::Generic &problem, NOX::StatusTest::CheckType checkType)=0
	Test the stopping criterion More...

virtual std::ostream &	print (std::ostream &stream, int indent=0) const =0
	Output formatted description of stopping test to output stream.

Protected Attributes
double	rtol
	Relative tolerance for convergence.

double	atol
	Absolute tolerance for convergence.

double	tol
	Tolerance for convergence.

double	normWRMS
	Current value of the parameter update norm.

NOX::StatusTest::StatusType	status
	Status of entire status test

Detailed Description

A convergence test based on the weighted root-mean-square norm of the update to the null vector component for turning point location.

Let $n$ be the approximation to the null vector for turning point tracking (see LOCA::Bifurcation::TPBord::ExtendedGroup). This convergence test defines convergence for the null vector when the following is true

$\sqrt{\frac{1}{N}\sum_{i=1}^N\left(\frac{|n_i-(n_0)_i|}{\epsilon_r|n_i| + \epsilon_a}\right)} < \tau$

where $n_0$ is the previous approximation to the null vector, $N$ is the length of $n$ , $\epsilon_r$ is the relative tolerance, $\epsilon_a$ is the absolute tolerance, and $\tau$ is an overall scale factor (typically $\tau = 1$ ).

Note that this status test deals only with the null vector component of the turning point equations. This status test should be combined with other status tests for the solution and parameter components (using NOX::StatusTest::Combo and LOCA::StatusTest::Wrapper) to build a composite status test for the entire system.

Also note that if the group returned by the getSolutionGroup() method of the solver supplied in checkStatus() is not a turning point group (i.e., not derived from LOCA::Bifurcation::TPBord::ExtendedGroup), checkStatus() returns NOX::StatusTest::Converged. This allows the status test to be used in situations other than turning point tracking, e.g., steady-state solves, without raising error conditions.

Constructor & Destructor Documentation

LOCA::Bifurcation::TPBord::StatusTest::NullVectorNormWRMS::NullVectorNormWRMS	(	double	rtol,
		double	atol,
		double	tol
	)

Constructor.

rtol is the relative tolerance $\epsilon_r$ , atol is the absolute tolerance $\epsilon_a$ , and tol is the overall scale factor $\tau$ defined above.

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

LOCA_Bifurcation_TPBord_StatusTest_NullVectorNormWRMS.H
LOCA_Bifurcation_TPBord_StatusTest_NullVectorNormWRMS.C

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation