Optimizing Hard Thresholding for Sparse Model Discovery