#include <LOCA_Pitchfork_MooreSpence_PhippsBordering.H>
Public Member Functions | |
PhippsBordering (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RCP< Teuchos::ParameterList > &solverParams) | |
Constructor. | |
virtual | ~PhippsBordering () |
Destructor. | |
virtual void | setBlocks (const Teuchos::RCP< LOCA::Pitchfork::MooreSpence::AbstractGroup > &group, const Teuchos::RCP< LOCA::Pitchfork::MooreSpence::ExtendedGroup > &pfGroup, const Teuchos::RCP< const NOX::Abstract::MultiVector > &asymMultiVector, const Teuchos::RCP< const NOX::Abstract::Vector > &nullVector, const Teuchos::RCP< const NOX::Abstract::Vector > &JnVector, const Teuchos::RCP< const NOX::Abstract::Vector > &dfdp, const Teuchos::RCP< const NOX::Abstract::Vector > &dJndp) |
Set blocks in extended linear system. | |
virtual NOX::Abstract::Group::ReturnType | solve (Teuchos::ParameterList ¶ms, const LOCA::Pitchfork::MooreSpence::ExtendedMultiVector &input, LOCA::Pitchfork::MooreSpence::ExtendedMultiVector &result) const |
Solves the extended system as defined above. | |
Protected Member Functions | |
NOX::Abstract::Group::ReturnType | solveContiguous (Teuchos::ParameterList ¶ms, const NOX::Abstract::MultiVector &input_x, const NOX::Abstract::MultiVector &input_null, const NOX::Abstract::MultiVector::DenseMatrix &input_slack, const NOX::Abstract::MultiVector::DenseMatrix &input_param, NOX::Abstract::MultiVector &result_x, NOX::Abstract::MultiVector &result_null, NOX::Abstract::MultiVector::DenseMatrix &result_slack, NOX::Abstract::MultiVector::DenseMatrix &result_param) 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::Pitchfork::MooreSpence::AbstractGroup > | group |
Underlying group. | |
Teuchos::RCP < LOCA::Pitchfork::MooreSpence::ExtendedGroup > | pfGroup |
Turning point group. | |
Teuchos::RCP< const NOX::Abstract::MultiVector > | asymMultiVector |
Antisymmetric vector (psi). | |
Teuchos::RCP< const NOX::Abstract::Vector > | asymVector |
Antisymmetric vector (psi). | |
Teuchos::RCP< const NOX::Abstract::Vector > | nullVector |
Null vector. | |
Teuchos::RCP< const NOX::Abstract::Vector > | JnVector |
Jacobian times null vector. | |
Teuchos::RCP< const NOX::Abstract::Vector > | dfdp |
df/dp | |
Teuchos::RCP< const NOX::Abstract::Vector > | dJndp |
d(Jn)/dp | |
Teuchos::RCP < LOCA::BorderedSolver::AbstractStrategy > | borderedSolver |
Bordered solver for solving (n+1)x(n+1) sets of equations. | |
Teuchos::RCP < NOX::Abstract::MultiVector > | nullMultiVector |
Null vector as a multivector. | |
Teuchos::RCP < NOX::Abstract::MultiVector > | JnMultiVector |
Jacobian times null vector as a multivector. | |
double | sigma |
Norm of JnVector. | |
Private Member Functions | |
PhippsBordering (const PhippsBordering &) | |
Private to prohibit copying. | |
PhippsBordering & | operator= (const PhippsBordering &) |
Private to prohibit copying. |
This class solves the Moore-Spence pitchfork Newton equations:
via the following modified block elimination scheme:
where and
. Each bordered solve is implemented by a LOCA::BorderedSolver::AbstractStrategy strategy object.
Definition at line 145 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
LOCA::Pitchfork::MooreSpence::PhippsBordering::PhippsBordering | ( | const Teuchos::RCP< LOCA::GlobalData > & | global_data, | |
const Teuchos::RCP< LOCA::Parameter::SublistParser > & | topParams, | |||
const Teuchos::RCP< Teuchos::ParameterList > & | solverParams | |||
) |
Constructor.
global_data | [in] Global data object | |
topParams | [in] Parsed top-level parameter list | |
solverParams | [in] Bordered solver parameters. Instantiates a bordered solver for solving the bordeded systems described above. See LOCA::BorderedSolver::Factory for a description of available solvers. |
Definition at line 52 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.C.
References borderedSolver, globalData, and solverParams.
LOCA::Pitchfork::MooreSpence::PhippsBordering::~PhippsBordering | ( | ) | [virtual] |
LOCA::Pitchfork::MooreSpence::PhippsBordering::PhippsBordering | ( | const PhippsBordering & | ) | [private] |
Private to prohibit copying.
void LOCA::Pitchfork::MooreSpence::PhippsBordering::setBlocks | ( | const Teuchos::RCP< LOCA::Pitchfork::MooreSpence::AbstractGroup > & | group, | |
const Teuchos::RCP< LOCA::Pitchfork::MooreSpence::ExtendedGroup > & | pfGroup, | |||
const Teuchos::RCP< const NOX::Abstract::MultiVector > & | asymMultiVector, | |||
const Teuchos::RCP< const NOX::Abstract::Vector > & | nullVector, | |||
const Teuchos::RCP< const NOX::Abstract::Vector > & | JnVector, | |||
const Teuchos::RCP< const NOX::Abstract::Vector > & | dfdp, | |||
const Teuchos::RCP< const NOX::Abstract::Vector > & | dJndp | |||
) | [virtual] |
Set blocks in extended linear system.
group | [in] Underlying group representing J | |
pfGroup | [in] Pitchfork group representing the pitchfork equations. | |
asymMultiVector | [in] Multivector representing the asymmetric vector | |
nullVector | [in] Vector representing v | |
JnVector | [in] Vector representing Jv | |
dfdp | [in] Vector representing df/dp | |
dJndp | [in] Vector representing d(Jv)/dp |
Implements LOCA::Pitchfork::MooreSpence::SolverStrategy.
Definition at line 81 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.C.
References asymMultiVector, asymVector, borderedSolver, NOX::DeepCopy, dfdp, dJndp, globalData, group, JnMultiVector, JnVector, nullMultiVector, nullVector, pfGroup, sigma, and NOX::Abstract::Vector::TwoNorm.
NOX::Abstract::Group::ReturnType LOCA::Pitchfork::MooreSpence::PhippsBordering::solve | ( | Teuchos::ParameterList & | params, | |
const LOCA::Pitchfork::MooreSpence::ExtendedMultiVector & | input, | |||
LOCA::Pitchfork::MooreSpence::ExtendedMultiVector & | result | |||
) | const [virtual] |
Solves the extended system as defined above.
The params argument is the linear solver parameters.
Implements LOCA::Pitchfork::MooreSpence::SolverStrategy.
Definition at line 118 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.C.
References asymVector, dfdp, dJndp, LOCA::Pitchfork::MooreSpence::ExtendedMultiVector::getBifParams(), LOCA::Pitchfork::MooreSpence::ExtendedMultiVector::getNullMultiVec(), LOCA::Pitchfork::MooreSpence::ExtendedMultiVector::getSlacks(), LOCA::Pitchfork::MooreSpence::ExtendedMultiVector::getXMultiVec(), LOCA::Extended::MultiVector::numVectors(), and solveContiguous().
NOX::Abstract::Group::ReturnType LOCA::Pitchfork::MooreSpence::PhippsBordering::solveContiguous | ( | Teuchos::ParameterList & | params, | |
const NOX::Abstract::MultiVector & | input_x, | |||
const NOX::Abstract::MultiVector & | input_null, | |||
const NOX::Abstract::MultiVector::DenseMatrix & | input_slack, | |||
const NOX::Abstract::MultiVector::DenseMatrix & | input_param, | |||
NOX::Abstract::MultiVector & | result_x, | |||
NOX::Abstract::MultiVector & | result_null, | |||
NOX::Abstract::MultiVector::DenseMatrix & | result_slack, | |||
NOX::Abstract::MultiVector::DenseMatrix & | result_param | |||
) | const [protected] |
Solves equations with contiguous arguments.
Definition at line 219 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.C.
References asymMultiVector, asymVector, borderedSolver, NOX::Abstract::MultiVector::clone(), NOX::Abstract::Group::Failed, globalData, group, JnVector, nullMultiVector, nullVector, NOX::Abstract::MultiVector::numVectors(), NOX::Abstract::Group::Ok, pfGroup, NOX::ShapeCopy, sigma, and NOX::Abstract::MultiVector::subView().
Referenced by solve().
PhippsBordering& LOCA::Pitchfork::MooreSpence::PhippsBordering::operator= | ( | const PhippsBordering & | ) | [private] |
Private to prohibit copying.
Teuchos::RCP<LOCA::GlobalData> LOCA::Pitchfork::MooreSpence::PhippsBordering::globalData [protected] |
Global data object.
Definition at line 226 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by PhippsBordering(), setBlocks(), and solveContiguous().
Teuchos::RCP<Teuchos::ParameterList> LOCA::Pitchfork::MooreSpence::PhippsBordering::solverParams [protected] |
Solver parameters.
Definition at line 229 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by PhippsBordering().
Teuchos::RCP<LOCA::Pitchfork::MooreSpence::AbstractGroup> LOCA::Pitchfork::MooreSpence::PhippsBordering::group [protected] |
Underlying group.
Definition at line 232 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<LOCA::Pitchfork::MooreSpence::ExtendedGroup> LOCA::Pitchfork::MooreSpence::PhippsBordering::pfGroup [protected] |
Turning point group.
Definition at line 235 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<const NOX::Abstract::MultiVector> LOCA::Pitchfork::MooreSpence::PhippsBordering::asymMultiVector [protected] |
Antisymmetric vector (psi).
Definition at line 238 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Pitchfork::MooreSpence::PhippsBordering::asymVector [protected] |
Antisymmetric vector (psi).
Definition at line 241 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), solve(), and solveContiguous().
Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Pitchfork::MooreSpence::PhippsBordering::nullVector [protected] |
Null vector.
Definition at line 244 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Pitchfork::MooreSpence::PhippsBordering::JnVector [protected] |
Jacobian times null vector.
Definition at line 247 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Pitchfork::MooreSpence::PhippsBordering::dfdp [protected] |
df/dp
Definition at line 250 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solve().
Teuchos::RCP<const NOX::Abstract::Vector> LOCA::Pitchfork::MooreSpence::PhippsBordering::dJndp [protected] |
d(Jn)/dp
Definition at line 253 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solve().
Teuchos::RCP<LOCA::BorderedSolver::AbstractStrategy> LOCA::Pitchfork::MooreSpence::PhippsBordering::borderedSolver [protected] |
Bordered solver for solving (n+1)x(n+1) sets of equations.
Definition at line 256 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by PhippsBordering(), setBlocks(), and solveContiguous().
Teuchos::RCP<NOX::Abstract::MultiVector> LOCA::Pitchfork::MooreSpence::PhippsBordering::nullMultiVector [protected] |
Null vector as a multivector.
Definition at line 259 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().
Teuchos::RCP<NOX::Abstract::MultiVector> LOCA::Pitchfork::MooreSpence::PhippsBordering::JnMultiVector [protected] |
Jacobian times null vector as a multivector.
Definition at line 262 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks().
double LOCA::Pitchfork::MooreSpence::PhippsBordering::sigma [protected] |
Norm of JnVector.
Definition at line 265 of file LOCA_Pitchfork_MooreSpence_PhippsBordering.H.
Referenced by setBlocks(), and solveContiguous().