Class GeneralSolver provides an abstract interface for the implementation of linear or nonlinear solvers for systems of equations defined on an AMR patch hierarchy.

#include <ibtk/GeneralSolver.h>

Inheritance diagram for IBTK::GeneralSolver:

[legend]

Public Member Functions
	GeneralSolver ()=default
	Constructor. More...

virtual	~GeneralSolver ()=default
	Empty virtual destructor. More...

General-purpose solver functionality.
const std::string &	getName () const
	Return the object name. More...

virtual bool	getIsInitialized () const
	Return whether the operator is initialized. More...

virtual void	setHomogeneousBc (bool homogeneous_bc)
	Set whether the solver should use homogeneous boundary conditions. More...

virtual bool	getHomogeneousBc () const
	Return whether the solver is using homogeneous boundary conditions. More...

virtual void	setSolutionTime (double solution_time)
	Set the time at which the solution is to be evaluated. More...

virtual double	getSolutionTime () const
	Get the time at which the solution is being evaluated. More...

virtual void	setTimeInterval (double current_time, double new_time)
	Set the current time interval. More...

virtual std::pair< double, double >	getTimeInterval () const
	Get the current time interval. More...

virtual double	getDt () const
	Get the current time step size. More...

virtual void	setHierarchyMathOps (SAMRAI::tbox::Pointer< HierarchyMathOps > hier_math_ops)
	Set the HierarchyMathOps object used by the solver. More...

virtual SAMRAI::tbox::Pointer< HierarchyMathOps >	getHierarchyMathOps () const
	Get the HierarchyMathOps object used by the solver. More...

virtual bool	solveSystem (SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &b)=0
	Solve the system of equations. More...

virtual void	initializeSolverState (const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x, const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &b)
	Compute hierarchy dependent data required for solving \(F[x]=b\). More...

virtual void	deallocateSolverState ()
	Remove all hierarchy dependent data allocated by initializeSolverState(). More...

Functions to access solver parameters.
virtual void	setMaxIterations (int max_iterations)
	Set the maximum number of nonlinear iterations to use per solve. More...

virtual int	getMaxIterations () const
	Get the maximum number of nonlinear iterations to use per solve. More...

virtual void	setAbsoluteTolerance (double abs_residual_tol)
	Set the absolute residual tolerance for convergence. More...

virtual double	getAbsoluteTolerance () const
	Get the absolute residual tolerance for convergence. More...

virtual void	setRelativeTolerance (double rel_residual_tol)
	Set the relative residual tolerance for convergence. More...

virtual double	getRelativeTolerance () const
	Get the relative residual tolerance for convergence. More...

Logging functions.
std::string	d_object_name = "unitialized"

bool	d_is_initialized = false

bool	d_homogeneous_bc = false

double	d_solution_time = std::numeric_limits<double>::quiet_NaN()

double	d_current_time = std::numeric_limits<double>::quiet_NaN()

double	d_new_time = std::numeric_limits<double>::quiet_NaN()

double	d_rel_residual_tol = 0.0

double	d_abs_residual_tol = 0.0

int	d_max_iterations = 100

int	d_current_iterations = 0

double	d_current_residual_norm = std::numeric_limits<double>::quiet_NaN()

SAMRAI::tbox::Pointer< HierarchyMathOps >	d_hier_math_ops

bool	d_hier_math_ops_external = false

bool	d_enable_logging = false

virtual void	setLoggingEnabled (bool enable_logging=true)
	Enable or disable logging. More...

virtual bool	getLoggingEnabled () const
	Determine whether logging is enabled or disabled. More...

virtual void	printClassData (std::ostream &stream)
	Print class data to stream. More...

void	init (const std::string &object_name, bool homogeneous_bc)

virtual void	initSpecialized (const std::string &object_name, bool homogeneous_bc)

	GeneralSolver (const GeneralSolver &from)=delete
	Copy constructor. More...

GeneralSolver &	operator= (const GeneralSolver &that)=delete
	Assignment operator. More...

Constructor & Destructor Documentation

◆ GeneralSolver() [1/2]

IBTK::GeneralSolver::GeneralSolver ( )

default

◆ ~GeneralSolver()

virtual IBTK::GeneralSolver::~GeneralSolver ( )

virtualdefault

◆ GeneralSolver() [2/2]

IBTK::GeneralSolver::GeneralSolver ( const GeneralSolver & from )

privatedelete

Note: This constructor is not implemented and should not be used.

