dune-common 2.11
Loading...
Searching...
No Matches
Dune::DynamicMatrix< K > Class Template Reference

Construct a matrix with a dynamic size. More...

#include <dune/common/dynmatrix.hh>

Inheritance diagram for Dune::DynamicMatrix< K >:
Inheritance graph

Public Types

typedef Base::size_type size_type
 
typedef Base::value_type value_type
 
typedef Base::row_type row_type
 
typedef Traits::derived_type derived_type
 type of derived matrix class
 
typedef Traits::value_type field_type
 export the type representing the field
 
typedef Traits::value_type block_type
 export the type representing the components
 
typedef Traits::row_reference row_reference
 The type used to represent a reference to a row (usually row_type &)
 
typedef Traits::const_row_reference const_row_reference
 The type used to represent a reference to a constant row (usually const row_type &)
 
typedef DenseIterator< DenseMatrix, row_type, row_referenceIterator
 Iterator class for sequential access.
 
typedef Iterator iterator
 typedef for stl compliant access
 
typedef Iterator RowIterator
 rename the iterators for easier access
 
typedef std::remove_reference< row_reference >::type::Iterator ColIterator
 rename the iterators for easier access
 
typedef DenseIterator< const DenseMatrix, const row_type, const_row_referenceConstIterator
 Iterator class for sequential access.
 
typedef ConstIterator const_iterator
 typedef for stl compliant access
 
typedef ConstIterator ConstRowIterator
 rename the iterators for easier access
 
typedef std::remove_reference< const_row_reference >::type::ConstIterator ConstColIterator
 rename the iterators for easier access
 

Public Member Functions

 DynamicMatrix ()
 Default constructor.
 
 DynamicMatrix (size_type r, size_type c, value_type v=value_type())
 Constructor initializing the whole matrix with a scalar.
 
 DynamicMatrix (std::initializer_list< DynamicVector< K > > const &ll)
 Constructor initializing the matrix from a list of vector.
 
template<class T , typename = std::enable_if_t<!Dune::IsNumber<T>::value && HasDenseMatrixAssigner<DynamicMatrix, T>::value>>
 DynamicMatrix (T const &rhs)
 
void resize (size_type r, size_type c, value_type v=value_type())
 resize matrix to r × c
 
template<typename T , typename = std::enable_if_t<!Dune::IsNumber<T>::value>>
DynamicMatrixoperator= (T const &rhs)
 
template<typename T , typename = std::enable_if_t<Dune::IsNumber<T>::value>>
DynamicMatrixoperator= (T scalar)
 
DynamicMatrix transposed () const
 Return transposed of the matrix as DynamicMatrix.
 
size_type mat_rows () const
 
size_type mat_cols () const
 
row_typemat_access (size_type i)
 
const row_typemat_access (size_type i) const
 
constexpr row_reference operator[] (size_type i)
 random access
 
constexpr const_row_reference operator[] (size_type i) const
 
constexpr size_type size () const
 size method (number of rows)
 
constexpr Iterator begin ()
 begin iterator
 
constexpr ConstIterator begin () const
 begin iterator
 
constexpr Iterator end ()
 end iterator
 
constexpr ConstIterator end () const
 end iterator
 
constexpr Iterator beforeEnd ()
 
constexpr ConstIterator beforeEnd () const
 
constexpr Iterator beforeBegin ()
 
constexpr ConstIterator beforeBegin () const
 
constexpr derived_typeoperator+= (const DenseMatrix< Other > &x)
 vector space addition
 
constexpr derived_type operator- () const
 Matrix negation.
 
constexpr derived_typeoperator-= (const DenseMatrix< Other > &x)
 vector space subtraction
 
constexpr derived_typeoperator*= (const field_type &k)
 vector space multiplication with scalar
 
constexpr derived_typeoperator/= (const field_type &k)
 vector space division by scalar
 
constexpr derived_typeaxpy (const field_type &a, const DenseMatrix< Other > &x)
 vector space axpy operation (*this += a x)
 
constexpr bool operator== (const DenseMatrix< Other > &x) const
 Binary matrix comparison.
 
constexpr bool operator!= (const DenseMatrix< Other > &x) const
 Binary matrix incomparison.
 
constexpr void mv (const X &x, Y &y) const
 y = A x
 
constexpr void mtv (const X &x, Y &y) const
 y = A^T x
 
constexpr void umv (const X &x, Y &y) const
 y += A x
 
