From 05fdbbf6c198babe132f2bc1fc2e410602f8c96e Mon Sep 17 00:00:00 2001
From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Sun, 12 Mar 2017 19:59:38 +0100
Subject: [PATCH] Test that norms on tensors of VectorizedArray work

---
 tests/base/vectorization_08.cc     | 38 ++++++++++++++++++++++++++++++
 tests/base/vectorization_08.output |  4 ++++
 2 files changed, 42 insertions(+)
 create mode 100644 tests/base/vectorization_08.cc
 create mode 100644 tests/base/vectorization_08.output

diff --git a/tests/base/vectorization_08.cc b/tests/base/vectorization_08.cc
new file mode 100644
index 0000000000..71d5099d41
--- /dev/null
+++ b/tests/base/vectorization_08.cc
@@ -0,0 +1,38 @@
+// ---------------------------------------------------------------------
+//
+// 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.
+//
+// ---------------------------------------------------------------------
+
+// check that operations on tensors of vectorized arrays are properly
+// supported
+
+#include "../tests.h"
+
+#include <deal.II/base/point.h>
+#include <deal.II/base/symmetric_tensor.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/base/vectorization.h>
+
+int main()
+{
+  initlog();
+
+  Tensor<1,2,VectorizedArray<double> > t;
+  deallog << "Tensor norm: " << t.norm()[0] << std::endl;
+  Point<3,VectorizedArray<float> > p;
+  deallog << "Point norm: " << p.norm()[0] << std::endl;
+  SymmetricTensor<2,3,VectorizedArray<double> > st2;
+  SymmetricTensor<4,2,VectorizedArray<double> > st4;
+  deallog << "Symmetric tensor norm: " << st2.norm()[0]
+          << " " << st4.norm()[0] << std::endl;
+}
diff --git a/tests/base/vectorization_08.output b/tests/base/vectorization_08.output
new file mode 100644
index 0000000000..f85ab0ceae
--- /dev/null
+++ b/tests/base/vectorization_08.output
@@ -0,0 +1,4 @@
+
+DEAL::Tensor norm: 0.00000
+DEAL::Point norm: 0
+DEAL::Symmetric tensor norm: 0.00000 0.00000
-- 
2.39.5