xc
|
Sparse nonsymmetric matrix linear system of equations. More...
#include <SparseGenColLinSOE.h>
Public Member Functions | |
virtual int | setSize (Graph &theGraph) |
Sets the size of the system from the number of vertices in the graph. More... | |
virtual int | addA (const Matrix &, const ID &, double fact=1.0) |
Assemblies the product fact*m into the system matrix. More... | |
virtual int | sendSelf (Communicator &) |
Send the object. More... | |
virtual int | recvSelf (const Communicator &) |
Receive the object. 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 | 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 Communicator &) |
void | receiveX (const Communicator &) |
void | receiveBX (const Communicator &) |
void | sendB (Communicator &) const |
void | sendBX (Communicator &) const |
![]() | |
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. | |
virtual double | getRCond (const char &norm='1') |
Returns the reciprocal of the condition number. More... | |
LinearSOESolver * | getSolver (void) |
Returns a pointer to the solver. | |
LinearSOESolver & | newSolver (const std::string &) |
![]() | |
int | checkSize (Graph &theGraph) const |
Check number of DOFs in the graph. | |
virtual void | save (void) const |
Stores the matrices and vectors in a temporal file. | |
virtual void | restore (void) |
Restore the matrices and vectors from a temporal file. | |
![]() | |
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 (Communicator &) |
Return the object identifier in the database. | |
void | setDbTag (int dbTag) |
Sets the object identifier in the database. | |
void | setDbTag (Communicator &) |
Sets the object identifier in the database if not already set. More... | |
virtual int | setParameter (const std::vector< std::string > &argv, Parameter ¶m) |
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 . | |
![]() | |
DistributedBase (void) | |
Constructor. | |
virtual | ~DistributedBase (void) |
Destructor. | |
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... | |
SparseGenColLinSOE (SolutionStrategy *) | |
Constructor. More... | |
SparseGenColLinSOE (SolutionStrategy *, int classTag) | |
Constructor. More... | |
SystemOfEqn * | getCopy (void) const |
Virtual constructor. | |
![]() | |
SparseGenSOEBase (SolutionStrategy *, int classTag, int N=0, int NNZ=0) | |
1d array containing coefficient of A More... | |
![]() | |
SparseSOEBase (SolutionStrategy *, int classTag, int N=0, int NNZ=0) | |
Constructor. More... | |
![]() | |
FactoredSOEBase (SolutionStrategy *, int classTag, int N=0) | |
Constructor. More... | |
![]() | |
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 (SolutionStrategy *, int classTag, int N=0) | |
Constructor. More... | |
![]() | |
int | setSolverSize (void) |
invoke setSize() on the Solver | |
LinearSOE (SolutionStrategy *, int classTag) | |
Constructor. More... | |
![]() | |
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 (SolutionStrategy *, int classTag) | |
Constructor. More... | |
Protected Attributes | |
ID | rowA |
ID | colStartA |
int arrays containing info about coeficientss in A | |
![]() | |
Vector | A |
![]() | |
int | nnz |
int | Bsize |
number of non-zeros in A | |
![]() | |
bool | factored |
True if the system is factored. | |
![]() | |
int | size |
Vector | B |
order of A | |
Vector | X |
![]() | |
std::string | tmpFileName |
File name to store matrices. | |
Friends | |
class | SolutionStrategy |
class | FEM_ObjectBroker |
class | SuperLU |
Sparse nonsymmetric matrix linear system of equations.
SparseGenColLinSOE is class which is used to store the matrix equation \(Ax=b\) of order \(size\) using a sparse column-compacted storage scheme for \(A\). The \(A\) matrix is stored in a 1d double array with \(nnz\) elements, where \(nnz\) is the number of non-zeroes in the matrix \(A\). Two additional 1d integer arrays \(rowA\) and \(colStartA\) are used to store information about the location of the coefficients, with \(colStartA(i)\) storing the location in the 1d double array of the start of column \(i\) and \(rowA(j)\) identifying the row in \(A\) to which the \(j'th\) component in the 1d double array. \(colStartA\) is of dimension \(size+1\) and \(rowA\) of dimension \(nnz\). For example
\(\) [ {array}{ccccc} a_{0,0} & 0 & a_{0,2} & a_{0,3} & 0 a_{1,0} & a_{1,1} & 0 & 0 & 0 0 & a_{2,1} & a_{2,2} & 0 & 0 0 & 0 & 0 & a_{3,3} & a_{3,4} a_{4,0} & a_{4,1} & 0 & 0 & a_{4,4} {array} ] \(\)
is stored in:
\(\) [ {array}{cccccccccccccc} a_{0,0} & a_{1,0} & a_{4,0} & a_{1,1} & a_{2,1} & a_{4,1} & a_{0,2} & a_{2,2} & a_{0,3} & a_{3,3} & a_{3,4} & a_{4,4} {array} ] \(\)
with
\(\) colStartA = [ {array}{cccccccccccccc} 0 & 3 & 6 & 8 & 10 & 12 {array} ] \(\)
and
\(\) rowA = [ {array}{cccccccccccccc} 0 & 1 & 4 & 1 & 2 & 4 & 0 & 2 & 0 & 3 & 3 & 4 {array} ] \(\) The \(x\) and \(b\) vectors are stored in 1d double arrays of length \(n\).
|
protected |
Constructor.
owr | analysis aggregation that owns this object. |
|
protected |
Constructor.
owr | analysis aggregation that owns this object. |
classTag | class identifier. |
Assemblies the product fact*m into the system matrix.
First tests that loc
and M
are of compatible 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::DistributedSparseGenColLinSOE.
|
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::DistributedSparseGenColLinSOE.
|
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::DistributedSparseGenColLinSOE.
|
virtual |
Sets the size of the system from the number of vertices in the graph.
The size of the system is determined from the Graph object {theGraph}, which must contain size
vertices labelled \(0\) through \(size-1\), the adjacency list for each vertex containing the associated vertices in a column of the matrix \(A\). The size is determined by invoking getNumVertex() on theGraph
and the number of non-zeros is determined by looking at the size of the adjacenecy list of each vertex in the graph, allowing space for the diagonal term. If the old space allocated for the 1d arrays is not big enough, it the old space is returned to the heap and new space is allocated from the heap. Prints a warning message, sets size to \(0\) and returns a \(-1\), if not enough memory is available on the heap for the 1d arrays. 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. The \(colStartA\) and \(rowA\) are then determined by looping through the vertices, setting \(colStartA(i) //! = colStartA(i-1) + 1 + \) size of Vertices \(i\) adjacency list and placing the contents of \(i\) and the adjacency list into \(rowA\) in ascending order. Finally, the result of invoking setSize() on the associated Solver object is returned.
Implements XC::LinearSOE.
Reimplemented in XC::DistributedSparseGenColLinSOE.
|
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 \(-1\) 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.