Class SCPoissonHypreLevelSolver is a concrete LinearSolver for solving elliptic equations of the form $ \mbox{$L u$} = \mbox{$(C I + \nabla \cdot D \nabla) u$} = f $ on a single SAMRAI::hier::PatchLevel using hypre. More...

#include <ibtk/SCPoissonHypreLevelSolver.h>

Inheritance diagram for IBTK::SCPoissonHypreLevelSolver:

[legend]

Public Member Functions
	SCPoissonHypreLevelSolver (const std::string &object_name, SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, const std::string &default_options_prefix)
	Constructor. More...

	~SCPoissonHypreLevelSolver ()
	Destructor. More...

virtual void	setPoissonSpecifications (const SAMRAI::solv::PoissonSpecifications &poisson_spec)
	Set the SAMRAI::solv::PoissonSpecifications object used to specify the coefficients for the scalar-valued or vector-valued Laplace operator. More...

virtual void	setPhysicalBcCoef (SAMRAI::solv::RobinBcCoefStrategy< NDIM > *bc_coef)
	Set the SAMRAI::solv::RobinBcCoefStrategy object used to specify physical boundary conditions. More...

virtual void	setPhysicalBcCoefs (const std::vector< SAMRAI::solv::RobinBcCoefStrategy< NDIM > * > &bc_coefs)
	Set the SAMRAI::solv::RobinBcCoefStrategy objects used to specify physical boundary conditions. More...

Static Public Member Functions
static SAMRAI::tbox::Pointer< PoissonSolver >	allocate_solver (const std::string &object_name, SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, const std::string &default_options_prefix)
	Static function to construct a SCPoissonHypreLevelSolver. More...

Protected Member Functions
void	initSpecialized (const std::string &object_name, bool homogeneous_bc) override

Protected Attributes
SAMRAI::solv::PoissonSpecifications	d_poisson_spec = SAMRAI::solv::PoissonSpecifications("")

std::unique_ptr< SAMRAI::solv::RobinBcCoefStrategy< NDIM > >	d_default_bc_coef

std::vector< SAMRAI::solv::RobinBcCoefStrategy< NDIM > * >	d_bc_coefs

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.
virtual void	setLoggingEnabled (bool enable_logging=true)
	Enable or disable logging. More...

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

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

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

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

Functions to access solver parameters.
virtual void	setInitialGuessNonzero (bool initial_guess_nonzero=true)
	Set whether the initial guess is non-zero. More...

virtual bool	getInitialGuessNonzero () const
	Get whether the initial guess is non-zero. More...

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

bool	d_initial_guess_nonzero = true

bool	d_nullspace_contains_constant_vec = false

std::vector< SAMRAI::tbox::Pointer< SAMRAI::solv::SAMRAIVectorReal< NDIM, double > > >	d_nullspace_basis_vecs

Linear solver functionality.
SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > >	d_hierarchy
	Associated hierarchy. More...

int	d_level_num = IBTK::invalid_level_number
	Associated patch level and C-F boundary (for level numbers > 0). More...

SAMRAI::tbox::Pointer< SAMRAI::hier::PatchLevel< NDIM > >	d_level

SAMRAI::tbox::Pointer< SAMRAI::hier::CoarseFineBoundary< NDIM > >	d_cf_boundary

bool	solveSystem (SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &b) override
	Solve the linear system of equations $Ax=b$ for $x$. More...

void	initializeSolverState (const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x, const SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &b) override
	Compute hierarchy dependent data required for solving $Ax=b$. More...

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

	SCPoissonHypreLevelSolver ()=delete
	Default constructor. More...

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

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

void	allocateHypreData ()
	Functions to allocate, initialize, access, and deallocate hypre data structures. More...

void	setMatrixCoefficients ()

void	setupHypreSolver ()

bool	solveSystem (int x_idx, int b_idx)

void	destroyHypreSolver ()

void	deallocateHypreData ()

hypre objects.
HYPRE_SStructGrid	d_grid = nullptr

