+ SaveFacet savedFacet( myCurFace );
+
+ // get an origin not lying in plane of any facet
+ int faceIndex = std::max( myCurFace.myIndex, 0 );
+ setFace( faceIndex );
+ double minProj, maxProj;
+ XYZ normal, origin;
+ projectNodesToNormal( faceIndex, minProj, maxProj, normal.data());
+ GetFaceBaryCenter( faceIndex, origin.x, origin.y, origin.z );
+ origin.x += normal.x * ( maxProj - minProj ) * 1e-1;
+ origin.y += normal.y * ( maxProj - minProj ) * 1e-2;
+ origin.z += normal.z * ( maxProj - minProj ) * 1e-3;
+