constexpr void umtv (const X &x, Y &y) const
 y += A^T x
 
constexpr void umhv (const X &x, Y &y) const
 y += A^H x
 
constexpr void mmv (const X &x, Y &y) const
 y -= A x
 
constexpr void mmtv (const X &x, Y &y) const
 y -= A^T x
 
constexpr void mmhv (const X &x, Y &y) const
 y -= A^H x
 
constexpr void usmv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const
 y += alpha A x
 
constexpr void usmtv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const
 y += alpha A^T x
 
constexpr void usmhv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const
 y += alpha A^H x
 
constexpr FieldTraits< value_type >::real_type frobenius_norm () const
 frobenius norm: sqrt(sum over squared values of entries)
 
constexpr FieldTraits< value_type >::real_type frobenius_norm2 () const
 square of frobenius norm, need for block recursion
 
constexpr FieldTraits< vt >::real_type infinity_norm () const
 infinity norm (row sum norm, how to generalize for blocks?)
 
constexpr FieldTraits< vt >::real_type infinity_norm () const
 infinity norm (row sum norm, how to generalize for blocks?)
 
constexpr FieldTraits< vt >::real_type infinity_norm_real () const
 simplified infinity norm (uses Manhattan norm for complex values)
 
constexpr FieldTraits< vt >::real_type infinity_norm_real () const
 simplified infinity norm (uses Manhattan norm for complex values)
 
void solve (V1 &x, const V2 &b, bool doPivoting=true) const
 Solve system A x = b.
 
void invert (bool doPivoting=true)
 Compute inverse.
 
field_type determinant (bool doPivoting=true) const
 calculates the determinant of this matrix
 
DynamicMatrix< K > & leftmultiply (const DenseMatrix< M2 > &M)
 Multiplies M from the left to this matrix.
 
DynamicMatrix< K > & rightmultiply (const DenseMatrix< M2 > &M)
 Multiplies M from the right to this matrix.
 
constexpr size_type N () const
 number of rows
 
constexpr size_type M () const
 number of columns
 
constexpr size_type rows () const
 number of rows
 
constexpr size_type cols () const
 number of columns
 
constexpr bool exists (size_type i, size_type j) const
 return true when (i,j) is in pattern
 

Static Public Attributes

static constexpr int blocklevel
 The number of block levels we contain. This is the leaf, that is, 1.
 

Static Protected Member Functions

static void luDecomposition (DenseMatrix< DynamicMatrix< K > > &A, Func func, Mask &nonsingularLanes, bool throwEarly, bool doPivoting)
 do an LU-Decomposition on matrix A
 

Detailed Description

template<class K>
class Dune::DynamicMatrix< K >

Construct a matrix with a dynamic size.

Template Parameters
Kis the field type (use float, double, complex, etc)

Member Typedef Documentation

◆ block_type

typedef Traits::value_type Dune::DenseMatrix< DynamicMatrix< K > >::block_type
inherited

export the type representing the components

◆ ColIterator

typedef std::remove_reference<row_reference>::type::Iterator Dune::DenseMatrix< DynamicMatrix< K > >::ColIterator
inherited

rename the iterators for easier access

◆ const_iterator

typedef ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::const_iterator
inherited

typedef for stl compliant access

◆ const_row_reference

typedef Traits::const_row_reference Dune::DenseMatrix< DynamicMatrix< K > >::const_row_reference
inherited

The type used to represent a reference to a constant row (usually const row_type &)

◆ ConstColIterator

typedef std::remove_reference<const_row_reference>::type::ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::ConstColIterator
inherited

rename the iterators for easier access

◆ ConstIterator

typedef DenseIterator<const DenseMatrix,const row_type,const_row_reference> Dune::DenseMatrix< DynamicMatrix< K > >::ConstIterator
inherited

Iterator class for sequential access.

◆ ConstRowIterator

typedef ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::ConstRowIterator
inherited

rename the iterators for easier access

◆ derived_type

typedef Traits::derived_type Dune::DenseMatrix< DynamicMatrix< K > >::derived_type
inherited

type of derived matrix class

◆ field_type

typedef Traits::value_type Dune::DenseMatrix< DynamicMatrix< K > >::field_type
inherited

export the type representing the field

◆ Iterator

Iterator class for sequential access.

◆ iterator

typedef Iterator Dune::DenseMatrix< DynamicMatrix< K > >::iterator
inherited