Parameters

from	The value to copy to this object.

Member Function Documentation

◆ getName()

const std::string& IBTK::GeneralSolver::getName ( ) const

◆ getIsInitialized()

virtual bool IBTK::GeneralSolver::getIsInitialized ( ) const

virtual

◆ setHomogeneousBc()

virtual void IBTK::GeneralSolver::setHomogeneousBc ( bool homogeneous_bc )

virtual

Reimplemented in IBTK::FACPreconditioner, IBTK::NewtonKrylovSolver, and IBTK::KrylovLinearSolver.

◆ getHomogeneousBc()

virtual bool IBTK::GeneralSolver::getHomogeneousBc ( ) const

virtual

◆ setSolutionTime()

virtual void IBTK::GeneralSolver::setSolutionTime ( double solution_time )

virtual

Reimplemented in IBAMR::CIBStaggeredStokesSolver, IBTK::FACPreconditioner, IBTK::NewtonKrylovSolver, and IBTK::KrylovLinearSolver.

◆ getSolutionTime()

virtual double IBTK::GeneralSolver::getSolutionTime ( ) const

virtual

◆ setTimeInterval()

virtual void IBTK::GeneralSolver::setTimeInterval	(	double	current_time,
		double	new_time
	)

virtual

Reimplemented in IBAMR::CIBStaggeredStokesSolver, IBTK::FACPreconditioner, IBTK::NewtonKrylovSolver, and IBTK::KrylovLinearSolver.

◆ getTimeInterval()

virtual std::pair<double, double> IBTK::GeneralSolver::getTimeInterval ( ) const

virtual

◆ getDt()

virtual double IBTK::GeneralSolver::getDt ( ) const

virtual

◆ setHierarchyMathOps()

virtual void IBTK::GeneralSolver::setHierarchyMathOps ( SAMRAI::tbox::Pointer< HierarchyMathOps > hier_math_ops )

virtual

Reimplemented in IBTK::NewtonKrylovSolver, and IBTK::KrylovLinearSolver.

◆ getHierarchyMathOps()

virtual SAMRAI::tbox::Pointer<HierarchyMathOps> IBTK::GeneralSolver::getHierarchyMathOps ( ) const

virtual

◆ solveSystem()

virtual bool IBTK::GeneralSolver::solveSystem	(	SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	x,
		SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	b
	)

pure virtual

Before calling solveSystem(), the form of the solution x and right-hand-side b vectors must be set properly by the user on all patch interiors on the specified range of levels in the patch hierarchy. The user is responsible for all data management for the quantities associated with the solution and right-hand-side vectors. In particular, patch data in these vectors must be allocated prior to calling this method.

Parameters

x	solution vector
b	right-hand-side vector

Conditions on Parameters:

vectors x and b must have same patch hierarchy
vectors x and b must have same structure, depth, etc.

Note: Subclasses must be implemented so that the vector arguments for solveSystem() need not match those for initializeSolverState(). However, they are allowed to require a certain degree of similarity, including:

hierarchy configuration (hierarchy pointer and range of levels)
number, type and alignment of vector component data
ghost cell widths of data in the solution x and right-hand-side b vectors

Note: Subclasses are required to be implemented so that the solver does not need to be initialized prior to calling solveSystem(); however, see initializeSolverState() and deallocateSolverState() for opportunities to save overhead when performing multiple consecutive solves.

See also: initializeSolverState; deallocateSolverState

Returns: true if the solver converged to the specified tolerances, false otherwise

Implemented in IBTK::PETScKrylovLinearSolver, IBTK::PETScNewtonKrylovSolver, IBTK::CCPoissonHypreLevelSolver, IBTK::SCPoissonHypreLevelSolver, IBTK::BGaussSeidelPreconditioner, IBTK::PETScLevelSolver, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBTK::PETScPCLSWrapper, IBAMR::StaggeredStokesProjectionPreconditioner, IBAMR::VCStaggeredStokesProjectionPreconditioner, IBAMR::CIBStaggeredStokesSolver, and IBAMR::IBImplicitStaggeredHierarchyIntegrator::IBImplicitStaggeredStokesSolver.

◆ initializeSolverState()

virtual void IBTK::GeneralSolver::initializeSolverState	(	const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	x,
		const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	b
	)

virtual

In a typical implementation, the solveSystem() method will compute some required hierarchy dependent data before the solve, and then remove that data after the solve. For multiple solves that use the same hierarchy configuration, it is more efficient to:

