Add comprehensive test for symmetric tensor constructor and initializer

diff --git a/tests/base/symmetric_tensor_42.cc b/tests/base/symmetric_tensor_42.cc
new file mode 100644 (file)
index 0000000..9e95cd8
--- /dev/null
+++ b/tests/base/symmetric_tensor_42.cc
@@ -0,0 +1,68 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+// Verify that symmetric tensor construction and initialization
+// works as expected for non-trivial number types
+
+#include "../tests.h"
+#include <deal.II/base/symmetric_tensor.h>
+
+#include<complex>
+
+template <int dim, int rank,  typename NumberType>
+void test_constructor ()
+{
+  const SymmetricTensor<rank,dim,NumberType> t1;
+}
+
+template <int dim, int rank, typename NumberType, typename OtherNumberType>
+void test_initializer ()
+{
+  const SymmetricTensor<rank,dim,OtherNumberType> t1;
+  const SymmetricTensor<rank,dim,NumberType> t2 (t1);
+}
+
+template <int rank, typename NumberType, typename OtherNumberType>
+void test ()
+{
+  deallog << "constructor" << std::endl;
+  test_constructor<1,rank,NumberType> ();
+  test_constructor<2,rank,NumberType> ();
+  test_constructor<3,rank,NumberType> ();
+
+  deallog << "initializer" << std::endl;
+  test_initializer<1,rank,NumberType,OtherNumberType> ();
+  test_initializer<2,rank,NumberType,OtherNumberType> ();
+  test_initializer<3,rank,NumberType,OtherNumberType> ();
+}
+
+
+int main ()
+{
+  initlog();
+
+  deallog << "rank 2, double" << std::endl;
+  test<2, double, float >();
+  deallog << "rank 2, std::complex<double>" << std::endl;
+  test<2, std::complex<double>, std::complex<float> >();
+
+  deallog << "rank 4, double" << std::endl;
+  test<4, double, float >();
+  deallog << "rank 4, std::complex<double>" << std::endl;
+  test<4, std::complex<double>, std::complex<float> >();
+
+  deallog << "OK" << std::endl;
+}

diff --git a/tests/base/symmetric_tensor_42.output b/tests/base/symmetric_tensor_42.output
new file mode 100644 (file)
index 0000000..19ec964
--- /dev/null
+++ b/tests/base/symmetric_tensor_42.output
@@ -0,0 +1,14 @@
+
+DEAL::rank 2, double
+DEAL::constructor
+DEAL::initializer
+DEAL::rank 2, std::complex<double>
+DEAL::constructor
+DEAL::initializer
+DEAL::rank 4, double
+DEAL::constructor
+DEAL::initializer
+DEAL::rank 4, std::complex<double>
+DEAL::constructor
+DEAL::initializer
+DEAL::OK