]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
bug fix: symmetrize the tensor of derivatives, not an intermediate one
authorleicht <leicht@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 3 Aug 2006 07:14:20 +0000 (07:14 +0000)
committerleicht <leicht@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 3 Aug 2006 07:14:20 +0000 (07:14 +0000)
git-svn-id: https://svn.dealii.org/trunk@13581 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/numerics/derivative_approximation.cc

index d5188042c498038757d0416e17eefd3f057463b6..f33de9df38a5d9767173db905274a435be2ec98f 100644 (file)
@@ -748,14 +748,14 @@ DerivativeApproximation::approximate (const Mapping<dim>    &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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.