HYPRE_SStructStencil	d_stencil [NVARS]

HYPRE_SStructGraph	d_graph = nullptr

HYPRE_SStructMatrix	d_matrix = nullptr

HYPRE_SStructVector	d_rhs_vec = nullptr

HYPRE_SStructVector	d_sol_vec = nullptr

HYPRE_SStructSolver	d_solver = nullptr

HYPRE_SStructSolver	d_precond = nullptr

std::vector< SAMRAI::hier::Index< NDIM > >	d_stencil_offsets

std::string	d_solver_type = "Split"

std::string	d_precond_type = "none"

std::string	d_split_solver_type = "PFMG"

int	d_rel_change = 0

int	d_num_pre_relax_steps = 1

int	d_num_post_relax_steps = 1

int	d_relax_type

int	d_skip_relax = 1

int	d_two_norm = 1

static const int	PART = 0

static const int	NPARTS = 1

static const int	NVARS = NDIM

static const int	X_VAR = 0

static const int	Y_VAR = 1

static const int	Z_VAR = 2

Detailed Description

This solver class uses the hypre library to solve linear equations of the form $ (C I + \nabla \cdot D \nabla ) u = f $, where $C$ and $D$ are constants, and $u$ and $f$ are side-centered arrays. The discretization is second-order accurate.

Robin boundary conditions may be specified through the interface class SAMRAI::solv::RobinBcCoefStrategy, but boundary conditions must be of either Dirichlet or Neumann type. In particular, mixed boundary conditions are not presently supported.

The user must perform the following steps to use class SCPoissonHypreLevelSolver:

Create a SCPoissonHypreLevelSolver object.
Set the problem specification via setPoissonSpecifications(), setPhysicalBcCoef(), and setHomogeneousBc().
Initialize SCPoissonHypreLevelSolver object using the function initializeSolverState().
Solve the linear system using the member function solveSystem(), passing in SAMRAI::solv::SAMRAIVectorReal objects corresponding to $u$ and $f$.

Sample parameters for initialization from database (and their default values):

enable_logging = FALSE         // see setLoggingEnabled()
solver_type = "Split"          // choices are: "Split", "SysPFMG", "PCG", "GMRES", "FlexGMRES"
,
"LGMRES", "BiCGSTAB"
precond_type = "none"          // choices are: "Split", "SysPFMG"
split_solver_type = "PFMG"     // choices are: "PFMG", "SMG", "Jacobi"
max_iterations = 25            // see setMaxIterations()
abs_residual_tol = 1.e-50      // see setAbsoluteTolerance() (only used by hypre Krylov
solvers)
rel_residual_tol = 1.0e-5      // see setRelativeTolerance()
initial_guess_nonzero = FALSE  // see setInitialGuessNonzero()
rel_change = 0                 // see hypre User's Manual (only used by SysPFMG or PCG solver)
num_pre_relax_steps = 1        // number of pre-sweeps (only used by SysPFMG solver)
num_post_relax_steps = 1       // number of post-sweeps (only used by SysPFMG solver)
relax_type = 1                 // see hypre User's Manual (only used by SysPFMG solver or
preconditioner)
skip_relax = 1                 // see hypre User's Manual (only used by SysPFMG solver or
preconditioner)
two_norm = 1                   // see hypre User's Manual (only used by PCG solver)

hypre is developed in the Center for Applied Scientific Computing (CASC) at Lawrence Livermore National Laboratory (LLNL). For more information about hypre, see https://computation.llnl.gov/casc/linear_solvers/sls_hypre.html.

Constructor & Destructor Documentation

◆ SCPoissonHypreLevelSolver() [1/3]

IBTK::SCPoissonHypreLevelSolver::SCPoissonHypreLevelSolver	(	const std::string &	object_name,
		SAMRAI::tbox::Pointer< SAMRAI::tbox::Database >	input_db,
		const std::string &	default_options_prefix
	)

◆ ~SCPoissonHypreLevelSolver()

