#include <grins_solver.h>

Inheritance diagram for GRINS::Solver:

Public Member Functions
	Solver (const GetPot &input)

virtual	~Solver ()

virtual void	initialize (const GetPot &input, std::tr1::shared_ptr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system)

virtual void	solve (SolverContext &context)=0

virtual void	adjoint_qoi_parameter_sensitivity (SolverContext &, const libMesh::QoISet &, const libMesh::ParameterVector &, libMesh::SensitivityData &) const

virtual void	forward_qoi_parameter_sensitivity (SolverContext &, const libMesh::QoISet &, const libMesh::ParameterVector &, libMesh::SensitivityData &) const

void	steady_adjoint_solve (SolverContext &context)
	Do steady version of adjoint solve. More...

Protected Member Functions
void	set_solver_options (libMesh::DiffSolver &solver)

virtual void	init_time_solver (GRINS::MultiphysicsSystem *system)=0

Protected Attributes
unsigned int	_max_nonlinear_iterations

double	_relative_step_tolerance

double	_absolute_step_tolerance

double	_relative_residual_tolerance

double	_absolute_residual_tolerance

double	_initial_linear_tolerance

double	_minimum_linear_tolerance

unsigned int	_max_linear_iterations

bool	_continue_after_backtrack_failure

bool	_continue_after_max_iterations

bool	_solver_quiet

bool	_solver_verbose

std::map< std::string, GRINS::NBCContainer >	_neumann_bc_funcs

Detailed Description

Definition at line 58 of file grins_solver.h.

Constructor & Destructor Documentation

GRINS::Solver::Solver ( const GetPot & input )

Definition at line 44 of file grins_solver.C.

     : _max_nonlinear_iterations( input("linear-nonlinear-solver/max_nonlinear_iterations", 10 ) ),
       _relative_step_tolerance( input("linear-nonlinear-solver/relative_step_tolerance", 1.e-6 ) ),
       _absolute_step_tolerance( input("linear-nonlinear-solver/absolute_step_tolerance", 0.0 ) ),
       _relative_residual_tolerance( input("linear-nonlinear-solver/relative_residual_tolerance", 1.e-15 ) ),
       _absolute_residual_tolerance( input("linear-nonlinear-solver/absolute_residual_tolerance", 0.0 ) ),
       _initial_linear_tolerance( input("linear-nonlinear-solver/initial_linear_tolerance", 1.e-3 ) ),
       _minimum_linear_tolerance( input("linear-nonlinear-solver/minimum_linear_tolerance", 1.e-3 ) ),
       _max_linear_iterations( input("linear-nonlinear-solver/max_linear_iterations", 500 ) ),
       _continue_after_backtrack_failure( input("linear-nonlinear-solver/continue_after_backtrack_failure", false ) ),
       _continue_after_max_iterations( input("linear-nonlinear-solver/continue_after_max_iterations", false ) ),
       _solver_quiet( input("screen-options/solver_quiet", false ) ),
       _solver_verbose( input("screen-options/solver_verbose", false ) )
   {
     return;
   }

GRINS::Solver::~Solver ( )

virtual

Definition at line 62 of file grins_solver.C.

   {
     return;
   }

Member Function Documentation

virtual void GRINS::Solver::adjoint_qoi_parameter_sensitivity	(	SolverContext &	,
		const libMesh::QoISet &	,
		const libMesh::ParameterVector &	,
		libMesh::SensitivityData &
	)		const

inlinevirtual

Reimplemented in GRINS::SteadyMeshAdaptiveSolver, and GRINS::SteadySolver.

Definition at line 71 of file grins_solver.h.

76 { libmesh_not_implemented(); }

virtual void GRINS::Solver::forward_qoi_parameter_sensitivity	(	SolverContext &	,
		const libMesh::QoISet &	,
		const libMesh::ParameterVector &	,
		libMesh::SensitivityData &
	)		const

inlinevirtual

Reimplemented in GRINS::SteadyMeshAdaptiveSolver, and GRINS::SteadySolver.

Definition at line 79 of file grins_solver.h.

84 { libmesh_not_implemented(); }

virtual void GRINS::Solver::init_time_solver ( GRINS::MultiphysicsSystem * system )

protectedpure virtual

Implemented in GRINS::SteadyMeshAdaptiveSolver, GRINS::SteadySolver, and GRINS::UnsteadySolver.

Referenced by initialize().

void GRINS::Solver::initialize	(	const GetPot &	input,
		std::tr1::shared_ptr< libMesh::EquationSystems >	equation_system,
		GRINS::MultiphysicsSystem *	system
	)

virtual

Reimplemented in GRINS::DisplacementContinuationSolver.

Definition at line 67 of file grins_solver.C.

References init_time_solver(), and set_solver_options().

Referenced by GRINS::DisplacementContinuationSolver::initialize().

   {
  
     // Defined in subclasses depending on the solver used.
     this->init_time_solver(system);
 
     // Initialize the system
     equation_system->init();
 
     // Get diff solver to set options
     libMesh::DiffSolver &solver = *(system->time_solver->diff_solver().get());
 
     // Set linear/nonlinear solver options
     this->set_solver_options( solver );
 
     return;
   }

void GRINS::Solver::set_solver_options ( libMesh::DiffSolver & solver )

protected

Definition at line 87 of file grins_solver.C.

References _absolute_residual_tolerance, _absolute_step_tolerance, _continue_after_backtrack_failure, _continue_after_max_iterations, _initial_linear_tolerance, _max_linear_iterations, _max_nonlinear_iterations, _minimum_linear_tolerance, _relative_residual_tolerance, _relative_step_tolerance, _solver_quiet, and _solver_verbose.

Referenced by initialize().

   {
     solver.quiet                       = this->_solver_quiet;
     solver.verbose                     = this->_solver_verbose;
     solver.max_nonlinear_iterations    = this->_max_nonlinear_iterations;
     solver.relative_step_tolerance     = this->_relative_step_tolerance;
     solver.absolute_step_tolerance     = this->_absolute_step_tolerance;
     solver.relative_residual_tolerance = this->_relative_residual_tolerance;
     solver.absolute_residual_tolerance = this->_absolute_residual_tolerance;
     solver.max_linear_iterations       = this->_max_linear_iterations;
     solver.continue_after_backtrack_failure = this->_continue_after_backtrack_failure;
     solver.initial_linear_tolerance    = this->_initial_linear_tolerance;
     solver.minimum_linear_tolerance    = this->_minimum_linear_tolerance;
     solver.continue_after_max_iterations    = this->_continue_after_max_iterations;
 
     return;
   }

virtual void GRINS::Solver::solve ( SolverContext & context )

pure virtual

Implemented in GRINS::SteadyMeshAdaptiveSolver, GRINS::DisplacementContinuationSolver, GRINS::UnsteadySolver, and GRINS::SteadySolver.

void GRINS::Solver::steady_adjoint_solve ( SolverContext & context )

Do steady version of adjoint solve.

We put this here since we may want to reuse this in multiple different steady solves.

Definition at line 105 of file grins_solver.C.

References GRINS::SolverContext::system.

Referenced by GRINS::SteadySolver::solve(), and GRINS::SteadyMeshAdaptiveSolver::solve().

   {
     libMesh::out << "==========================================================" << std::endl
                  << "Solving adjoint problem." << std::endl
                  << "==========================================================" << std::endl;
 
     context.system->adjoint_solve();
     context.system->set_adjoint_already_solved(true);
   }