Class PETScLevelSolver is an abstract LinearSolver for solving systems of linear equations on a single SAMRAI::hier::PatchLevel using PETSc. More...

#include <ibtk/PETScLevelSolver.h>

Inheritance diagram for IBTK::PETScLevelSolver:

[legend]

Public Member Functions
	PETScLevelSolver ()
	Default constructor. More...

	~PETScLevelSolver ()
	Destructor. More...

void	setKSPType (const std::string &ksp_type)
	Set the KSP type. More...

void	setOptionsPrefix (const std::string &options_prefix)
	Set the options prefix used by this PETSc solver object. More...

const KSP &	getPETScKSP () const
	Get the PETSc KSP object. More...

void	getASMSubdomains (std::vector< IS > nonoverlapping_subdomains, std::vector< IS > overlapping_subdomains)
	Get ASM subdomains. More...

Linear solver functionality.
virtual bool	getNullSpaceContainsConstantVector () const
	Get whether the nullspace of the linear system contains th constant vector. More...

virtual const std::vector< SAMRAI::tbox::Pointer< SAMRAI::solv::SAMRAIVectorReal< NDIM, double > > > &	getNullSpaceBasisVectors () const
	Get the basis vectors for the nullspace of the linear system. More...

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...

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...

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...

Functions to access data on the most recent solve.
virtual int	getNumIterations () const
	Return the iteration count from the most recent solve. More...

virtual double	getResidualNorm () const
	Return the residual norm from the most recent iteration. More...

Protected Attributes
PETSc objects.
std::string	d_ksp_type = KSPGMRES

std::string	d_pc_type = PCILU

std::string	d_shell_pc_type

std::string	d_options_prefix

KSP	d_petsc_ksp = nullptr

Mat	d_petsc_mat = nullptr

Mat	d_petsc_pc = nullptr

MatNullSpace	d_petsc_nullsp

Vec	d_petsc_x = nullptr

Vec	d_petsc_b = nullptr

Support for additive and multiplicative Schwarz preconditioners.
Vec	d_local_x

Vec	d_local_y

SAMRAI::hier::IntVector< NDIM >	d_box_size

SAMRAI::hier::IntVector< NDIM >	d_overlap_size

int	d_n_local_subdomains

int	d_n_subdomains_max

std::vector< IS >	d_overlap_is

std::vector< IS >	d_nonoverlap_is

std::vector< IS >	d_local_overlap_is

std::vector< IS >	d_local_nonoverlap_is

std::vector< VecScatter >	d_restriction

std::vector< VecScatter >	d_prolongation

std::vector< KSP >	d_sub_ksp

Mat *	d_sub_mat

Mat *	d_sub_bc_mat

std::vector< Vec >	d_sub_x

std::vector< Vec >	d_sub_y

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

SAMRAIDataCache	d_cached_eulerian_data
	Scratch data. More...

void	setNullSpace (bool 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 >>>()) override
	Set the nullspace of the linear system. More...

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...

void	init (SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, const std::string &default_options_prefix)
	Basic initialization. More...

virtual void	generateASMSubdomains (std::vector< std::set< int >> &overlap_is, std::vector< std::set< int >> &nonoverlap_is)
	Generate IS/subdomains for Schwartz type preconditioners. More...

virtual void	generateFieldSplitSubdomains (std::vector< std::string > &field_names, std::vector< std::set< int >> &field_is)
	Generate IS/subdomains for fieldsplit type preconditioners. More...

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

virtual void	deallocateSolverStateSpecialized ()=0
	Remove all hierarchy dependent data allocated by initializeSolverStateSpecialized(). More...

virtual void	copyToPETScVec (Vec &petsc_x, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x)=0
	Copy a generic vector to the PETSc representation. More...

virtual void	copyFromPETScVec (Vec &petsc_x, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x)=0
	Copy a generic vector from the PETSc representation. More...

virtual void	setupKSPVecs (Vec &petsc_x, Vec &petsc_b, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &x, SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &b)=0
	Copy solution and right-hand-side data to the PETSc representation, including any modifications to account for boundary conditions. More...

virtual void	setupNullSpace ()
	Setup the solver nullspace (if any). More...

Field split preconditioner.
std::vector< std::string >	d_field_name

std::vector< IS >	d_field_is

	PETScLevelSolver (const PETScLevelSolver &from)
	Copy constructor. More...

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

