#include <Albany_Application.hpp>

Collaboration diagram for Albany::Application:

Public Types
enum	SolutionMethod { Steady, Transient, Continuation, Eigensolve }
Public Member Functions
	Application (const Teuchos::RCP< const Epetra_Comm > &comm, const Teuchos::RCP< Teuchos::ParameterList > &params, const Teuchos::RCP< const Epetra_Vector > &initial_guess=Teuchos::null)
	Constructor.
	~Application ()
	Destructor.
Teuchos::RCP < Albany::AbstractDiscretization >	getDiscretization () const
	Get underlying abstract discretization.
Teuchos::RCP < Albany::AbstractProblem >	getProblem () const
	Get problem object.
Teuchos::RCP< const Epetra_Comm >	getComm () const
	Get communicator.
Teuchos::RCP< const Epetra_Map >	getMap () const
	Get DOF map.
Teuchos::RCP< const Epetra_CrsGraph >	getJacobianGraph () const
	Get Jacobian graph.
Teuchos::RCP< Epetra_Operator >	getPreconditioner ()
	Get Preconditioner Operator.
Teuchos::RCP < AAdapt::AdaptiveSolutionManager >	getAdaptSolMgr ()
	Get the solution memory manager.
Teuchos::RCP< ParamLib >	getParamLib ()
	Get parameter library.
SolutionMethod	getSolutionMethod () const
	Get solution method.
int	getNumResponses () const
	Get number of responses.
Teuchos::RCP < AbstractResponseFunction >	getResponse (int i) const
	Get response function.
bool	suppliesPreconditioner () const
	Return whether problem wants to use its own preconditioner.
Teuchos::RCP < Stokhos::OrthogPolyExpansion < int, double > >	getStochasticExpansion ()
	Get stochastic expansion.
void	computeGlobalResidual (const double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, Epetra_Vector &f)
	Intialize stochastic Galerkin method.
void	computeGlobalJacobian (const double alpha, const double beta, const double omega, const double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, Epetra_Vector *f, Epetra_CrsMatrix &jac)
	Compute global Jacobian.
void	computeGlobalPreconditioner (const Teuchos::RCP< Epetra_CrsMatrix > &jac, const Teuchos::RCP< Epetra_Operator > &prec)
	Compute global Preconditioner.
void	computeGlobalTangent (const double alpha, const double beta, const double omega, const double current_time, bool sum_derivs, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, ParamVec deriv_p, const Epetra_MultiVector Vx, const Epetra_MultiVector Vxdot, const Epetra_MultiVector Vxdotdot, const Epetra_MultiVector Vp, Epetra_Vector f, Epetra_MultiVector JV, Epetra_MultiVector fp)
	Compute global Tangent.
void	evaluateResponse (int response_index, const double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, Epetra_Vector &g)
	Evaluate response functions.
void	evaluateResponseTangent (int response_index, const double alpha, const double beta, const double omega, const double current_time, bool sum_derivs, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, ParamVec deriv_p, const Epetra_MultiVector Vxdot, const Epetra_MultiVector Vxdotdot, const Epetra_MultiVector Vx, const Epetra_MultiVector Vp, Epetra_Vector g, Epetra_MultiVector gx, Epetra_MultiVector gp)
	Evaluate tangent = alphadg/dxVx + betadg/dxdotVxdot + dg/dp*Vp.
void	evaluateResponseDerivative (int response_index, const double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x, const Teuchos::Array< ParamVec > &p, ParamVec deriv_p, Epetra_Vector g, const EpetraExt::ModelEvaluator::Derivative &dg_dx, const EpetraExt::ModelEvaluator::Derivative &dg_dxdot, const EpetraExt::ModelEvaluator::Derivative &dg_dxdotdot, const EpetraExt::ModelEvaluator::Derivative &dg_dp)
	Evaluate gradient = dg/dx, dg/dxdot, dg/dp.
PHAL::AlbanyTraits::Residual::ScalarT &	getValue (const std::string &n)
	Provide access to shapeParameters -- no AD.
StateManager &	getStateMgr ()
	Class to manage state variables (a.k.a. history).
void	evaluateStateFieldManager (const double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector &x)
	Evaluate state field manager.
int	getNumWorksets ()
	Access to number of worksets - needed for working with StateManager.
template<typename EvalT >
void	loadWorksetBucketInfo (PHAL::Workset &workset, const int &ws)
	Routine to get workset (bucket) sized info needed by all Evaluation types.
void	loadBasicWorksetInfo (PHAL::Workset &workset, double current_time)
	Routine to load some basic workset info needed by many Evaluation types.
void	loadWorksetJacobianInfo (PHAL::Workset &workset, const double &alpha, const double &beta, const double &omega)
void	loadWorksetNodesetInfo (PHAL::Workset &workset)
	Routine to load common nodeset info into workset.
