Dune::FMatrixHelp Namespace Reference

Functions
template<typename K >
static K	invertMatrix (const FieldMatrix< K, 1, 1 > &matrix, FieldMatrix< K, 1, 1 > &inverse)
	invert scalar without changing the original matrix
template<typename K >
static K	invertMatrix_retTransposed (const FieldMatrix< K, 1, 1 > &matrix, FieldMatrix< K, 1, 1 > &inverse)
	invert scalar without changing the original matrix
template<typename K >
static K	invertMatrix (const FieldMatrix< K, 2, 2 > &matrix, FieldMatrix< K, 2, 2 > &inverse)
	invert 2x2 Matrix without changing the original matrix
template<typename K >
static K	invertMatrix_retTransposed (const FieldMatrix< K, 2, 2 > &matrix, FieldMatrix< K, 2, 2 > &inverse)
template<typename K >
static K	invertMatrix (const FieldMatrix< K, 3, 3 > &matrix, FieldMatrix< K, 3, 3 > &inverse)
	invert 3x3 Matrix without changing the original matrix
template<typename K >
static K	invertMatrix_retTransposed (const FieldMatrix< K, 3, 3 > &matrix, FieldMatrix< K, 3, 3 > &inverse)
	invert 3x3 Matrix without changing the original matrix
template<class K , int m, int n, int p>
static void	multMatrix (const FieldMatrix< K, m, n > &A, const FieldMatrix< K, n, p > &B, FieldMatrix< K, m, p > &ret)
	calculates ret = A * B
template<typename K , int rows, int cols>
static void	multTransposedMatrix (const FieldMatrix< K, rows, cols > &matrix, FieldMatrix< K, cols, cols > &ret)
	calculates ret= A_t*A
template<typename K , int rows, int cols>
static void	multAssignTransposed (const FieldMatrix< K, rows, cols > &matrix, const FieldVector< K, rows > &x, FieldVector< K, cols > &ret)
	calculates ret = matrix^T * x
template<typename K , int rows, int cols>
static FieldVector< K, rows >	mult (const FieldMatrix< K, rows, cols > &matrix, const FieldVector< K, cols > &x)
	calculates ret = matrix * x
template<typename K , int rows, int cols>
static FieldVector< K, cols >	multTransposed (const FieldMatrix< K, rows, cols > &matrix, const FieldVector< K, rows > &x)
	calculates ret = matrix^T * x
void	eigenValuesLapackCall (const char *jobz, const char *uplo, const long int *n, double *a, const long int *lda, double *w, double *work, const long int *lwork, long int *info)
template<typename K >
static void	eigenValues (const FieldMatrix< K, 1, 1 > &matrix, FieldVector< K, 1 > &eigenvalues)
	calculates the eigenvalues of a symetric field matrix
template<typename K >
static void	eigenValues (const FieldMatrix< K, 2, 2 > &matrix, FieldVector< K, 2 > &eigenvalues)
	calculates the eigenvalues of a symetric field matrix
template<int dim, typename K >
static void	eigenValues (const FieldMatrix< K, dim, dim > &matrix, FieldVector< K, dim > &eigenvalues)
	calculates the eigenvalues of a symetric field matrix

---

Function Documentation

template<typename K >

static void Dune::FMatrixHelp::eigenValues	(	const FieldMatrix< K, 1, 1 > &	matrix,
		FieldVector< K, 1 > &	eigenvalues
	)		[static]

calculates the eigenvalues of a symetric field matrix

Parameters:

[in]	matrix	matrix eigenvalues are calculated for
[out]	eigenvalues	FieldVector that contains eigenvalues in ascending order

template<typename K >

static void Dune::FMatrixHelp::eigenValues	(	const FieldMatrix< K, 2, 2 > &	matrix,
		FieldVector< K, 2 > &	eigenvalues
	)		[static]

calculates the eigenvalues of a symetric field matrix

Parameters:

[in]	matrix	matrix eigenvalues are calculated for
[out]	eigenvalues	FieldVector that contains eigenvalues in ascending order

template<int dim, typename K >

