HLIBpro 3.1
Loading...
Searching...
No Matches
TLinearOperator< T_value > Class Template Referenceabstract

Base class for all linear operators mapping vectors to vectors. More...

#include <TLinearOperator.hh>

Inheritance diagram for TLinearOperator< T_value >:
TTypeInfo TMatrix< value_t > TFacInvMatrix< T_value > TGaussSeidel< T_value > TJacobi< T_value > TMatrix< T_value > TMatrixProduct< T_value > TMatrixSum< T_value > TNearfieldMulVec< T_value > TPermMatrix< T_value > TSOR< T_value >

Public Member Functions

virtual bool is_complex () const
 return true, if field type is complex valued
 
virtual bool is_real () const
 return true, if field type is real valued
 
virtual bool is_self_adjoint () const =0
 return true, of operator is self adjoint
 
virtual void apply (const TVector< value_t > *x, TVector< value_t > *y, const matop_t op=apply_normal) const =0
 
virtual void apply_add (const value_t alpha, const TVector< value_t > *x, TVector< value_t > *y, const matop_t op=apply_normal) const =0
 
virtual void apply_add (const value_t alpha, const BLAS::Vector< value_t > &x, BLAS::Vector< value_t > &y, const matop_t op=apply_normal) const =0
 
virtual size_t domain_dim () const =0
 return dimension of domain
 
virtual size_t range_dim () const =0
 return dimension of range
 
virtual auto domain_vector () const -> std::unique_ptr< TVector< value_t > >=0
 return vector in domain space
 
virtual auto range_vector () const -> std::unique_ptr< TVector< value_t > >=0
 return vector in range space
 
- Public Member Functions inherited from TTypeInfo
virtual typeid_t type () const =0
 return type ID of object
 
virtual bool is_type (const typeid_t t) const
 return true if local object is of given type ID t
 
virtual std::string typestr () const
 return string representation of type
 

Detailed Description

template<typename T_value>
class Hpro::TLinearOperator< T_value >
     Many standard arithmetic operations only depend upon a linear operator
     providing the mapping between vectors, e.g. iterativ solvers. Instead of
     requiring a full matrix and hence the need for an implementation of the
     full matrix algebra, an object of type TLinearOperator is fully sufficient
     in such cases.

Member Function Documentation

◆ apply()

template<typename T_value >
virtual void apply ( const TVector< value_t > *  x,
TVector< value_t > *  y,
const matop_t  op = apply_normal 
) const
pure virtual

mapping function of linear operator \(A\), e.g. \( y := A(x)\). Depending on op, either \(A\), \(A^T\) or \(A^H\) is applied.

Implemented in TFacInvMatrix< T_value >, TJacobi< T_value >, TMatrix< T_value >, TMatrixProduct< T_value >, TMatrixSum< T_value >, TNearfieldMulVec< T_value >, TPermMatrix< T_value >, TSOR< T_value >, and TGaussSeidel< T_value >.

◆ apply_add() [1/2]

template<typename T_value >
virtual void apply_add ( const value_t  alpha,
const BLAS::Vector< value_t > &  x,
BLAS::Vector< value_t > &  y,
const matop_t  op = apply_normal 
) const
pure virtual

◆ apply_add() [2/2]

template<typename T_value >
virtual void apply_add ( const value_t  alpha,
const TVector< value_t > *  x,
TVector< value_t > *  y,
const matop_t  op = apply_normal 
) const
pure virtual

mapping function with update: \( y := y + \alpha A(x)\). Depending on op, either \(A\), \(A^T\) or \(A^H\) is applied.

Implemented in TFacInvMatrix< T_value >, TJacobi< T_value >, TMatrix< T_value >, TMatrixProduct< T_value >, TMatrixSum< T_value >, TNearfieldMulVec< T_value >, TPermMatrix< T_value >, TSOR< T_value >, and TGaussSeidel< T_value >.