progressive Rutishauser QD algorithm