LOCA::Hopf::MooreSpence::SalingerBordering Class Reference

Moore-Spence Hopf solver strategy based on "Salinger" bordering. This is the classic 5-solve Hopf bordering method. More...

#include <LOCA_Hopf_MooreSpence_SalingerBordering.H>

Inheritance diagram for LOCA::Hopf::MooreSpence::SalingerBordering:

Inheritance graph
[legend]
Collaboration diagram for LOCA::Hopf::MooreSpence::SalingerBordering:

Collaboration graph
[legend]

List of all members.

Public Member Functions

 SalingerBordering (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RCP< Teuchos::ParameterList > &solverParams)
 Constructor.
virtual ~SalingerBordering ()
 Destructor.
virtual void setBlocks (const Teuchos::RCP< LOCA::Hopf::MooreSpence::AbstractGroup > &group, const Teuchos::RCP< LOCA::Hopf::MooreSpence::ExtendedGroup > &hopfGroup, const Teuchos::RCP< const NOX::Abstract::Vector > &yVector, const Teuchos::RCP< const NOX::Abstract::Vector > &zVector, const Teuchos::RCP< const NOX::Abstract::Vector > &CeRealVector, const Teuchos::RCP< const NOX::Abstract::Vector > &CeImagVector, const Teuchos::RCP< const NOX::Abstract::Vector > &dfdp, const Teuchos::RCP< const NOX::Abstract::Vector > &dCedpReal, const Teuchos::RCP< const NOX::Abstract::Vector > &dCedpImag, const Teuchos::RCP< const NOX::Abstract::Vector > &ByVector, const Teuchos::RCP< const NOX::Abstract::Vector > &mBzVector, double w)
 Set blocks in extended linear system.
virtual
NOX::Abstract::Group::ReturnType 		solve (Teuchos::ParameterList &params, const LOCA::Hopf::MooreSpence::ExtendedMultiVector &input, LOCA::Hopf::MooreSpence::ExtendedMultiVector &result) const
 Solves the extended system as defined above.

Protected Member Functions

NOX::Abstract::Group::ReturnType solveContiguous (Teuchos::ParameterList &params, const NOX::Abstract::MultiVector &input_x, const NOX::Abstract::MultiVector &input_y, const NOX::Abstract::MultiVector &input_z, const NOX::Abstract::MultiVector::DenseMatrix &input_w, const NOX::Abstract::MultiVector::DenseMatrix &input_p, NOX::Abstract::MultiVector &result_x, NOX::Abstract::MultiVector &result_y, NOX::Abstract::MultiVector &result_z, NOX::Abstract::MultiVector::DenseMatrix &result_w, NOX::Abstract::MultiVector::DenseMatrix &result_p) const
 Solves equations with contiguous arguments.

Protected Attributes

Teuchos::RCP< LOCA::GlobalData > globalData
 Global data object.
Teuchos::RCP
< Teuchos::ParameterList > 		solverParams
 Solver parameters.
Teuchos::RCP
< LOCA::Hopf::MooreSpence::AbstractGroup > 		group
 Underlying group.
Teuchos::RCP
< LOCA::Hopf::MooreSpence::ExtendedGroup > 		hopfGroup
 Hopf group.
Teuchos::RCP< const
NOX::Abstract::Vector > 		yVector
 Real component of eigenvector.
Teuchos::RCP< const
NOX::Abstract::Vector > 		zVector
 Imaginary component of eigenvector.
Teuchos::RCP< const
NOX::Abstract::Vector > 		CeRealVector
 Jy-wBz vector.
Teuchos::RCP< const
NOX::Abstract::Vector > 		CeImagVector
 Jz+wBy vector.
Teuchos::RCP< const
NOX::Abstract::Vector > 		dfdp
 df/dp vector
Teuchos::RCP< const
NOX::Abstract::Vector > 		dCedpReal
 d(Jy-wBz)/dp vector
Teuchos::RCP< const
NOX::Abstract::Vector > 		dCedpImag
 d(Jz+wBy)/dp vector
Teuchos::RCP< const
NOX::Abstract::Vector > 		ByVector
 By vector.
Teuchos::RCP< const
NOX::Abstract::Vector > 		minusBzVector
 -Bz vector
double w
 Frequency w.

Private Member Functions

 SalingerBordering (const SalingerBordering &)
 Private to prohibit copying.
SalingerBordering & operator= (const SalingerBordering &)
 Private to prohibit copying.

Detailed Description

Moore-Spence Hopf solver strategy based on "Salinger" bordering. This is the classic 5-solve Hopf bordering method.

This class solves the Moore-Spence Hopf Newton equations:

