Reference documentation for deal.II version 9.6.0
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
TridiagonalMatrix< number > Class Template Reference

#include <deal.II/lac/tridiagonal_matrix.h>

Public Types

Constructors
using size_type = types::global_dof_index
 

Public Member Functions

Constructors and initialization
 TridiagonalMatrix (size_type n=0, bool symmetric=false)
 
void reinit (size_type n, bool symmetric=false)
 
Non-modifying operators
size_type m () const
 
size_type n () const
 
bool all_zero () const
 
Element access
number operator() (size_type i, size_type j) const
 
number & operator() (size_type i, size_type j)
 
Multiplications with vectors
void vmult (Vector< number > &w, const Vector< number > &v, const bool adding=false) const
 
void vmult_add (Vector< number > &w, const Vector< number > &v) const
 
void Tvmult (Vector< number > &w, const Vector< number > &v, const bool adding=false) const
 
void Tvmult_add (Vector< number > &w, const Vector< number > &v) const
 
number matrix_scalar_product (const Vector< number > &u, const Vector< number > &v) const
 
number matrix_norm_square (const Vector< number > &v) const
 
LAPACK operations
void compute_eigenvalues ()
 
number eigenvalue (const size_type i) const
 
Miscellanea
template<class OutputStream >
void print (OutputStream &s, const unsigned int width=5, const unsigned int precision=2) const
 

Private Attributes

std::vector< number > diagonal
 
std::vector< number > left
 
std::vector< number > right
 
bool is_symmetric
 
LAPACKSupport::State state
 

Detailed Description

template<typename number>
class TridiagonalMatrix< number >

A quadratic tridiagonal matrix. That is, a matrix where all entries are zero, except the diagonal and the entries left and right of it.

The matrix has an additional symmetric mode, in which case only the upper triangle of the matrix is stored and mirrored to the lower one for matrix vector operations.

Definition at line 52 of file tridiagonal_matrix.h.

Member Typedef Documentation

◆ size_type

template<typename number >
using TridiagonalMatrix< number >::size_type = types::global_dof_index

Declare type for container size.

Definition at line 60 of file tridiagonal_matrix.h.

Constructor & Destructor Documentation

◆ TridiagonalMatrix()

template<typename number >
TridiagonalMatrix< number >::TridiagonalMatrix ( size_type n = 0,
bool symmetric = false )

Constructor generating an empty matrix of dimension n.

Definition at line 27 of file tridiagonal_matrix.cc.

Member Function Documentation

◆ reinit()

template<typename number >
void TridiagonalMatrix< number >::reinit ( size_type n,
bool symmetric = false )

Reinitialize the matrix to a new size and reset all entries to zero. The symmetry properties may be set as well.

Definition at line 39 of file tridiagonal_matrix.cc.

◆ m()

template<typename number >
size_type TridiagonalMatrix< number >::m ( ) const

Number of rows of this matrix. Note that the matrix is an m x m matrix.

◆ n()

template<typename number >
size_type TridiagonalMatrix< number >::n ( ) const

Number of columns of this matrix. Note that the matrix is an n x n matrix.

◆ all_zero()

template<typename number >
bool TridiagonalMatrix< number >::all_zero ( ) const

Return whether the matrix contains only elements with value zero. This function is mainly for internal consistency checks and should seldom be used when not in debug mode since it uses quite some time.

Definition at line 52 of file tridiagonal_matrix.cc.

◆ operator()() [1/2]

template<typename number >
number TridiagonalMatrix< number >::operator() ( size_type i,
size_type j ) const

Read-only access to a value. This is restricted to the case where |i-j| <= 1.

◆ operator()() [2/2]

template<typename number >
number & TridiagonalMatrix< number >::operator() ( size_type i,
size_type j )

Read-write access to a value. This is restricted to the case where |i-j| <= 1.

Note
In case of symmetric storage technique, the entries (i,j) and (j,i) are identified and both exist. This must be taken into account if adding up is used for matrix assembling in order not to obtain doubled entries.

◆ vmult()

template<typename number >
void TridiagonalMatrix< number >::vmult ( Vector< number > & w,
const Vector< number > & v,
const bool adding = false ) const

Matrix-vector-multiplication. Multiplies v from the right and stores the result in w.

If the optional parameter adding is true, the result is added to w.

Source and destination must not be the same vector.

Definition at line 80 of file tridiagonal_matrix.cc.