IBTK::SCPoissonHypreLevelSolver::~SCPoissonHypreLevelSolver ( )

◆ SCPoissonHypreLevelSolver() [2/3]

IBTK::SCPoissonHypreLevelSolver::SCPoissonHypreLevelSolver ( )

privatedelete

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

◆ SCPoissonHypreLevelSolver() [3/3]

IBTK::SCPoissonHypreLevelSolver::SCPoissonHypreLevelSolver ( const SCPoissonHypreLevelSolver & 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

◆ allocate_solver()

static SAMRAI::tbox::Pointer<PoissonSolver> IBTK::SCPoissonHypreLevelSolver::allocate_solver	(	const std::string &	object_name,
		SAMRAI::tbox::Pointer< SAMRAI::tbox::Database >	input_db,
		const std::string &	default_options_prefix
	)

inlinestatic

◆ solveSystem() [1/2]

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

overridevirtual

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: The vector arguments for solveSystem() need not match those for initializeSolverState(). However, there must be 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: The solver need not 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

Implements IBTK::GeneralSolver.

◆ initializeSolverState()

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

overridevirtual

By default, the solveSystem() method computes some required hierarchy dependent data before solving and removes 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: The vector arguments for solveSystem() need not match those for initializeSolverState(). However, there must be 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: It is safe to call initializeSolverState() when the state is already initialized. In this case, the solver state is first deallocated and then reinitialized.

See also: deallocateSolverState

Reimplemented from IBTK::GeneralSolver.

◆ deallocateSolverState()

void IBTK::SCPoissonHypreLevelSolver::deallocateSolverState ( )

overridevirtual

Note: It is safe to call deallocateSolverState() when the solver state is already deallocated.

See also: initializeSolverState

Reimplemented from IBTK::GeneralSolver.

◆ operator=()

SCPoissonHypreLevelSolver& IBTK::SCPoissonHypreLevelSolver::operator= ( const SCPoissonHypreLevelSolver & 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.

◆ allocateHypreData()

void IBTK::SCPoissonHypreLevelSolver::allocateHypreData ( )

private

◆ setMatrixCoefficients()

void IBTK::SCPoissonHypreLevelSolver::setMatrixCoefficients ( )

private

◆ setupHypreSolver()

void IBTK::SCPoissonHypreLevelSolver::setupHypreSolver ( )

private

◆ solveSystem() [2/2]

bool IBTK::SCPoissonHypreLevelSolver::solveSystem	(	int	x_idx,
		int	b_idx
	)

private

◆ destroyHypreSolver()

void IBTK::SCPoissonHypreLevelSolver::destroyHypreSolver ( )

private

◆ deallocateHypreData()

void IBTK::SCPoissonHypreLevelSolver::deallocateHypreData ( )

private

◆ setNullSpace()

virtual void IBTK::LinearSolver::setNullSpace	(	bool	nullspace_contains_constant_vec,
		const std::vector< SAMRAI::tbox::Pointer< SAMRAI::solv::SAMRAIVectorReal< NDIM, double > > > &	nullspace_basis_vecs = `std::vector< SAMRAI::tbox::Pointer< SAMRAI::solv::SAMRAIVectorReal< NDIM, double > > >()`
	)

virtualinherited

Implementations can require the nullspace basis vectors to be orthogonal but should not assume the basis vectors to be orthonormal. If the basis vectors are not orthonormal, the solver may normalize them in place.

Reimplemented in IBTK::PETScKrylovLinearSolver, and IBTK::PETScLevelSolver.

◆ getNullSpaceContainsConstantVector()

virtual bool IBTK::LinearSolver::getNullSpaceContainsConstantVector ( ) const

virtualinherited

◆ getNullSpaceBasisVectors()

virtual const std::vector<SAMRAI::tbox::Pointer<SAMRAI::solv::SAMRAIVectorReal<NDIM, double> > >& IBTK::LinearSolver::getNullSpaceBasisVectors ( ) const

virtualinherited

◆ setInitialGuessNonzero()

virtual void IBTK::LinearSolver::setInitialGuessNonzero ( bool initial_guess_nonzero = true )

virtualinherited

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

◆ getInitialGuessNonzero()

virtual bool IBTK::LinearSolver::getInitialGuessNonzero ( ) const

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ printClassData()

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

overridevirtualinherited

Reimplemented from IBTK::GeneralSolver.

◆ getName()

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

inherited

◆ getIsInitialized()

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

virtualinherited

◆ setHomogeneousBc()

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

virtualinherited

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

◆ getHomogeneousBc()

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

virtualinherited

◆ setSolutionTime()

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

virtualinherited

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

◆ getSolutionTime()

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

virtualinherited

◆ setTimeInterval()

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

virtualinherited

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

◆ getTimeInterval()

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

virtualinherited

◆ getDt()

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

virtualinherited

◆ setHierarchyMathOps()

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

virtualinherited

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

◆ getHierarchyMathOps()

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

virtualinherited

◆ setMaxIterations()

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

virtualinherited

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

◆ getMaxIterations()

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

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ setAbsoluteTolerance()

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

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getAbsoluteTolerance()

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

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ setRelativeTolerance()

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

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getRelativeTolerance()

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

virtualinherited

Reimplemented in IBTK::PETScPCLSWrapper.

◆ getNumIterations()

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

virtualinherited

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

◆ getResidualNorm()

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

virtualinherited

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

◆ setLoggingEnabled()

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

virtualinherited

◆ getLoggingEnabled()

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

virtualinherited

◆ init()

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

protectedinherited

◆ initSpecialized() [1/2]

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

protectedvirtualinherited

Reimplemented in IBTK::PoissonSolver.

◆ setPoissonSpecifications()

virtual void IBTK::PoissonSolver::setPoissonSpecifications ( const SAMRAI::solv::PoissonSpecifications & poisson_spec )

virtualinherited

Reimplemented in IBTK::PoissonFACPreconditioner, and IBTK::KrylovLinearSolverPoissonSolverInterface.

◆ setPhysicalBcCoef()

virtual void IBTK::PoissonSolver::setPhysicalBcCoef ( SAMRAI::solv::RobinBcCoefStrategy< NDIM > * bc_coef )

virtualinherited

Note: bc_coef may be nullptr. In this case, default boundary conditions (as supplied to the class constructor) are employed.

Parameters

bc_coef Pointer to an object that can set the Robin boundary condition coefficients

Reimplemented in IBTK::PoissonFACPreconditioner, and IBTK::KrylovLinearSolverPoissonSolverInterface.

◆ setPhysicalBcCoefs()

virtual void IBTK::PoissonSolver::setPhysicalBcCoefs ( const std::vector< SAMRAI::solv::RobinBcCoefStrategy< NDIM > * > & bc_coefs )

virtualinherited

Note: Any of the elements of bc_coefs may be nullptr. In this case, default boundary conditions (as supplied to the class constructor) are employed for that data depth.

Parameters

bc_coefs Vector of pointers to objects that can set the Robin boundary condition coefficients

Reimplemented in IBTK::PoissonFACPreconditioner, and IBTK::KrylovLinearSolverPoissonSolverInterface.

◆ initSpecialized() [2/2]

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

overrideprotectedvirtualinherited

Reimplemented from IBTK::GeneralSolver.

Member Data Documentation

◆ d_hierarchy

SAMRAI::tbox::Pointer<SAMRAI::hier::PatchHierarchy<NDIM> > IBTK::SCPoissonHypreLevelSolver::d_hierarchy

private

◆ d_level_num

int IBTK::SCPoissonHypreLevelSolver::d_level_num = IBTK::invalid_level_number

private

◆ d_level

SAMRAI::tbox::Pointer<SAMRAI::hier::PatchLevel<NDIM> > IBTK::SCPoissonHypreLevelSolver::d_level

private

◆ d_cf_boundary

SAMRAI::tbox::Pointer<SAMRAI::hier::CoarseFineBoundary<NDIM> > IBTK::SCPoissonHypreLevelSolver::d_cf_boundary

private

◆ PART

const int IBTK::SCPoissonHypreLevelSolver::PART = 0

staticprivate

◆ NPARTS

const int IBTK::SCPoissonHypreLevelSolver::NPARTS = 1

staticprivate

◆ NVARS

const int IBTK::SCPoissonHypreLevelSolver::NVARS = NDIM

staticprivate

◆ X_VAR

const int IBTK::SCPoissonHypreLevelSolver::X_VAR = 0

staticprivate

◆ Y_VAR

const int IBTK::SCPoissonHypreLevelSolver::Y_VAR = 1

staticprivate

◆ Z_VAR

const int IBTK::SCPoissonHypreLevelSolver::Z_VAR = 2

staticprivate

◆ d_grid

HYPRE_SStructGrid IBTK::SCPoissonHypreLevelSolver::d_grid = nullptr

private

◆ d_stencil

HYPRE_SStructStencil IBTK::SCPoissonHypreLevelSolver::d_stencil[NVARS]

private

◆ d_graph

HYPRE_SStructGraph IBTK::SCPoissonHypreLevelSolver::d_graph = nullptr

private

◆ d_matrix

HYPRE_SStructMatrix IBTK::SCPoissonHypreLevelSolver::d_matrix = nullptr

private

◆ d_rhs_vec

HYPRE_SStructVector IBTK::SCPoissonHypreLevelSolver::d_rhs_vec = nullptr

private

◆ d_sol_vec

HYPRE_SStructVector IBTK::SCPoissonHypreLevelSolver::d_sol_vec = nullptr

private

◆ d_solver

HYPRE_SStructSolver IBTK::SCPoissonHypreLevelSolver::d_solver = nullptr

private

◆ d_precond

HYPRE_SStructSolver IBTK::SCPoissonHypreLevelSolver::d_precond = nullptr

private

◆ d_stencil_offsets

std::vector<SAMRAI::hier::Index<NDIM> > IBTK::SCPoissonHypreLevelSolver::d_stencil_offsets

private

◆ d_solver_type

std::string IBTK::SCPoissonHypreLevelSolver::d_solver_type = "Split"

private

◆ d_precond_type

std::string IBTK::SCPoissonHypreLevelSolver::d_precond_type = "none"

private

◆ d_split_solver_type

std::string IBTK::SCPoissonHypreLevelSolver::d_split_solver_type = "PFMG"

private

◆ d_rel_change

int IBTK::SCPoissonHypreLevelSolver::d_rel_change = 0

private

◆ d_num_pre_relax_steps

int IBTK::SCPoissonHypreLevelSolver::d_num_pre_relax_steps = 1

private

◆ d_num_post_relax_steps

int IBTK::SCPoissonHypreLevelSolver::d_num_post_relax_steps = 1

private

◆ d_relax_type

int IBTK::SCPoissonHypreLevelSolver::d_relax_type

private

◆ d_skip_relax

int IBTK::SCPoissonHypreLevelSolver::d_skip_relax = 1

private

◆ d_two_norm

int IBTK::SCPoissonHypreLevelSolver::d_two_norm = 1

private

◆ d_initial_guess_nonzero

bool IBTK::LinearSolver::d_initial_guess_nonzero = true

protectedinherited

◆ d_nullspace_contains_constant_vec

bool IBTK::LinearSolver::d_nullspace_contains_constant_vec = false

protectedinherited

◆ d_nullspace_basis_vecs

std::vector<SAMRAI::tbox::Pointer<SAMRAI::solv::SAMRAIVectorReal<NDIM, double> > > IBTK::LinearSolver::d_nullspace_basis_vecs

protectedinherited

◆ d_object_name

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

protectedinherited

◆ d_is_initialized

bool IBTK::GeneralSolver::d_is_initialized = false

protectedinherited

◆ d_homogeneous_bc

bool IBTK::GeneralSolver::d_homogeneous_bc = false

protectedinherited

◆ d_solution_time

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

protectedinherited

◆ d_current_time

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

protectedinherited

◆ d_new_time

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

protectedinherited

◆ d_rel_residual_tol

double IBTK::GeneralSolver::d_rel_residual_tol = 0.0

protectedinherited

◆ d_abs_residual_tol

double IBTK::GeneralSolver::d_abs_residual_tol = 0.0

protectedinherited

◆ d_max_iterations

int IBTK::GeneralSolver::d_max_iterations = 100

protectedinherited

◆ d_current_iterations

int IBTK::GeneralSolver::d_current_iterations = 0

protectedinherited

◆ d_current_residual_norm

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

protectedinherited

◆ d_hier_math_ops

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

protectedinherited

◆ d_hier_math_ops_external

bool IBTK::GeneralSolver::d_hier_math_ops_external = false

protectedinherited

◆ d_enable_logging

bool IBTK::GeneralSolver::d_enable_logging = false

protectedinherited

◆ d_poisson_spec

SAMRAI::solv::PoissonSpecifications IBTK::PoissonSolver::d_poisson_spec = SAMRAI::solv::PoissonSpecifications("")

protectedinherited

◆ d_default_bc_coef

std::unique_ptr<SAMRAI::solv::RobinBcCoefStrategy<NDIM> > IBTK::PoissonSolver::d_default_bc_coef

protectedinherited

◆ d_bc_coefs

std::vector<SAMRAI::solv::RobinBcCoefStrategy<NDIM>*> IBTK::PoissonSolver::d_bc_coefs

protectedinherited

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

ibtk/include/ibtk/SCPoissonHypreLevelSolver.h

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Functions to access solver parameters.

Logging functions.

Functions to access solver parameters.

Logging functions.

Linear solver functionality.

hypre objects.

Detailed Description

Constructor & Destructor Documentation

◆ SCPoissonHypreLevelSolver() [1/3]

◆ ~SCPoissonHypreLevelSolver()

◆ SCPoissonHypreLevelSolver() [2/3]

◆ SCPoissonHypreLevelSolver() [3/3]

Member Function Documentation

◆ allocate_solver()

◆ solveSystem() [1/2]

◆ initializeSolverState()

◆ deallocateSolverState()

◆ operator=()

◆ allocateHypreData()

◆ setMatrixCoefficients()

◆ setupHypreSolver()

◆ solveSystem() [2/2]

◆ destroyHypreSolver()

◆ deallocateHypreData()

◆ setNullSpace()

◆ getNullSpaceContainsConstantVector()

◆ getNullSpaceBasisVectors()

◆ setInitialGuessNonzero()

◆ getInitialGuessNonzero()

◆ printClassData()

◆ getName()

◆ getIsInitialized()

◆ setHomogeneousBc()

◆ getHomogeneousBc()

◆ setSolutionTime()

◆ getSolutionTime()

◆ setTimeInterval()

◆ getTimeInterval()

◆ getDt()

◆ setHierarchyMathOps()

◆ getHierarchyMathOps()

◆ setMaxIterations()

◆ getMaxIterations()

◆ setAbsoluteTolerance()

◆ getAbsoluteTolerance()

◆ setRelativeTolerance()

◆ getRelativeTolerance()

◆ getNumIterations()

◆ getResidualNorm()

◆ setLoggingEnabled()

◆ getLoggingEnabled()

◆ init()

◆ initSpecialized() [1/2]

◆ setPoissonSpecifications()

◆ setPhysicalBcCoef()

◆ setPhysicalBcCoefs()

◆ initSpecialized() [2/2]

Member Data Documentation

◆ d_hierarchy

◆ d_level_num

◆ d_level

◆ d_cf_boundary

◆ PART

◆ NPARTS

◆ NVARS

◆ X_VAR

◆ Y_VAR

◆ Z_VAR

◆ d_grid

◆ d_stencil

◆ d_graph

◆ d_matrix

◆ d_rhs_vec

◆ d_sol_vec

◆ d_solver

◆ d_precond