Methods for numerical computation of characteristic roots for delay differential equations: experimental comparison