concepts::Pardiso Class Reference

#include <pardiso.hh>

Inheritance diagram for concepts::Pardiso:

Public Types
enum	MATRIX_TYPE { STRUCT_SYMM = 1 , SPD = 2 , SYMM_INDEF = -2 , UNSYMM = 11 }
	Type of matrix to solve. More...
typedef Real	type
	Type of data, e.g. matrix entries.
typedef Realtype< Real >::type	r_type
	Real type of data type.
typedef Cmplxtype< Real >::type	c_type
	Real type of data type.

Public Member Functions
	Pardiso (const SparseMatrix< Real > &matrix, enum MATRIX_TYPE type)
virtual void	operator() (const Function< Real > &fncY, Function< Real > &fncX)
void	operator() (const Vector< Real > &fncY, Vector< Real > &fncX)
virtual const Space< Real > &	spaceX () const
virtual const Space< Real > &	spaceY () const
virtual void	operator() (const Function< r_type > &fncY, Function< Real > &fncX)
virtual void	operator() (const Function< c_type > &fncY, Function< c_type > &fncX)
virtual void	operator() ()
virtual const uint	dimX () const
virtual const uint	dimY () const
virtual void	show_messages ()

Protected Member Functions
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.

Protected Attributes
uint	dimX_
	Dimension of image space and the source space.
uint	dimY_

Detailed Description

Sparse direct solver for symmetric and unsymmetric matrices. This class interfaces to PARDISO [1]. Use the environment variable OMP_NUM_THREADS to change the number of threas used (default is 1).

The package PARDISO is a thread-safe, high-performance, robust, memory efficient and easy to use software for solving large sparse symmetric and unsymmetric linear systems of equations on shared memory multiprocessors. The solver uses a combination of left- and right-looking Level-3 BLAS supernode techniques to exploit pipelining parallelism and memory hierarchies. In order to improve sequential and parallel sparse numerical factorization performance, the algorithms are based on a Level-3 BLAS update and pipelining parallelism is exploited with a combination of left- and right-looking supernode techniques.

See also

[1] Olaf Schenk, Klaus Gärtner, Wolfgang Fichtner, Andreas Stricker, PARDISO: a high-performance serial and parallel sparse linear Linear Solvers in semiconductor device simulation, Future Generation Computer System 18 (2001), 69-78.

[2] Pardiso Homepage

Author

Philipp Frauenfelder, 2003

Test:

test::PardisoTest

Examples

hpFEM2d.cc.

Definition at line 43 of file pardiso.hh.

Member Typedef Documentation

c_type

typedef Cmplxtype<Real >::type concepts::Operator< Real >::c_type

inherited

Real type of data type.

Definition at line 49 of file compositions.hh.

r_type

typedef Realtype<Real >::type concepts::Operator< Real >::r_type

inherited

Real type of data type.

Definition at line 47 of file compositions.hh.

type

typedef Real concepts::Operator< Real >::type

inherited

Type of data, e.g. matrix entries.

Definition at line 45 of file compositions.hh.

Member Enumeration Documentation

MATRIX_TYPE

enum concepts::Pardiso::MATRIX_TYPE

Type of matrix to solve.

Definition at line 46 of file pardiso.hh.

Constructor & Destructor Documentation

Pardiso()

concepts::Pardiso::Pardiso	(	const SparseMatrix< Real > &	matrix,
		enum MATRIX_TYPE	type
	)

Constructor.

Parameters

matrix	Matrix
type	Type of the matrix

Member Function Documentation

dimX()

virtual const uint concepts::Operator< Real >::dimX ( ) const

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.

dimY()

virtual const uint concepts::Operator< Real >::dimY ( ) const

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.

info()

virtual std::ostream & concepts::Pardiso::info ( std::ostream &  os ) const

protectedvirtual

Returns information in an output stream.

Reimplemented from concepts::Operator< Real >.

operator()() [1/3]

virtual void concepts::Operator< Real >::operator() ( )

virtualinherited

Application operator without argument

operator()() [2/3]

virtual void concepts::Operator< Real >::operator() ( const Function< c_type > &  fncY, Function< c_type > &  fncX  )

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.

operator()() [3/3]

virtual void concepts::Operator< Real >::operator() ( const Function< r_type > &  fncY, Function< Real > &  fncX  )

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.

show_messages()

virtual void concepts::Operator< Real >::show_messages ( )

inlinevirtualinherited

Definition at line 100 of file compositions.hh.

spaceX()

virtual const Space< Real > & concepts::Pardiso::spaceX ( ) const

inlinevirtual

Definition at line 64 of file pardiso.hh.

spaceY()

virtual const Space< Real > & concepts::Pardiso::spaceY ( ) const

inlinevirtual

Definition at line 65 of file pardiso.hh.

Member Data Documentation

dimX_

uint concepts::Operator< Real >::dimX_

protectedinherited

Dimension of image space and the source space.

Definition at line 104 of file compositions.hh.

dimY_

uint concepts::Operator< Real >::dimY_

protectedinherited

Definition at line 104 of file compositions.hh.

---

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

• operator/pardiso.hh