Fix SymmetricTensor::deviator for VectorizedArray

diff --git a/include/deal.II/base/symmetric_tensor.h b/include/deal.II/base/symmetric_tensor.h
index 4a8397a6d0c1b8957a7d1b2d2bf7f84f9f6712a2..a9c7da529ba84a990a6fa5ac38477eb54c5f6301 100644 (file)
--- a/include/deal.II/base/symmetric_tensor.h
+++ b/include/deal.II/base/symmetric_tensor.h
@@ -3311,7 +3311,7 @@ deviator(const SymmetricTensor<2, dim, Number> &t)
   SymmetricTensor<2, dim, Number> tmp = t;
 
   // subtract scaled trace from the diagonal
-  const Number tr = trace(t) / dim;
+  const Number tr = trace(t) * internal::NumberType<Number>::value(1.0 / dim);
   for (unsigned int i = 0; i < dim; ++i)
     tmp.data[i] -= tr;