#include <fetisolver.h>

Inheritance diagram for oofem::FETISolver:

Collaboration diagram for oofem::FETISolver:

Public Types
enum	{ FETISolverZeroTag , NumberOfRBMMsg , RBMMessage , QQMessage , SolutionMessage , ResidualMessage , DirectionVectorMessage , PPVectorMessage , GammasMessage , FETISolverIterationContinue , FETISolverIterationBreak }

Public Member Functions
	FETISolver (Domain d, EngngModel m)
virtual	~FETISolver ()
ConvergedReason	solve (SparseMtrx &A, FloatArray &b, FloatArray &x) override
int	estimateMaxPackSize (IntArray &, DataStream &, int &)
void	setUpCommunicationMaps ()
	Sets up the communication maps.
void	initializeFrom (InputRecord &ir) override
const char *	giveClassName () const override
LinSystSolverType	giveLinSystSolverType () const override
void	projection (FloatArray &v, FloatMatrix &l, FloatMatrix &l1)
int	packRBM (ProcessCommunicator &processComm)
int	masterUnpackRBM (ProcessCommunicator &processComm)
int	packQQProducts (ProcessCommunicator &processComm)
int	masterUnpackQQProduct (ProcessCommunicator &processComm)
int	packSolution (ProcessCommunicator &processComm)
int	unpackSolution (ProcessCommunicator &processComm)
int	packResiduals (ProcessCommunicator &processComm)
int	unpackResiduals (ProcessCommunicator &processComm)
int	packDirectionVector (ProcessCommunicator &processComm)
int	unpackDirectionVector (ProcessCommunicator &processComm)
int	packPPVector (ProcessCommunicator &processComm)
int	unpackPPVector (ProcessCommunicator &processComm)
int	packGammas (ProcessCommunicator &processComm)
int	unpackGammas (ProcessCommunicator &processComm)
int	masterMapRBM ()
int	masterMapQQProduct ()
int	masterMapSolution ()
int	masterMapResiduals ()
int	masterMapDirectionVector ()
int	masterMapPPVector ()
int	masterMapGammas ()
virtual SparseMtrxType	giveRecommendedMatrix (bool symmetric) const
Public Member Functions inherited from oofem::SparseLinearSystemNM
	SparseLinearSystemNM (Domain d, EngngModel m)
	Constructor.
virtual	~SparseLinearSystemNM ()
	Destructor.
std::string	errorInfo (const char *func)
virtual ConvergedReason	solve (SparseMtrx &A, FloatMatrix &B, FloatMatrix &X)
Public Member Functions inherited from oofem::NumericalMethod
	NumericalMethod (Domain d, EngngModel m)
virtual	~NumericalMethod ()
	Destructor.
EngngModel *	giveEngngModel ()
virtual void	reinitialize ()
virtual void	setDomain (Domain *d)
virtual void	saveContext (DataStream &stream, ContextMode mode)
virtual void	restoreContext (DataStream &stream, ContextMode mode)

Private Attributes
int	nse
int	ni
	Max number of iterations.
double	err
	Max allowed error.
double	limit
	Linear dep./indep. trigger.
FloatMatrix	rbm
	Rigid body motions.
FloatMatrix	l
	Rigid body motions of all partitions. On master only.
IntArray	se
	Indices of singular equations.
IntArray	rbmAddr
	Addresses of initial partition contribution to rbm matrix.
FloatArray	w
	Primary unknowns.
FloatArray	qq
FloatArray	q
FloatArray	dd
FloatArray	g
FloatArray	d
FloatArray	p
FloatArray	pp
FloatArray	gamma
FloatArray	localGammas
IntArray	nsem
ProcessCommunicatorBuff	pcbuff
ProcessCommunicator	processCommunicator
CommunicatorBuff *	commBuff
	Common Communicator buffer.
FETICommunicator *	masterCommunicator
IntArray	masterCommMap
	List of local nodes (at master) participating in communication (list of boundary dof managers).
int	energyNorm_comput_flag
	Flag indicating computation of energy norm.

Additional Inherited Members
Protected Attributes inherited from oofem::NumericalMethod
Domain *	domain
	Pointer to domain.
EngngModel *	engngModel
	Pointer to engineering model.

Detailed Description

This class implements the class NumericalMethod instance FETI linear algebraic equation parallel solver.

Definition at line 68 of file fetisolver.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

Enumerator
FETISolverZeroTag
NumberOfRBMMsg
RBMMessage
QQMessage
SolutionMessage
ResidualMessage
DirectionVectorMessage
PPVectorMessage
GammasMessage
FETISolverIterationContinue
FETISolverIterationBreak

