From 1b82a7353fb6534c3e8af9d4ebfd4fbbe2404ef9 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Fri, 11 Feb 2011 03:40:45 +0000 Subject: [PATCH] New test. git-svn-id: https://svn.dealii.org/trunk@23328 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/base/symmetric_tensor_25.cc | 74 ++++++++++++++++++++++ tests/base/symmetric_tensor_25/cmp/generic | 52 +++++++++++++++ 2 files changed, 126 insertions(+) create mode 100644 tests/base/symmetric_tensor_25.cc create mode 100644 tests/base/symmetric_tensor_25/cmp/generic diff --git a/tests/base/symmetric_tensor_25.cc b/tests/base/symmetric_tensor_25.cc new file mode 100644 index 0000000000..487f89ee0c --- /dev/null +++ b/tests/base/symmetric_tensor_25.cc @@ -0,0 +1,74 @@ +//---------------------------- symmetric_tensor_25.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2008, 2009, 2011 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------- symmetric_tensor_25.cc --------------------------- + +// Test the invariants of tensors using the Cayley-Hamilton theorem, +// see http://en.wikipedia.org/wiki/Invariants_of_tensors + +#include "../tests.h" +#include +#include +#include +#include +#include + + +template +void check () +{ + for (unsigned int round=0; round<10; ++round) + { + SymmetricTensor<2,dim> S; + for (unsigned int i=0; i (); + + deallog << "OK" << std::endl; +} diff --git a/tests/base/symmetric_tensor_25/cmp/generic b/tests/base/symmetric_tensor_25/cmp/generic new file mode 100644 index 0000000000..5aec41ba54 --- /dev/null +++ b/tests/base/symmetric_tensor_25/cmp/generic @@ -0,0 +1,52 @@ + +DEAL::S = 3.00 5.00 3.00 5.00 6.00 5.00 3.00 5.00 7.00 +DEAL::first invariant = 16.0 +DEAL::second invariant = 22.0 +DEAL::third invariant = -28.0 +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 6.00 1.00 2.00 1.00 2.00 7.00 2.00 7.00 9.00 +DEAL::first invariant = 17.0 +DEAL::second invariant = 30.0 +DEAL::third invariant = -175. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 0.00 6.00 0.00 6.00 9.00 6.00 0.00 6.00 3.00 +DEAL::first invariant = 12.0 +DEAL::second invariant = -45.0 +DEAL::third invariant = -108. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 2.00 8.00 7.00 8.00 6.00 9.00 7.00 9.00 1.00 +DEAL::first invariant = 9.00 +DEAL::second invariant = -174. +DEAL::third invariant = 500. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 2.00 3.00 7.00 3.00 0.00 5.00 7.00 5.00 2.00 +DEAL::first invariant = 4.00 +DEAL::second invariant = -79.0 +DEAL::third invariant = 142. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 9.00 8.00 9.00 8.00 2.00 7.00 9.00 7.00 2.00 +DEAL::first invariant = 13.0 +DEAL::second invariant = -154. +DEAL::third invariant = 313. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 3.00 2.00 9.00 2.00 6.00 3.00 9.00 3.00 1.00 +DEAL::first invariant = 10.0 +DEAL::second invariant = -67.0 +DEAL::third invariant = -391. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 1.00 7.00 8.00 7.00 9.00 4.00 8.00 4.00 4.00 +DEAL::first invariant = 14.0 +DEAL::second invariant = -80.0 +DEAL::third invariant = -304. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 5.00 6.00 1.00 6.00 0.00 0.00 1.00 0.00 3.00 +DEAL::first invariant = 8.00 +DEAL::second invariant = -22.0 +DEAL::third invariant = -108. +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::S = 6.00 0.00 6.00 0.00 3.00 1.00 6.00 1.00 2.00 +DEAL::first invariant = 11.0 +DEAL::second invariant = -1.00 +DEAL::third invariant = -78.0 +DEAL::0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 +DEAL::OK -- 2.39.5