typedef for stl compliant access

◆ row_reference

typedef Traits::row_reference Dune::DenseMatrix< DynamicMatrix< K > >::row_reference
inherited

The type used to represent a reference to a row (usually row_type &)

◆ row_type

template<class K >
typedef Base::row_type Dune::DynamicMatrix< K >::row_type

◆ RowIterator

typedef Iterator Dune::DenseMatrix< DynamicMatrix< K > >::RowIterator
inherited

rename the iterators for easier access

◆ size_type

template<class K >
typedef Base::size_type Dune::DynamicMatrix< K >::size_type

◆ value_type

template<class K >
typedef Base::value_type Dune::DynamicMatrix< K >::value_type

Constructor & Destructor Documentation

◆ DynamicMatrix() [1/4]

template<class K >
Dune::DynamicMatrix< K >::DynamicMatrix ( )
inline

Default constructor.

◆ DynamicMatrix() [2/4]

template<class K >
Dune::DynamicMatrix< K >::DynamicMatrix ( size_type  r,
size_type  c,
value_type  v = value_type() 
)
inline

Constructor initializing the whole matrix with a scalar.

◆ DynamicMatrix() [3/4]

template<class K >
Dune::DynamicMatrix< K >::DynamicMatrix ( std::initializer_list< DynamicVector< K > > const &  ll)
inline

Constructor initializing the matrix from a list of vector.

◆ DynamicMatrix() [4/4]

template<class K >
template<class T , typename = std::enable_if_t<!Dune::IsNumber<T>::value && HasDenseMatrixAssigner<DynamicMatrix, T>::value>>
Dune::DynamicMatrix< K >::DynamicMatrix ( T const &  rhs)
inline

Member Function Documentation

◆ axpy()

constexpr derived_type & Dune::DenseMatrix< DynamicMatrix< K > >::axpy ( const field_type a,
const DenseMatrix< Other > &  x 
)
inlineconstexprinherited

vector space axpy operation (*this += a x)

◆ beforeBegin() [1/2]

constexpr Iterator Dune::DenseMatrix< DynamicMatrix< K > >::beforeBegin ( )
inlineconstexprinherited
Returns
an iterator that is positioned before the first entry of the vector.

◆ beforeBegin() [2/2]

constexpr ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::beforeBegin ( ) const
inlineconstexprinherited
Returns
an iterator that is positioned before the first entry of the vector.

◆ beforeEnd() [1/2]

constexpr Iterator Dune::DenseMatrix< DynamicMatrix< K > >::beforeEnd ( )
inlineconstexprinherited
Returns
an iterator that is positioned before the end iterator of the vector, i.e. at the last entry.

◆ beforeEnd() [2/2]

constexpr ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::beforeEnd ( ) const
inlineconstexprinherited
Returns
an iterator that is positioned before the end iterator of the vector. i.e. at the last element

◆ begin() [1/2]

constexpr Iterator Dune::DenseMatrix< DynamicMatrix< K > >::begin ( )
inlineconstexprinherited

begin iterator

◆ begin() [2/2]

constexpr ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::begin ( ) const
inlineconstexprinherited

begin iterator

◆ cols()

constexpr size_type Dune::DenseMatrix< DynamicMatrix< K > >::cols ( ) const
inlineconstexprinherited

number of columns

◆ determinant()

field_type Dune::DenseMatrix< DynamicMatrix< K > >::determinant ( bool  doPivoting = true) const
inherited

calculates the determinant of this matrix

◆ end() [1/2]

constexpr Iterator Dune::DenseMatrix< DynamicMatrix< K > >::end ( )
inlineconstexprinherited

end iterator

◆ end() [2/2]

constexpr ConstIterator Dune::DenseMatrix< DynamicMatrix< K > >::end ( ) const
inlineconstexprinherited

end iterator

◆ exists()

constexpr bool Dune::DenseMatrix< DynamicMatrix< K > >::exists ( size_type  i,
size_type  j 
) const
inlineconstexprinherited

return true when (i,j) is in pattern

◆ frobenius_norm()

constexpr FieldTraits< value_type >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::frobenius_norm ( ) const
inlineconstexprinherited

frobenius norm: sqrt(sum over squared values of entries)

◆ frobenius_norm2()

constexpr FieldTraits< value_type >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::frobenius_norm2 ( ) const
inlineconstexprinherited

