Is normalize (getNormal mesh vi) going to give accurate results on a boundary edge between smoothing groups, where one welded vertex has two normals?