Profile matrix system of equations. More...

#include <ProfileSPDLinSOE.h>

Inheritance diagram for XC::ProfileSPDLinSOE:

Public Member Functions
virtual int	setSize (Graph &theGraph)
	Sets the system size. More...

virtual int	addA (const Matrix &, const ID &, double fact=1.0)
	Assembles the product of m by fact into A. More...

virtual void	zeroA (void)
	Zeros the entries in the 1d array for \(A\) and marks the system as not having been factored. More...

virtual int	sendSelf (CommParameters &)
	Send the object. More...

virtual int	recvSelf (const CommParameters &)
	Receive the object. More...

Public Member Functions inherited from XC::LinearSOEData
virtual int	getNumEqn (void) const
	Return the current size of the system.

virtual void	zeroB (void)
	Zeros the entries in the 1d array for \(b\).

virtual void	zeroX (void)
	Zeros the entries in the 1d array for \(x\).

virtual void	zero (void)
	Zeros the entries in the 1d array for \(x\) and \(b\).

virtual int	addB (const Vector &, const double &fact=1.0)
	add v multiplied by factor to vector B. More...

virtual int	addB (const Vector &, const ID &, const double &fact=1.0)
	assemblies v multiplied by factor into vector B. More...

virtual int	setB (const XC::Vector &v, const double &fact=1.0)
	Sets B= v*fact. More...

virtual void	setX (int loc, double value)
	If `loc` is within the range of \(x\), sets \(x(loc) = value\).

virtual void	setX (const Vector &x)
	Assigns \(x\).

virtual const Vector &	getX (void) const
	Returns the Vector object created for \(x\).

virtual Vector &	getX (void)
	Returns the Vector object created for \(x\).

virtual const Vector &	getB (void) const
	Returns the Vector object created for \(b\).

virtual Vector &	getB (void)
	Returns the Vector object created for \(b\).

virtual double	normRHS (void) const
	Returns the 2-norm of the vector \(x\).

void	receiveB (const CommParameters &)

void	receiveX (const CommParameters &)

void	receiveBX (const CommParameters &)

void	sendB (CommParameters &) const

void	sendBX (CommParameters &) const

Public Member Functions inherited from XC::LinearSOE
virtual	~LinearSOE (void)
	Destructor.

virtual int	solve (void)
	Computes the solution of the system of equations. More...

virtual double	getDeterminant (void)
	Returns the determinant of the system matrix.

LinearSOESolver *	getSolver (void)
	Returns a pointer to the solver.

LinearSOESolver &	newSolver (const std::string &)

Public Member Functions inherited from XC::SystemOfEqn
int	checkSize (Graph &theGraph) const
	Check number of DOFs in the graph.

Public Member Functions inherited from XC::MovableObject
	MovableObject (int classTag, int dbTag)
	Constructor. More...

	MovableObject (int classTag)
	Constructor. More...

	MovableObject (const MovableObject &)
	Copy constructor. Doesn't copy the dbTag.

MovableObject &	operator= (const MovableObject &)
	Assignment operator. Doesn't copy the dbTag.

int	getClassTag (void) const
	Return the class identifier.

int	getDbTag (void) const
	Return the object identifier in the database.

int	getDbTag (CommParameters &)
	Return the object identifier in the database.

void	setDbTag (int dbTag)
	Sets the object identifier in the database.

void	setDbTag (CommParameters &)
	Sets the object identifier in the database if not already set. More...

virtual int	setParameter (const std::vector< std::string > &argv, Parameter &param)
	Sets the value `param` to the parameter `argv`.

virtual int	updateParameter (int parameterID, Information &info)
	Updates the parameter identified by `parameterID` with `info`.

virtual int	activateParameter (int parameterID)
	Activates the parameter identified by `parameterID`.

virtual int	setVariable (const std::string &variable, Information &)
	Set the value of the variable idenfified by `var`.

virtual int	getVariable (const std::string &variable, Information &)
	Return the value of the variable idenfified by `var`.

Public Member Functions inherited from XC::DistributedBase
	DistributedBase (void)
	Constructor.

virtual DbTagData &	getDbTagData (void) const
	Returns a vector to store class dbTags.

const int &	getDbTagDataPos (const int &i) const
	Returns the data at the i-th position.

void	setDbTagDataPos (const int &i, const int &v)
	Sets the data at the i-th position.

void	inicComm (const int &dataSize) const
	Initializes communication.

Protected Member Functions
virtual bool	setSolver (LinearSOESolver *)
	Set the solver to use. More...

	ProfileSPDLinSOE (AnalysisAggregation *)
	Constructor. More...

	ProfileSPDLinSOE (AnalysisAggregation *, int classTag)
	Constructor. More...

	ProfileSPDLinSOE (AnalysisAggregation , int N, int iLoc, ProfileSPDLinSolver *theSolver=nullptr)
	Constructor. More...

