|
dune-common 2.10
|
A dense n x m matrix. More...
#include <dune/common/fmatrix.hh>

Public Types | |
| typedef Base::size_type | size_type |
| typedef Base::row_type | row_type |
| typedef Base::row_reference | row_reference |
| typedef Base::const_row_reference | const_row_reference |
| typedef Traits::derived_type | derived_type |
| type of derived matrix class | |
| typedef Traits::value_type | value_type |
| export the type representing the field | |
| 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 DenseIterator< DenseMatrix, row_type, row_reference > | Iterator |
| 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_reference > | ConstIterator |
| 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 | |
| constexpr | FieldMatrix ()=default |
| Default constructor. | |
| constexpr | FieldMatrix (std::initializer_list< Dune::FieldVector< K, cols > > const &l) |
| Constructor initializing the matrix from a list of vector. | |
| template<class T , typename = std::enable_if_t<HasDenseMatrixAssigner<FieldMatrix, T>::value>> | |
| FieldMatrix (T const &rhs) | |
| FieldMatrix & | operator= (const FieldMatrix &)=default |
| copy assignment operator | |
| template<typename T > | |
| FieldMatrix & | operator= (const FieldMatrix< T, ROWS, COLS > &x) |
| copy assignment from FieldMatrix over a different field | |
| template<typename T , int rows, int cols> | |
| FieldMatrix & | operator= (FieldMatrix< T, rows, cols > const &)=delete |
| no copy assignment from FieldMatrix of different size | |
| FieldMatrix< K, COLS, ROWS > | transposed () const |
| Return transposed of the matrix as FieldMatrix. | |
| template<int l> | |
| FieldMatrix< K, l, cols > | leftmultiplyany (const FieldMatrix< K, l, rows > &M) const |
| Multiplies M from the left to this matrix, this matrix is not modified. | |
| template<int r, int c> | |
| FieldMatrix & | rightmultiply (const FieldMatrix< K, r, c > &M) |
| Multiplies M from the right to this matrix. | |
| template<int l> | |
| FieldMatrix< K, rows, l > | rightmultiplyany (const FieldMatrix< K, cols, l > &M) const |
| Multiplies M from the right to this matrix, this matrix is not modified. | |
| row_reference | mat_access (size_type i) |
| const_row_reference | mat_access (size_type i) const |
| row_reference | operator[] (size_type i) |
| random access | |
| const_row_reference | operator[] (size_type i) const |
| size_type | size () const |
| size method (number of rows) | |
| Iterator | begin () |
| begin iterator | |
| ConstIterator | begin () const |
| begin iterator | |
| Iterator | end () |
| end iterator | |
| ConstIterator | end () const |
| end iterator | |
| Iterator | beforeEnd () |
| ConstIterator | beforeEnd () const |
| Iterator | beforeBegin () |
| ConstIterator | beforeBegin () const |
| derived_type & | operator+= (const DenseMatrix< Other > &x) |
| vector space addition | |
| derived_type | operator- () const |
| Matrix negation. | |
| derived_type & | operator-= (const DenseMatrix< Other > &x) |
| vector space subtraction | |
| derived_type & | operator*= (const field_type &k) |
| vector space multiplication with scalar | |
| derived_type & | operator/= (const field_type &k) |
| vector space division by scalar | |
| derived_type & | axpy (const field_type &a, const DenseMatrix< Other > &x) |
| vector space axpy operation (*this += a x) | |
| bool | operator== (const DenseMatrix< Other > &x) const |
| Binary matrix comparison. | |
| bool | operator!= (const DenseMatrix< Other > &x) const |
| Binary matrix incomparison. | |
| void | mv (const X &x, Y &y) const |
| y = A x | |
| void | mtv (const X &x, Y &y) const |
| y = A^T x | |
| void | umv (const X &x, Y &y) const |
| y += A x | |
| void | umtv (const X &x, Y &y) const |
| y += A^T x | |
| void | umhv (const X &x, Y &y) const |
| y += A^H x | |
| void | mmv (const X &x, Y &y) const |
| y -= A x | |
| void | mmtv (const X &x, Y &y) const |
| y -= A^T x | |
| void | mmhv (const X &x, Y &y) const |
| y -= A^H x | |
| void | usmv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const |
| y += alpha A x | |
| void | usmtv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const |
| y += alpha A^T x | |
| void | usmhv (const typename FieldTraits< Y >::field_type &alpha, const X &x, Y &y) const |
| y += alpha A^H x | |
| FieldTraits< value_type >::real_type | frobenius_norm () const |
| frobenius norm: sqrt(sum over squared values of entries) | |
| FieldTraits< value_type >::real_type | frobenius_norm2 () const |
| square of frobenius norm, need for block recursion | |
| FieldTraits< vt >::real_type | infinity_norm () const |
| infinity norm (row sum norm, how to generalize for blocks?) | |
| FieldTraits< vt >::real_type | infinity_norm () const |
| infinity norm (row sum norm, how to generalize for blocks?) | |
| FieldTraits< vt >::real_type | infinity_norm_real () const |
| simplified infinity norm (uses Manhattan norm for complex values) | |
| 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 | |
| FieldMatrix< K, ROWS, COLS > & | leftmultiply (const DenseMatrix< M2 > &M) |
| Multiplies M from the left to this matrix. | |
| FieldMatrix< K, ROWS, COLS > & | 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 | |
| bool | exists (size_type i, size_type j) const |
| return true when (i,j) is in pattern | |
Static Public Member Functions | |
| static constexpr size_type | mat_rows () |
| static constexpr size_type | mat_cols () |
Static Public Attributes | |
| static constexpr int | rows = ROWS |
| The number of rows. | |
| static constexpr int | cols = COLS |
| The number of columns. | |
| 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< FieldMatrix< K, ROWS, COLS > > &A, Func func, Mask &nonsingularLanes, bool throwEarly, bool doPivoting) |
| do an LU-Decomposition on matrix A | |
Detailed Description
class Dune::FieldMatrix< K, ROWS, COLS >
A dense n x m matrix.
Matrices represent linear maps from a vector space V to a vector space W. This class represents such a linear map by storing a two-dimensional array of numbers of a given field type K. The number of rows and columns is given at compile time.
Member Typedef Documentation
◆ block_type
|
inherited |
export the type representing the components
◆ ColIterator
|
inherited |
rename the iterators for easier access
◆ const_iterator
|
inherited |
typedef for stl compliant access
◆ const_row_reference
| typedef Base::const_row_reference Dune::FieldMatrix< K, ROWS, COLS >::const_row_reference |
◆ ConstColIterator
|
inherited |
rename the iterators for easier access
◆ ConstIterator
|
inherited |
Iterator class for sequential access.
◆ ConstRowIterator
|
inherited |
rename the iterators for easier access
◆ derived_type
|
inherited |
type of derived matrix class
◆ field_type
|
inherited |
export the type representing the field
◆ Iterator
|
inherited |
Iterator class for sequential access.
◆ iterator
|
inherited |
typedef for stl compliant access
◆ row_reference
| typedef Base::row_reference Dune::FieldMatrix< K, ROWS, COLS >::row_reference |
◆ row_type
| typedef Base::row_type Dune::FieldMatrix< K, ROWS, COLS >::row_type |
◆ RowIterator
|
inherited |
rename the iterators for easier access
◆ size_type
| typedef Base::size_type Dune::FieldMatrix< K, ROWS, COLS >::size_type |
◆ value_type
|
inherited |
export the type representing the field
Constructor & Destructor Documentation
◆ FieldMatrix() [1/3]
|
constexprdefault |
Default constructor.
◆ FieldMatrix() [2/3]
|
inlineconstexpr |
Constructor initializing the matrix from a list of vector.
◆ FieldMatrix() [3/3]
|
inline |
Member Function Documentation
◆ axpy()
|
inlineinherited |
vector space axpy operation (*this += a x)
◆ beforeBegin() [1/2]
|
inlineinherited |
- Returns
- an iterator that is positioned before the first entry of the vector.
◆ beforeBegin() [2/2]
|
inlineinherited |
- Returns
- an iterator that is positioned before the first entry of the vector.
◆ beforeEnd() [1/2]
|
inlineinherited |
- Returns
- an iterator that is positioned before the end iterator of the vector, i.e. at the last entry.
◆ beforeEnd() [2/2]
|
inlineinherited |
- Returns
- an iterator that is positioned before the end iterator of the vector. i.e. at the last element
◆ begin() [1/2]
|
inlineinherited |
begin iterator
◆ begin() [2/2]
|
inlineinherited |
begin iterator
◆ cols()
|
inlineconstexprinherited |
number of columns
◆ determinant()
|
inherited |
calculates the determinant of this matrix
◆ end() [1/2]
|
inlineinherited |
end iterator
◆ end() [2/2]
|
inlineinherited |
end iterator
◆ exists()
|
inlineinherited |
return true when (i,j) is in pattern
◆ frobenius_norm()
|
inlineinherited |
frobenius norm: sqrt(sum over squared values of entries)
◆ frobenius_norm2()
|
inlineinherited |
square of frobenius norm, need for block recursion
◆ infinity_norm() [1/2]
|
inlineinherited |
infinity norm (row sum norm, how to generalize for blocks?)
◆ infinity_norm() [2/2]
|
inlineinherited |
infinity norm (row sum norm, how to generalize for blocks?)
◆ infinity_norm_real() [1/2]
|
inlineinherited |
simplified infinity norm (uses Manhattan norm for complex values)
◆ infinity_norm_real() [2/2]
|
inlineinherited |
simplified infinity norm (uses Manhattan norm for complex values)
◆ invert()
|
inherited |
Compute inverse.
- Exceptions
-
FMatrixError if the matrix is singular
◆ leftmultiply()
|
inlineinherited |
Multiplies M from the left to this matrix.
◆ leftmultiplyany()
|
inline |
Multiplies M from the left to this matrix, this matrix is not modified.
◆ luDecomposition()
|
staticprotectedinherited |
do an LU-Decomposition on matrix A
- Parameters
-
A The matrix to decompose, and to store the result in. func Functor 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.nonsingularLanes SimdMask of lanes that are nonsingular. throwEarly Whether to throw an FMatrixErrorimmediately as soon as one lane is discovered to be singular. Iffalse, do not throw, instead continue until finished or all lanes are singular, and exit via return in both cases.doPivoting Enable 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)andthrowEarly==trueshould hold. After early termination, the contents ofAshould be considered bogus, andnonsingularLaneshas the lane(s) that triggered the early termination unset. There may be more singular lanes than the one reported innonsingularLanes, which just haven't been discovered yet; so the value ofnonsingularLanesis 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,nonsingularLanesmay have any value andthrowEarly==falseshould 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,nonsingularLanescontains the map of lanes that are valid inA.
◆ M()
|
inlineconstexprinherited |
number of columns
◆ mat_access() [1/2]
|
inline |
◆ mat_access() [2/2]
|
inline |
◆ mat_cols()
|
inlinestaticconstexpr |
◆ mat_rows()
|
inlinestaticconstexpr |
◆ mmhv()
|
inlineinherited |
y -= A^H x
◆ mmtv()
|
inlineinherited |
y -= A^T x
◆ mmv()
|
inlineinherited |
y -= A x
◆ mtv()
|
inlineinherited |
y = A^T x
◆ mv()
|
inlineinherited |
y = A x
◆ N()
|
inlineconstexprinherited |
number of rows
◆ operator!=()
|
inlineinherited |
Binary matrix incomparison.
◆ operator*=()
|
inlineinherited |
vector space multiplication with scalar
◆ operator+=()
|
inlineinherited |
vector space addition
◆ operator-()
|
inlineinherited |
Matrix negation.
◆ operator-=()
|
inlineinherited |
vector space subtraction
◆ operator/=()
|
inlineinherited |
vector space division by scalar
◆ operator=() [1/3]
|
default |
copy assignment operator
◆ operator=() [2/3]
|
inline |
copy assignment from FieldMatrix over a different field
◆ operator=() [3/3]
|
delete |
no copy assignment from FieldMatrix of different size
◆ operator==()
|
inlineinherited |
Binary matrix comparison.
◆ operator[]() [1/2]
|
inlineinherited |
random access
◆ operator[]() [2/2]
|
inlineinherited |
◆ rightmultiply() [1/2]
|
inlineinherited |
Multiplies M from the right to this matrix.
◆ rightmultiply() [2/2]
|
inline |
Multiplies M from the right to this matrix.
◆ rightmultiplyany()
|
inline |
Multiplies M from the right to this matrix, this matrix is not modified.
◆ rows()
|
inlineconstexprinherited |
number of rows
◆ size()
|
inlineinherited |
size method (number of rows)
◆ solve()
|
inherited |
Solve system A x = b.
- Exceptions
-
FMatrixError if the matrix is singular
◆ transposed()
|
inline |
Return transposed of the matrix as FieldMatrix.
◆ umhv()
|
inlineinherited |
y += A^H x
◆ umtv()
|
inlineinherited |
y += A^T x
◆ umv()
|
inlineinherited |
y += A x
◆ usmhv()
|
inlineinherited |
y += alpha A^H x
◆ usmtv()
|
inlineinherited |
y += alpha A^T x
◆ usmv()
|
inlineinherited |
y += alpha A x
Member Data Documentation
◆ blocklevel
|
staticconstexprinherited |
The number of block levels we contain. This is the leaf, that is, 1.
◆ cols
|
staticconstexpr |
The number of columns.
◆ rows
|
staticconstexpr |
The number of rows.
The documentation for this class was generated from the following files:
Legal Statements / Impressum | Hosted by TU Dresden & Uni Heidelberg | Generated by
1.9.8