Use larger loop for OpenMP in ph-ph interaction calculation in C