Optimizing two-electron repulsion integral calculation on knights landing architecture

In this paper, we introduced some optimization methods we used to optimize two-electron repulsion integral calculation on Knights Landing architecture. We developed a schedule for parallelism and vectorization, and we compared two different methods for calculating lower incomplete gamma function. Our optimization achieved 1.7 speedup on KNL than CPU platform.