square of frobenius norm, need for block recursion

◆ infinity_norm() [1/2]

constexpr FieldTraits< vt >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::infinity_norm ( ) const
inlineconstexprinherited

infinity norm (row sum norm, how to generalize for blocks?)

◆ infinity_norm() [2/2]

constexpr FieldTraits< vt >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::infinity_norm ( ) const
inlineconstexprinherited

infinity norm (row sum norm, how to generalize for blocks?)

◆ infinity_norm_real() [1/2]

constexpr FieldTraits< vt >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::infinity_norm_real ( ) const
inlineconstexprinherited

simplified infinity norm (uses Manhattan norm for complex values)

◆ infinity_norm_real() [2/2]

constexpr FieldTraits< vt >::real_type Dune::DenseMatrix< DynamicMatrix< K > >::infinity_norm_real ( ) const
inlineconstexprinherited

simplified infinity norm (uses Manhattan norm for complex values)

◆ invert()

void Dune::DenseMatrix< DynamicMatrix< K > >::invert ( bool  doPivoting = true)
inherited

Compute inverse.

Exceptions
FMatrixErrorif the matrix is singular

◆ leftmultiply()

DynamicMatrix< K > & Dune::DenseMatrix< DynamicMatrix< K > >::leftmultiply ( const DenseMatrix< M2 > &  M)
inlineinherited

Multiplies M from the left to this matrix.

◆ luDecomposition()

static void Dune::DenseMatrix< DynamicMatrix< K > >::luDecomposition ( DenseMatrix< DynamicMatrix< K > > &  A,
Func  func,
Mask &  nonsingularLanes,
bool  throwEarly,
bool  doPivoting 
)
staticprotectedinherited

do an LU-Decomposition on matrix A

Parameters
AThe matrix to decompose, and to store the result in.
funcFunctor used for swapping lanes and to conduct the elimination. Depending on the functor, luDecomposition() can be used for solving, for inverting, or to compute the determinant.
nonsingularLanesSimdMask of lanes that are nonsingular.
throwEarlyWhether to throw an FMatrixError immediately as soon as one lane is discovered to be singular. If false, do not throw, instead continue until finished or all lanes are singular, and exit via return in both cases.
doPivotingEnable pivoting.

There are two modes of operation:

  • Terminate as soon as one lane is discovered to be singular. Early termination is done by throwing an FMatrixError. On entry, Simd::allTrue(nonsingularLanes) and throwEarly==true should hold. After early termination, the contents of A should be considered bogus, and nonsingularLanes has the lane(s) that triggered the early termination unset. There may be more singular lanes than the one reported in nonsingularLanes, which just haven't been discovered yet; so the value of nonsingularLanes is mostly useful for diagnostics.
  • Terminate only when all lanes are discovered to be singular. Use this when you want to apply special postprocessing in singular lines (e.g. setting the determinant of singular lanes to 0 in determinant()). On entry, nonsingularLanes may have any value and throwEarly==false should hold. The function will not throw an exception if some lanes are discovered to be singular, instead it will continue running until all lanes are singular or until finished, and terminate only via normal return. On exit, nonsingularLanes contains the map of lanes that are valid in A.

◆ M()

constexpr size_type Dune::DenseMatrix< DynamicMatrix< K > >::M ( ) const
inlineconstexprinherited

number of columns

◆ mat_access() [1/2]

template<class K >
row_type & Dune::DynamicMatrix< K >::mat_access ( size_type  i)
inline

◆ mat_access() [2/2]

template<class K >
const row_type & Dune::DynamicMatrix< K >::mat_access ( size_type  i) const
inline

◆ mat_cols()

template<class K >
size_type Dune::DynamicMatrix< K >::mat_cols ( ) const
inline

◆ mat_rows()

template<class K >
size_type Dune::DynamicMatrix< K >::mat_rows ( ) const
inline

◆ mmhv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::mmhv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y -= A^H x

◆ mmtv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::mmtv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y -= A^T x

◆ mmv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::mmv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y -= A x

◆ mtv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::mtv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y = A^T x

◆ mv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::mv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y = A x

◆ N()

constexpr size_type Dune::DenseMatrix< DynamicMatrix< K > >::N ( ) const
inlineconstexprinherited

number of rows

◆ operator!=()

constexpr bool Dune::DenseMatrix< DynamicMatrix< K > >::operator!= ( const DenseMatrix< Other > &  x) const
inlineconstexprinherited

