diff --git a/include/dem/find_cell_neighbors.h b/include/dem/find_cell_neighbors.h
index 9035a6da59..adf9c41ab4 100644
--- a/include/dem/find_cell_neighbors.h
+++ b/include/dem/find_cell_neighbors.h
@@ -25,15 +25,6 @@
 
 using namespace dealii;
 
-/**
- * @brief Finds the neighbors lists of all the active cells in the input triangulation.
- * find_cell_neighbors() is written to avoid any repetition, for instance if
- * cell B is recognized as the neighbor of cell A once, cell A will not appear
- * in the neighbor list of cell B again. On the other hand,
- * find_full_cell_neighbors() function finds the neighbors list with this
- * repetition.
- */
-
 /**
  * @brief Finds the neighbor list (without repetition) of all the active
  * cells in the triangulation. It gets the vertices of the cells to get lists
@@ -117,7 +108,6 @@ find_full_cell_neighbors(
   typename DEM::dem_data_structures<dim>::cells_total_neighbor_list
     &cells_total_neighbor_list);
 
-
 /**
  * @brief Generate a periodic neighbor cells list of the cell. With the
  * coinciding_vertex_groups and the vertex_to_coinciding_vertex_group, we can
diff --git a/source/dem/find_cell_neighbors.cc b/source/dem/find_cell_neighbors.cc
index d6df7b97b5..c2286038e2 100644
--- a/source/dem/find_cell_neighbors.cc
+++ b/source/dem/find_cell_neighbors.cc
@@ -391,7 +391,6 @@ find_cell_periodic_neighbors<3>(
   typename dem_data_structures<3>::cells_neighbor_list
     &cells_ghost_local_periodic_neighbor_list);
 
-
 template void
 find_full_cell_neighbors<2>(
   const parallel::distributed::Triangulation<2> &triangulation,
@@ -403,23 +402,3 @@ find_full_cell_neighbors<3>(
   const parallel::distributed::Triangulation<3> &triangulation,
   typename dem_data_structures<3>::cells_total_neighbor_list
     &cells_total_neighbor_list);
-
-template void
-get_periodic_neighbor_list<2>(
-  const typename Triangulation<2>::active_cell_iterator &cell,
-  const std::map<unsigned int, std::vector<unsigned int>>
-                                             &coinciding_vertex_groups,
-  const std::map<unsigned int, unsigned int> &vertex_to_coinciding_vertex_group,
-  const std::vector<std::set<typename Triangulation<2>::active_cell_iterator>>
-                                               &v_to_c,
-  typename dem_data_structures<2>::cell_vector &periodic_neighbor_list);
-
-template void
-get_periodic_neighbor_list<3>(
-  const typename Triangulation<3>::active_cell_iterator &cell,
-  const std::map<unsigned int, std::vector<unsigned int>>
-                                             &coinciding_vertex_groups,
-  const std::map<unsigned int, unsigned int> &vertex_to_coinciding_vertex_group,
-  const std::vector<std::set<typename Triangulation<3>::active_cell_iterator>>
-                                               &v_to_c,
-  typename dem_data_structures<3>::cell_vector &periodic_neighbor_list);
diff --git a/source/dem/lagrangian_post_processing.cc b/source/dem/lagrangian_post_processing.cc
index 128d993ab0..3fe71ecb0b 100644
--- a/source/dem/lagrangian_post_processing.cc
+++ b/source/dem/lagrangian_post_processing.cc
@@ -277,43 +277,3 @@ write_post_processing_results<3>(
   const unsigned int                             step_number,
   const MPI_Comm                                &mpi_communicator,
   AdaptiveSparseContacts<3>                     &sparse_contacts_object);
-
-template void
-calculate_average_particles_velocity<2>(
-  const parallel::distributed::Triangulation<2> &triangulation,
-  const Particles::ParticleHandler<2>           &particle_handler,
-  Vector<double>                                &velocity_average_x,
-  Vector<double>                                &velocity_average_y,
-  Vector<double>                                &velocity_average_z,
-  Vector<double>                                &velocity_average_magnitude);
-
-template void
-calculate_average_particles_velocity<3>(
-  const parallel::distributed::Triangulation<3> &triangulation,
-  const Particles::ParticleHandler<3>           &particle_handler,
-  Vector<double>                                &velocity_average_x,
-  Vector<double>                                &velocity_average_y,
-  Vector<double>                                &velocity_average_z,
-  Vector<double>                                &velocity_average_magnitude);
-
-template void
-calculate_average_granular_temperature<2>(
-  const parallel::distributed::Triangulation<2> &triangulation,
-  const Particles::ParticleHandler<2>           &particle_handler,
-  Vector<double>                                &granular_temperature_average);
-
-template void
-calculate_average_granular_temperature<3>(
-  const parallel::distributed::Triangulation<3> &triangulation,
-  const Particles::ParticleHandler<3>           &particle_handler,
-  Vector<double>                                &granular_temperature_average);
-
-template Tensor<1, 2>
-calculate_cell_average_particles_velocity<2>(
-  const typename parallel::distributed::Triangulation<2>::cell_iterator &cell,
-  const Particles::ParticleHandler<2> &particle_handler);
-
-template Tensor<1, 3>
-calculate_cell_average_particles_velocity<3>(
-  const typename parallel::distributed::Triangulation<3>::cell_iterator &cell,
-  const Particles::ParticleHandler<3> &particle_handler);