add test.sage, compare rho/pmin1/ecm