Least-Squared Optimized Polynomials for Smeared Link Actions
  We introduce a numerical method for generating the approximating polynomials used in fermionic calculations with smeared link actions. We investigate the stability of the algorithm and determine the optimal weight function and the optimal type of discretization. The achievable order of polynomial approximation reaches several thousands allowing fermionic calculations using the Hypercubic Smeared Link action even with physical quark masses. 
S.D. Katz, B.C. Toth