diff --git a/Mesh/mwis.hpp b/Mesh/mwis.hpp
index fe8ee1f887de9395f86ee368dd200dd175737827..68d0b7ca2508be2f831004a79bac867be1208b8a 100644
--- a/Mesh/mwis.hpp
+++ b/Mesh/mwis.hpp
@@ -528,6 +528,8 @@ public:
       selected_vertex.push_back(_vertex);
       new_weight += get(state.weight_map, _vertex);
 
+      removed.push_back(_vertex);
+
       std::pair<out_edge_iterator, out_edge_iterator> edges =
         out_edges(_vertex, state.graph);
       for (out_edge_iterator eit = edges.first; eit != edges.second; eit++) {