From 9bd869e7d6d36a2a1c98728ad9b7a4e994e5f02d Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 16 Sep 2015 09:02:40 -0500 Subject: [PATCH] Use the new variants of contract(). --- source/fe/mapping_q_generic.cc | 9 ++------- 1 file changed, 2 insertions(+), 7 deletions(-) diff --git a/source/fe/mapping_q_generic.cc b/source/fe/mapping_q_generic.cc index 5429fc3b4f..2cc850d8eb 100644 --- a/source/fe/mapping_q_generic.cc +++ b/source/fe/mapping_q_generic.cc @@ -1185,9 +1185,8 @@ namespace } // Solve [f'(x)]d=f(x) - Tensor<1,spacedim> delta; Tensor<2,spacedim> df_inverse = invert(df); - contract (delta, df_inverse, static_cast&>(f)); + const Tensor<1,spacedim> delta = df_inverse * static_cast&>(f); #ifdef DEBUG_TRANSFORM_REAL_TO_UNIT_CELL std::cout << " delta=" << delta << std::endl; @@ -1326,11 +1325,7 @@ namespace while (f.norm()>eps && loop++ d; - Tensor<2,dim> df_1; - - df_1 = invert(df); - contract (d, df_1, static_cast&>(f)); + const Tensor<1,dim> d = invert(df) * static_cast&>(f); p_unit -= d; for (unsigned int j=0; j