Variable selection is central to high-dimensional data analysis, and various algorithms have been developed. Ideally, a variable algorithm shall be flexible, scalable, with theoretical guarantee, yet most existing cannot attain these properties at the same time. In this article, three-step developed, involving kernel-based estimation of regression function its gradient functions as well hard th...