removed separate threshold tolerance, now the criterion is for |f| < ftol for omega replaced with Re[omega]