Definition at line 152 of file fetisolver.h.

Constructor & Destructor Documentation

◆ FETISolver()

oofem::FETISolver::FETISolver	(	Domain *	d,
		EngngModel *	m )

Definition at line 49 of file fetisolver.C.

References oofem::CBT_static, d, energyNorm_comput_flag, err, ni, pcbuff, processCommunicator, and oofem::SparseLinearSystemNM::SparseLinearSystemNM().

◆ ~FETISolver()

oofem::FETISolver::~FETISolver ( )

virtual

Definition at line 59 of file fetisolver.C.

References commBuff, oofem::NumericalMethod::engngModel, and masterCommunicator.

Member Function Documentation

◆ estimateMaxPackSize()

int oofem::FETISolver::estimateMaxPackSize	(	IntArray &	map,
		DataStream &	buff,
		int &	packUnpackType )

Definition at line 69 of file fetisolver.C.

References oofem::NumericalMethod::domain, FETISOLVER_MAX_RBM, oofem::DataStream::givePackSizeOfDouble(), and masterCommunicator.

◆ giveClassName()

const char * oofem::FETISolver::giveClassName ( ) const

inlineoverridevirtual

Implements oofem::SparseLinearSystemNM.

Definition at line 124 of file fetisolver.h.

◆ giveLinSystSolverType()

LinSystSolverType oofem::FETISolver::giveLinSystSolverType ( ) const

inlineoverridevirtual

Returns: LinSystSolverType value, corresponding to receiver.

Implements oofem::SparseLinearSystemNM.

Definition at line 125 of file fetisolver.h.

References oofem::ST_Feti.

◆ giveRecommendedMatrix()

virtual SparseMtrxType oofem::FETISolver::giveRecommendedMatrix ( bool symmetric ) const

inlinevirtual

Returns the recommended sparse matrix type for this solver.

Implements oofem::SparseLinearSystemNM.

Definition at line 153 of file fetisolver.h.

References oofem::SMT_Skyline, and oofem::SMT_SkylineU.

◆ initializeFrom()

void oofem::FETISolver::initializeFrom ( InputRecord & ir )

overridevirtual

Reimplemented from oofem::NumericalMethod.

Definition at line 95 of file fetisolver.C.

References _IFT_FETISolver_energynormflag, _IFT_FETISolver_limit, _IFT_FETISolver_maxerr, _IFT_FETISolver_maxiter, energyNorm_comput_flag, err, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, limit, and ni.

◆ masterMapDirectionVector()

int oofem::FETISolver::masterMapDirectionVector ( )

Definition at line 774 of file fetisolver.C.

References oofem::IntArray::at(), d, dd, oofem::NumericalMethod::domain, oofem::IntArray::giveSize(), masterCommMap, masterCommunicator, and OOFEM_ERROR.

Referenced by solve().

◆ masterMapGammas()

int oofem::FETISolver::masterMapGammas ( )

Definition at line 995 of file fetisolver.C.

References gamma, localGammas, nse, and rbmAddr.

Referenced by solve().

◆ masterMapPPVector()

int oofem::FETISolver::masterMapPPVector ( )

Definition at line 903 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, oofem::IntArray::giveSize(), masterCommMap, masterCommunicator, OOFEM_ERROR, p, and pp.

Referenced by solve().

◆ masterMapQQProduct()

int oofem::FETISolver::masterMapQQProduct ( )

Definition at line 422 of file fetisolver.C.

References nsem, q, qq, and rbmAddr.

Referenced by solve().

◆ masterMapRBM()

int oofem::FETISolver::masterMapRBM ( )

Definition at line 322 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, oofem::IntArray::giveSize(), l, masterCommMap, masterCommunicator, nsem, OOFEM_ERROR, rbm, and rbmAddr.

Referenced by solve().

◆ masterMapResiduals()

int oofem::FETISolver::masterMapResiduals ( )

Definition at line 641 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, g, oofem::IntArray::giveSize(), masterCommMap, masterCommunicator, OOFEM_ERROR, and pp.

Referenced by solve().

◆ masterMapSolution()

int oofem::FETISolver::masterMapSolution ( )

Definition at line 511 of file fetisolver.C.

References oofem::IntArray::at(), dd, oofem::NumericalMethod::domain, oofem::IntArray::giveSize(), masterCommMap, masterCommunicator, OOFEM_ERROR, and w.

Referenced by solve().

◆ masterUnpackQQProduct()

int oofem::FETISolver::masterUnpackQQProduct ( ProcessCommunicator & processComm )

