Commit 06384382 authored by akbib's avatar akbib Committed by akbib@FU-BERLIN.DE
Browse files

bugfix in linearisedDefDet: index was computed wrong

[[Imported from SVN: r11291]]
parent b8fb1fe2
...@@ -5,6 +5,8 @@ ...@@ -5,6 +5,8 @@
#include <dune/common/fmatrix.hh> #include <dune/common/fmatrix.hh>
#include <dune/fufem/symmetrictensor.hh>
namespace Dune { namespace Dune {
namespace Elasticity { namespace Elasticity {
...@@ -133,7 +135,7 @@ namespace Dune { ...@@ -133,7 +135,7 @@ namespace Dune {
linDet = 0; linDet = 0;
for (int i=0; i<2; i++) for (int i=0; i<2; i++)
for (int j=(i+1)%3; j<3; j++) { for (int j=(i+1)%3; j<3; j++) {
int k=(j+1)%3 + (j+1)/3; int k=(-(i+j))%3;
linDet[i][j] = u[j][k]*u[k][i] - u[j][i]*(1+u[k][k]); linDet[i][j] = u[j][k]*u[k][i] - u[j][i]*(1+u[k][k]);
linDet[j][i] = u[k][j]*u[i][k] - u[i][j]*(1+u[k][k]); linDet[j][i] = u[k][j]*u[i][k] - u[i][j]*(1+u[k][k]);
} }
...@@ -213,7 +215,7 @@ namespace Dune { ...@@ -213,7 +215,7 @@ namespace Dune {
// 1 1 1 0 --> 0 - 0 // 1 1 1 0 --> 0 - 0
// 1 1 1 1 --> 0 - 0 // 1 1 1 1 --> 0 - 0
} }
double det_val(const Dune::FieldMatrix<double,3,3>& u) { double det_val(const Dune::FieldMatrix<double,3,3>& u) {
return (1+u[0][0])*(1+u[1][1])*(1+u[2][2]) return (1+u[0][0])*(1+u[1][1])*(1+u[2][2])
+ u[0][1]*u[1][2]*u[2][0] + u[0][1]*u[1][2]*u[2][0]
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment