From 8c2bb3ee6dd7efaa639703c314e109d538ae61b0 Mon Sep 17 00:00:00 2001
From: Martin Kronbichler <martin.kronbichler@rub.de>
Date: Mon, 19 Aug 2024 21:28:49 +0200
Subject: [PATCH] Test case

---
 .../matrix_free/fe_face_evaluation_vector.cc  | 56 +++++++++++++++++++
 .../fe_face_evaluation_vector.output          |  2 +
 2 files changed, 58 insertions(+)
 create mode 100644 tests/matrix_free/fe_face_evaluation_vector.cc
 create mode 100644 tests/matrix_free/fe_face_evaluation_vector.output

diff --git a/tests/matrix_free/fe_face_evaluation_vector.cc b/tests/matrix_free/fe_face_evaluation_vector.cc
new file mode 100644
index 0000000000..bbe14430bb
--- /dev/null
+++ b/tests/matrix_free/fe_face_evaluation_vector.cc
@@ -0,0 +1,56 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2024 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+
+
+// Ensure that FEFaceEvaluation can be collected in a vector
+
+#include <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+
+#include <deal.II/matrix_free/fe_evaluation.h>
+
+#include "../tests.h"
+
+
+int
+main()
+{
+  initlog();
+
+  constexpr int      dim = 2;
+  Triangulation<dim> tria;
+  GridGenerator::hyper_cube(tria);
+  FE_Q<dim>       fe(1);
+  DoFHandler<dim> dof_handler(tria);
+  dof_handler.distribute_dofs(fe);
+
+  MappingQ<dim> mapping(1);
+
+  MatrixFree<dim, double>                 matrix_free;
+  MatrixFree<dim, double>::AdditionalData mf_data;
+  mf_data.mapping_update_flags_boundary_faces = update_values;
+  matrix_free.reinit(
+    mapping, dof_handler, AffineConstraints<double>(), QGauss<1>(2), mf_data);
+  FEFaceEvaluation<dim, -1, 0, 1, double>              evaluator(matrix_free);
+  std::vector<FEFaceEvaluation<dim, -1, 0, 1, double>> evaluators(2, evaluator);
+  evaluators[0].reinit(0);
+  evaluators[1].reinit(0);
+
+  deallog << "OK" << std::endl;
+}
diff --git a/tests/matrix_free/fe_face_evaluation_vector.output b/tests/matrix_free/fe_face_evaluation_vector.output
new file mode 100644
index 0000000000..0fd8fc12f0
--- /dev/null
+++ b/tests/matrix_free/fe_face_evaluation_vector.output
@@ -0,0 +1,2 @@
+
+DEAL::OK
-- 
2.39.5