#include <mumpsoverlap.hh>
Public Types | |
typedef Realtype< F >::type | r_type |
Real type of data type. | |
typedef Cmplxtype< F >::type | c_type |
Real type of data type. | |
typedef F | type |
Type of data, e.g. matrix entries. | |
Public Member Functions | |
MumpsOverlap (const int NGlob, const Operator< F > &A, const int *LocToGLob, const int rank, const int comm=-987654, const bool CConvention=true) | |
MumpsOverlap (const Operator< F > &A, const int rank, const int comm=-987654, const bool CConvention=true) | |
virtual | ~MumpsOverlap () |
virtual void | show_messages () |
virtual void | hide_messages () |
virtual void | apply_ (const Vector< F > &fncY, Vector< F > &fncX) |
virtual void | apply_ () |
virtual void | operator() (const Function< r_type > &fncY, Function< F > &fncX) |
virtual void | operator() (const Function< c_type > &fncY, Function< c_type > &fncX) |
virtual void | operator() (const Vector< r_type > &fncY, Vector< F > &fncX) |
virtual void | operator() (const Vector< c_type > &fncY, Vector< c_type > &fncX) |
void | operator() (const Matrix< r_type > &mX, Matrix< F > &mY) |
Application method to real matrices. Calls function apply() | |
void | operator() (const Matrix< c_type > &mX, Matrix< c_type > &mY) |
Application method to complex matrices. Calls apply_() | |
void | operator() () |
Application method without second argument. Used for parallel solvers. | |
virtual void | operator() (const Function< r_type > &fncY, Function< F > &fncX) |
virtual const uint | dimX () const |
virtual const uint | dimY () const |
Protected Member Functions | |
virtual std::ostream & | info (std::ostream &) const |
Returns information in an output stream. | |
Protected Attributes | |
uint | dimX_ |
Dimension of image space and the source space. | |
uint | dimY_ |
MUMPS : MUltifrontal Massively Parallel sparse direct Solver
MUMPS is a general purpose library for the direct solution of large, sparse, nonsymmetric systems of linear equations on high perormances machines, for both real and complex matrices. The library routines will perform an partial factorization and Schur complement matrix (centralized or 2D block-cyclic).
MUMPS offers different ordering options for a column permutation .
MUMPS has also parallel implementation.
Definition at line 73 of file mumpsoverlap.hh.
|
inherited |
Real type of data type.
Definition at line 120 of file compositions.hh.
|
inherited |
Real type of data type.
Definition at line 118 of file compositions.hh.
|
inherited |
Type of data, e.g. matrix entries.
Definition at line 45 of file compositions.hh.
concepts::MumpsOverlap< F >::MumpsOverlap | ( | const int | NGlob, |
const Operator< F > & | A, | ||
const int * | LocToGLob, | ||
const int | rank, | ||
const int | comm = -987654 , |
||
const bool | CConvention = true |
||
) |
Constructor
NGlob | global number of degrees of freedom |
A | Local part of the operator which should be inverted. Everything which is needed from A is copied in the constructor. |
LocToGLob | Vector that contains local to global numbering |
rank | Current thread rank. |
comm | MUMPS internal communicator |
CConvention | Define if |
concepts::MumpsOverlap< F >::MumpsOverlap | ( | const Operator< F > & | A, |
const int | rank, | ||
const int | comm = -987654 , |
||
const bool | CConvention = true |
||
) |
Constructor
A | Local part of the operator which should be inverted (of size NGlob times NGlob). Everything which is needed from A is copied in the constructor. We suppose that A is already scaled to the global matrix problem. |
rank | Current thread rank. |
comm | MUMPS internal communicator |
CConvention | Define if |
|
virtual |
Destructor
|
virtual |
Solver (call this solver on all other threads)
Implements concepts::VecOperator< F >.
|
virtual |
Solver (call this solver on thread rank=0 (master))
Implements concepts::VecOperator< F >.
|
inlinevirtualinherited |
Returns the size of the image space of the operator (number of rows of the corresponding matrix)
Definition at line 93 of file compositions.hh.
|
inlinevirtualinherited |
Returns the size of the source space of the operator (number of columns of the corresponding matrix)
Definition at line 98 of file compositions.hh.
|
virtual |
Function to hide solver messages
|
protectedvirtual |
Returns information in an output stream.
Reimplemented from concepts::VecOperator< F >.
|
virtualinherited |
Application method without second argument. Used for parallel solvers.
Reimplemented from concepts::Operator< F >.
|
virtualinherited |
Application operator for complex function fncY
.
Computes fncX
= A(fncY
) where A is this operator. fncX
becomes complex.
In derived classes its enough to implement the operator() for complex Operator's. If a real counterpart is not implemented, the function fncY
is splitted into real and imaginary part and the application operator for real functions is called for each. Then the result is combined.
If in a derived class the operator() for complex Operator's is not implemented, a exception is thrown from here.
Reimplemented from concepts::Operator< F >.
|
virtualinherited |
Application operator for real function fncY
.
Computes fncX
= A(fncY
) where A is this operator.
fncX
becomes the type of the operator, for real data it becomes real, for complex data it becomes complex.
In derived classes its enough to implement the operator() for real Operator's. If a complex counterpart is not implemented, the function fncY
is transformed to a complex function and then the application operator for complex functions is called.
If in a derived class the operator() for real Operator's is not implemented, a exception is thrown from here.
Reimplemented in aglowav::C2W< F >, aglowav2::C2W< F >, aglowav::W2C< F >, aglowav2::W2C< F >, aglowav::C2_tl2< F >, aglowav::C2tl2< F >, aglowav::CGt2< F >, aglowav::ComposeN< F >, aglowav2::Operator00< F >, bem::D< F >, bem::D_1< F >, concepts::TrivExtendRestrict< F >, sparseqr::GivensRotations< F >, vectorial::BlockOperator< F >, concepts::TrivExtendRestrict< Real >, sparseqr::GivensRotations< Real >, concepts::AfterIteration< F >, concepts::Compose< F, H >, concepts::DDSolver< F, G >, concepts::Multiple< F >, concepts::LiCoI< F >, concepts::LiCo< F >, concepts::DenseMatrix< F >, concepts::DiagonalMatrix< F >, concepts::Permutation< F >, concepts::Matrix< F >, concepts::Matrix< F::type >, and concepts::DiagonalSolver< F >.
|
virtualinherited |
Application operator for complex function fncY
.
Computes fncX
= A(fncY
) where A is this operator. fncX
becomes complex.
In derived classes its enough to implement the operator() for complex Operator's. If a real counterpart is not implemented, the vector fncY
is splitted into real and imaginary part and the application operator for real vectors is called for each. Then the result is combined
If in a derived class the operator() for complex Operator's i not implemented, a exception is thrown from here.
|
virtualinherited |
Application operator for real vector fncY
.
Computes fncX
= A(fncY
) where A is this operator.
Type of fncX
becomes that of the operator, for real data it becomes real, for complex data it becomes complex.
In derived classes its enough to implement the operator() for real Operator's. If a complex counterpart is not implemented, the vector fncY
is transformed to a complex vector and then the application for complex vectors is called.
If in a derived class the operator() for real Operator's is not implemented, a exception is thrown from here.
|
virtual |
Function to show solver messages
Reimplemented from concepts::Operator< F >.
|
protectedinherited |
Dimension of image space and the source space.
Definition at line 104 of file compositions.hh.
|
protectedinherited |
Definition at line 104 of file compositions.hh.