From: wolf Date: Sat, 28 May 2005 19:05:24 +0000 (+0000) Subject: Fix a bug in the deviator. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7e19d700094fa3dc72156a4362a6efbd639604f3;p=dealii-svn.git Fix a bug in the deviator. git-svn-id: https://svn.dealii.org/trunk@10771 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/include/base/symmetric_tensor.h b/deal.II/base/include/base/symmetric_tensor.h index 97fd8c4ec8..68f93150e2 100644 --- a/deal.II/base/include/base/symmetric_tensor.h +++ b/deal.II/base/include/base/symmetric_tensor.h @@ -2050,12 +2050,22 @@ SymmetricTensor<4,dim> deviator_tensor () { SymmetricTensor<4,dim> tmp; + + // fill the elements treating the diagonal for (unsigned int i=0; i::n_rank2_components; ++i) - tmp.data[i][i] = 1; + i + ::n_rank2_components; + ++i) + tmp.data[i][i] = 0.5; return tmp; }