Definition at line 401 of file fetisolver.C.

References oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicator::giveRank(), oofem::ProcessCommunicatorBuff::giveRecvBuff(), nsem, q, rbmAddr, and oofem::CommunicationBuffer::read().

◆ masterUnpackRBM()

int oofem::FETISolver::masterUnpackRBM ( ProcessCommunicator & processComm )

Definition at line 277 of file fetisolver.C.

References oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicator::giveRank(), oofem::ProcessCommunicatorBuff::giveRecvBuff(), oofem::ProcessCommunicator::giveToRecvMap(), l, masterCommunicator, nsem, rbmAddr, and oofem::CommunicationBuffer::read().

◆ packDirectionVector()

int oofem::FETISolver::packDirectionVector ( ProcessCommunicator & processComm )

Definition at line 703 of file fetisolver.C.

References oofem::IntArray::at(), d, oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicator::giveRank(), oofem::ProcessCommunicatorBuff::giveSendBuff(), oofem::ProcessCommunicator::giveToSendMap(), masterCommunicator, and oofem::CommunicationBuffer::write().

◆ packGammas()

int oofem::FETISolver::packGammas ( ProcessCommunicator & processComm )

Definition at line 963 of file fetisolver.C.

References gamma, oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicator::giveRank(), oofem::ProcessCommunicatorBuff::giveSendBuff(), nsem, rbmAddr, and oofem::CommunicationBuffer::write().

◆ packPPVector()

int oofem::FETISolver::packPPVector ( ProcessCommunicator & processComm )

Definition at line 832 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicatorBuff::giveSendBuff(), oofem::IntArray::giveSize(), oofem::ProcessCommunicator::giveToSendMap(), pp, and oofem::CommunicationBuffer::write().

◆ packQQProducts()

int oofem::FETISolver::packQQProducts ( ProcessCommunicator & processComm )

Definition at line 386 of file fetisolver.C.

References oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicatorBuff::giveSendBuff(), nse, qq, and oofem::CommunicationBuffer::write().

◆ packRBM()

int oofem::FETISolver::packRBM ( ProcessCommunicator & processComm )

Definition at line 247 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicatorBuff::giveSendBuff(), oofem::IntArray::giveSize(), oofem::ProcessCommunicator::giveToSendMap(), nse, rbm, and oofem::CommunicationBuffer::write().

◆ packResiduals()

int oofem::FETISolver::packResiduals ( ProcessCommunicator & processComm )

Definition at line 570 of file fetisolver.C.

References oofem::IntArray::at(), oofem::NumericalMethod::domain, oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicatorBuff::giveSendBuff(), oofem::IntArray::giveSize(), oofem::ProcessCommunicator::giveToSendMap(), pp, and oofem::CommunicationBuffer::write().

◆ packSolution()

int oofem::FETISolver::packSolution ( ProcessCommunicator & processComm )

Definition at line 436 of file fetisolver.C.

References oofem::IntArray::at(), oofem::ProcessCommunicator::giveProcessCommunicatorBuff(), oofem::ProcessCommunicator::giveRank(), oofem::ProcessCommunicatorBuff::giveSendBuff(), oofem::ProcessCommunicator::giveToSendMap(), masterCommunicator, w, and oofem::CommunicationBuffer::write().

◆ projection()

void oofem::FETISolver::projection	(	FloatArray &	v,
		FloatMatrix &	l,
		FloatMatrix &	l1 )

Definition at line 219 of file fetisolver.C.

References oofem::FloatArray::beProductOf(), oofem::FloatArray::beTProductOf(), l, and oofem::FloatArray::subtract().

Referenced by solve().

◆ setUpCommunicationMaps()

void oofem::FETISolver::setUpCommunicationMaps ( )

Sets up the communication maps.

Definition at line 116 of file fetisolver.C.

References oofem::IntArray::at(), commBuff, oofem::DofManager_shared, oofem::NumericalMethod::domain, oofem::NumericalMethod::giveEngngModel(), oofem::IntArray::giveSize(), masterCommMap, OOFEM_LOG_DEBUG, processCommunicator, and oofem::IntArray::resize().

Referenced by solve().

◆ solve()

ConvergedReason oofem::FETISolver::solve	(	SparseMtrx &	A,
		FloatArray &	b,
		FloatArray &	x )

overridevirtual

Solves the given linear system by LDL^T factorization. Implementation rely on factorization support provided by mapped sparse matrix. It calls Lhs->factorized()->backSubstitutionWith(*solutionArray). Sets solved flag to 1 if o.k.

Parameters