SystemOfEqn *	getCopy (void) const
	Virtual constructor.

Protected Member Functions inherited from XC::FactoredSOEBase
	FactoredSOEBase (AnalysisAggregation *, int classTag, int N=0)
	Constructor. More...

Protected Member Functions inherited from XC::LinearSOEData
void	inic (const size_t &)
	1d arrays containing coefficients of B and X More...

const double &	getB (const size_t &i) const

double &	getB (const size_t &i)

const double &	getX (const size_t &i) const

double &	getX (const size_t &i)

double *	getPtrB (void)

const double *	getPtrB (void) const

double *	getPtrX (void)

const double *	getPtrX (void) const

	LinearSOEData (AnalysisAggregation *, int classTag, int N=0)
	Constructor. More...

Protected Member Functions inherited from XC::LinearSOE
int	setSolverSize (void)
	invoke setSize() on the Solver

	LinearSOE (AnalysisAggregation *, int classTag)
	Constructor. More...

Protected Member Functions inherited from XC::SystemOfEqn
virtual AnalysisModel *	getAnalysisModelPtr (void)
	Returns a pointer to the analysis model.

virtual const AnalysisModel *	getAnalysisModelPtr (void) const
	Returns a const pointer to the analysis model.

	SystemOfEqn (AnalysisAggregation *, int classTag)
	Constructor. More...

Protected Attributes
int	profileSize

Vector	A

ID	iDiagLoc

bool	isAcondensed

int	numInt

Protected Attributes inherited from XC::FactoredSOEBase
bool	factored
	True if the system is factored.

Protected Attributes inherited from XC::LinearSOEData
int	size

Vector	B
	order of A

Vector	X

Friends
class	AnalysisAggregation

class	FEM_ObjectBroker

class	ProfileSPDLinSolver

class	ProfileSPDLinDirectSolver

class	ProfileSPDLinDirectBlockSolver

class	ProfileSPDLinDirectThreadSolver

class	ProfileSPDLinDirectSkypackSolver

class	ProfileSPDLinSubstrSolver

class	ProfileSPDLinSubstrThreadSolver

Detailed Description

Profile matrix system of equations.

ProfileSPDLinSOE is class which is used to store a symmetric system of equations using a profile storage scheme. The upper triangular part of \(A\) is stored in a 1d double array with the diagonals of \(A\) located at positions given by an integer array \(iLoc\). For example when \(n=5\) and \(A\) as shown below:

\(\) [ {array}{ccccc} a_{0,0} & a_{0,1} & 0 & 0 & a_{0,4} a_{1,0} & a_{1,1} & a_{1,2} & a_{1,3} & 0 a_{2,0} & a_{2,1} & a_{2,2} & a_{2,3} & a_{2,4} 0 & a_{3,1} & a_{3,2} & a_{3,3} & a_{3,4} 0 & 0 & a_{4,2} & a_{4,3} & a_{4,4} {array} ] \(\)

this is stored using:

\(\) A = [ {array}{cccccccccccccccccccc} a_{0,0} & a_{0,1} & a_{1,1} & a_{1,2} & a_{2,2} & a_{1,3} & a_{2,3} & a_{3,3} & a_{0,4} & 0 & a_{2,4} & a_{3,4} & a_{4,4} {array} ] \(\)

and

\(\) iLoc = [ {array}{cccccccccccccccccccc} 1 & 3 & 5 & 8 & 13 {array} ] \(\) Note \(iLoc\) stores the diagonal locations using Fortran indexing. This is to facilitate calls to Fortran libraries, e.g. Digital's DXML. The \(x\) and \(b\) vectors are stored in 1d double arrays of length \(N\).

Constructor & Destructor Documentation

§ ProfileSPDLinSOE() [1/3]

XC::ProfileSPDLinSOE::ProfileSPDLinSOE ( AnalysisAggregation * owr )

protected

Constructor.

Parameters

owr	analysis aggregation that owns this object.

§ ProfileSPDLinSOE() [2/3]

XC::ProfileSPDLinSOE::ProfileSPDLinSOE	(	AnalysisAggregation *	owr,
		int	classTag
	)

protected

Constructor.

Parameters

owr	analysis aggregation that owns this object.
classTag	identifier of the class.

§ ProfileSPDLinSOE() [3/3]

XC::ProfileSPDLinSOE::ProfileSPDLinSOE	(	AnalysisAggregation *	owr,
		int	N,
		int *	iLoc,
		ProfileSPDLinSolver *	the_Solver = `nullptr`
	)

protected

Constructor.

