CFD Online Discussion Forums - How to implement PISO algorithm based on the colocated grid

I have written my lid-driven code with SIMPLE on MATLAB following the instruction from Ferziger's book. And I want to change it to get a PISO case, however, it occured to me that due to previous work based on the colocated grid, the $\tilde{u_{f}'}=-\sum_{}^{}A_{l}u_{l}'/A_{P}$ cant be obtained directly on the colocated grid. So I think that I need to interpolate the nodal value to get that.

The question is :
Can $\tilde{u_{f}'}$ be calculated by simple interpolation or something like momentum interpolation? Or PISO requires a different way to implement on the colocated grid, which is not same with the description in Ferziger's book?

I will appreciate it if you can provide any suggetion or tutorial.
Thank you in advance!