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

#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 ()

virtual NM_Status	solve (SparseMtrx &A, FloatArray &b, FloatArray &x)
	Solves the given linear system by LDL^T factorization. More...

int	estimateMaxPackSize (IntArray &, DataStream &, int &)

void	setUpCommunicationMaps ()
	Sets up the communication maps. More...

virtual IRResultType	initializeFrom (InputRecord *ir)

virtual const char *	giveClassName () const

virtual LinSystSolverType	giveLinSystSolverType () const

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
	Returns the recommended sparse matrix type for this solver. More...

Public Member Functions inherited from oofem::SparseLinearSystemNM
	SparseLinearSystemNM (Domain d, EngngModel m)
	Constructor. More...

virtual	~SparseLinearSystemNM ()
	Destructor. More...

std::string	errorInfo (const char *func)

virtual NM_Status	solve (SparseMtrx &A, FloatMatrix &B, FloatMatrix &X)
	Solves the given sparse linear system of equations $A\cdot X=B$ . More...

Public Member Functions inherited from oofem::NumericalMethod
	NumericalMethod (Domain d, EngngModel m)
	Constructor. More...

virtual	~NumericalMethod ()
	Destructor. More...

EngngModel *	giveEngngModel ()

virtual void	reinitialize ()
	Reinitializes the receiver. More...

virtual void	setDomain (Domain *d)

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)

Private Attributes
int	nse

int	ni
	Max number of iterations. More...

double	err
	Max allowed error. More...

double	limit
	Linear dep./indep. trigger. More...

FloatMatrix	rbm
	Rigid body motions. More...

FloatMatrix	l
	Rigid body motions of all partitions. On master only. More...

IntArray	se
	Indices of singular equations. More...

IntArray	rbmAddr
	Addresses of initial partition contribution to rbm matrix. More...

FloatArray	w
	Primary unknowns. More...

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

FETICommunicator *	masterCommunicator

IntArray	masterCommMap
	List of local nodes (at master) participating in communication (list of boundary dof managers). More...

int	energyNorm_comput_flag
	Flag indicating computation of energy norm. More...

Additional Inherited Members
Protected Attributes inherited from oofem::NumericalMethod
Domain *	domain
	Pointer to domain. More...

EngngModel *	engngModel
	Pointer to engineering model. More...

Detailed Description

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

Definition at line 67 of file fetisolver.h.

Member Enumeration Documentation

anonymous enum

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

Definition at line 151 of file fetisolver.h.

Constructor & Destructor Documentation

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

Definition at line 49 of file fetisolver.C.

References energyNorm_comput_flag, err, and ni.

oofem::FETISolver::~FETISolver ( )

virtual

Definition at line 57 of file fetisolver.C.

References commBuff, oofem::NumericalMethod::engngModel, oofem::EngngModel::giveRank(), and masterCommunicator.

Member Function Documentation

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

Definition at line 67 of file fetisolver.C.

References oofem::NumericalMethod::domain, FETISOLVER_MAX_RBM, oofem::FETICommunicator::giveDofManager(), oofem::Domain::giveDofManager(), oofem::Domain::giveEngngModel(), oofem::FETIBoundaryDofManager::giveNumberOfDofs(), oofem::DataStream::givePackSizeOfDouble(), oofem::EngngModel::giveRank(), and masterCommunicator.

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

inlinevirtual

Implements oofem::SparseLinearSystemNM.

Definition at line 123 of file fetisolver.h.

virtual LinSystSolverType oofem::FETISolver::giveLinSystSolverType ( ) const

inlinevirtual

Returns: LinSystSolverType value, corresponding to receiver.

Implements oofem::SparseLinearSystemNM.

Definition at line 124 of file fetisolver.h.

References masterMapDirectionVector(), masterMapGammas(), masterMapPPVector(), masterMapQQProduct(), masterMapRBM(), masterMapResiduals(), masterMapSolution(), masterUnpackQQProduct(), masterUnpackRBM(), packDirectionVector(), packGammas(), packPPVector(), packQQProducts(), packRBM(), packResiduals(), packSolution(), projection(), oofem::ST_Feti, unpackDirectionVector(), unpackGammas(), unpackPPVector(), unpackResiduals(), and unpackSolution().

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

inlinevirtual

Returns the recommended sparse matrix type for this solver.

Implements oofem::SparseLinearSystemNM.

Definition at line 152 of file fetisolver.h.

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

IRResultType oofem::FETISolver::initializeFrom ( InputRecord * ir )

virtual

Reimplemented from oofem::NumericalMethod.

Definition at line 93 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, oofem::IRRT_OK, limit, and ni.

int oofem::FETISolver::masterMapDirectionVector ( )

int oofem::FETISolver::masterMapGammas ( )

Definition at line 1037 of file fetisolver.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), gamma, localGammas, nse, rbmAddr, and oofem::FloatArray::resize().

Referenced by giveLinSystSolverType(), and solve().

int oofem::FETISolver::masterMapPPVector ( )

int oofem::FETISolver::masterMapQQProduct ( )

Definition at line 432 of file fetisolver.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), nsem, q, qq, and rbmAddr.

Referenced by giveLinSystSolverType(), and solve().

int oofem::FETISolver::masterMapRBM ( )

int oofem::FETISolver::masterMapResiduals ( )

int oofem::FETISolver::masterMapSolution ( )

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

Definition at line 410 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

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

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

Definition at line 1005 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

Definition at line 867 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

Definition at line 395 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

Definition at line 249 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

Definition at line 589 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

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

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

Definition at line 221 of file fetisolver.C.

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

Referenced by giveLinSystSolverType(), and solve().

void oofem::FETISolver::setUpCommunicationMaps ( )

Sets up the communication maps.

Definition at line 118 of file fetisolver.C.

References oofem::IntArray::at(), oofem::FETICommunicator::BoundaryDofManagersRecMsg, commBuff, oofem::DofManager_shared, oofem::NumericalMethod::domain, oofem::DofManager::giveCompleteLocationArray(), oofem::Domain::giveDofManager(), oofem::NumericalMethod::giveEngngModel(), oofem::Domain::giveEngngModel(), oofem::DofManager::giveGlobalNumber(), oofem::Domain::giveNumberOfDofManagers(), oofem::CommunicationBuffer::givePackSizeOfInt(), oofem::DofManager::giveParallelMode(), oofem::EngngModel::giveRank(), oofem::IntArray::giveSize(), oofem::StaticCommunicationBuffer::init(), oofem::StaticCommunicationBuffer::iSend(), masterCommMap, oofem::FETICommunicator::NumberOfBoundaryDofManagersMsg, OOFEM_LOG_DEBUG, processCommunicator, oofem::IntArray::resize(), oofem::StaticCommunicationBuffer::resize(), oofem::ProcessCommunicator::setToRecvArry(), oofem::ProcessCommunicator::setToSendArry(), and oofem::StaticCommunicationBuffer::write().

Referenced by solve().

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

virtual

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