\[ \begin{bmatrix} J & 0 & 0 & 0 & f_p \\ (Jy-wBz)_x & J & -wB & -Bz & (Jy-wBz)_p \\ (Jz+wBy)_x & wB & J & By & (Jz+wBy)_p \\ 0 & \phi^T & 0 & 0 & 0 \\ 0 & 0 & \phi^T & 0 & 0 \end{bmatrix} \begin{bmatrix} X \\ Y \\ Z \\ \omega \\ \lambda \end{bmatrix} = \begin{bmatrix} F \\ G \\ H \\ u \\ v \end{bmatrix}. \]

via the following block elimination scheme:

\[ \begin{split} J [A \; b] &= [F \; f_p] \\ \begin{bmatrix} J & -wB \\ wB & J \end{bmatrix} \begin{bmatrix} C & e & g \\ D & f & h \end{bmatrix} &= \begin{bmatrix} G - (Jy-wBz)_x A & (Jy-wBz)_p - (Jy-wBz)_x b & -Bz \\ H - (Jz+wBy)_x A & (Jz+wBy)_p - (Jy-wBz)_x b & By \end{bmatrix} \\ \lambda &= \frac{(\phi^T h)(\phi^T C-u)-(\phi^T g)(\phi^T D-v)} {(\phi^T h)(\phi^T e)-(\phi^T g)(\phi^T f)} \\ \omega &= \frac{\phi^T D - v - (\phi^T f)\lambda}{\phi^T h} \\ X &= A - b \lambda \\ \begin{bmatrix} Y \\ Z \end{bmatrix} &= \begin{bmatrix} C \\ D \end{bmatrix} - \begin{bmatrix} e \\ f \end{bmatrix}\lambda - \begin{bmatrix} g \\ h \end{bmatrix}\omega \end{split} \]

Definition at line 131 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Constructor & Destructor Documentation

LOCA::Hopf::MooreSpence::SalingerBordering::SalingerBordering ( const Teuchos::RCP< LOCA::GlobalData > &  global_data,
const Teuchos::RCP< LOCA::Parameter::SublistParser > &  topParams,
const Teuchos::RCP< Teuchos::ParameterList > &  solverParams 
)

Constructor.

Parameters:
global_data [in] Global data object
topParams [in] Parsed top-level parameter list
solverParams [in] Bordered solver parameters. Currently none are referenced.

Definition at line 48 of file LOCA_Hopf_MooreSpence_SalingerBordering.C.

LOCA::Hopf::MooreSpence::SalingerBordering::~SalingerBordering (  )  [virtual]

Destructor.

Definition at line 68 of file LOCA_Hopf_MooreSpence_SalingerBordering.C.

LOCA::Hopf::MooreSpence::SalingerBordering::SalingerBordering ( const SalingerBordering &   )  [private]

Private to prohibit copying.

Member Function Documentation

void LOCA::Hopf::MooreSpence::SalingerBordering::setBlocks ( const Teuchos::RCP< LOCA::Hopf::MooreSpence::AbstractGroup > &  group,
const Teuchos::RCP< LOCA::Hopf::MooreSpence::ExtendedGroup > &  hopfGroup,
const Teuchos::RCP< const NOX::Abstract::Vector > &  yVector,
const Teuchos::RCP< const NOX::Abstract::Vector > &  zVector,
const Teuchos::RCP< const NOX::Abstract::Vector > &  CeRealVector,
const Teuchos::RCP< const NOX::Abstract::Vector > &  CeImagVector,
const Teuchos::RCP< const NOX::Abstract::Vector > &  dfdp,
const Teuchos::RCP< const NOX::Abstract::Vector > &  dCedpReal,
const Teuchos::RCP< const NOX::Abstract::Vector > &  dCedpImag,
const Teuchos::RCP< const NOX::Abstract::Vector > &  ByVector,
const Teuchos::RCP< const NOX::Abstract::Vector > &  mBzVector,
double  w 
) [virtual]

Set blocks in extended linear system.

Parameters:
group [in] Underlying group representing J
hopfGroup [in] Hopf group representing the Hopf equations.
yVector [in] Vector representing y
zVector [in] Vector representing z
CeRealVector [in] Vector representing Jy-wBz
CeImagVector [in] Vector representing Jz+wBy
dfdp [in] Vector representing df/dp
dCedpReal [in] Vector representing d(Jy-wBz)/dp
dCedpImag [in] Vector representing d(Jz+wBy)/dp
ByVector [in] Vector representing By
mBzVector [in] Vector representing -Bz
w [in] Bifurcation frequency w

Implements LOCA::Hopf::MooreSpence::SolverStrategy.

Definition at line 73 of file LOCA_Hopf_MooreSpence_SalingerBordering.C.

References ByVector, CeImagVector, CeRealVector, dCedpImag, dCedpReal, dfdp, group, hopfGroup, minusBzVector, w, yVector, and zVector.

