Patch by k3daevin: Implement at least compute the area of planar faces of 3d cells.

Is this robust w.r.t large dimensions of the domain? Maybe we should try to make this one relative too, o take into account possible cancellation errors