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

All sugestions pass the numerical test, although this variation needs no square roots:
const double v03normal = v03 * normal;
if (std::abs(v03normal * v03normal / ((normal * normal) * (v03 * v03))) >= 1e-12)