◆ vmult_add()

template<typename number >
void TridiagonalMatrix< number >::vmult_add ( Vector< number > & w,
const Vector< number > & v ) const

Adding Matrix-vector-multiplication. Same as vmult() with parameter adding=true, but widely used in deal.II classes.

Source and destination must not be the same vector.

Definition at line 131 of file tridiagonal_matrix.cc.

◆ Tvmult()

template<typename number >
void TridiagonalMatrix< number >::Tvmult ( Vector< number > & w,
const Vector< number > & v,
const bool adding = false ) const

Transpose matrix-vector-multiplication. Multiplies vT from the left and stores the result in w.

If the optional parameter adding is true, the result is added to w.

Source and destination must not be the same vector.

Definition at line 141 of file tridiagonal_matrix.cc.

◆ Tvmult_add()

template<typename number >
void TridiagonalMatrix< number >::Tvmult_add ( Vector< number > & w,
const Vector< number > & v ) const

Adding transpose matrix-vector-multiplication. Same as Tvmult() with parameter adding=true, but widely used in deal.II classes.

Source and destination must not be the same vector.

Definition at line 186 of file tridiagonal_matrix.cc.

◆ matrix_scalar_product()

template<typename number >
number TridiagonalMatrix< number >::matrix_scalar_product ( const Vector< number > & u,
const Vector< number > & v ) const

Build the matrix scalar product u^T M v. This function is mostly useful when building the cellwise scalar product of two functions in the finite element context.

Definition at line 196 of file tridiagonal_matrix.cc.

◆ matrix_norm_square()

template<typename number >
number TridiagonalMatrix< number >::matrix_norm_square ( const Vector< number > & v) const

Return the square of the norm of the vector v with respect to the norm induced by this matrix, i.e. (v,Mv). This is useful, e.g. in the finite element context, where the L2 norm of a function equals the matrix norm with respect to the mass matrix of the vector representing the nodal values of the finite element function.

Obviously, the matrix needs to be quadratic for this operation.

Definition at line 223 of file tridiagonal_matrix.cc.

◆ compute_eigenvalues()

template<typename number >
void TridiagonalMatrix< number >::compute_eigenvalues ( )

Compute the eigenvalues of the symmetric tridiagonal matrix.

Note
This function requires configuration of deal.II with LAPACK support. Additionally, the matrix must use symmetric storage technique.

Definition at line 232 of file tridiagonal_matrix.cc.

◆ eigenvalue()

template<typename number >
number TridiagonalMatrix< number >::eigenvalue ( const size_type i) const

After calling compute_eigenvalues(), you can access each eigenvalue here.

Definition at line 260 of file tridiagonal_matrix.cc.

◆ print()

template<typename number >
template<class OutputStream >
void TridiagonalMatrix< number >::print ( OutputStream & s,
const unsigned int width = 5,
const unsigned int precision = 2 ) const

Output of the matrix in user-defined format.

Member Data Documentation

◆ diagonal

template<typename number >
std::vector<number> TridiagonalMatrix< number >::diagonal
private

The diagonal entries.

Definition at line 235 of file tridiagonal_matrix.h.

◆ left

template<typename number >
std::vector<number> TridiagonalMatrix< number >::left
private

The entries left of the diagonal. The entry with index zero is always zero, since the first row has no entry left of the diagonal. Therefore, the length of this vector is the same as that of diagonal.

The length of this vector is zero for symmetric storage. In this case, the second element of left is identified with the first element of right.

Definition at line 245 of file tridiagonal_matrix.h.

◆ right

template<typename number >
std::vector<number> TridiagonalMatrix< number >::right
private

The entries right of the diagonal. The last entry is always zero, since the last row has no entry right of the diagonal. Therefore, the length of this vector is the same as that of diagonal.

Definition at line 251 of file tridiagonal_matrix.h.

◆ is_symmetric

template<typename number >
bool TridiagonalMatrix< number >::is_symmetric
private

If this flag is true, only the entries to the right of the diagonal are stored and the matrix is assumed symmetric.

Definition at line 257 of file tridiagonal_matrix.h.

◆ state

template<typename number >
LAPACKSupport::State TridiagonalMatrix< number >::state
private

The state of the matrix. Normally, the state is LAPACKSupport::matrix, indicating that the object can be used for regular matrix operations.

See explanation of this data type for details.

Definition at line 265 of file tridiagonal_matrix.h.


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