From faec3374ac323ab6d3f54ddc8ab4e757b56196f1 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Fri, 20 Feb 2015 17:38:13 -0600
Subject: [PATCH] Add a testcase.

---
 tests/fe/fe_system_from_list.cc     | 151 ++++++++++++++++++++++++++++
 tests/fe/fe_system_from_list.output |  19 ++++
 2 files changed, 170 insertions(+)
 create mode 100644 tests/fe/fe_system_from_list.cc
 create mode 100644 tests/fe/fe_system_from_list.output

diff --git a/tests/fe/fe_system_from_list.cc b/tests/fe/fe_system_from_list.cc
new file mode 100644
index 0000000000..ab903ec9ba
--- /dev/null
+++ b/tests/fe/fe_system_from_list.cc
@@ -0,0 +1,151 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 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 the idea discussed in the multiple-element constructor FESystem
+// really works
+
+#include "../tests.h"
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+#include <deal.II/fe/fe_nedelec.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+
+#include <vector>
+#include <fstream>
+#include <string>
+
+#define PRECISION 2
+
+template <int dim>
+struct MyFESystem
+{
+  MyFESystem (const std::vector<const FiniteElement<dim>*> &fes,
+	      const std::vector<unsigned int>              &multiplicities)
+    {
+      deallog << "Constructing FESystem from list." << std::endl;
+    }
+};
+
+  
+
+template <int dim>
+class MySimulator {
+public:
+  MySimulator (const unsigned int polynomial_degree);
+
+private:
+  MyFESystem<dim> fe;
+
+  struct VectorElementDestroyer {
+    const std::vector<const FiniteElement<dim>*> data;
+    VectorElementDestroyer (const std::vector<const FiniteElement<dim>*> &pointers);
+    ~VectorElementDestroyer (); // destructor to delete the pointers
+    const std::vector<const FiniteElement<dim>*> & get_data () const;
+  };
+
+  static std::vector<const FiniteElement<dim>*>
+  create_fe_list (const unsigned int polynomial_degree);
+
+  static std::vector<unsigned int>
+  create_fe_multiplicities ();
+};
+
+
+template <int dim>
+std::vector<const FiniteElement<dim>*>
+MySimulator<dim>::create_fe_list (const unsigned int polynomial_degree)
+{
+  std::vector<const FiniteElement<dim>*> fe_list;
+  fe_list.push_back (new FE_Q<dim>(1));
+  fe_list.push_back (new FE_Q<dim>(2));
+  fe_list.push_back (new FE_Q<dim>(3));
+  fe_list.push_back (new FE_Q<dim>(4));
+  deallog << "Created list of elements" << std::endl;
+  return fe_list;
+}
+
+template <int dim>
+std::vector<unsigned int>
+MySimulator<dim>::create_fe_multiplicities ()
+{
+  std::vector<unsigned int> multiplicities;
+  multiplicities.push_back (1);
+  multiplicities.push_back (2);
+  multiplicities.push_back (3);
+  multiplicities.push_back (4);
+  return multiplicities;
+}
+
+template <int dim>
+MySimulator<dim>::VectorElementDestroyer::
+VectorElementDestroyer (const std::vector<const FiniteElement<dim>*> &pointers)
+: data(pointers)
+{}
+
+template <int dim>
+MySimulator<dim>::VectorElementDestroyer::
+~VectorElementDestroyer ()
+{
+  for (unsigned int i=0; i<data.size(); ++i)
+    {
+      deallog << "Destroying element " << data[i]->get_name() << std::endl;
+      delete data[i];
+    }
+}
+
+template <int dim>
+const std::vector<const FiniteElement<dim>*> &
+MySimulator<dim>::VectorElementDestroyer::
+get_data () const
+{
+  return data;
+}
+
+
+template <int dim>
+MySimulator<dim>::MySimulator (const unsigned int polynomial_degree)
+:
+fe (VectorElementDestroyer(create_fe_list (polynomial_degree)).get_data(),
+    create_fe_multiplicities ())
+{}
+
+
+int
+main()
+{
+  std::ofstream logfile ("output");
+  deallog << std::setprecision(PRECISION);
+  deallog << std::fixed;
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  MySimulator<1>(1);
+  MySimulator<2>(1);
+  MySimulator<3>(1);
+
+  return 0;
+}
+
+
+
diff --git a/tests/fe/fe_system_from_list.output b/tests/fe/fe_system_from_list.output
new file mode 100644
index 0000000000..ac100ab274
--- /dev/null
+++ b/tests/fe/fe_system_from_list.output
@@ -0,0 +1,19 @@
+
+DEAL::Created list of elements
+DEAL::Constructing FESystem from list.
+DEAL::Destroying element FE_Q<1>(1)
+DEAL::Destroying element FE_Q<1>(2)
+DEAL::Destroying element FE_Q<1>(3)
+DEAL::Destroying element FE_Q<1>(4)
+DEAL::Created list of elements
+DEAL::Constructing FESystem from list.
+DEAL::Destroying element FE_Q<2>(1)
+DEAL::Destroying element FE_Q<2>(2)
+DEAL::Destroying element FE_Q<2>(3)
+DEAL::Destroying element FE_Q<2>(4)
+DEAL::Created list of elements
+DEAL::Constructing FESystem from list.
+DEAL::Destroying element FE_Q<3>(1)
+DEAL::Destroying element FE_Q<3>(2)
+DEAL::Destroying element FE_Q<3>(3)
+DEAL::Destroying element FE_Q<3>(4)
-- 
2.39.5