From: leicht Date: Thu, 3 Aug 2006 07:14:20 +0000 (+0000) Subject: bug fix: symmetrize the tensor of derivatives, not an intermediate one X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3fc6504771b1ebfca1e48080b65bed0983038107;p=dealii-svn.git 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 --- 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);