Skip to content
Snippets Groups Projects
Commit f235807a authored by Klaus Böhnlein's avatar Klaus Böhnlein
Browse files

remove sym-function as its already in dunge-gfe/linearalgebra

parent 48e0139a
No related branches found
No related tags found
1 merge request!15Move to new Hermite-basis from dune-functions.
Pipeline #17198 failed
......@@ -319,26 +319,26 @@ auto compare(LocalDiscreteKirchhoffFunction& kirchhoffFunction,
/**
* @brief Symmetric part of a matrix.
*
* @param M
* @return FieldMatrix<RT,2,2>
*/
template<class RT>
static Dune::FieldMatrix<RT,2,2> sym (Dune::FieldMatrix<RT,2,2> M) { // 1/2 (M^T + M)
Dune::FieldMatrix<RT,2,2> ret(0);
for (int i=0; i<2; i++)
{
ret[i][i] = M[i][i];
for (int j=i+1; j<2; j++)
{
ret[i][j] = 0.5*(M[i][j] + M[j][i]);
ret[j][i] = ret[i][j];
}
}
return ret;
}
// /**
// * @brief Symmetric part of a matrix.
// *
// * @param M
// * @return FieldMatrix<RT,2,2>
// */
// template<class RT>
// static Dune::FieldMatrix<RT,2,2> sym (Dune::FieldMatrix<RT,2,2> M) { // 1/2 (M^T + M)
// Dune::FieldMatrix<RT,2,2> ret(0);
// for (int i=0; i<2; i++)
// {
// ret[i][i] = M[i][i];
// for (int j=i+1; j<2; j++)
// {
// ret[i][j] = 0.5*(M[i][j] + M[j][i]);
// ret[j][i] = ret[i][j];
// }
// }
// return ret;
// }
......
......@@ -9,6 +9,7 @@
#include <dune/gfe/assemblers/localenergy.hh>
#include <dune/gfe/bendingisometryhelper.hh>
#include <dune/gfe/tensor3.hh>
#include <dune/gfe/linearalgebra.hh>
#include <dune/localfunctions/lagrange/lagrangesimplex.hh>
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment