From be2b9873dd5400396a66c4c693226a63db124f90 Mon Sep 17 00:00:00 2001 From: wolf Date: Sat, 28 May 2005 17:01:48 +0000 Subject: [PATCH] Fix two problems that made this not compile. git-svn-id: https://svn.dealii.org/trunk@10769 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/include/base/symmetric_tensor.h | 29 ++++++++++++++++++++ 1 file changed, 29 insertions(+) diff --git a/deal.II/base/include/base/symmetric_tensor.h b/deal.II/base/include/base/symmetric_tensor.h index 60fc732ed1..97fd8c4ec8 100644 --- a/deal.II/base/include/base/symmetric_tensor.h +++ b/deal.II/base/include/base/symmetric_tensor.h @@ -21,6 +21,7 @@ template class SymmetricTensor; template SymmetricTensor<2,dim> unit_symmetric_tensor (); +template SymmetricTensor<4,dim> deviator_tensor (); namespace internal { @@ -878,6 +879,9 @@ class SymmetricTensor template friend SymmetricTensor<2,dim2> unit_symmetric_tensor (); + + template + friend SymmetricTensor<4,dim2> deviator_tensor (); }; @@ -2033,6 +2037,31 @@ unit_symmetric_tensor<3> () + +/** + * Return a unit symmetric tensor of rank 2 and dimension 2. + * + * @relates SymmetricTensor + * @author Wolfgang Bangerth, 2005 + */ +template +inline +SymmetricTensor<4,dim> +deviator_tensor () +{ + SymmetricTensor<4,dim> tmp; + for (unsigned int i=0; i::n_rank2_components; ++i) + tmp.data[i][i] = 1; + + return tmp; +} + + + /** * Return the symmetrized version of the full rank-2 tensor, * i.e. (t+transpose(t))/2, as a symmetric rank-2 tensor. -- 2.39.5