NOX::Abstract::Group::ReturnType LOCA::Hopf::MooreSpence::SalingerBordering::solve ( Teuchos::ParameterList &  params,
const LOCA::Hopf::MooreSpence::ExtendedMultiVector &  input,
LOCA::Hopf::MooreSpence::ExtendedMultiVector &  result 
) const [virtual]

Solves the extended system as defined above.

The params argument is the linear solver parameters.

Implements LOCA::Hopf::MooreSpence::SolverStrategy.

Definition at line 102 of file LOCA_Hopf_MooreSpence_SalingerBordering.C.

References ByVector, dCedpImag, dCedpReal, dfdp, LOCA::Hopf::MooreSpence::ExtendedMultiVector::getBifParams(), LOCA::Hopf::MooreSpence::ExtendedMultiVector::getFrequencies(), LOCA::Hopf::MooreSpence::ExtendedMultiVector::getImagEigenMultiVec(), LOCA::Hopf::MooreSpence::ExtendedMultiVector::getRealEigenMultiVec(), LOCA::Hopf::MooreSpence::ExtendedMultiVector::getXMultiVec(), minusBzVector, LOCA::Extended::MultiVector::numVectors(), and solveContiguous().

NOX::Abstract::Group::ReturnType LOCA::Hopf::MooreSpence::SalingerBordering::solveContiguous ( Teuchos::ParameterList &  params,
const NOX::Abstract::MultiVector &  input_x,
const NOX::Abstract::MultiVector &  input_y,
const NOX::Abstract::MultiVector &  input_z,
const NOX::Abstract::MultiVector::DenseMatrix &  input_w,
const NOX::Abstract::MultiVector::DenseMatrix &  input_p,
NOX::Abstract::MultiVector &  result_x,
NOX::Abstract::MultiVector &  result_y,
NOX::Abstract::MultiVector &  result_z,
NOX::Abstract::MultiVector::DenseMatrix &  result_w,
NOX::Abstract::MultiVector::DenseMatrix &  result_p 
) const [protected]

Solves equations with contiguous arguments.

Definition at line 218 of file LOCA_Hopf_MooreSpence_SalingerBordering.C.

References CeImagVector, CeRealVector, NOX::Abstract::MultiVector::clone(), globalData, group, hopfGroup, NOX::Abstract::MultiVector::numVectors(), NOX::Abstract::Group::Ok, NOX::ShapeCopy, NOX::Abstract::MultiVector::subView(), w, yVector, and zVector.

Referenced by solve().

SalingerBordering& LOCA::Hopf::MooreSpence::SalingerBordering::operator= ( const SalingerBordering &   )  [private]

Private to prohibit copying.

Member Data Documentation

Teuchos::RCP<LOCA::GlobalData> LOCA::Hopf::MooreSpence::SalingerBordering::globalData [protected]

Global data object.

Definition at line 220 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by solveContiguous().

Teuchos::RCP<Teuchos::ParameterList> LOCA::Hopf::MooreSpence::SalingerBordering::solverParams [protected]

Solver parameters.

Definition at line 223 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Teuchos::RCP<LOCA::Hopf::MooreSpence::AbstractGroup> LOCA::Hopf::MooreSpence::SalingerBordering::group [protected]

Underlying group.

Definition at line 226 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<LOCA::Hopf::MooreSpence::ExtendedGroup> LOCA::Hopf::MooreSpence::SalingerBordering::hopfGroup [protected]

Hopf group.

Definition at line 229 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::yVector [protected]

Real component of eigenvector.

Definition at line 232 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::zVector [protected]

Imaginary component of eigenvector.

Definition at line 235 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::CeRealVector [protected]

Jy-wBz vector.

Definition at line 238 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::CeImagVector [protected]

Jz+wBy vector.

Definition at line 241 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::dfdp [protected]

df/dp vector

Definition at line 244 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solve().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::dCedpReal [protected]

d(Jy-wBz)/dp vector

Definition at line 247 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solve().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::dCedpImag [protected]

d(Jz+wBy)/dp vector

Definition at line 250 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solve().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::ByVector [protected]

By vector.

Definition at line 253 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solve().

Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Hopf::MooreSpence::SalingerBordering::minusBzVector [protected]

-Bz vector

Definition at line 256 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solve().

double LOCA::Hopf::MooreSpence::SalingerBordering::w [protected]

Frequency w.

Definition at line 259 of file LOCA_Hopf_MooreSpence_SalingerBordering.H.

Referenced by setBlocks(), and solveContiguous().

The documentation for this class was generated from the following files:
Generated on Wed Oct 21 14:28:34 2009 for Nonlinear Solver Project by  doxygen 1.5.9