Remove Quadrature(unsigned int) from the public interface.

diff --git a/include/deal.II/base/quadrature.h b/include/deal.II/base/quadrature.h
index d71f739850ff12ea61cff9903524131994fb832d..c680be4c5ff5d4e4fb72da35dd6e48a1d3b728d4 100644 (file)
--- a/include/deal.II/base/quadrature.h
+++ b/include/deal.II/base/quadrature.h
@@ -130,14 +130,9 @@ public:
   using SubQuadrature = Quadrature<dim == 0 ? 0 : dim - 1>;
 
   /**
-   * Constructor.
-   *
-   * This constructor is marked as explicit to avoid involuntary accidents
-   * like in <code>hp::QCollection@<dim@> q_collection(3)</code> where
-   * <code>hp::QCollection@<dim@> q_collection(QGauss@<dim@>(3))</code> was
-   * meant.
+   * Default constructor.
    */
-  explicit Quadrature(const unsigned int n_quadrature_points = 0);
+  Quadrature();
 
   /**
    * Build this quadrature formula as the tensor product of a formula in a
@@ -336,6 +331,21 @@ public:
   get_tensor_basis() const;
 
 protected:
+  /**
+   * Constructor.
+   *
+   * This constructor is marked as explicit to avoid involuntary accidents
+   * like in <code>hp::QCollection@<dim@> q_collection(3)</code> where
+   * <code>hp::QCollection@<dim@> q_collection(QGauss@<dim@>(3))</code> was
+   * meant. Nonetheless, it is easy to accidentally write
+   * @code
+   *   Quadrature<dim> quadrature(3);
+   * @endcode
+   * where QGauss was meant. As a consequence, this constructor is `protected`
+   * and so only available to derived classes initializing their base class.
+   */
+  explicit Quadrature(const unsigned int n_quadrature_points);
+
   /**
    * List of quadrature points. To be filled by the constructors of derived
    * classes.

diff --git a/source/base/quadrature.cc b/source/base/quadrature.cc
index 0aade1577def7eb9caa5f91c1f4bb5f1ca430d17..d39f88611e7a3804b8694cc15659d98dda6d40c3 100644 (file)
--- a/source/base/quadrature.cc
+++ b/source/base/quadrature.cc
@@ -27,6 +27,11 @@ DEAL_II_NAMESPACE_OPEN
 
 
 #ifndef DOXYGEN
+template <>
+Quadrature<0>::Quadrature()
+  : is_tensor_product_flag(false)
+{}
+
 template <>
 Quadrature<0>::Quadrature(const unsigned int n_q)
   : quadrature_points(n_q)
@@ -37,6 +42,13 @@ Quadrature<0>::Quadrature(const unsigned int n_q)
 
 
 
+template <int dim>
+Quadrature<dim>::Quadrature()
+  : is_tensor_product_flag(dim == 1)
+{}
+
+
+
 template <int dim>
 Quadrature<dim>::Quadrature(const unsigned int n_q)
   : quadrature_points(n_q, Point<dim>())