#include <grins_steady_solver.h>

Inheritance diagram for GRINS::SteadySolver:

Collaboration diagram for GRINS::SteadySolver:

Public Member Functions
	SteadySolver (const GetPot &input)

virtual	~SteadySolver ()

virtual void	solve (SolverContext &context)

virtual void	adjoint_qoi_parameter_sensitivity (SolverContext &context, const libMesh::QoISet &qoi_indices, const libMesh::ParameterVector &parameters_in, libMesh::SensitivityData &sensitivities) const

virtual void	forward_qoi_parameter_sensitivity (SolverContext &context, const libMesh::QoISet &qoi_indices, const libMesh::ParameterVector &parameters_in, libMesh::SensitivityData &sensitivities) const

Public Member Functions inherited from GRINS::Solver
	Solver (const GetPot &input)

virtual	~Solver ()

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

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

void	print_scalar_vars (SolverContext &context)

void	print_qoi (SolverContext &context, std::ostream &output)

Protected Member Functions
virtual void	init_time_solver (GRINS::MultiphysicsSystem *system)

Protected Member Functions inherited from GRINS::Solver
void	set_solver_options (libMesh::DiffSolver &solver)

Additional Inherited Members
Protected Attributes inherited from GRINS::Solver
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	_require_residual_reduction

bool	_solver_quiet

bool	_solver_verbose

Detailed Description

Definition at line 34 of file grins_steady_solver.h.

Constructor & Destructor Documentation

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

Definition at line 43 of file grins_steady_solver.C.

     : Solver( input )
   {
     return;
   }

GRINS::SteadySolver::~SteadySolver ( )

virtual

Definition at line 49 of file grins_steady_solver.C.

   {
     return;
   }

Member Function Documentation

void GRINS::SteadySolver::adjoint_qoi_parameter_sensitivity	(	SolverContext &	context,
		const libMesh::QoISet &	qoi_indices,
		const libMesh::ParameterVector &	parameters_in,
		libMesh::SensitivityData &	sensitivities
	)		const

virtual

Reimplemented from GRINS::Solver.

Definition at line 108 of file grins_steady_solver.C.

References GRINS::SolverContext::system.

   {
     // Get the linear solver
     libMesh::LinearSolver<libMesh::Number> *linear_solver = context.system->get_linear_solver();
 
     // Set ourselves to reuse the preconditioner
     linear_solver->reuse_preconditioner(true);
 
     context.system->adjoint_qoi_parameter_sensitivity
       (qoi_indices, parameters_in, sensitivities);
 
     // Go back to not reusing the preconditioner
     linear_solver->reuse_preconditioner(false);
   }

void GRINS::SteadySolver::forward_qoi_parameter_sensitivity	(	SolverContext &	context,
		const libMesh::QoISet &	qoi_indices,
		const libMesh::ParameterVector &	parameters_in,
		libMesh::SensitivityData &	sensitivities
	)		const

virtual

Reimplemented from GRINS::Solver.

Definition at line 127 of file grins_steady_solver.C.

References GRINS::SolverContext::equation_system, GRINS::SolverContext::output_residual_sensitivities, GRINS::SolverContext::output_solution_sensitivities, GRINS::SolverContext::system, and GRINS::SolverContext::vis.

   {
     context.system->forward_qoi_parameter_sensitivity
       (qoi_indices, parameters_in, sensitivities);
 
     if( context.output_residual_sensitivities )
       context.vis->output_residual_sensitivities
         ( context.equation_system, context.system, parameters_in );
 
     if( context.output_solution_sensitivities )
       context.vis->output_solution_sensitivities
         ( context.equation_system, context.system, parameters_in );
   }

void GRINS::SteadySolver::init_time_solver ( GRINS::MultiphysicsSystem * system )

protectedvirtual

Implements GRINS::Solver.

Definition at line 54 of file grins_steady_solver.C.

   {
     libMesh::SteadySolver* time_solver = new libMesh::SteadySolver( *(system) );
 
     system->time_solver = libMesh::UniquePtr<libMesh::TimeSolver>(time_solver);
     return;
   }

void GRINS::SteadySolver::solve ( SolverContext & context )

virtual

Implements GRINS::Solver.

Definition at line 62 of file grins_steady_solver.C.

References GRINS::SolverContext::do_adjoint_solve, GRINS::SolverContext::equation_system, GRINS::SolverContext::output_adjoint, GRINS::SolverContext::output_residual, GRINS::SolverContext::output_vis, GRINS::SolverContext::postprocessing, GRINS::Solver::print_scalar_vars(), GRINS::SolverContext::print_scalars, GRINS::Solver::steady_adjoint_solve(), GRINS::SolverContext::system, and GRINS::SolverContext::vis.

   {
     libmesh_assert( context.system );
 
     if( context.output_vis ) 
       {
         context.postprocessing->update_quantities( *(context.equation_system) );
         context.vis->output( context.equation_system );
       }
 
     // GRVY timers contained in here (if enabled)
     context.system->solve();
 
     if ( context.print_scalars )
       this->print_scalar_vars(context);
 
     if( context.do_adjoint_solve )
       {
         // Get the linear solver
         libMesh::LinearSolver<libMesh::Number> *linear_solver = context.system->get_linear_solver();
 
         // Set ourselves to reuse the preconditioner
         linear_solver->reuse_preconditioner(true);
 
         // Solve the adjoint problem
         this->steady_adjoint_solve(context);
 
         // Go back to not reusing the preconditioner
         linear_solver->reuse_preconditioner(false);
       }
 
     if( context.output_adjoint )
       context.vis->output_adjoint( context.equation_system, context.system );
 
     if( context.output_vis )
       {
         context.postprocessing->update_quantities( *(context.equation_system) );
         context.vis->output( context.equation_system );
       }
 
     if( context.output_residual ) context.vis->output_residual( context.equation_system, context.system );
 
     return;
   }

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

src/solver/include/grins/grins_steady_solver.h
src/solver/src/grins_steady_solver.C

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation