Find the root of pi in the outer loop and update it on the fly in the inner loop. Saves having to call fof_find at every iteration of the inner loop.