REF: add is_sorted option to fdr to avoid duplicate calculations
     delete intermediate arrays (all work except Hommel)