initialize the hierarchy-dependent data required by the solver via initializeSolverState(),
solve the system one or more times via solveSystem(), and
remove the hierarchy-dependent data via deallocateSolverState().

Note that it is generally necessary to reinitialize the solver state when the hierarchy configuration changes.

Parameters

x	solution vector
b	right-hand-side vector

Conditions on Parameters:

vectors x and b must have same patch hierarchy
vectors x and b must have same structure, depth, etc.

Note: Subclasses must be implemented so that the vector arguments for solveSystem() need not match those for initializeSolverState(). However, they are allowed to require a certain degree of similarity, including:

hierarchy configuration (hierarchy pointer and range of levels)
number, type and alignment of vector component data
ghost cell widths of data in the solution x and right-hand-side b vectors

Note: Subclasses are required to be implemented so that it is safe to call initializeSolverState() when the solver state is already initialized. In this case, the solver state should be first deallocated and then reinitialized.; Subclasses are required to be implemented so that when any operator objects have been registered with the solver via setOperator() or setJacobian(), they are also initialized by initializeSolverState().

See also: deallocateSolverState

Reimplemented in IBTK::PETScKrylovLinearSolver, IBTK::PETScNewtonKrylovSolver, IBTK::CCPoissonHypreLevelSolver, IBTK::SCPoissonHypreLevelSolver, IBTK::BGaussSeidelPreconditioner, IBTK::PETScLevelSolver, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, IBAMR::StaggeredStokesBlockPreconditioner, IBTK::PETScPCLSWrapper, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBAMR::StaggeredStokesProjectionPreconditioner, IBAMR::VCStaggeredStokesProjectionPreconditioner, and IBAMR::CIBStaggeredStokesSolver.

◆ deallocateSolverState()

virtual void IBTK::GeneralSolver::deallocateSolverState ( )

virtual

Note: Subclasses are required to be implemented so that it is safe to call deallocateSolverState() when the solver state is already deallocated.; Subclasses are required to be implemented so that when any operator objects have been registered with the solver via setOperator() or setJacobian(), they are also deallocated by deallocateSolverState().

See also: initializeSolverState

Reimplemented in IBTK::PETScKrylovLinearSolver, IBTK::PETScNewtonKrylovSolver, IBTK::CCPoissonHypreLevelSolver, IBTK::SCPoissonHypreLevelSolver, IBTK::BGaussSeidelPreconditioner, IBTK::PETScLevelSolver, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, IBAMR::StaggeredStokesBlockPreconditioner, IBTK::PETScPCLSWrapper, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBAMR::StaggeredStokesProjectionPreconditioner, IBAMR::VCStaggeredStokesProjectionPreconditioner, and IBAMR::CIBStaggeredStokesSolver.

◆ setMaxIterations()

virtual void IBTK::GeneralSolver::setMaxIterations ( int max_iterations )

virtual

Reimplemented in IBTK::BGaussSeidelPreconditioner, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, IBTK::PETScPCLSWrapper, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBAMR::StaggeredStokesProjectionPreconditioner, and IBAMR::VCStaggeredStokesProjectionPreconditioner.

◆ getMaxIterations()

virtual int IBTK::GeneralSolver::getMaxIterations ( ) const

virtual

Reimplemented in IBTK::PETScPCLSWrapper.

◆ setAbsoluteTolerance()

virtual void IBTK::GeneralSolver::setAbsoluteTolerance ( double abs_residual_tol )

virtual

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getAbsoluteTolerance()

virtual double IBTK::GeneralSolver::getAbsoluteTolerance ( ) const

virtual

Reimplemented in IBTK::PETScPCLSWrapper.

◆ setRelativeTolerance()

virtual void IBTK::GeneralSolver::setRelativeTolerance ( double rel_residual_tol )

virtual

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getRelativeTolerance()

virtual double IBTK::GeneralSolver::getRelativeTolerance ( ) const

virtual

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getNumIterations()

virtual int IBTK::GeneralSolver::getNumIterations ( ) const

virtual

Reimplemented in IBTK::BGaussSeidelPreconditioner, IBTK::PETScPCLSWrapper, and IBTK::BJacobiPreconditioner.

◆ getResidualNorm()

virtual double IBTK::GeneralSolver::getResidualNorm ( ) const

virtual

Reimplemented in IBTK::BGaussSeidelPreconditioner, IBTK::PETScPCLSWrapper, and IBTK::BJacobiPreconditioner.

