Skip to content
Snippets Groups Projects
Commit 3ad0e2c6 authored by Lisa Julia Nebel's avatar Lisa Julia Nebel
Browse files

Fix error in nonplanarcosseratshellenergy

The contravariant base vectors were not calculated correctly.
The contravariant base vectors are the *columns* of the inverse of the covariant matrix, not the rows.
To fix this, take the rows of the transpose of inverse of the covariant matrix.
parent 4bebd0ef
Branches
Tags
1 merge request!80Fix error in nonplanarcosseratshellenergy
......@@ -230,6 +230,9 @@ energy(const typename Basis::LocalView& localView,
auto aContravariant = aCovariant;
aContravariant.invert();
// The contravariant base vectors are the *columns* of the inverse of the covariant matrix
// To get an easier access to the columns, we use the transpose of the contravariant matrix
aContravariant = Dune::GFE::transpose(aContravariant);
Dune::FieldMatrix<double,3,3> a(0);
for (int alpha=0; alpha<gridDim; alpha++)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment