From: Matthias Maier Date: Thu, 28 Mar 2024 01:32:13 +0000 (-0500) Subject: base/tensor.h: make a statement easier to read X-Git-Tag: v9.6.0-rc1~442^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7f4eee490ae1b5c666fba1251e4994ab02bba654;p=dealii.git base/tensor.h: make a statement easier to read --- diff --git a/include/deal.II/base/tensor.h b/include/deal.II/base/tensor.h index 07e53f92c8..d1d8b5ad35 100644 --- a/include/deal.II/base/tensor.h +++ b/include/deal.II/base/tensor.h @@ -2887,34 +2887,29 @@ constexpr inline DEAL_II_ALWAYS_INLINE Tensor<2, 2, Number> return return_tensor; } - template constexpr inline DEAL_II_ALWAYS_INLINE Tensor<2, 3, Number> invert(const Tensor<2, 3, Number> &t) { Tensor<2, 3, Number> return_tensor; - return_tensor[0][0] = internal::NumberType::value(t[1][1] * t[2][2]) - - internal::NumberType::value(t[1][2] * t[2][1]); - return_tensor[0][1] = internal::NumberType::value(t[0][2] * t[2][1]) - - internal::NumberType::value(t[0][1] * t[2][2]); - return_tensor[0][2] = internal::NumberType::value(t[0][1] * t[1][2]) - - internal::NumberType::value(t[0][2] * t[1][1]); - return_tensor[1][0] = internal::NumberType::value(t[1][2] * t[2][0]) - - internal::NumberType::value(t[1][0] * t[2][2]); - return_tensor[1][1] = internal::NumberType::value(t[0][0] * t[2][2]) - - internal::NumberType::value(t[0][2] * t[2][0]); - return_tensor[1][2] = internal::NumberType::value(t[0][2] * t[1][0]) - - internal::NumberType::value(t[0][0] * t[1][2]); - return_tensor[2][0] = internal::NumberType::value(t[1][0] * t[2][1]) - - internal::NumberType::value(t[1][1] * t[2][0]); - return_tensor[2][1] = internal::NumberType::value(t[0][1] * t[2][0]) - - internal::NumberType::value(t[0][0] * t[2][1]); - return_tensor[2][2] = internal::NumberType::value(t[0][0] * t[1][1]) - - internal::NumberType::value(t[0][1] * t[1][0]); - const Number inv_det_t = internal::NumberType::value( - 1.0 / (t[0][0] * return_tensor[0][0] + t[0][1] * return_tensor[1][0] + - t[0][2] * return_tensor[2][0])); + const auto value = [](const auto &t) { + return internal::NumberType::value(t); + }; + + return_tensor[0][0] = value(t[1][1] * t[2][2]) - value(t[1][2] * t[2][1]); + return_tensor[0][1] = value(t[0][2] * t[2][1]) - value(t[0][1] * t[2][2]); + return_tensor[0][2] = value(t[0][1] * t[1][2]) - value(t[0][2] * t[1][1]); + return_tensor[1][0] = value(t[1][2] * t[2][0]) - value(t[1][0] * t[2][2]); + return_tensor[1][1] = value(t[0][0] * t[2][2]) - value(t[0][2] * t[2][0]); + return_tensor[1][2] = value(t[0][2] * t[1][0]) - value(t[0][0] * t[1][2]); + return_tensor[2][0] = value(t[1][0] * t[2][1]) - value(t[1][1] * t[2][0]); + return_tensor[2][1] = value(t[0][1] * t[2][0]) - value(t[0][0] * t[2][1]); + return_tensor[2][2] = value(t[0][0] * t[1][1]) - value(t[0][1] * t[1][0]); + + const Number inv_det_t = + value(1.0 / (t[0][0] * return_tensor[0][0] + t[0][1] * return_tensor[1][0] + + t[0][2] * return_tensor[2][0])); return_tensor *= inv_det_t; return return_tensor;