static PetscErrorCode	PCApply_Additive (PC pc, Vec x, Vec y)
	Apply the preconditioner to x and store the result in y. More...

static PetscErrorCode	PCApply_Multiplicative (PC pc, Vec x, Vec y)
	Apply the preconditioner to x and store the result in y. More...

static PetscErrorCode	PCApply_RedBlackMultiplicative (PC pc, Vec x, Vec y)
	Apply the preconditioner to x and store the result in y. More...

Logging functions.
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

virtual void	printClassData (std::ostream &stream) override
	Print class data to stream. 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...

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

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

Detailed Description

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

options_prefix = ""           // see setOptionsPrefix()
ksp_type = "gmres"            // see setKSPType()
initial_guess_nonzero = TRUE  // see setInitialGuessNonzero()
rel_residual_tol = 1.0e-5     // see setRelativeTolerance()
abs_residual_tol = 1.0e-50    // see setAbsoluteTolerance()
max_iterations = 10000        // see setMaxIterations()
enable_logging = FALSE        // see setLoggingEnabled()

PETSc is developed at the Argonne National Laboratory Mathematics and Computer Science Division. For more information about PETSc, see http://www.mcs.anl.gov/petsc.

Constructor & Destructor Documentation

◆ PETScLevelSolver() [1/2]

IBTK::PETScLevelSolver::PETScLevelSolver ( )

◆ ~PETScLevelSolver()

IBTK::PETScLevelSolver::~PETScLevelSolver ( )

◆ PETScLevelSolver() [2/2]

IBTK::PETScLevelSolver::PETScLevelSolver ( const PETScLevelSolver & from )

private

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

Parameters

from	The value to copy to this object.

Member Function Documentation

◆ setKSPType()

void IBTK::PETScLevelSolver::setKSPType ( const std::string & ksp_type )

◆ setOptionsPrefix()

void IBTK::PETScLevelSolver::setOptionsPrefix ( const std::string & options_prefix )

◆ getPETScKSP()

const KSP& IBTK::PETScLevelSolver::getPETScKSP ( ) const

◆ getASMSubdomains()

void IBTK::PETScLevelSolver::getASMSubdomains	(	std::vector< IS > **	nonoverlapping_subdomains,
		std::vector< IS > **	overlapping_subdomains
	)

◆ setNullSpace()

