Remove remaining deprecations

doc/news/changes/incompatibilities/20220629Arndt

@@ -0,0 +1,4 @@
+Removed: The deprecated QTrapezoid, ParticleHandler::locally_relevant_ids(),
+and Arkode::reinit_vector have been removed.
+<br>
+(Daniel Arndt, 2022/06/29)

include/deal.II/base/quadrature_lib.h

@@ -127,22 +127,6 @@ class QTrapezoid : public Quadrature<dim>
 };
 
 
-/**
- * An alias for QTrapezoid available for historic reasons. This name is
- * deprecated.
- *
- * The class was originally named QTrapez, a poorly named choice since the
- * proper name of the quadrature formula
- * is "trapezoidal rule", or sometimes also called the "trapezoid rule". The
- * misnomer resulted from the fact that its original authors' poor English
- * language skills led them to translate the name incorrectly from the German
- * "Trapezregel".
- */
-template <int dim>
-using QTrapez DEAL_II_DEPRECATED = QTrapezoid<dim>;
-
-
-
 /**
  * The Milne rule for numerical quadrature formula. The Milne rule is a closed
  * Newton-Cotes formula and is exact for polynomials of degree 5.

include/deal.II/fe/fe_tools.templates.h

@@ -2489,8 +2489,7 @@ namespace FETools
                     // find sub-quadrature
                     position = name.find('(');
                     const std::string subquadrature_name(name, 0, position);
-                    AssertThrow(subquadrature_name == "QTrapez" ||
-                                  subquadrature_name == "QTrapezoid",
+                    AssertThrow(subquadrature_name == "QTrapezoid",
                                 ExcNotImplemented(
                                   "Could not detect quadrature of name " +
                                   subquadrature_name));

include/deal.II/grid/reference_cell.h

@@ -236,7 +236,7 @@ class ReferenceCell
    *
    * @note The weights of the quadrature object are left unfilled and
    *   consequently the object cannot usefully be used for actually
-   *   computing integrals. This is in contrast to, for example, the QTrapez
+   *   computing integrals. This is in contrast to, for example, the QTrapezoid
    *   class that correctly sets quadrature weights.
    */
   template <int dim>

include/deal.II/particles/particle_handler.h

@@ -654,27 +654,6 @@ namespace Particles
     types::particle_index
     get_next_free_particle_index() const;
 
-    /**
-     * Extract an IndexSet with global dimensions equal to
-     * get_next_free_particle_index(), containing the locally owned
-     * particle indices.
-     *
-     * This function can be used to construct distributed vectors and matrices
-     * to manipulate particles using linear algebra operations.
-     *
-     * Notice that it is the user's responsibility to guarantee that particle
-     * indices are unique, and no check is performed to verify that this is the
-     * case, nor that the union of all IndexSet objects on each mpi process is
-     * complete.
-     *
-     * @return An IndexSet of size get_next_free_particle_index(), containing
-     * n_locally_owned_particle() indices.
-     *
-     * @deprecated Use locally_owned_particle_ids() instead.
-     */
-    DEAL_II_DEPRECATED IndexSet
-    locally_relevant_ids() const;
-
     /**
      * Extract an IndexSet with global dimensions equal to
      * get_next_free_particle_index(), containing the locally owned
@@ -1403,15 +1382,6 @@ namespace Particles
       output_vector.compress(VectorOperation::insert);
   }
 
-
-
-  template <int dim, int spacedim>
-  inline IndexSet
-  ParticleHandler<dim, spacedim>::locally_relevant_ids() const
-  {
-    return this->locally_owned_particle_ids();
-  }
-
 } // namespace Particles
 
 DEAL_II_NAMESPACE_CLOSE

include/deal.II/sundials/arkode.h

@@ -569,17 +569,6 @@ namespace SUNDIALS
     void *
     get_arkode_memory() const;
 
-    /**
-     * A function object that was used to `reinit` the given vector. Setting
-     * this field does no longer have any effect and all auxiliary vectors are
-     * reinit-ed automatically based on the user-supplied vector in solve_ode().
-     *
-     * @deprecated This function is no longer used and can be safely removed in
-     *   user code.
-     */
-    DEAL_II_DEPRECATED
-    std::function<void(VectorType &)> reinit_vector;
-
     /**
      * A function object that users may supply and that is intended to compute
      * the explicit part of the IVP right hand side. Sets $explicit_f = f_E(t,

source/grid/manifold_lib.cc

@@ -1689,10 +1689,10 @@ TransfiniteInterpolationManifold<dim, spacedim>::initialize(
   //
   // In the co-dimension one case (meaning  dim < spacedim) we have to fall
   // back to a simple GridTools::affine_cell_approximation<dim>() which
-  // requires 2^dim points, instead. Thus, initialize the QIteraded
+  // requires 2^dim points, instead. Thus, initialize the QIterated
   // quadrature with no subdivisions.
   std::vector<Point<dim>> unit_points =
-    QIterated<dim>(QTrapez<1>(), (dim == spacedim ? 2 : 1)).get_points();
+    QIterated<dim>(QTrapezoid<1>(), (dim == spacedim ? 2 : 1)).get_points();
   std::vector<Point<spacedim>> real_points(unit_points.size());
 
   for (const auto &cell : triangulation.active_cell_iterators())

source/sundials/arkode.cc

@@ -679,10 +679,6 @@ namespace SUNDIALS
   void
   ARKode<VectorType>::set_functions_to_trigger_an_assert()
   {
-    reinit_vector = [](VectorType &) {
-      AssertThrow(false, ExcFunctionNotProvided("reinit_vector"));
-    };
-
     solver_should_restart = [](const double, VectorType &) -> bool {
       return false;
     };

tests/bits/point_value_03.cc

@@ -124,7 +124,7 @@ check()
   parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
   make_mesh(tria);
 
-  FE_Q<dim>       element(QIterated<1>(QTrapez<1>(), 3));
+  FE_Q<dim>       element(QIterated<1>(QTrapezoid<1>(), 3));
   DoFHandler<dim> dof(tria);
   dof.distribute_dofs(element);
 

tests/matrix_free/tensor_product_evaluate_01.cc

@@ -65,7 +65,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;

tests/matrix_free/tensor_product_evaluate_02.cc

@@ -60,7 +60,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;

tests/matrix_free/tensor_product_evaluate_03.cc

@@ -61,7 +61,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;

tests/matrix_free/tensor_product_evaluate_04.cc

@@ -58,7 +58,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;

tests/matrix_free/tensor_product_evaluate_05.cc

@@ -59,7 +59,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;

tests/matrix_free/tensor_product_evaluate_06.cc

@@ -55,7 +55,7 @@ test(const unsigned int degree)
 
   const std::vector<Point<dim>> evaluation_points =
     dim == 3 ? QGauss<dim>(2).get_points() :
-               QIterated<dim>(QTrapez<1>(), 3).get_points();
+               QIterated<dim>(QTrapezoid<1>(), 3).get_points();
 
   deallog << "Evaluate in " << dim << "d with polynomial degree " << degree
           << std::endl;