We propose a practical iterative remeshing algorithm for multi-material tetrahedral meshes which is solely based on simple local topological operations, such as edge collapse, flip, split and vertex smoothing. To do so, we exploit an intermediate implicit feature complex which reconstructs piecewise smooth multi-material boundaries made of surface patches, feature edges and corner vertices. Furthermore, we design specific feature-aware local remeshing rules which, combined with a moving least square projection, result in high quality isotropic meshes representing the input mesh at a user defined resolution while preserving important features. Our algorithm uses only topology-aware local operations, which allows us to process difficult input meshes such as self-intersecting ones. We evaluate our approach on a collection of examples and experimentally show that it is fast and scales well. Graphical abstractDisplay Omitted HighlightsRobust tetrahedral remeshing: can process self-intersecting and poor quality meshes.Multi-material tetrahedral meshes: high-quality segemented meshes at any resolution.Using local operators: edge collapse, flip, split and vertex smoothing.Feature preserving: feature-aware local rules with a moving least square projection.