◆ setLoggingEnabled()

virtual void IBTK::GeneralSolver::setLoggingEnabled ( bool enable_logging = true )

virtual

◆ getLoggingEnabled()

virtual bool IBTK::GeneralSolver::getLoggingEnabled ( ) const

virtual

◆ printClassData()

virtual void IBTK::GeneralSolver::printClassData ( std::ostream & stream )

virtual

Reimplemented in IBTK::LinearSolver.

◆ init()

void IBTK::GeneralSolver::init	(	const std::string &	object_name,
		bool	homogeneous_bc
	)

protected

◆ initSpecialized()

virtual void IBTK::GeneralSolver::initSpecialized	(	const std::string &	object_name,
		bool	homogeneous_bc
	)

protectedvirtual

Reimplemented in IBTK::PoissonSolver.

◆ operator=()

GeneralSolver& IBTK::GeneralSolver::operator= ( const GeneralSolver & that )

privatedelete

Note: This operator is not implemented and should not be used.

Parameters

that	The value to assign to this object.

Returns: A reference to this object.

Member Data Documentation

◆ d_object_name

std::string IBTK::GeneralSolver::d_object_name = "unitialized"

protected

◆ d_is_initialized

bool IBTK::GeneralSolver::d_is_initialized = false

protected

◆ d_homogeneous_bc

bool IBTK::GeneralSolver::d_homogeneous_bc = false

protected

◆ d_solution_time

double IBTK::GeneralSolver::d_solution_time = std::numeric_limits<double>::quiet_NaN()

protected

◆ d_current_time

double IBTK::GeneralSolver::d_current_time = std::numeric_limits<double>::quiet_NaN()

protected

◆ d_new_time

double IBTK::GeneralSolver::d_new_time = std::numeric_limits<double>::quiet_NaN()

protected

◆ d_rel_residual_tol

double IBTK::GeneralSolver::d_rel_residual_tol = 0.0

protected

◆ d_abs_residual_tol

double IBTK::GeneralSolver::d_abs_residual_tol = 0.0

protected

◆ d_max_iterations

int IBTK::GeneralSolver::d_max_iterations = 100

protected

◆ d_current_iterations

int IBTK::GeneralSolver::d_current_iterations = 0

protected

◆ d_current_residual_norm

double IBTK::GeneralSolver::d_current_residual_norm = std::numeric_limits<double>::quiet_NaN()

protected

◆ d_hier_math_ops

SAMRAI::tbox::Pointer<HierarchyMathOps> IBTK::GeneralSolver::d_hier_math_ops

protected

◆ d_hier_math_ops_external

bool IBTK::GeneralSolver::d_hier_math_ops_external = false

protected

◆ d_enable_logging

bool IBTK::GeneralSolver::d_enable_logging = false

protected

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

ibtk/include/ibtk/GeneralSolver.h

Public Member Functions

General-purpose solver functionality.

Functions to access solver parameters.

Logging functions.

Constructor & Destructor Documentation

◆ GeneralSolver() [1/2]

◆ ~GeneralSolver()

◆ GeneralSolver() [2/2]

Member Function Documentation

◆ getName()

◆ getIsInitialized()

◆ setHomogeneousBc()

◆ getHomogeneousBc()

◆ setSolutionTime()

◆ getSolutionTime()

◆ setTimeInterval()

◆ getTimeInterval()

◆ getDt()

◆ setHierarchyMathOps()

◆ getHierarchyMathOps()

◆ solveSystem()

◆ initializeSolverState()

◆ deallocateSolverState()

◆ setMaxIterations()

◆ getMaxIterations()

◆ setAbsoluteTolerance()

◆ getAbsoluteTolerance()

◆ setRelativeTolerance()

◆ getRelativeTolerance()

◆ getNumIterations()

◆ getResidualNorm()

◆ setLoggingEnabled()

◆ getLoggingEnabled()

◆ printClassData()

◆ init()

◆ initSpecialized()

◆ operator=()

Member Data Documentation

◆ d_object_name

◆ d_is_initialized

◆ d_homogeneous_bc

◆ d_solution_time

◆ d_current_time

◆ d_new_time

◆ d_rel_residual_tol

◆ d_abs_residual_tol

◆ d_max_iterations

◆ d_current_iterations

◆ d_current_residual_norm

◆ d_hier_math_ops

◆ d_hier_math_ops_external

◆ d_enable_logging