Gradient calculation in unstructured mesh.
 

I've been writing my own research code using unstructured, pressure-based FV approach. I uses node centered scheme, edge-based data structure, and Gauss theorem. As I know, Gauss Theorem, Least square, or weighted least square method are generally used to calculate cell gradients. Can anyone give me some advice and information about differences of these methods?