The kernel polynomial method based on Jacobi polynomials

The kernel polynomial method based on Jacobi polynomials Pn(α,β)(x) is proposed. The optimal-resolution positivity-preserving kernels and the corresponding damping factors are calculated. The results provide a generalization of the Jackson damping factors for arbitrary Jacobi polynomials. For α=±1/2, β=±1/2 (Chebyshev polynomials of the first to fourth kinds), explicit trigonometric expressions for the damping factors are obtained. The resulting algorithm can be easily introduced into existing implementations of the kernel polynomial method.