From: Jean-Paul Pelteret <jppelteret@gmail.com>
Date: Tue, 18 Sep 2018 12:59:32 +0000 (+0200)
Subject: Improve in-code documentation of SymmetricTensor<2,3> inverse.
X-Git-Tag: v9.1.0-rc1~705^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F7209%2Fhead;p=dealii.git

Improve in-code documentation of SymmetricTensor<2,3> inverse.
---

diff --git a/include/deal.II/base/symmetric_tensor.h b/include/deal.II/base/symmetric_tensor.h
index 634c36ceac..64e4105f82 100644
--- a/include/deal.II/base/symmetric_tensor.h
+++ b/include/deal.II/base/symmetric_tensor.h
@@ -1227,21 +1227,38 @@ namespace internal
 
         // Sympy result: ([
         // [  (t11*t22 - t12**2)/(t00*t11*t22 - t00*t12**2 - t01**2*t22 +
-        // 2*t01*t02*t12 - t02**2*t11),
+        //                        2*t01*t02*t12 - t02**2*t11),
         //    (-t01*t22 + t02*t12)/(t00*t11*t22 - t00*t12**2 - t01**2*t22 +
-        //    2*t01*t02*t12 - t02**2*t11), (t01*t12 - t02*t11)/(t00*t11*t22 -
-        //    t00*t12**2 - t01**2*t22 + 2*t01*t02*t12 - t02**2*t11)],
+        //                          2*t01*t02*t12 - t02**2*t11),
+        //    (t01*t12 - t02*t11)/(t00*t11*t22 -  t00*t12**2 - t01**2*t22 +
+        //                         2*t01*t02*t12 - t02**2*t11)],
         // [  (-t01*t22 + t02*t12)/(t00*t11*t22 - t00*t12**2 - t01**2*t22 +
-        // 2*t01*t02*t12 - t02**2*t11),
+        //                          2*t01*t02*t12 - t02**2*t11),
         //    (t00*t22 - t02**2)/(t00*t11*t22 - t00*t12**2 - t01**2*t22 +
-        //    2*t01*t02*t12 - t02**2*t11), (t00*t12 - t01*t02)/(-t00*t11*t22 +
-        //    t00*t12**2 + t01**2*t22 - 2*t01*t02*t12 + t02**2*t11)],
+        //                        2*t01*t02*t12 - t02**2*t11),
+        //    (t00*t12 - t01*t02)/(-t00*t11*t22 + t00*t12**2 + t01**2*t22 -
+        //                         2*t01*t02*t12 + t02**2*t11)],
         // [  (t01*t12 - t02*t11)/(t00*t11*t22 - t00*t12**2 - t01**2*t22 +
-        // 2*t01*t02*t12 - t02**2*t11),
+        //                         2*t01*t02*t12 - t02**2*t11),
         //    (t00*t12 - t01*t02)/(-t00*t11*t22 + t00*t12**2 + t01**2*t22 -
-        //    2*t01*t02*t12 + t02**2*t11),
+        //                         2*t01*t02*t12 + t02**2*t11),
         //    (-t00*t11 + t01**2)/(-t00*t11*t22 + t00*t12**2 + t01**2*t22 -
-        //    2*t01*t02*t12 + t02**2*t11)]  ])
+        //                         2*t01*t02*t12 + t02**2*t11)]  ])
+        //
+        // =
+        //
+        // [  (t11*t22 - t12**2)/det_t,
+        //    (-t01*t22 + t02*t12)/det_t,
+        //    (t01*t12 - t02*t11)/det_t],
+        // [  (-t01*t22 + t02*t12)/det_t,
+        //    (t00*t22 - t02**2)/det_t,
+        //    (-t00*t12 + t01*t02)/det_t],
+        // [  (t01*t12 - t02*t11)/det_t,
+        //    (-t00*t12 + t01*t02)/det_t,
+        //    (t00*t11 - t01**2)/det_t]   ])
+        //
+        // with det_t = (t00*t11*t22 - t00*t12**2 - t01**2*t22 +
+        //               2*t01*t02*t12 - t02**2*t11)
         const TableIndices<2> idx_00(0, 0);
         const TableIndices<2> idx_01(0, 1);
         const TableIndices<2> idx_02(0, 2);