From b6396a20a39aa4152e8b0c225e4293aa3c6463ce Mon Sep 17 00:00:00 2001 From: wolf Date: Wed, 30 Mar 2005 17:30:10 +0000 Subject: [PATCH] operator* doesn't do complete contractions of tensors of rank 4 any more. git-svn-id: https://svn.dealii.org/trunk@10321 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/base/symmetric_tensor_03.cc | 9 ++++++++- tests/base/symmetric_tensor_04.cc | 8 +++++++- 2 files changed, 15 insertions(+), 2 deletions(-) diff --git a/tests/base/symmetric_tensor_03.cc b/tests/base/symmetric_tensor_03.cc index 72d35a9b45..5c43cd5fc1 100644 --- a/tests/base/symmetric_tensor_03.cc +++ b/tests/base/symmetric_tensor_03.cc @@ -48,7 +48,14 @@ int main () // make sure norm is induced by scalar // product - Assert (std::fabs (t.norm()*t.norm() - t*t) < 1e-14, + double norm_sqr = 0; + for (unsigned int i=0; i<2; ++i) + for (unsigned int j=0; j<2; ++j) + for (unsigned int k=0; k<2; ++k) + for (unsigned int l=0; l<2; ++l) + norm_sqr += t[i][j][k][l] * t[i][j][k][l]; + + Assert (std::fabs (t.norm()*t.norm() - norm_sqr) < 1e-14, ExcInternalError()); deallog << "OK" << std::endl; diff --git a/tests/base/symmetric_tensor_04.cc b/tests/base/symmetric_tensor_04.cc index 9a6af12cd7..93e16f6fe2 100644 --- a/tests/base/symmetric_tensor_04.cc +++ b/tests/base/symmetric_tensor_04.cc @@ -48,7 +48,13 @@ int main () // make sure norm is induced by scalar // product - Assert (std::fabs (t.norm()*t.norm() - t*t) < 1e-14, + double norm_sqr = 0; + for (unsigned int i=0; i<2; ++i) + for (unsigned int j=0; j<2; ++j) + for (unsigned int k=0; k<2; ++k) + for (unsigned int l=0; l<2; ++l) + norm_sqr += t[i][j][k][l] * t[i][j][k][l]; + Assert (std::fabs (t.norm()*t.norm() - norm_sqr) < 1e-14, ExcInternalError()); deallog << "OK" << std::endl; -- 2.39.5