tridiaginv.h File Reference

#include <boost/math/special_functions.hpp>
#include "lapacke.h"
#include "dg/algorithm.h"
#include "functors.h"

Go to the source code of this file.

Classes
class	dg::mat::TridiagInvHMGTI< real_type >
	Compute the inverse of a general tridiagonal matrix. More...
class	dg::mat::TridiagInvDF< real_type >
	USE THIS ONE Compute the inverse of a general tridiagonal matrix. The algorithm does not rely on the determinant. More...
class	dg::mat::TridiagInvD< real_type >
	Compute the inverse of a general tridiagonal matrix. More...

Namespaces
namespace	dg
	Classes for Krylov space approximations of a Matrix-Vector product.
namespace	dg::mat

Functions
template<class value_type >
value_type	dg::mat::compute_Tinv_m1 (const dg::TriDiagonal< thrust::host_vector< value_type > > &T, unsigned size)
	Computes the value of $(T^{-1}{)}_{m1}=⟨{\overrightarrow{e}}_m,T^{-1}{\overrightarrow{e}}_1⟩$ via a Thomas algorithm.
template<class value_type >
void	dg::mat::compute_Tinv_y (const dg::TriDiagonal< thrust::host_vector< value_type > > &T, thrust::host_vector< value_type > &x, const thrust::host_vector< value_type > &y, value_type a=1., value_type d=0.)
	Computes the value of $x=((aT+dI{)}^{-1})y$ via Thomas algorithm.
template<class value_type >
void	dg::mat::invert (const dg::TriDiagonal< thrust::host_vector< value_type > > &T, dg::SquareMatrix< value_type > &Tinv)
	Invert a tridiagonal matrix.
template<class value_type >
dg::SquareMatrix< value_type >	dg::mat::invert (const dg::TriDiagonal< thrust::host_vector< value_type > > &T)
	Invert a tridiagonal matrix.
template<class value_type >
std::array< value_type, 2 >	dg::mat::compute_extreme_EV (const dg::TriDiagonal< thrust::host_vector< value_type > > &T)
	Compute extreme Eigenvalues of a symmetric tridiangular matrix.