an implementation of the iterative hard thresholding method for cardinality constrained least squares