#include <cl_SolverPETSC.hpp>

Inheritance diagram for belfem::solver::PETSC:

Collaboration diagram for belfem::solver::PETSC:

Public Member Functions
	PETSC (const SolverParameters *aParams)
	~PETSC () override
void	set (const Preconditioner aPreconditioner, const KrylovMethod aKrylovMethod, const real aEpsilon=1e-8)
void	solve (SpMatrix &aMatrix, Vector< real > &aLHS, Vector< real > &aRHS) override
void	free () override
Public Member Functions inherited from belfem::solver::Wrapper
	Wrapper (const string &aLabel, const bool aUsesMPI)
virtual	~Wrapper ()
const string &	label () const
	returns the name of the solver as string
void	set_soft_fail (const bool aFlag)
	arm/disarm the soft-fail contract ( see the member note )
bool	failed () const
	true if the last factorization/solve failed softly
void	clear_failure ()
	clear a recorded soft failure before the next solve
virtual void	solve (SpMatrix &aMatrix, Matrix< real > &aLHS, Matrix< real > &aRHS)
bool	is_initialized () const
	tells if the initialitzation routine has already been called
proc_t	rank () const
	returns the communication rank of this proc
proc_t	comm_size () const
	returns the communication size of this proc
bool	uses_mpi () const
void	free () override
virtual real	get_determinant () const
	returns the determinant, if supported by the solver and computation was requested
virtual real	get_cond1 () const
	returns the conditioning numbner, if supported by the solver and computation was requested
virtual real	get_cond2 () const
	returns the conditioning numbner, if supported by the solver and computation was requested
virtual real	get_forward_error () const
	estimated FORWARD error of the last solve, \|\| dx \|\| / \|\| x \|\|, if supported by the solver and computation was requested.
virtual real	get_backward_error () const
	componentwise BACKWARD error of the last solve, if supported by the solver and computation was requested: the smallest relative perturbation of A and b for which the computed x is exact.
virtual real	get_omega2 () const
	the SECOND half of that backward error on its own, if supported.
virtual bool	supports_factorization_reuse () const
	true if this wrapper can solve against an EXISTING factorization rather than rebuilding it.
virtual void	freeze_factorization (const SpMatrix &aMatrix)
	arm the frozen state: subsequent solves reuse the factorization that the LAST SUCCESSFUL solve built, instead of rebuilding it.
virtual void	unfreeze_factorization ()
	drop the frozen state and return to ordinary re-factorizing solves.
virtual bool	factorization_is_frozen () const
	true while the frozen scope is armed
Public Member Functions inherited from belfem::CommunicationObject
	CommunicationObject ()
	~CommunicationObject ()
index_t	index () const

Protected Member Functions
void	initialize (SpMatrix &aMatrix, const SymmetryMode aSymmetryMode=SymmetryMode::Unsymmetric, const int_t aNumRhsColumns=1) override
Protected Member Functions inherited from belfem::solver::Wrapper
bool	soft_fail () const
	subclass hooks for the soft-fail contract
void	flag_failure ()
virtual void	initialize ()
void	hatch_turtle ()
	print a warning ( the turtle ) when the MPI ranks sharing a node ask for more OpenMP threads than the node has PHYSICAL cores.
void	mat2vec (const Matrix< real > &aM, Vector< real > &aV)
void	vec2mat (const Vector< real > &aV, Matrix< real > &aM)
SpMatrix *	matrix ()
Vector< real > &	x ()
Vector< real > &	y ()

Additional Inherited Members
Protected Attributes inherited from belfem::solver::Wrapper
SpMatrix *	mMatrix = nullptr
Vector< real > *	mX = nullptr
Vector< real > *	mY = nullptr

Constructor & Destructor Documentation

◆ PETSC()

belfem::solver::PETSC::PETSC ( const SolverParameters * aParams )

◆ ~PETSC()

belfem::solver::PETSC::~PETSC ( )

override

Member Function Documentation

◆ free()

void belfem::solver::PETSC::free ( )

overridevirtual

Reimplemented from belfem::CommunicationObject.

◆ initialize()

void belfem::solver::PETSC::initialize	(	SpMatrix &	aMatrix,
		const SymmetryMode	aSymmetryMode = SymmetryMode::Unsymmetric,
		const int_t	aNumRhsColumns = 1 )

overrideprotectedvirtual

Reimplemented from belfem::solver::Wrapper.

◆ set()

void belfem::solver::PETSC::set	(	const Preconditioner	aPreconditioner,
		const KrylovMethod	aKrylovMethod,
		const real	aEpsilon = 1e-8 )

◆ solve()

void belfem::solver::PETSC::solve	(	SpMatrix &	aMatrix,
		Vector< real > &	aLHS,
		Vector< real > &	aRHS )