Linear solver using Lapack subroutine DPOSV. More...

#include <lapack.hh>

Inheritance diagram for concepts::LapackChol:

Public Types
typedef Cmplxtype< Real >::type	c_type
	Real type of data type. More...

typedef Realtype< Real >::type	r_type
	Real type of data type. More...

typedef Real	type
	Type of data, e.g. matrix entries. More...

enum	UpLo { U = 0, L = 1 }

Public Member Functions
virtual void	apply (DiagonalMatrix< Real > &Dinv, const Vector< Real > &fncY, Vector< Real > &fncX)

virtual const uint	dimX () const
	Returns the size of the image space of the operator (number of rows of the corresponding matrix) More...

virtual const uint	dimY () const
	Returns the size of the source space of the operator (number of columns of the corresponding matrix) More...

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

void	operator() ()
	Application method without second argument. Used for parallel solvers. More...

virtual void	operator() (const Function< c_type > &fncY, Function< c_type > &fncX)
	Application operator for complex function `fncY`. More...

virtual void	operator() (const Function< r_type > &fncY, Function< Real > &fncX)
	Application operator for real function `fncY`. More...

void	operator() (const Matrix< c_type > &mX, Matrix< c_type > &mY)
	Application method to complex matrices. Calls apply_() More...

void	operator() (const Matrix< r_type > &mX, Matrix< Real > &mY)
	Application method to real matrices. Calls function apply() More...

virtual void	operator() (const Vector< c_type > &fncY, Vector< c_type > &fncX)
	Application operator for complex function `fncY`. More...

virtual void	operator() (const Vector< r_type > &fncY, Vector< Real > &fncX)
	Application operator for real vector `fncY`. More...

virtual void	show_messages ()

Protected Member Functions
virtual void	apply_ ()
	Intrinsic application method without argument. More...

virtual void	apply_ (const Vector< Real > &fncY, Vector< Real > &fncX)
	Intrinsic application method, i.e. More...

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

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

uint	dimY_

Private Attributes
ElementMatrix< Real >	A_
	Storage for factorised matrix. More...

bool	factorized_
	Flag if the matrix is factorized. More...

char	uplo_
	Flag, if the upper or the lower triangular part is stored in the given matrix. More...

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.

Author: Kersten Schmidt, 2010

Definition at line 30 of file lapack.hh.

Member Typedef Documentation

◆ c_type

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

inherited

Real type of data type.

Definition at line 120 of file compositions.hh.

◆ r_type

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

inherited

Real type of data type.

Definition at line 118 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

◆ UpLo

enum concepts::LapackChol::UpLo

Enumerator
U
L

Definition at line 33 of file lapack.hh.

Constructor & Destructor Documentation

◆ LapackChol()

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

Member Function Documentation

◆ apply()

virtual void concepts::LapackChol::apply	(	DiagonalMatrix< Real > &	Dinv,
		const Vector< Real > &	fncY,
		Vector< Real > &	fncX
	)

virtual

◆ apply_() [1/2]

virtual void concepts::LapackChol::apply_ ( )

protectedvirtual

Intrinsic application method without argument.

Implements concepts::VecOperator< Real >.

◆ apply_() [2/2]

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

protectedvirtual

Intrinsic application method, i.e.

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

Implements concepts::VecOperator< Real >.

◆ 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::LapackChol::info ( std::ostream & os ) const

protectedvirtual

Reimplemented from concepts::VecOperator< Real >.

◆ operator()() [1/7]

void concepts::VecOperator< Real >::operator()

virtualinherited

Application method without second argument. Used for parallel solvers.

Reimplemented from concepts::Operator< Real >.

◆ operator()() [2/7]

virtual void concepts::VecOperator< 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.

Reimplemented from concepts::Operator< Real >.

◆ operator()() [3/7]

virtual void concepts::VecOperator< 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.

Reimplemented from concepts::Operator< Real >.

◆ operator()() [4/7]

void concepts::VecOperator< Real >::operator()	(	const Matrix< c_type > &	mX,
		Matrix< c_type > &	mY
	)

inherited

Application method to complex matrices. Calls apply_()

◆ operator()() [5/7]

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

inherited

Application method to real matrices. Calls function apply()

◆ operator()() [6/7]

virtual void concepts::VecOperator< Real >::operator()	(	const Vector< c_type > &	fncY,
		Vector< 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 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.

◆ operator()() [7/7]

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

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.

◆ show_messages()

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

inlinevirtualinherited

Definition at line 100 of file compositions.hh.

Member Data Documentation

◆ A_

ElementMatrix<Real> concepts::LapackChol::A_

mutableprivate

Storage for factorised matrix.

Definition at line 53 of file lapack.hh.

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

◆ factorized_

bool concepts::LapackChol::factorized_

mutableprivate

Flag if the matrix is factorized.

Definition at line 60 of file lapack.hh.

◆ uplo_

char concepts::LapackChol::uplo_

private

Flag, if the upper or the lower triangular part is stored in the given matrix.

Definition at line 58 of file lapack.hh.

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

operator/lapack.hh

concepts::LapackChol Class Reference

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

◆ c_type

◆ r_type

◆ type

Member Enumeration Documentation

◆ UpLo

Constructor & Destructor Documentation

◆ LapackChol()

Member Function Documentation

◆ apply()

◆ apply_() [1/2]

◆ apply_() [2/2]

◆ dimX()

◆ dimY()

◆ info()

◆ operator()() [1/7]

◆ operator()() [2/7]

◆ operator()() [3/7]

◆ operator()() [4/7]

◆ operator()() [5/7]

◆ operator()() [6/7]

◆ operator()() [7/7]

◆ show_messages()

Member Data Documentation

◆ A_

◆ dimX_

◆ dimY_

◆ factorized_

◆ uplo_