diff --git a/src/Meshing/Mage.py b/src/Meshing/Mage.py
index 43d5cdb6a953f2e97f13c24f0d43a40633758d33..5a32f6149613742305de2040c43ce4e9b2e19cbe 100644
--- a/src/Meshing/Mage.py
+++ b/src/Meshing/Mage.py
@@ -49,6 +49,7 @@ class MeshingWithMage(AMeshingTool):
         self._init_c_set_bief_name()
         self._init_c_st_to_m_compl()
         self._init_c_interpolate_profils_pas_transversal()
+        self._init_c_purge()
         self._init_c_output_bief()
 
     def _lib_path(self):
@@ -104,6 +105,11 @@ class MeshingWithMage(AMeshingTool):
         ]
         self._c_interpolate_profils_pas_transversal.restype = None
 
+    def _init_c_purge(self):
+        self._c_purge = getattr(self._bief_lib, 'c_purge')
+        self._c_purge.argtypes = None
+        self._c_purge.restype = None
+
     def _init_c_output_bief(self):
         self._c_output_bief = getattr(self._bief_lib, 'c_output_bief')
         self._c_output_bief.argtypes = [c_char_p]
@@ -184,6 +190,7 @@ class MeshingWithMage(AMeshingTool):
             gl = reach.compute_guidelines()
             self.complete_cross_section(gl[0])
             self.interpolate_cross_section(ns, step)
+            self.purge()
 
             self.export_to_m(m_file)
 
@@ -226,6 +233,9 @@ class MeshingWithMage(AMeshingTool):
             step
         )
 
+    def purge(self):
+        self._c_purge()
+
     def export_to_m(self, m):
         self.output_bief(m)