The size of \(A\) is given by \(newIloc(N-1)\), if this is not a valid address in newIloc a segmentation fault or erronious results will result. The contents of \(iLoc\) are set equal to those of newIloc. Invokes {setLinearSOE(*this)} and setSize() on solver, printing a warning message if setSize() returns a negative number. Also creates Vector objects for \(x\) and \(b\).

Parameters

owr	analysis aggregation that owns this object.
N	size of the system.
the_Solver	pointer to the solver to use.

Member Function Documentation

§ addA()

int XC::ProfileSPDLinSOE::addA	(	const Matrix &	m,
		const ID &	id,
		double	fact = `1.0`
	)

virtual

Assembles the product of m by fact into A.

First tests that loc and M are of compatable sizes; if not a warning message is printed and a \(-1\) is returned. The LinearSOE object then assembles fact times the Matrix {M} into the matrix \(A\). The Matrix is assembled into \(A\) at the locations given by the ID object loc, i.e. \(a_{loc(i),loc(j)} += //! fact * M(i,j)\). If the location specified is outside the range, i.e. \((-1,-1)\) the corrseponding entry in M is not added to \(A\). If fact is equal to \(0.0\) or \(1.0\), more efficient steps are performed. Returns \(0\).

Implements XC::LinearSOE.

Reimplemented in XC::DistributedProfileSPDLinSOE.

§ recvSelf()

int XC::ProfileSPDLinSOE::recvSelf ( const CommParameters & cp )

virtual

Receive the object.

Each object has to receive the data needed to be able to recreate itself in the new process after it has been sent through cp. If the object is an aggregation containing other objects, new objects of the correct type can be constructed using #theBroker. To return 0 if successful or a -1 if not.

Implements XC::MovableObject.

Reimplemented in XC::DistributedProfileSPDLinSOE.

§ sendSelf()

int XC::ProfileSPDLinSOE::sendSelf ( CommParameters & cp )

virtual

Send the object.

Each object has to send the data needed to be able to reproduce that object in a remote process. The object uses the methods provided by cp object to send the data to another channel at the remote actor, the address of the channel is set before this method is called. An object of similar type at the remote actor is invoked with a receiveSelf() to receive the data. Returns 0 if successful (successful in that the data got to the channel), or a - if no data was sent.

Implements XC::MovableObject.

Reimplemented in XC::DistributedProfileSPDLinSOE.

§ setSize()

int XC::ProfileSPDLinSOE::setSize ( Graph & theGraph )

virtual

Sets the system size.

The size of the system is determined by looking at the adjacency ID of each Vertex in the Graph object G. This is done by first determining the column height for each Vertex \(i\) in G, done by setting \(iLoc(i)\) equal to \(0\) and then checking for each Vertex in G, whether any of the vertex tags in the Vertices adjacency ID results in \(iLoc(i)\) being increased. Knowing the col height of each column, the values of iLoc can be determined. Knowing {iLoc} and the size of the system (the number of Vertices in G, a check to see if the previously allocated 1d arrays for \(A\), \(x\) and \(b\) are large enough. If the memory portions allocated for the 1d arrays are not big enough, the old space is returned to the heap and new space is allocated from the heap. Printins a warning message if not enough memory is available on the heap for the 1d arrays and returns a \(-1\). If memory is available, the components of the arrays are zeroed and \(A\) is marked as being unfactored. If the system size has increased, new Vector objects for \(x\) and \(b\) using the {(double *,int)} Vector constructor are created. Finally, the result of invoking setSize() on the associated Solver object is returned.

Implements XC::LinearSOE.

Reimplemented in XC::DistributedProfileSPDLinSOE.

§ setSolver()

bool XC::ProfileSPDLinSOE::setSolver ( LinearSOESolver * newSolver )

protectedvirtual

Set the solver to use.

Invokes {setLinearSOE(*this)} on newSolver. If the system size is not equal to \(0\), it also invokes setSize() on newSolver, printing a warning and returning the returned value if this method returns a number less than \(0\). Finally it returns the result of invoking the LinearSOE classes setSolver() method.

Reimplemented from XC::LinearSOE.

§ zeroA()

void XC::ProfileSPDLinSOE::zeroA ( void )

virtual

Zeros the entries in the 1d array for \(A\) and marks the system as not having been factored.

Implements XC::LinearSOE.

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

src/solution/system_of_eqn/linearSOE/profileSPD/ProfileSPDLinSOE.h
src/solution/system_of_eqn/linearSOE/profileSPD/ProfileSPDLinSOE.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

§ ProfileSPDLinSOE() [1/3]

§ ProfileSPDLinSOE() [2/3]

§ ProfileSPDLinSOE() [3/3]

Member Function Documentation

§ addA()

§ recvSelf()

§ sendSelf()

§ setSize()

§ setSolver()

§ zeroA()