From cbd4f2f01a0a1f15096ba20f61225ec6d8a86b43 Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Wed, 16 Aug 2017 13:35:24 -0600 Subject: [PATCH] Generalised symmetric tensor operator* (tensor contrac.) on number type --- include/deal.II/base/symmetric_tensor.h | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/include/deal.II/base/symmetric_tensor.h b/include/deal.II/base/symmetric_tensor.h index 4238592d42..98fde9de2b 100644 --- a/include/deal.II/base/symmetric_tensor.h +++ b/include/deal.II/base/symmetric_tensor.h @@ -3895,12 +3895,12 @@ double_contract (SymmetricTensor<2,3,typename ProductType:: * @relates SymmetricTensor * @author Wolfgang Bangerth, 2005 */ -template -Tensor<1,dim,Number> +template +Tensor<1,dim,typename ProductType::type> operator * (const SymmetricTensor<2,dim,Number> &src1, - const Tensor<1,dim,Number> &src2) + const Tensor<1,dim,OtherNumber> &src2) { - Tensor<1,dim,Number> dest; + Tensor<1,dim,typename ProductType::type> dest; for (unsigned int i=0; i &src1, * @relates SymmetricTensor * @author Wolfgang Bangerth, 2005 */ -template -Tensor<1,dim,Number> -operator * (const Tensor<1,dim,Number> &src1, - const SymmetricTensor<2,dim,Number> &src2) +template +Tensor<1,dim,typename ProductType::type> +operator * (const Tensor<1,dim,Number> &src1, + const SymmetricTensor<2,dim,OtherNumber> &src2) { // this is easy for symmetric tensors: return src2 * src1; -- 2.39.5