f_max(const void *a, const void *b)
{
        return *(double *)a > *(double *)b ? 1 :
                        (*(double *)a < *(double *)b ? -1 : 0);
}