sor_heap() now returns the number of bisections under which the calculation was carried out as well.