Crash on flat after phong surface

Closes #17738 (closed)

The present fix considers a change in 'surface interpolation' relevant for recalculating polydata mapper internals. Previously, only 'surface representation' was taken into account.

The branch contains a unit test to demonstrates the problem.

If you consider integrating this, please check/verify that my CMake code is ok, I do not know how your "NO_DATA" markers work..