Revised so that the output of sor_bisect... also includes the aggregation parameter psi.