find_poly_roots
nroots = find_poly_roots(bch, 1, bch->elp, errloc);
cnt += find_poly_roots(bch, k+1, f1, roots);
cnt += find_poly_roots(bch, k+1, f2, roots+cnt);