void	loadWorksetSidesetInfo (PHAL::Workset &workset, const int ws)
	Routine to load common sideset info into workset.
void	setupBasicWorksetInfo (PHAL::Workset &workset, double current_time, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector *x, const Teuchos::Array< ParamVec > &p)
void	setupTangentWorksetInfo (PHAL::Workset &workset, double current_time, bool sum_derivs, const Epetra_Vector xdot, const Epetra_Vector xdotdot, const Epetra_Vector x, const Teuchos::Array< ParamVec > &p, ParamVec deriv_p, const Epetra_MultiVector Vxdot, const Epetra_MultiVector Vxdotdot, const Epetra_MultiVector Vx, const Epetra_MultiVector Vp)
void	postRegSetup (std::string eval)
Public Attributes
bool	is_adjoint
Protected Member Functions
void	determinePiroSolver (const Teuchos::RCP< Teuchos::ParameterList > &topLevelParams)
Protected Attributes
Teuchos::RCP< const Epetra_Comm >	comm
	Communicator.
Teuchos::RCP < Teuchos::FancyOStream >	out
	Output stream, defaults to pronting just Proc 0.
Teuchos::RCP < Albany::AbstractDiscretization >	disc
	Element discretization.
Teuchos::RCP < Albany::AbstractProblem >	problem
	Problem class.
Teuchos::RCP< ParamLib >	paramLib
	Parameter library.
Teuchos::RCP < AAdapt::AdaptiveSolutionManager >	solMgr
	Solution memory manager.
Teuchos::Array< Teuchos::RCP < Albany::AbstractResponseFunction > >	responses
	Response functions.
Teuchos::ArrayRCP < Teuchos::RCP < PHX::FieldManager < PHAL::AlbanyTraits > > >	fm
	Phalanx Field Manager for volumetric fills.
Teuchos::RCP < PHX::FieldManager < PHAL::AlbanyTraits > >	dfm
	Phalanx Field Manager for Dirichlet Conditions.
Teuchos::ArrayRCP < Teuchos::RCP < PHX::FieldManager < PHAL::AlbanyTraits > > >	nfm
	Phalanx Field Manager for Neumann Conditions.
Teuchos::Array< Teuchos::RCP < PHX::FieldManager < PHAL::AlbanyTraits > > >	sfm
	Phalanx Field Manager for states.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > >	sg_basis
	Stochastic Galerkin basis.
Teuchos::RCP< const Stokhos::Quadrature< int, double > >	sg_quad
	Stochastic Galerkin quadrature.
Teuchos::RCP < Stokhos::OrthogPolyExpansion < int, double > >	sg_expansion
	Stochastic Galerkin expansion.
Teuchos::RCP< const EpetraExt::MultiComm >	product_comm
	Product multi-comm.
Teuchos::RCP< const Epetra_BlockMap >	sg_overlap_map
	Overlap stochastic map.
Teuchos::RCP < Stokhos::EpetraVectorOrthogPoly >	sg_overlapped_x
	SG overlapped solution vectors.
Teuchos::RCP < Stokhos::EpetraVectorOrthogPoly >	sg_overlapped_xdot
	SG overlapped time derivative vectors.
Teuchos::RCP < Stokhos::EpetraVectorOrthogPoly >	sg_overlapped_xdotdot
Teuchos::RCP < Stokhos::EpetraVectorOrthogPoly >	sg_overlapped_f
	SG overlapped residual vectors.
Teuchos::RCP < Stokhos::VectorOrthogPoly < Epetra_CrsMatrix > >	sg_overlapped_jac
	Overlapped Jacobian matrixs.
Teuchos::RCP < Stokhos::ProductEpetraVector >	mp_overlapped_x
	MP overlapped solution vectors.
Teuchos::RCP < Stokhos::ProductEpetraVector >	mp_overlapped_xdot
	MP overlapped time derivative vectors.
Teuchos::RCP < Stokhos::ProductEpetraVector >	mp_overlapped_xdotdot
Teuchos::RCP < Stokhos::ProductEpetraVector >	mp_overlapped_f
	MP overlapped residual vectors.
Teuchos::RCP < Stokhos::ProductContainer < Epetra_CrsMatrix > >	mp_overlapped_jac
	Overlapped Jacobian matrixs.
bool	physicsBasedPreconditioner
	Data for Physics-Based Preconditioners.
Teuchos::RCP < Teuchos::ParameterList >	tekoParams
SolutionMethod	solMethod
	Type of solution method.
int	writeToMatrixMarketJac
	Integer specifying whether user wants to write Jacobian to MatrixMarket file.
