Linear solver using Lapack subroutine DPOSV. More...

#include <lapack.hh>

Inheritance diagram for concepts::LapackChol:

Collaboration diagram for concepts::LapackChol:

Public Types
typedef Cmplxtype< Real >::type	c_type
	Real type of data type.
typedef Realtype< Real >::type	r_type
	Real type of data type.
typedef Real	type
	Type of data, e.g. matrix entries.
enum	UpLo { U = 0, L = 1 }
Public Member Functions
virtual const uint	dimX () const
	Returns the size of the image space of the operator (number of rows of the corresponding matrix)
virtual const uint	dimY () const
	Returns the size of the source space of the operator (number of columns of the corresponding matrix)
	LapackChol (const ElementMatrix< Real > A, enum UpLo uplo=U)
void	operator() (const Matrix< c_type > &mX, Matrix< c_type > &mY)
	Application method to complex matrices. Calls apply_()
void	operator() (const Matrix< r_type > &mX, Matrix< Real > &mY)
	Application method to real matrices. Calls function apply()
virtual void	operator() (const Vector< c_type > &fncY, Vector< c_type > &fncX)
	Application operator for complex function `fncY`.
virtual void	operator() (const Function< r_type > &fncY, Function< Real > &fncX)
	Application operator for real function `fncY`.
virtual void	operator() (const Function< c_type > &fncY, Function< c_type > &fncX)
	Application operator for complex function `fncY`.
virtual void	operator() (const Vector< r_type > &fncY, Vector< Real > &fncX)
	Application operator for real vector `fncY`.
Protected Member Functions
virtual void	apply_ (const Vector< Real > &fncY, Vector< Real > &fncX)
	Intrinsic application method, i.e.
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_
Private Attributes
ElementMatrix< Real >	A_
	Storage for factorised matrix.
bool	factorized_
	Flag if the matrix is factorized.
char	uplo_
	Flag, if the upper or the lower triangular part is stored in the given matrix.

Detailed Description

Linear solver using Lapack subroutine DPOSV.

Solves for a symmetric positive definite (dense) matrix.

The application operator can be called with the same objects for right hand side and solution.

Definition at line 23 of file lapack.hh.

Member Typedef Documentation

typedef Cmplxtype<Real >::type concepts::VecOperator< Real >::c_type [inherited]

Real type of data type.

Reimplemented from concepts::Operator< Real >.

Definition at line 112 of file compositions.hh.

typedef Realtype<Real >::type concepts::VecOperator< Real >::r_type [inherited]

Real type of data type.

Reimplemented from concepts::Operator< Real >.

Definition at line 110 of file compositions.hh.

typedef Real concepts::Operator< Real >::type [inherited]

Type of data, e.g. matrix entries.

Definition at line 43 of file compositions.hh.

Member Enumeration Documentation

enum concepts::LapackChol::UpLo

Enumerator:

U
L

Definition at line 26 of file lapack.hh.

Constructor & Destructor Documentation

concepts::LapackChol::LapackChol	(	const ElementMatrix< Real >	A,
		enum UpLo	uplo = `U`
	)

Member Function Documentation

virtual void concepts::LapackChol::apply_	(	const Vector< Real > &	fncY,
		Vector< Real > &	fncX
	)		`[protected, virtual]`

Intrinsic application method, i.e.

real Operator and real Vector or complex Operator and real Vector.

Implements concepts::VecOperator< Real >.

virtual const uint concepts::Operator< Real >::dimX ( ) const [inline, virtual, inherited]

Returns the size of the image space of the operator (number of rows of the corresponding matrix)

Examples:: hpFEM2d-simple.cc, and hpFEM2d.cc.

Definition at line 87 of file compositions.hh.

virtual const uint concepts::Operator< Real >::dimY ( ) const [inline, virtual, inherited]

Returns the size of the source space of the operator (number of columns of the corresponding matrix)

Definition at line 92 of file compositions.hh.

virtual std::ostream& concepts::LapackChol::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from concepts::VecOperator< Real >.

void concepts::VecOperator< Real >::operator()	(	const Matrix< r_type > &	mX,
		Matrix< Real > &	mY
	)		`[inherited]`

Application method to real matrices. Calls function apply()

virtual void concepts::VecOperator< Real >::operator()	(	const Vector< c_type > &	fncY,
		Vector< c_type > &	fncX
	)		`[virtual, inherited]`

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.

virtual void concepts::VecOperator< Real >::operator()	(	const Vector< r_type > &	fncY,
		Vector< Real > &	fncX
	)		`[virtual, inherited]`

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 void concepts::VecOperator< Real >::operator()	(	const Function< c_type > &	fncY,
		Function< c_type > &	fncX
	)		`[virtual, inherited]`

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

virtual void concepts::VecOperator< Real >::operator()	(	const Function< r_type > &	fncY,
		Function< Real > &	fncX
	)		`[virtual, inherited]`

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 from concepts::Operator< Real >.