Some code for potential accurate calculation of cell center

(the code is commented out)