int	writeToMatrixMarketRes
int	writeToCoutJac
	Integer specifying whether user wants to write Jacobian and residual to Standard output (cout).
int	writeToCoutRes
bool	shapeParamsHaveBeenReset
	Shape Optimization data.
std::vector< RealType >	shapeParams
std::vector< std::string >	shapeParamNames
unsigned int	neq
Teuchos::RCP < Teko::InverseLibrary >	inverseLib
	Teko stuff.
Teuchos::RCP < Teko::InverseFactory >	inverseFac
Teuchos::RCP< Epetra_Operator >	wrappedJac
std::vector< int >	blockDecomp
std::set< std::string >	setupSet
int	phxGraphVisDetail
int	stateGraphVisDetail
StateManager	stateMgr
bool	morphFromInit
bool	ignore_residual_in_jacobian
double	perturbBetaForDirichlets
	To prevent a singular mass matrix associated with Dirichlet.
Private Member Functions
	Application (const Application &)
	Private to prohibit copying.
Application &	operator= (const Application &)
	Private to prohibit copying.
Teuchos::RCP< Epetra_Operator >	buildWrappedOperator (const Teuchos::RCP< Epetra_Operator > &Jac, const Teuchos::RCP< Epetra_Operator > &wrapInput, bool reorder=false) const
	Call to Teko to build strided block operator.
void	registerShapeParameters ()
	Utility function to set up ShapeParameters through Sacado.
void	defineTimers ()

Detailed Description

Definition at line 60 of file Albany_Application.hpp.

Member Enumeration Documentation

enum Albany::Application::SolutionMethod

Enumerator:

Steady
Transient
Continuation
Eigensolve

Definition at line 64 of file Albany_Application.hpp.

Constructor & Destructor Documentation

Albany::Application::Application	(	const Teuchos::RCP< const Epetra_Comm > &	comm,
		const Teuchos::RCP< Teuchos::ParameterList > &	params,
		const Teuchos::RCP< const Epetra_Vector > &	initial_guess = `Teuchos::null`
	)

Constructor.

Albany::Application::~Application ( )

Destructor.

Definition at line 285 of file Albany_Application.cpp.

Albany::Application::Application ( const Application & ) [private]

Private to prohibit copying.

Member Function Documentation

RCP< Albany::AbstractDiscretization > Albany::Application::getDiscretization ( ) const

Get underlying abstract discretization.

Definition at line 294 of file Albany_Application.cpp.

RCP< Albany::AbstractProblem > Albany::Application::getProblem ( ) const

Get problem object.

Definition at line 301 of file Albany_Application.cpp.

RCP< const Epetra_Comm > Albany::Application::getComm ( ) const

Get communicator.

Definition at line 308 of file Albany_Application.cpp.

RCP< const Epetra_Map > Albany::Application::getMap ( ) const

Get DOF map.

Definition at line 315 of file Albany_Application.cpp.

RCP< const Epetra_CrsGraph > Albany::Application::getJacobianGraph ( ) const

Get Jacobian graph.

Definition at line 322 of file Albany_Application.cpp.

RCP< Epetra_Operator > Albany::Application::getPreconditioner ( )

Get Preconditioner Operator.

Definition at line 330 of file Albany_Application.cpp.

Teuchos::RCP<AAdapt::AdaptiveSolutionManager> Albany::Application::getAdaptSolMgr ( ) [inline]

Get the solution memory manager.

Definition at line 94 of file Albany_Application.hpp.

RCP< ParamLib > Albany::Application::getParamLib ( )

Get parameter library.

Definition at line 357 of file Albany_Application.cpp.

SolutionMethod Albany::Application::getSolutionMethod ( ) const [inline]

Get solution method.

Definition at line 100 of file Albany_Application.hpp.

int Albany::Application::getNumResponses ( ) const

Get number of responses.

Definition at line 364 of file Albany_Application.cpp.

Teuchos::RCP< Albany::AbstractResponseFunction > Albany::Application::getResponse ( int i ) const

Get response function.

Definition at line 370 of file Albany_Application.cpp.

bool Albany::Application::suppliesPreconditioner ( ) const

Return whether problem wants to use its own preconditioner.

Definition at line 377 of file Albany_Application.cpp.

RCP< Stokhos::OrthogPolyExpansion< int, double > > Albany::Application::getStochasticExpansion ( )

Get stochastic expansion.

Definition at line 384 of file Albany_Application.cpp.

void Albany::Application::computeGlobalResidual	(	const double	current_time,
		const Epetra_Vector *	xdot,
		const Epetra_Vector *	xdotdot,
		const Epetra_Vector &	x,
		const Teuchos::Array< ParamVec > &	p,
		Epetra_Vector &	f
	)

Intialize stochastic Galerkin method.

Compute global residual