core: Sanitize script interface shared pointers

src/core/accumulators.cpp

@@ -62,9 +62,17 @@ int auto_update_next_update() {
 
 void auto_update_add(AccumulatorBase *acc) {
   assert(acc);
+  assert(std::find_if(auto_update_accumulators.begin(),
+                      auto_update_accumulators.end(), [acc](auto const &a) {
+                        return a.acc == acc;
+                      }) == auto_update_accumulators.end());
   auto_update_accumulators.emplace_back(acc);
 }
 void auto_update_remove(AccumulatorBase *acc) {
+  assert(std::find_if(auto_update_accumulators.begin(),
+                      auto_update_accumulators.end(), [acc](auto const &a) {
+                        return a.acc == acc;
+                      }) != auto_update_accumulators.end());
   auto_update_accumulators.erase(
       boost::remove_if(
           auto_update_accumulators,

src/core/constraints/Constraints.hpp

@@ -23,6 +23,7 @@
 #include "grid.hpp"
 #include "statistics.hpp"
 
+#include <algorithm>
 #include <memory>
 #include <vector>
 
@@ -51,11 +52,15 @@ template <class ParticleRange, class Constraint> class Constraints {
     if (not c->fits_in_box(box_geo.length())) {
       throw std::runtime_error("Constraint not compatible with box size.");
     }
+    assert(std::find(m_constraints.begin(), m_constraints.end(), c) ==
+           m_constraints.end());
 
     m_constraints.emplace_back(c);
     on_constraint_change();
   }
   void remove(std::shared_ptr<Constraint> const &c) {
+    assert(std::find(m_constraints.begin(), m_constraints.end(), c) !=
+           m_constraints.end());
     m_constraints.erase(
         std::remove(m_constraints.begin(), m_constraints.end(), c),
         m_constraints.end());

src/core/grid_based_algorithms/lb_boundaries.cpp

@@ -59,12 +59,16 @@ std::vector<std::shared_ptr<LBBoundary>> lbboundaries;
 #if defined(LB_BOUNDARIES) || defined(LB_BOUNDARIES_GPU)
 
 void add(const std::shared_ptr<LBBoundary> &b) {
+  assert(std::find(lbboundaries.begin(), lbboundaries.end(), b) ==
+         lbboundaries.end());
   lbboundaries.emplace_back(b);
 
   on_lbboundary_change();
 }
 
 void remove(const std::shared_ptr<LBBoundary> &b) {
+  assert(std::find(lbboundaries.begin(), lbboundaries.end(), b) !=
+         lbboundaries.end());
   lbboundaries.erase(std::remove(lbboundaries.begin(), lbboundaries.end(), b),
                      lbboundaries.end());