Comparing the normal will fail with small features or close to degenerate faces.
To match what SketchUp to consider two faces coplanar you’d have to get best fit plane for all the vertices in the two faces and then check if they are all on the plane.
Something like this: Bitbucket