From: Jean-Paul Pelteret 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=ec76803879462bf0b1f05a8dd8cec36165c44f98;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);