static void Dune::FMatrixHelp::eigenValues	(	const FieldMatrix< K, dim, dim > &	matrix,
		FieldVector< K, dim > &	eigenvalues
	)		[static]

calculates the eigenvalues of a symetric field matrix

Parameters:

[in]	matrix	matrix eigenvalues are calculated for
[out]	eigenvalues	FieldVector that contains eigenvalues in ascending order

Note:

LAPACK::dsyev is used to calculate the eigen values

References DUNE_THROW, and eigenValuesLapackCall().

void Dune::FMatrixHelp::eigenValuesLapackCall	(	const char *	jobz,
		const char *	uplo,
		const long int *	n,
		double *	a,
		const long int *	lda,
		double *	w,
		double *	work,
		const long int *	lwork,
		long int *	info
	)

References DUNE_THROW.

Referenced by eigenValues().

template<typename K >

static K Dune::FMatrixHelp::invertMatrix	(	const FieldMatrix< K, 1, 1 > &	matrix,
		FieldMatrix< K, 1, 1 > &	inverse
	)		[inline, static]

invert scalar without changing the original matrix

Referenced by invertMatrix_retTransposed().

template<typename K >

static K Dune::FMatrixHelp::invertMatrix	(	const FieldMatrix< K, 2, 2 > &	matrix,
		FieldMatrix< K, 2, 2 > &	inverse
	)		[inline, static]

invert 2x2 Matrix without changing the original matrix

template<typename K >

static K Dune::FMatrixHelp::invertMatrix	(	const FieldMatrix< K, 3, 3 > &	matrix,
		FieldMatrix< K, 3, 3 > &	inverse
	)		[inline, static]

invert 3x3 Matrix without changing the original matrix

template<typename K >

static K Dune::FMatrixHelp::invertMatrix_retTransposed	(	const FieldMatrix< K, 1, 1 > &	matrix,
		FieldMatrix< K, 1, 1 > &	inverse
	)		[inline, static]

invert scalar without changing the original matrix

References invertMatrix().

template<typename K >

static K Dune::FMatrixHelp::invertMatrix_retTransposed	(	const FieldMatrix< K, 2, 2 > &	matrix,
		FieldMatrix< K, 2, 2 > &	inverse
	)		[inline, static]

invert 2x2 Matrix without changing the original matrix return transposed matrix

template<typename K >

static K Dune::FMatrixHelp::invertMatrix_retTransposed	(	const FieldMatrix< K, 3, 3 > &	matrix,
		FieldMatrix< K, 3, 3 > &	inverse
	)		[inline, static]

invert 3x3 Matrix without changing the original matrix

template<typename K , int rows, int cols>

static FieldVector<K,rows> Dune::FMatrixHelp::mult	(	const FieldMatrix< K, rows, cols > &	matrix,
		const FieldVector< K, cols > &	x
	)		[inline, static]

calculates ret = matrix * x

References Dune::DenseMatrixHelp::multAssign().

template<typename K , int rows, int cols>

static void Dune::FMatrixHelp::multAssignTransposed	(	const FieldMatrix< K, rows, cols > &	matrix,
		const FieldVector< K, rows > &	x,
		FieldVector< K, cols > &	ret
	)		[inline, static]

calculates ret = matrix^T * x

Referenced by multTransposed().

template<class K , int m, int n, int p>

static void Dune::FMatrixHelp::multMatrix	(	const FieldMatrix< K, m, n > &	A,
		const FieldMatrix< K, n, p > &	B,
		FieldMatrix< K, m, p > &	ret
	)		[inline, static]

calculates ret = A * B

template<typename K , int rows, int cols>

static FieldVector<K,cols> Dune::FMatrixHelp::multTransposed	(	const FieldMatrix< K, rows, cols > &	matrix,
		const FieldVector< K, rows > &	x
	)		[inline, static]

calculates ret = matrix^T * x

References multAssignTransposed().

template<typename K , int rows, int cols>

static void Dune::FMatrixHelp::multTransposedMatrix	(	const FieldMatrix< K, rows, cols > &	matrix,
		FieldMatrix< K, cols, cols > &	ret
	)		[inline, static]

calculates ret= A_t*A