“…Moreover, for a particular order of the reduced model, a reduced model in terms of delay differential equations has in general the potential to be more accurate than a finite-dimensional approximation of the same order [9]. Therefore, delay-structure preserving methods, i.e., methods that preserve the infinite-dimensional nature of the time delay system during model reduction, *This research has been carried out in the HYDRA project, which has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No 675731. have gained much attention [9], [10], [11], [12], [13]. In many cases, however, there is a need to preserve additional properties during model reduction.…”