From 3fc6504771b1ebfca1e48080b65bed0983038107 Mon Sep 17 00:00:00 2001 From: leicht Date: Thu, 3 Aug 2006 07:14:20 +0000 Subject: [PATCH] bug fix: symmetrize the tensor of derivatives, not an intermediate one git-svn-id: https://svn.dealii.org/trunk@13581 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/derivative_approximation.cc | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/deal.II/deal.II/source/numerics/derivative_approximation.cc b/deal.II/deal.II/source/numerics/derivative_approximation.cc index d5188042c4..f33de9df38 100644 --- a/deal.II/deal.II/source/numerics/derivative_approximation.cc +++ b/deal.II/deal.II/source/numerics/derivative_approximation.cc @@ -748,14 +748,14 @@ DerivativeApproximation::approximate (const Mapping &mapping, AssertThrow (determinant(Y) != 0, ExcInsufficientDirections()); - // first symmetrize g - DerivativeDescription::symmetrize (projected_derivative); - // compute Y^-1 g const Tensor<2,dim> Y_inverse = invert(Y); typename DerivativeDescription::Derivative derivative; contract (derivative, Y_inverse, projected_derivative); + + // finally symmetrize the derivative + DerivativeDescription::symmetrize (derivative); *derivative_norm_on_this_cell = DerivativeDescription::derivative_norm (derivative); -- 2.39.5