diff --git a/dune/common/dynmatrixev.hh b/dune/common/dynmatrixev.hh index 5ad91d34712ce57e3dfd7741698cbd107ceda71c..f5e63269f76553d25b5c9cd1b96360b6c8ba3789 100644 --- a/dune/common/dynmatrixev.hh +++ b/dune/common/dynmatrixev.hh @@ -16,6 +16,13 @@ namespace Dune { const long int* ldvl, double* vr, const long int* ldvr, double* work, const long int* lwork, const long int* info); + /** \brief calculates the eigenvalues of a symetric field matrix + \param[in] matrix matrix eigenvalues are calculated for + \param[out] eigenvalues FieldVector that contains eigenvalues in + ascending order + + \note LAPACK::dgeev is used to calculate the eigen values + */ template <typename K, class C> static void eigenValuesNonSym(const DynamicMatrix<K>& matrix, DynamicVector<C>& eigenValues) diff --git a/dune/common/fmatrixev.hh b/dune/common/fmatrixev.hh index cee1aa190ae14820d68b436d73a4d33a45a98f68..5b39e74ac1ab3d63f9e87e1287758842327c9cdc 100644 --- a/dune/common/fmatrixev.hh +++ b/dune/common/fmatrixev.hh @@ -127,7 +127,13 @@ namespace Dune { } } } + /** \brief calculates the eigenvalues of a symetric field matrix + \param[in] matrix matrix eigenvalues are calculated for + \param[out] eigenvalues FieldVector that contains eigenvalues in + ascending order + \note LAPACK::dgeev is used to calculate the eigen values + */ template <int dim, typename K, class C> static void eigenValuesNonSym(const FieldMatrix<K, dim, dim>& matrix, FieldVector<C, dim>& eigenValues)