Solver for a complex system. More...

#include <solverConjugate.hh>

Inheritance diagram for concepts::SolverConjugate:
concepts::VecOperator< Cmplx > concepts::Operator< Cmplx >

Public Types

typedef Cmplxtype< Cmplx >::type c_type
 Real type of data type. More...
 
typedef Realtype< Cmplx >::type r_type
 Real type of data type. More...
 
typedef Cmplx type
 Type of data, e.g. matrix entries. More...
 

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) More...
 
virtual const uint dimY () const
 Returns the size of the source space of the operator (number of columns of the corresponding matrix) More...
 
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< Cmplx > &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< Cmplx > &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< Cmplx > &fncX)
 Application operator for real vector fncY. More...
 
virtual void show_messages ()
 
 SolverConjugate (Operator< Cmplx > &A, Operator< Cmplx > &B, SolverFabric< Real > &solverfabric)
 Constructor. More...
 

Protected Member Functions

virtual void apply_ ()
 Intrinsic application method without argument. More...
 
virtual void apply_ (const Vector< Cmplx > &fncY, Vector< Cmplx > &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

std::unique_ptr< SparseMatrix< Real > > S_
 Real system matrix. More...
 
std::unique_ptr< Operator< Real > > solver_
 Solver for the real system matrix. More...
 

Detailed Description

Solver for a complex system.

\[A x + B \overline{x} = b\]

with given complex matrices $A$ and $B$ and the complex vector of the right hand side $b$.

It solves indeed the double sized real system

\[\left(\begin{array}{cc}A_R + B_R & B_I - A_I\\A_I + B_I & A_R - B_R\end{array}\right) \left(\begin{array}{c}x_R \\ x_I\end{array}\right) = \left(\begin{array}{c}b_R \\ b_I\end{array}\right),\]

where the subscripts $R$ and $I$ mean the real and imaginary parts, respectively.

Author
Kersten Schmidt, 2014

Definition at line 37 of file solverConjugate.hh.

Member Typedef Documentation

◆ c_type

Real type of data type.

Definition at line 120 of file compositions.hh.

◆ r_type

Real type of data type.

Definition at line 118 of file compositions.hh.

◆ type

typedef Cmplx concepts::Operator< Cmplx >::type
inherited

Type of data, e.g. matrix entries.

Definition at line 45 of file compositions.hh.

Constructor & Destructor Documentation

◆ SolverConjugate()

concepts::SolverConjugate::SolverConjugate ( Operator< Cmplx > &  A,
Operator< Cmplx > &  B,
SolverFabric< Real > &  solverfabric 
)

Constructor.

Parameters
Acomplex matrix related to $x$
Bcomplex matrix related to $\overline{x}$
solverfabricfabric of the solver for the real system

Member Function Documentation

◆ apply_() [1/2]

virtual void concepts::SolverConjugate::apply_ ( )
protectedvirtual

Intrinsic application method without argument.

Implements concepts::VecOperator< Cmplx >.

◆ apply_() [2/2]

virtual void concepts::SolverConjugate::apply_ ( const Vector< Cmplx > &  fncY,
Vector< Cmplx > &  fncX 
)
protectedvirtual

Intrinsic application method, i.e.

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

Implements concepts::VecOperator< Cmplx >.

◆ dimX()

virtual const uint concepts::Operator< Cmplx >::dimX
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< Cmplx >::dimY
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::SolverConjugate::info ( std::ostream &  os) const
protectedvirtual

Reimplemented from concepts::VecOperator< Cmplx >.

◆ operator()() [1/7]

void concepts::VecOperator< Cmplx >::operator()
virtualinherited

Application method without second argument. Used for parallel solvers.

Reimplemented from concepts::Operator< Cmplx >.

◆ operator()() [2/7]

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

◆ operator()() [3/7]

virtual void concepts::VecOperator< Cmplx >::operator() ( const Function< r_type > &  fncY,
Function< Cmplx > &  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< Cmplx >.

◆ operator()() [4/7]

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

Application method to complex matrices. Calls apply_()

◆ operator()() [5/7]

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

Application method to real matrices. Calls function apply()

◆ operator()() [6/7]

virtual void concepts::VecOperator< Cmplx >::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< Cmplx >::operator() ( const Vector< r_type > &  fncY,
Vector< Cmplx > &  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< Cmplx >::show_messages
inlinevirtualinherited

Definition at line 100 of file compositions.hh.

Member Data Documentation

◆ dimX_

uint concepts::Operator< Cmplx >::dimX_
protectedinherited

Dimension of image space and the source space.

Definition at line 104 of file compositions.hh.

◆ dimY_

uint concepts::Operator< Cmplx >::dimY_
protectedinherited

Definition at line 104 of file compositions.hh.

◆ S_

std::unique_ptr<SparseMatrix<Real> > concepts::SolverConjugate::S_
private

Real system matrix.

Definition at line 58 of file solverConjugate.hh.

◆ solver_

std::unique_ptr<Operator<Real> > concepts::SolverConjugate::solver_
private

Solver for the real system matrix.

Definition at line 60 of file solverConjugate.hh.


The documentation for this class was generated from the following file:
Page URL: http://wiki.math.ethz.ch/bin/view/Concepts/WebHome
21 August 2020
© 2020 Eidgenössische Technische Hochschule Zürich