“…Consequently, a sufficiently small step size is required to account for the lack of synchronization with events. To prevent this, several methods, including those of NETOMAC and PSCAD, employ linear interpolation to rollback the solution to the moment of switching before executing the corresponding switch actions [4], [14], [17], [19], [24], [25], [29]- [31], [53]- [57]. However, as the approximation error of linear interpolation is O(h 2 ), the local O(h 3 ) accuracy of the trapezoidal method is still not matched in the respective simulation steps.…”