polynomial_test.exe: polynomial_test.g.$(OBJEXT) abort.g.$(OBJEXT) $(lib-base.g)
polynomial1d.exe : polynomial1d.g.$(OBJEXT) $(lib-base.g)
quadrature_test.exe: quadrature_test.g.$(OBJEXT) $(lib-base.g)
+quadrature_selector.exe: quadrature_selector.g.$(OBJEXT) $(lib-base.g)
reference.exe : reference.g.$(OBJEXT) abort.g.$(OBJEXT) $(lib-base.g)
table.exe : table.g.$(OBJEXT) $(lib-base.g)
tensor.exe : tensor.g.$(OBJEXT) $(lib-base.g)
anisotropic_2.exe : anisotropic_2.g.$(OBJEXT) $(lib-base.g)
hierarchical.exe : hierarchical.g.$(OBJEXT) $(lib-base.g)
-tests = logtest reference quadrature_test table tensor \
+tests = logtest reference quadrature_test quadrature_selector table tensor \
timer threads polynomial1d polynomial_test \
auto_derivative_function anisotropic_1 anisotropic_2 \
hierarchical data_out_base
--- /dev/null
+//---------------------------- quadrature_selector.cc ---------------------------
+// quadrature_selector.cc,v 1.18 2003/01/08 17:58:18 wolf Exp
+// Version:
+//
+// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//---------------------------- quadrature_selector.cc ---------------------------
+
+
+// make sure that the QuadratureSelector works for a selection of
+// arguments
+
+
+#include <fstream>
+
+#include <base/logstream.h>
+#include <base/quadrature_lib.h>
+#include <base/quadrature_selector.h>
+#include <string>
+#include <cmath>
+
+
+template <int dim>
+void check (const std::string &name,
+ const unsigned int order,
+ const Quadrature<dim> &q)
+{
+ Assert (QuadratureSelector<dim>(name, order).get_points() ==
+ q.get_points(),
+ ExcInternalError());
+ deallog << name << ' ' << order << " ok" << std::endl;
+}
+
+
+int main()
+{
+ std::ofstream logfile("quadrature_selector.output");
+ deallog.attach(logfile);
+ deallog.depth_console(0);
+
+ check ("gauss", 2, QGauss2<1>());
+ check ("gauss", 2, QGauss2<2>());
+ check ("gauss", 2, QGauss2<3>());
+
+ check ("gauss", 2, QGauss2<3>());
+ check ("gauss", 6, QGauss6<3>());
+ check ("gauss", 10, QGauss<3>(10));
+
+ check ("weddle", 0, QWeddle<2>());
+}
+
+