--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2019 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.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+
+// test compare_and_apply_mask
+
+#include <deal.II/base/vectorization.h>
+
+#include "../tests.h"
+
+
+template <typename VectorizedArrayType>
+void
+do_test()
+{
+ deallog << " test " << VectorizedArrayType::n_array_elements
+ << " array elements" << std::endl;
+
+ VectorizedArrayType left;
+ for (unsigned int i = 0; i < VectorizedArrayType::n_array_elements; i++)
+ left[i] = i + 1.;
+
+ VectorizedArrayType right(3.);
+ VectorizedArrayType true_values(1.);
+ VectorizedArrayType false_values(-1.);
+
+ auto result = compare_and_apply_mask<SIMDComparison::equal>(left,
+ right,
+ true_values,
+ false_values);
+ deallog << result << std::endl;
+
+ result = compare_and_apply_mask<SIMDComparison::not_equal>(left,
+ right,
+ true_values,
+ false_values);
+ deallog << result << std::endl;
+
+ result = compare_and_apply_mask<SIMDComparison::less_than>(left,
+ right,
+ true_values,
+ false_values);
+ deallog << result << std::endl;
+
+ result = compare_and_apply_mask<SIMDComparison::less_than_or_equal>(
+ left, right, true_values, false_values);
+ deallog << result << std::endl;
+
+ result = compare_and_apply_mask<SIMDComparison::greater_than>(left,
+ right,
+ true_values,
+ false_values);
+ deallog << result << std::endl;
+
+ result = compare_and_apply_mask<SIMDComparison::greater_than_or_equal>(
+ left, right, true_values, false_values);
+ deallog << result << std::endl;
+}
+
+
+int
+main()
+{
+ initlog();
+
+#if DEAL_II_COMPILER_VECTORIZATION_LEVEL >= 3 && defined(__AVX512F__)
+ do_test<VectorizedArray<double, 8>>();
+ do_test<VectorizedArray<float, 16>>();
+#endif
+
+#if DEAL_II_COMPILER_VECTORIZATION_LEVEL >= 2 && defined(__AVX__)
+ do_test<VectorizedArray<double, 4>>();
+ do_test<VectorizedArray<float, 8>>();
+#endif
+
+#if DEAL_II_COMPILER_VECTORIZATION_LEVEL >= 1 && defined(__SSE2__)
+ do_test<VectorizedArray<double, 2>>();
+ do_test<VectorizedArray<float, 4>>();
+#endif
+
+ do_test<VectorizedArray<double, 1>>();
+ do_test<VectorizedArray<float, 1>>();
+}
--- /dev/null
+
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
--- /dev/null
+
+DEAL:: test 4 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000
+DEAL:: test 8 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL:: test 2 array elements
+DEAL::-1.00000 -1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::-1.00000 -1.00000
+DEAL::-1.00000 -1.00000
+DEAL:: test 4 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
--- /dev/null
+
+DEAL:: test 8 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL:: test 16 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL:: test 4 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000
+DEAL:: test 8 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
+DEAL:: test 2 array elements
+DEAL::-1.00000 -1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::-1.00000 -1.00000
+DEAL::-1.00000 -1.00000
+DEAL:: test 4 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
--- /dev/null
+
+DEAL:: test 2 array elements
+DEAL::-1.00000 -1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000
+DEAL::-1.00000 -1.00000
+DEAL::-1.00000 -1.00000
+DEAL:: test 4 array elements
+DEAL::-1.00000 -1.00000 1.00000 -1.00000
+DEAL::1.00000 1.00000 -1.00000 1.00000
+DEAL::1.00000 1.00000 -1.00000 -1.00000
+DEAL::1.00000 1.00000 1.00000 -1.00000
+DEAL::-1.00000 -1.00000 -1.00000 1.00000
+DEAL::-1.00000 -1.00000 1.00000 1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000
+DEAL:: test 1 array elements
+DEAL::-1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::1.00000
+DEAL::-1.00000
+DEAL::-1.00000