Binary matrix incomparison.

◆ operator*=()

constexpr derived_type & Dune::DenseMatrix< DynamicMatrix< K > >::operator*= ( const field_type k)
inlineconstexprinherited

vector space multiplication with scalar

◆ operator+=()

constexpr derived_type & Dune::DenseMatrix< DynamicMatrix< K > >::operator+= ( const DenseMatrix< Other > &  x)
inlineconstexprinherited

vector space addition

◆ operator-()

constexpr derived_type Dune::DenseMatrix< DynamicMatrix< K > >::operator- ( ) const
inlineconstexprinherited

Matrix negation.

◆ operator-=()

constexpr derived_type & Dune::DenseMatrix< DynamicMatrix< K > >::operator-= ( const DenseMatrix< Other > &  x)
inlineconstexprinherited

vector space subtraction

◆ operator/=()

constexpr derived_type & Dune::DenseMatrix< DynamicMatrix< K > >::operator/= ( const field_type k)
inlineconstexprinherited

vector space division by scalar

◆ operator=() [1/2]

template<class K >
template<typename T , typename = std::enable_if_t<!Dune::IsNumber<T>::value>>
DynamicMatrix & Dune::DynamicMatrix< K >::operator= ( T const &  rhs)
inline

◆ operator=() [2/2]

template<class K >
template<typename T , typename = std::enable_if_t<Dune::IsNumber<T>::value>>
DynamicMatrix & Dune::DynamicMatrix< K >::operator= ( scalar)
inline

◆ operator==()

constexpr bool Dune::DenseMatrix< DynamicMatrix< K > >::operator== ( const DenseMatrix< Other > &  x) const
inlineconstexprinherited

Binary matrix comparison.

◆ operator[]() [1/2]

constexpr row_reference Dune::DenseMatrix< DynamicMatrix< K > >::operator[] ( size_type  i)
inlineconstexprinherited

random access

◆ operator[]() [2/2]

constexpr const_row_reference Dune::DenseMatrix< DynamicMatrix< K > >::operator[] ( size_type  i) const
inlineconstexprinherited

◆ resize()

template<class K >
void Dune::DynamicMatrix< K >::resize ( size_type  r,
size_type  c,
value_type  v = value_type() 
)
inline

resize matrix to r × c

Resize the matrix to r × c, using v as the value of all entries.

Warning
All previous entries are lost, even when the matrix was not actually resized.
Parameters
rnumber of rows
cnumber of columns
vvalue of matrix entries

◆ rightmultiply()

DynamicMatrix< K > & Dune::DenseMatrix< DynamicMatrix< K > >::rightmultiply ( const DenseMatrix< M2 > &  M)
inlineinherited

Multiplies M from the right to this matrix.

◆ rows()

constexpr size_type Dune::DenseMatrix< DynamicMatrix< K > >::rows ( ) const
inlineconstexprinherited

number of rows

◆ size()

constexpr size_type Dune::DenseMatrix< DynamicMatrix< K > >::size ( ) const
inlineconstexprinherited

size method (number of rows)

◆ solve()

void Dune::DenseMatrix< DynamicMatrix< K > >::solve ( V1 &  x,
const V2 &  b,
bool  doPivoting = true 
) const
inherited

Solve system A x = b.

Exceptions
FMatrixErrorif the matrix is singular

◆ transposed()

template<class K >
DynamicMatrix Dune::DynamicMatrix< K >::transposed ( ) const
inline

Return transposed of the matrix as DynamicMatrix.

◆ umhv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::umhv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += A^H x

◆ umtv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::umtv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += A^T x

◆ umv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::umv ( const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += A x

◆ usmhv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::usmhv ( const typename FieldTraits< Y >::field_type &  alpha,
const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += alpha A^H x

◆ usmtv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::usmtv ( const typename FieldTraits< Y >::field_type &  alpha,
const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += alpha A^T x

◆ usmv()

constexpr void Dune::DenseMatrix< DynamicMatrix< K > >::usmv ( const typename FieldTraits< Y >::field_type &  alpha,
const X &  x,
Y &  y 
) const
inlineconstexprinherited

y += alpha A x

Member Data Documentation

◆ blocklevel

constexpr int Dune::DenseMatrix< DynamicMatrix< K > >::blocklevel
staticconstexprinherited

The number of block levels we contain. This is the leaf, that is, 1.


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