void IBTK::PETScLevelSolver::setNullSpace	(	bool	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 >>>()`
	)

overridevirtual

Reimplemented from IBTK::LinearSolver.

◆ solveSystem()

bool IBTK::PETScLevelSolver::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::PETScLevelSolver::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.; Subclasses of class PETScLevelSolver should not override this method. Instead, they should override the protected method initializeSolverStateSpecialized().

See also: deallocateSolverState

Reimplemented from IBTK::GeneralSolver.

◆ deallocateSolverState()

void IBTK::PETScLevelSolver::deallocateSolverState ( )

overridevirtual

Note: It is safe to call deallocateSolverState() when the solver state is already deallocated.; Subclasses of class PETScLevelSolver should not override this method. Instead, they should override the protected method deallocatedSolverStateSpecialized().

See also: initializeSolverState

Reimplemented from IBTK::GeneralSolver.

◆ init() [1/2]

void IBTK::PETScLevelSolver::init	(	SAMRAI::tbox::Pointer< SAMRAI::tbox::Database >	input_db,
		const std::string &	default_options_prefix
	)

protected

◆ generateASMSubdomains()

virtual void IBTK::PETScLevelSolver::generateASMSubdomains	(	std::vector< std::set< int >> &	overlap_is,
		std::vector< std::set< int >> &	nonoverlap_is
	)

protectedvirtual

Reimplemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ generateFieldSplitSubdomains()

virtual void IBTK::PETScLevelSolver::generateFieldSplitSubdomains	(	std::vector< std::string > &	field_names,
		std::vector< std::set< int >> &	field_is
	)

protectedvirtual

Reimplemented in IBAMR::StaggeredStokesPETScLevelSolver.

◆ initializeSolverStateSpecialized()

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

protectedpure virtual

Implemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::VCSCViscousPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ deallocateSolverStateSpecialized()

virtual void IBTK::PETScLevelSolver::deallocateSolverStateSpecialized ( )

protectedpure virtual

Implemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ copyToPETScVec()

virtual void IBTK::PETScLevelSolver::copyToPETScVec	(	Vec &	petsc_x,
		SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	x
	)

protectedpure virtual

Implemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ copyFromPETScVec()

virtual void IBTK::PETScLevelSolver::copyFromPETScVec	(	Vec &	petsc_x,
		SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	x
	)

protectedpure virtual

Implemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ setupKSPVecs()

virtual void IBTK::PETScLevelSolver::setupKSPVecs	(	Vec &	petsc_x,
		Vec &	petsc_b,
		SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	x,
		SAMRAI::solv::SAMRAIVectorReal< NDIM, double > &	b
	)

protectedpure virtual

Implemented in IBAMR::StaggeredStokesPETScLevelSolver, IBTK::VCSCViscousPETScLevelSolver, IBTK::SCPoissonPETScLevelSolver, and IBTK::CCPoissonPETScLevelSolver.

◆ setupNullSpace()

virtual void IBTK::PETScLevelSolver::setupNullSpace ( )

protectedvirtual

◆ operator=()

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

◆ PCApply_Additive()

static PetscErrorCode IBTK::PETScLevelSolver::PCApply_Additive	(	PC	pc,
		Vec	x,
		Vec	y
	)

staticprivate

◆ PCApply_Multiplicative()

static PetscErrorCode IBTK::PETScLevelSolver::PCApply_Multiplicative	(	PC	pc,
		Vec	x,
		Vec	y
	)

staticprivate

◆ PCApply_RedBlackMultiplicative()

static PetscErrorCode IBTK::PETScLevelSolver::PCApply_RedBlackMultiplicative	(	PC	pc,
		Vec	x,
		Vec	y
	)

staticprivate

◆ 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 IBAMR::VCStaggeredStokesProjectionPreconditioner, IBAMR::StaggeredStokesProjectionPreconditioner, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBTK::PETScPCLSWrapper, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, and IBTK::BGaussSeidelPreconditioner.

◆ 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::NewtonKrylovSolver, IBTK::KrylovLinearSolver, and IBTK::FACPreconditioner.

◆ getHomogeneousBc()

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

virtualinherited

◆ setSolutionTime()

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

virtualinherited

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

◆ 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::NewtonKrylovSolver, IBTK::KrylovLinearSolver, and IBTK::FACPreconditioner.

◆ 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 IBAMR::VCStaggeredStokesProjectionPreconditioner, IBAMR::StaggeredStokesProjectionPreconditioner, IBAMR::StaggeredStokesBlockFactorizationPreconditioner, IBTK::PETScPCLSWrapper, IBTK::FACPreconditioner, IBTK::BJacobiPreconditioner, and IBTK::BGaussSeidelPreconditioner.

◆ 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::PETScPCLSWrapper, IBTK::BJacobiPreconditioner, and IBTK::BGaussSeidelPreconditioner.

◆ getResidualNorm()

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

virtualinherited

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

◆ setLoggingEnabled()

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

virtualinherited

◆ getLoggingEnabled()

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

virtualinherited

◆ init() [2/2]

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

protectedinherited

◆ initSpecialized()

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

protectedvirtualinherited

Reimplemented in IBTK::PoissonSolver.

Member Data Documentation

◆ d_hierarchy

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

protected

◆ d_level_num

int IBTK::PETScLevelSolver::d_level_num = IBTK::invalid_level_number

protected

◆ d_level

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

protected

◆ d_cf_boundary

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

protected

◆ d_cached_eulerian_data

SAMRAIDataCache IBTK::PETScLevelSolver::d_cached_eulerian_data

protected

◆ d_ksp_type

std::string IBTK::PETScLevelSolver::d_ksp_type = KSPGMRES

protected

◆ d_pc_type

std::string IBTK::PETScLevelSolver::d_pc_type = PCILU

protected

◆ d_shell_pc_type

std::string IBTK::PETScLevelSolver::d_shell_pc_type

protected

◆ d_options_prefix

std::string IBTK::PETScLevelSolver::d_options_prefix

protected

◆ d_petsc_ksp

KSP IBTK::PETScLevelSolver::d_petsc_ksp = nullptr

protected

◆ d_petsc_mat

Mat IBTK::PETScLevelSolver::d_petsc_mat = nullptr

protected

◆ d_petsc_pc

Mat IBTK::PETScLevelSolver::d_petsc_pc = nullptr

protected

◆ d_petsc_nullsp

MatNullSpace IBTK::PETScLevelSolver::d_petsc_nullsp

protected

◆ d_petsc_x

Vec IBTK::PETScLevelSolver::d_petsc_x = nullptr

protected

◆ d_petsc_b

Vec IBTK::PETScLevelSolver::d_petsc_b = nullptr

protected

◆ d_local_x

Vec IBTK::PETScLevelSolver::d_local_x

protected

◆ d_local_y

Vec IBTK::PETScLevelSolver::d_local_y

protected

◆ d_box_size

SAMRAI::hier::IntVector<NDIM> IBTK::PETScLevelSolver::d_box_size

protected

◆ d_overlap_size

SAMRAI::hier::IntVector<NDIM> IBTK::PETScLevelSolver::d_overlap_size

protected

◆ d_n_local_subdomains

int IBTK::PETScLevelSolver::d_n_local_subdomains

protected

◆ d_n_subdomains_max

int IBTK::PETScLevelSolver::d_n_subdomains_max

protected

◆ d_overlap_is

std::vector<IS> IBTK::PETScLevelSolver::d_overlap_is

protected

◆ d_nonoverlap_is

std::vector<IS> IBTK::PETScLevelSolver::d_nonoverlap_is

protected

◆ d_local_overlap_is

std::vector<IS> IBTK::PETScLevelSolver::d_local_overlap_is

protected

◆ d_local_nonoverlap_is

std::vector<IS> IBTK::PETScLevelSolver::d_local_nonoverlap_is

protected

◆ d_restriction

std::vector<VecScatter> IBTK::PETScLevelSolver::d_restriction

protected

◆ d_prolongation

std::vector<VecScatter> IBTK::PETScLevelSolver::d_prolongation

protected

◆ d_sub_ksp

std::vector<KSP> IBTK::PETScLevelSolver::d_sub_ksp

protected

◆ d_sub_mat

Mat* IBTK::PETScLevelSolver::d_sub_mat

protected

◆ d_sub_bc_mat

Mat * IBTK::PETScLevelSolver::d_sub_bc_mat

protected

◆ d_sub_x

std::vector<Vec> IBTK::PETScLevelSolver::d_sub_x

protected

◆ d_sub_y

std::vector<Vec> IBTK::PETScLevelSolver::d_sub_y

protected

◆ d_field_name

std::vector<std::string> IBTK::PETScLevelSolver::d_field_name

protected

◆ d_field_is

std::vector<IS> IBTK::PETScLevelSolver::d_field_is

protected

◆ 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

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

ibtk/include/ibtk/PETScLevelSolver.h

Public Member Functions

Protected Attributes

Linear solver functionality.

Field split preconditioner.

Logging functions.

Logging functions.

Detailed Description

Constructor & Destructor Documentation

◆ PETScLevelSolver() [1/2]

◆ ~PETScLevelSolver()

◆ PETScLevelSolver() [2/2]

Member Function Documentation

◆ setKSPType()

◆ setOptionsPrefix()

◆ getPETScKSP()

◆ getASMSubdomains()

◆ setNullSpace()

◆ solveSystem()

◆ initializeSolverState()

◆ deallocateSolverState()

◆ init() [1/2]

◆ generateASMSubdomains()

◆ generateFieldSplitSubdomains()

◆ initializeSolverStateSpecialized()

◆ deallocateSolverStateSpecialized()

◆ copyToPETScVec()

◆ copyFromPETScVec()

◆ setupKSPVecs()

◆ setupNullSpace()

◆ operator=()

◆ PCApply_Additive()

◆ PCApply_Multiplicative()

◆ PCApply_RedBlackMultiplicative()

◆ 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() [2/2]

◆ initSpecialized()

Member Data Documentation

◆ d_hierarchy

◆ d_level_num

◆ d_level

◆ d_cf_boundary

◆ d_cached_eulerian_data

◆ d_ksp_type

◆ d_pc_type

◆ d_shell_pc_type

◆ d_options_prefix

◆ d_petsc_ksp

◆ d_petsc_mat

◆ d_petsc_pc

◆ d_petsc_nullsp

◆ d_petsc_x

◆ d_petsc_b

◆ d_local_x

◆ d_local_y

◆ d_box_size