compute.PFN.par function

Parallelized PFN computation