Complex modes computation for the low-rank non-proportionally damped structures by the homotopy method